diff options
Diffstat (limited to 'Engineering_Heat_Transfer_by_W_S_Janna/CHAPTER8.ipynb')
| -rwxr-xr-x | Engineering_Heat_Transfer_by_W_S_Janna/CHAPTER8.ipynb | 481 |
1 files changed, 481 insertions, 0 deletions
diff --git a/Engineering_Heat_Transfer_by_W_S_Janna/CHAPTER8.ipynb b/Engineering_Heat_Transfer_by_W_S_Janna/CHAPTER8.ipynb new file mode 100755 index 00000000..759f39bd --- /dev/null +++ b/Engineering_Heat_Transfer_by_W_S_Janna/CHAPTER8.ipynb @@ -0,0 +1,481 @@ +{
+ "metadata": {
+ "name": ""
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter 8: Natural Convection System"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 8.1 Page No. 413"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "rou= 0.0551 # density in Ibm/cu.ft \n",
+ "cp=0.2420 # specific heat BTU/(lbm-degree Rankine) \n",
+ "v= 27.88e-5 # viscosity in sq.ft/s \n",
+ "kf = 0.01944 # thermal conductivity in BTU/(hr.ft.degree Rankine) \n",
+ "a = 1.457 # diffusivity in sq.ft/hr \n",
+ "Pr = 0.689 # Prandtl Number\n",
+ "T_inf=120.0+460.0 # wall temperature in degree R\n",
+ "Tw=400.0+460.0 # inside wall temperature in degree R\n",
+ "Beta=1/T_inf\n",
+ "\n",
+ "Beta_=0.00116\n",
+ "gc=32.2\n",
+ "L=1.0 # length of wall in ft\n",
+ "W=2.0 # width in ft\n",
+ "Gr=(gc*Beta_*(Tw-T_inf)*L**3)/v**2 # Grashof Number\n",
+ "temperature_slope=0.505 #temperature slope from table 8.1 \n",
+ "hL=(kf/L)*(4/3.0)*(Gr/4.0)**(1/4.0)*temperature_slope \n",
+ "A=L*W # cross sectional area in sq.ft\n",
+ "qw=hL*A*(Tw-T_inf)\n",
+ "\n",
+ "print\"The heat transferred is\",round(qw,0),\"BTU/hr\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat transferred is 558.0 BTU/hr\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 8.2 Page No. 414"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "rou1=1.295 # density in kg/cu.m\n",
+ "cp1=1005.5 # specific heat in J/(kg*K) \n",
+ "v1=12.59e-6 # viscosity in sq.m/s \n",
+ "Pr1=0.713 # Prandtl Number \n",
+ "kf1=0.02426 # thermal conductivity in W/(m.K)\n",
+ "a1=0.17661e-4 # diffusivity in sq.m/s \n",
+ "T_inf1=0 # inside and outside temperature in K\n",
+ "Beta1=1/(T_inf1+273.0) # volumetric thermal expansion coefficient at 295 K and 273 K\n",
+ "\n",
+ "rou2=1.177 # density in kg/cu.m\n",
+ "cp2=1005 # specific heat in J/(kg*K) \n",
+ "v2=15.68e-6 # viscosity in sq.m/s \n",
+ "Pr2=0.708 # Prandtl Number \n",
+ "kf2=0.02624 # thermal conductivity in W/(m.K)\n",
+ "a2=0.22160e-4 # diffusivity in sq.m/s \n",
+ "T_inf2=22.0 # inside and outside temperature in K\n",
+ "Beta2=1/(T_inf2+273.0) # volumetric thermal expansion coefficient at 295 K and 273 K\n",
+ "\n",
+ "g=9.81\n",
+ "t=0.005 # thickness of glass\n",
+ "L=0.60 # window length in m\n",
+ "k=0.81 # thermal conductivity of glass from appendix table B3\n",
+ "Tw1=18\n",
+ "Tw2=4\n",
+ "Ra1=(g*Beta1*(Tw2-T_inf1)*L**3)/(v1*a1)\n",
+ "hL1=(kf1/L)*(0.68+((0.67*((abs(Ra1)))**(1/4.0))/(1+(0.492/Pr1)**(9/16.0))**(4/9.0)))\n",
+ "Ra2=(g*Beta2*(Tw1-T_inf2)*L**3)/(v2*a2)\n",
+ "hL2=(kf2/L)*(0.68+((0.67*(abs(Ra2))**(1/4.0))/(1+(0.492/Pr2)**(9/16.0))**(4/9.0)))\n",
+ "q1=(T_inf1-T_inf2)/((1/hL2)+(t/k)+(1/hL1))\n",
+ "Tw2_=T_inf2-(q1/hL2)\n",
+ "Tw1_=q1/hL1+T_inf1\n",
+ "\n",
+ "Ra1_=3.7*10**8\n",
+ "hL1_=2.92\n",
+ "Ra2_=2.31*10**8\n",
+ "hL2_=2.80\n",
+ "q2=(T_inf2-T_inf1)/((1/hL2_)+(t/k)+(1/hL1_))\n",
+ "\n",
+ "Tw2final=q2-T_inf2\n",
+ "Tw1final=10.7\n",
+ "\n",
+ "print\"The heat loss is \",round(q2,1),\" W/sq.m\"\n",
+ "\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat loss is 31.2 W/sq.m\n"
+ ]
+ }
+ ],
+ "prompt_number": 15
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 8.3 Page No.419"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "rou= 0.0735 # density in Ibm/cu.ft \n",
+ "cp=0.240 # specific heat BTU/(lbm-degree Rankine) \n",
+ "v= 16.88e-5 # viscosity in sq.ft/s \n",
+ "kf = 0.01516 # thermal conductivity in BTU/(hr.ft.degree Rankine) \n",
+ "a = 0.859 # diffusivity in sq.ft/hr \n",
+ "Pr = 0.708 # Prandtl Number\n",
+ "Tw=90\n",
+ "T_inf=70\n",
+ "g=32.2\n",
+ "L=5.5 # length in ft\n",
+ "W=2+(4/12.0) # width in ft\n",
+ "Beta=1/(Tw+460.0) # volumetric thermal expansion coefficient in per degree Rankine\n",
+ "Ra=(g*Beta*(Tw-T_inf)*L**3)/(v*a/3600)\n",
+ "hc=(kf/L)*(0.825+((0.387*(Ra)**(1/6.0))/(1+(0.492/Pr)**(9/16.0))**(8/27.0)))**2\n",
+ "q=hc*L*W*(Tw-T_inf)\n",
+ "\n",
+ "print\"The heat gained is %d BTU/hr\",round(q,0),\"BTU/hr\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat gained is %d BTU/hr 142.0 BTU/hr\n"
+ ]
+ }
+ ],
+ "prompt_number": 30
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 8.4 Page no. 421"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "rou= 1123 # density in kg/m^3 \n",
+ "cp= 1006.7 # specific heat in J/(kg*K) \n",
+ "v= 17.204e-6 # vismath.cosity in m^2/s \n",
+ "Pr =0.703 # Prandtl Number \n",
+ "kf= 0.02738 # thermal conductivity in W/(m.K)\n",
+ "a = 0.2446e-4 # diffusivity in m^2/s \n",
+ "g=9.81\n",
+ "L=5.0\n",
+ "theta=45\n",
+ "T_inf=20.0 # ambient air temperature in degree C\n",
+ "Tw=65 # roof surface temperature in degree C\n",
+ "Beta=1/(T_inf+273.0) # volumetric thermal math.expansion coefficient in per K\n",
+ "\n",
+ "import math\n",
+ "x=((3e5*math.exp(0.1368*math.cos(90-theta))*v*a)/(g*math.cos(theta)*Beta*(Tw-T_inf)))**(1/3.0)\n",
+ "x=0.051\n",
+ "print\"The Laminar-turbulent transition length by Vliet equation is \",round(x,3),\"m\\n\\n\"\n",
+ "lists=[0.02,0.04,0.051,0.051,0.1,1.0,3,5]\n",
+ "Ra=[0,0,0,0,0,0,0,0]\n",
+ "hc=[0,0,0,0,0,0,0,0]\n",
+ "print\"_______________________________________________________\"\n",
+ "print\"x(m)\\t\\tRaL\\t\\t\\thc(W/[m.K])\"\n",
+ "print\"_______________________________________________________\"\n",
+ "for i in range(0,8):\n",
+ " if lists[i]<x:\n",
+ " # Laminar Flow regime exists\n",
+ " Ra[i]=(g*math.cos(math.pi*45.0/180.0)*Beta*(Tw-T_inf)*lists[i]**3)/(v*a)\n",
+ " hc[i]=(kf/lists[i])*(0.68+(0.670*Ra[i]**(1/4.0))/(1+(0.492/Pr)**(9/16.0))**(4.0/9.0))\n",
+ " print lists[i],\"\\t\\t%.4g\"%Ra[i],\"\\t\\t%.3g\"%hc[i]\n",
+ " else:\n",
+ " # Turbulent Flow regime exists\n",
+ " Ra[i]=(g*Beta*(Tw-T_inf)*lists[i]**3)/(v*a)\n",
+ " hc[i]=(0.02738/lists[i])*(0.825+0.324*Ra[i]**(1/6.0))**2\n",
+ " print lists[i],\"\\t\\t%.4g\"%Ra[i],\"\\t\\t%.3g\"%hc[i]\n",
+ " \n",
+ "print\"\\n\\nNOTE:\\nCalculation mistake in book in calculation of Ral and hc,when x=0.04(2nd step in loop)\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The Laminar-turbulent transition length by Vliet equation is 0.051 m\n",
+ "\n",
+ "\n",
+ "_______________________________________________________\n",
+ "x(m)\t\tRaL\t\t\thc(W/[m.K])\n",
+ "_______________________________________________________\n",
+ "0.02 \t\t2.025e+04 \t\t9.32\n",
+ "0.04 \t\t1.62e+05 \t\t7.52\n",
+ "0.051 \t\t4.749e+05 \t\t7.3\n",
+ "0.051 \t\t4.749e+05 \t\t7.3\n",
+ "0.1 \t\t3.58e+06 \t\t6.39\n",
+ "1.0 \t\t3.58e+09 \t\t4.99\n",
+ "3 \t\t9.667e+10 \t\t4.73\n",
+ "5 \t\t4.475e+11 \t\t4.66\n",
+ "\n",
+ "\n",
+ "NOTE:\n",
+ "Calculation mistake in book in calculation of Ral and hc,when x=0.04(2nd step in loop)\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 8.5 Page No. 424"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "rou= 0.0735 # density in lbm/cu.ft\n",
+ "cp=0.240 # specific heat BTU/(lbm-degree Rankine) \n",
+ "v= 16.88e-5 # viscosity in sq.ft/s \n",
+ "kf = 0.01516 # thermal conductivity in BTU/(hr.ft.degree Rankine) \n",
+ "a = 0.859 # diffusivity in sq.ft/hr \n",
+ "Pr = 0.708 # Prandtl Number\n",
+ "Tw=100.0 # temperature of outside surface temperature of oven in degree F\n",
+ "T_inf=60.0 # ambient temperature in degree F\n",
+ "g=32.2\n",
+ "L=2.0 # length in ft\n",
+ "W=2.0 # width in ft\n",
+ "\n",
+ "Beta=1/(T_inf+460.0) # volumetric thermal expansion coefficient in per degree Rankine\n",
+ "Ra=(g*Beta*(Tw-T_inf)*L**3)/(v*a/3600.0)\n",
+ "hc=(kf/L)*(0.68+(0.670*Ra**(0.25))/(1+(0.492/Pr)**(9/16.0))**(4/9.0))\n",
+ "q1side=hc*L*W*(Tw-T_inf)\n",
+ "Lc=0.5\n",
+ "Ra_L=(g*Beta*(Tw-T_inf)*Lc**3)/(v*a/3600.0) # Rayleigh number based on characteristic length\n",
+ "hc_L=(kf/Lc)*0.54*(Ra_L)**(1/4.0)\n",
+ "qtop=hc_L*L*W*(Tw-T_inf)\n",
+ "\n",
+ "print\"The heat transferred from one side is \",round(q1side,1),\"BTU/hr\"\n",
+ "print\"The heat transferred from top is \",round(qtop,0),\"BTU/hr\"\n",
+ "\n",
+ "if qtop < q1side:\n",
+ " \n",
+ " print\"More heat is transfered from side\" \n",
+ "else:\n",
+ " print \"More heat is transfered top side\"\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat transferred from one side is 93.7 BTU/hr\n",
+ "The heat transferred from top is 138.0 BTU/hr\n",
+ "More heat is transfered top side\n"
+ ]
+ }
+ ],
+ "prompt_number": 17
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 8.6 Page No.427"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "rou= 1.177 # density in kg/cu.m\n",
+ "cp= 1005.7 # specific heat in J/(kg*K) \n",
+ "v= 15.68e-6 # viscosity in sq.m/s \n",
+ "Pr =0.708 # Prandtl Number \n",
+ "kf=0.02624 # thermal conductivity in W/(m.K)\n",
+ "a=0.22160e-4 # diffusivity in sq.m/s \n",
+ "g=9.81\n",
+ "L=4.0 # length in m\n",
+ "D=15/100.0 # diameter in m\n",
+ "T_inf=5.0 # ambient air temperature in degree C\n",
+ "Tw=50.0 # outside surface temperature in degree C\n",
+ "\n",
+ "import math\n",
+ "Beta=1/(T_inf+273.0) # volumetric thermal expansion coefficient in per K\n",
+ "Ra=(g*Beta*(Tw-T_inf)*D**3)/(v*a)\n",
+ "hc_h=(kf/D)*(0.60+(0.387*Ra**(1/6.0))/(1+(0.559/Pr)**(9/16.0))**(8/27.0))**2\n",
+ "As=math.pi*D*L\n",
+ "q_hor=hc_h*As*(Tw-T_inf)\n",
+ "hc_v=(kf/D)*0.6*(Ra*(D/L))**(1/4.0)\n",
+ "q_ver=hc_v*As*(Tw-T_inf)\n",
+ "q=round(q_ver,0)+round(q_hor,0)\n",
+ "\n",
+ "print\"The heat transferred from the horizontal length of 4 m is \",round(q_hor,0),\"W\"\n",
+ "print\"The heat transferred from the vertical length of 4 m is \",round(q_ver,0),\"W\"\n",
+ "print\"nThe total heat lost from the pipe is \",round(q,2),\"W\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat transferred from the horizontal length of 4 m is 477.0 W\n",
+ "The heat transferred from the vertical length of 4 m is 246.0 W\n",
+ "nThe total heat lost from the pipe is 723.0 W\n"
+ ]
+ }
+ ],
+ "prompt_number": 52
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 8.7 Page No. 430"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "rou= 0.0809 # density in lbm/cu.ft \n",
+ "cp=0.240 # specific heat BTU/(lbm-degree Rankine) \n",
+ "v= 13.54e-5 # viscosity in sq.ft/s \n",
+ "kf = 0.01402 # thermal conductivity in BTU/(hr.ft.degree Rankine) \n",
+ "a = 0.685 # diffusivity in sq.ft/hr \n",
+ "Pr = 0.712 # Prandtl Number\n",
+ "Tw=0 # temperature of outside surface temperature of oven in degree F\n",
+ "T_inf=70.0 # ambient temperature in degree F\n",
+ "g=32.2\n",
+ "Beta=1/(T_inf+460.0) # volumetric thermal expansion coefficient in per degree Rankine\n",
+ "Lc=1/((1/1)+(1/1.2))\n",
+ "Ra=(g*Beta*abs(Tw-T_inf)*Lc**3)/(v*a/3600.0)\n",
+ "hc=(kf/Lc)*0.6*(Ra)**(1/4.0)\n",
+ "\n",
+ "print\"The value of convection coefficient is \",round(hc,2),\"BTU/(hr.sq.ft.degree R)\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The value of convection coefficient is 1.11 BTU/(hr.sq.ft.degree R)\n"
+ ]
+ }
+ ],
+ "prompt_number": 54
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 8.8 Page No.433"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "\n",
+ "rou= 0.998 # density in kg/cu.m\n",
+ "cp= 1009.0 # specific heat in J/(kg*K) \n",
+ "v= 20.76e-6 # viscosity in sq.m/s \n",
+ "Pr =0.697 # Prandtl Number \n",
+ "kf= 0.03003 # thermal conductivity in W/(m.K)\n",
+ "a = 0.2983e-4 # diffusivity in sq.m/s \n",
+ "g=9.81\n",
+ "T_inf=35 # ambient air temperature in degree C\n",
+ "Tw=100 # surface temperature in degree C\n",
+ "Beta=1/(T_inf+273.0) # volumetric thermal expansion coefficient in per K\n",
+ "rou_Al=2702 # density in kg/cu.m\n",
+ "k_Al=236 # thermal conductivity in W/(m.K)\n",
+ "cp_Al=896 # specific heat in J/(kg*K) \n",
+ "a_Al=97.5e-6 # diffusivity in sq.m/s \n",
+ "b=46/100.0\n",
+ "w=24/100.0\n",
+ "\n",
+ "zeta=((w*v**2)/(g*Beta*(Tw-T_inf)*Pr))**(1/4.0)\n",
+ "L=1.54*(k_Al/kf)**(1/2)*zeta\n",
+ "S=2.89*zeta\n",
+ "q=(b*w*(Tw-T_inf)*1.3*(k_Al*kf)**(1/2.0))/(6*zeta)\n",
+ "N=b/(2*S)\n",
+ "\n",
+ "print\"The heat transfer rate is \",round(q,0),\"W\"\n",
+ "print\"The number of fins can be atmost\",round(N,0)\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat transfer rate is 1423.0 W\n",
+ "The number of fins can be atmost 27.0\n"
+ ]
+ }
+ ],
+ "prompt_number": 59
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [],
+ "language": "python",
+ "metadata": {},
+ "outputs": []
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |