diff options
author | debashisdeb | 2014-06-21 00:52:25 +0530 |
---|---|---|
committer | debashisdeb | 2014-06-21 00:52:25 +0530 |
commit | 7c756fcc12d21693818e58f6936cab5b7c112868 (patch) | |
tree | 009cb02ec85f4a75ac7b64239751f15361df2bfe /Introduction_To_Chemical_Engineering/ch5.ipynb | |
parent | 83c1bfceb1b681b4bb7253b47491be2d8b2014a1 (diff) | |
download | Python-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.ipynb | 93 |
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": {}, |