summaryrefslogtreecommitdiff
path: root/Introduction_To_Chemical_Engineering/ch5.ipynb
diff options
context:
space:
mode:
authordebashisdeb2014-06-21 00:52:25 +0530
committerdebashisdeb2014-06-21 00:52:25 +0530
commit7c756fcc12d21693818e58f6936cab5b7c112868 (patch)
tree009cb02ec85f4a75ac7b64239751f15361df2bfe /Introduction_To_Chemical_Engineering/ch5.ipynb
parent83c1bfceb1b681b4bb7253b47491be2d8b2014a1 (diff)
downloadPython-Textbook-Companions-7c756fcc12d21693818e58f6936cab5b7c112868.tar.gz
Python-Textbook-Companions-7c756fcc12d21693818e58f6936cab5b7c112868.tar.bz2
Python-Textbook-Companions-7c756fcc12d21693818e58f6936cab5b7c112868.zip
Removed Problem Statements Completely
Diffstat (limited to 'Introduction_To_Chemical_Engineering/ch5.ipynb')
-rw-r--r--Introduction_To_Chemical_Engineering/ch5.ipynb93
1 files changed, 70 insertions, 23 deletions
diff --git a/Introduction_To_Chemical_Engineering/ch5.ipynb b/Introduction_To_Chemical_Engineering/ch5.ipynb
index 6af51e50..24e997c2 100644
--- a/Introduction_To_Chemical_Engineering/ch5.ipynb
+++ b/Introduction_To_Chemical_Engineering/ch5.ipynb
@@ -1,6 +1,7 @@
{
"metadata": {
- "name": ""
+ "name": "",
+ "signature": "sha256:a84d51472594c75f9afdf6cd721b95cb9d112fa36b409d534c3bc68aa928266a"
},
"nbformat": 3,
"nbformat_minor": 0,
@@ -27,18 +28,21 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"A=5.*4 #in m2\n",
"T1=100.; #in K\n",
"T2=30.; #in K\n",
"\n",
+ "# Calculations\n",
"delta_T=T1-T2;\n",
"\n",
"x=0.25 #in m\n",
"k=0.70 #in W/mK\n",
"Q=k*A*(delta_T/x);\n",
"\n",
+ "# Results\n",
"print \"rate of heat loss = %f W\"%(Q)\n"
],
"language": "python",
@@ -66,12 +70,14 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"d1=0.15 #in m\n",
"d2=0.16 #in m\n",
"l=1. #in m\n",
"\n",
+ "# Calculations\n",
"A1=3.14*d1*l;\n",
"A2=3.14*d2*l\n",
"Am=(A1-A2)/math.log (A1/A2);\n",
@@ -84,6 +90,7 @@
"k=50. #in W/mK\n",
"Q=k*Am*(delta_T/x);\n",
"\n",
+ "# Results\n",
"print \"rate of heat loss per unit length = %f W/m\"%(Q)\n"
],
"language": "python",
@@ -111,13 +118,15 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"ri=0.5 #in m\n",
"ro=0.6; #in m\n",
"A1=4*3.14*ri**2;\n",
"A2=4*3.14*ro**2;\n",
"\n",
+ "# Calculations\n",
"Am=(A1*A2)**0.5;\n",
"\n",
"Ti=140.; #in K\n",
@@ -128,6 +137,7 @@
"\n",
"Q=k*Am*(delta_T/x);\n",
"\n",
+ "# Results\n",
"print \"Heat loss through sphere = %f W\"%(Q)\n"
],
"language": "python",
@@ -155,13 +165,16 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"x1=0.250; #in m\n",
"k1=0.7; #in W/mK\n",
"A1=1.; #in m2\n",
"R1=x1/(k1*A1); #in K/W\n",
"\n",
+ "# Calculations and Results\n",
+ "#for the felt layer\n",
"x2=0.020; #in m\n",
"k2=0.046; #in W/mK\n",
"A2=1.; #in m2\n",
@@ -201,8 +214,9 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"d1=0.15 #in m\n",
"d2=0.16 #in m\n",
"l=1. #in m\n",
@@ -213,6 +227,8 @@
"k1=50. #in W/mK\n",
"R1=x1/(k1*Am1);\n",
"\n",
+ "# Calculations and Results\n",
+ "#resistance by insulation\n",
"d2=0.16 #in m\n",
"d3=0.26 #in m\n",
"l=1. #in m\n",
@@ -259,8 +275,9 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"x1=0.1; #in m\n",
"x2= 0.25; #in m\n",
"k_rb=0.93; #in W/mK\n",
@@ -268,6 +285,8 @@
"k_al=203.6 #in W/mK\n",
"A=0.1 #in m2\n",
"\n",
+ "# Calculations and Results\n",
+ "#to find resistance without rivets\n",
"R=(1/A)*((x1/k_rb)+(x2/k_ib));\n",
"T1=225 #in K\n",
"T2=37 #in K\n",
@@ -275,6 +294,7 @@
"Q=delta_T/R;\n",
"print \"heat transfer rate = %f W\"%(Q)\n",
"\n",
+ "#to find resistance with rivet\n",
"d=0.03 #in m\n",
"rivet_area= (3.14/4)*d**2;\n",
"R_r=(x1+x2)/(k_al*rivet_area);\n",
@@ -314,20 +334,17 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "'''\n",
- "calculate the heat tranfer coefficient based on \n",
- "i) the arithmetic mean difference between the temperatures of the water and the wall of the tube\n",
- "ii) the logarithmic mean difference between the temperatures of the water and the wall of the tube.\n",
- "'''\n",
- "\n",
+ " \n",
"import math\n",
"\n",
+ "# variables\n",
"Cp = 4.178 # kJ/kg K for water\n",
"q = 1838. # rate at which heat is transfered\n",
"A = .1005 # heat transfer area\n",
"dt1 = 80. - 24 # temperature diffference at hot end\n",
"dt2 = 36.-24 # temperature difference at cold end\n",
"\n",
+ "# Calculations and Results\n",
"dtm = (56 + 12)/2.0\n",
"h = q/(A*dtm)\n",
"print \"Heat transfer coefficient, h = %.0f W/m**2 K\"%h\n",
@@ -363,9 +380,10 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
"\n",
+ "# Variables\n",
"density=984.1 #in kg/cubic meter\n",
"v=3. #in m/s\n",
"viscosity=485*10**-6; #in Pa-s\n",
@@ -373,12 +391,15 @@
"cp=4178. #in J/kg K\n",
"d=0.016 #in m\n",
"\n",
+ "# Calculations and Results\n",
"Re=(density*v*d)/viscosity;\n",
"Pr=(cp*viscosity)/k;\n",
"\n",
+ "#dittus boelter equation\n",
"h=0.023*Re**0.8*Pr**0.3*(k/d);\n",
"print \"heat transfer coefficient = %f W/sq meter K\"%(h)\n",
"\n",
+ "#Sieder Tate equation\n",
"viscosity_w=920*10**-6.\n",
"h1=0.023*Re**0.8*Pr**(1./3)*(k/d)*(viscosity/viscosity_w)**0.14;\n",
"print \"heat transfer coefficient = %f W/sq meter K\"%(h1)\n"
@@ -409,14 +430,17 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"T_sun = 5973 #in degree C\n",
"d = 1.5*10**13 #in cm\n",
"R = 7.1*10**10; #in cm\n",
"\n",
+ "# Calculations\n",
"T_earth = ((R/(2*d))**0.5)*T_sun;\n",
"\n",
+ "# Results\n",
"print \"Temperature of earth = %f C\"%(T_earth-273) \n"
],
"language": "python",
@@ -444,14 +468,17 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"R=7*10**10; #in cm\n",
"Ts=6000; #in K\n",
"\n",
+ "# Calculations\n",
"l=1.5*10**13; #in m\n",
"To=((R**2/(4*l**2))**0.25)*Ts;\n",
"\n",
+ "# Results\n",
"print \"temperature of earth = %f K\"%(To)\n"
],
"language": "python",
@@ -479,14 +506,17 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"R=6.92*10**5 #in km\n",
"l=14.97*10**7 #in km\n",
"Ts=6200; #in K\n",
"\n",
+ "# Calculations\n",
"To=(R**2/l**2)**0.25*Ts;\n",
"\n",
+ "# Results\n",
"print \"Equilibrium temperature = %f K\"%(To)\n"
],
"language": "python",
@@ -514,15 +544,18 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"view_factor=0.5;\n",
"R=6.92*10**5 #in km\n",
"l=14.97*10**7 #in km\n",
"Ts=6200; #in K\n",
"\n",
+ "# Calculations\n",
"To=(view_factor*(R**2/l**2))**0.25*Ts;\n",
"\n",
+ "# Results\n",
"print \"Equilibrium temperature = %f K\"%(To)\n",
"\n",
"\n"
@@ -552,8 +585,9 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"view_factor=0.25;\n",
"R=7.1*10**10 #in cm\n",
"l=1.5*10**13 #in cm\n",
@@ -561,10 +595,12 @@
"alpha=0.2;\n",
"epsilon=0.1;\n",
"\n",
+ "# Calculations\n",
"ratio=alpha/epsilon;\n",
"To=(ratio*view_factor*(R**2/l**2))**0.25*Ts;\n",
"\n",
"\n",
+ "# Results\n",
"print \"Equilibrium temperature = %f K\"%(To)\n"
],
"language": "python",
@@ -592,15 +628,18 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"R=7*10**10; #in cm\n",
"l=1.5*10**13; #in cm\n",
"sigma=5.3*10**-5; #in erd/s(cm2)(K)4\n",
"T=6000; #in K\n",
"\n",
+ "# Calculations\n",
"S=(R/l)**2*(sigma)*(T**4)*60;\n",
"\n",
+ "# Results\n",
"print \"solar constant = %f J/sq cm min\"%(S/10**7)\n"
],
"language": "python",
@@ -628,12 +667,14 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"F = 5000. #in kg/hr\n",
"xF = 0.01\n",
"xL = 0.02;\n",
"\n",
+ "# Calculations and Results\n",
"L = F*xF/xL;\n",
"V = F-L;\n",
"print \"L = %f Kg/hr V = %f kg/hr\"%(L,V)\n",
@@ -682,15 +723,18 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
"from numpy import *\n",
+ "# Variables\n",
"b1 = 6000*125.79+3187.56*2691.5-3187.56*461.30; #data from previous problem\n",
"b2 = 6000;\n",
"A = array([[419.04, 2676.1],[1, 1]])\n",
"\n",
+ "# Calculations and Results\n",
"b = array([[b1],[b2]]);\n",
"x = linalg.solve(A,b)\n",
+ "#x = x*b\n",
"L = x[0];\n",
"V = x[1];\n",
"\n",
@@ -727,11 +771,13 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "\n",
+ " \n",
"import math \n",
+ "# Variables\n",
"Hv=2635.3 #kJ/kg\n",
"hL=313.93 #in kJ/kg\n",
"\n",
+ "# Calculations and Results\n",
"S=(2500*313.93+2500*2635.3-5000*125.79)/(2691.5-461.30);\n",
"print \"steam flow rate = %f kg steam/hr\"%(S)\n",
"\n",
@@ -744,6 +790,7 @@
"print \"Area = %f sq meter\"%(A)\n",
"print \"in this case a condensor and vaccum pump should be used\"\n",
"\n",
+ "# Note : there is mistake in calculation in Book. Please calculate manually."
],
"language": "python",
"metadata": {},