diff options
| author | hardythe1 | 2015-04-07 15:58:05 +0530 |
|---|---|---|
| committer | hardythe1 | 2015-04-07 15:58:05 +0530 |
| commit | 92cca121f959c6616e3da431c1e2d23c4fa5e886 (patch) | |
| tree | 205e68d0ce598ac5caca7de839a2934d746cce86 /Thermodynamics_An_Engineering_Approach/Chapter5.ipynb | |
| parent | b14c13fcc6bb6d01c468805d612acb353ec168ac (diff) | |
| download | Python-Textbook-Companions-92cca121f959c6616e3da431c1e2d23c4fa5e886.tar.gz Python-Textbook-Companions-92cca121f959c6616e3da431c1e2d23c4fa5e886.tar.bz2 Python-Textbook-Companions-92cca121f959c6616e3da431c1e2d23c4fa5e886.zip | |
added books
Diffstat (limited to 'Thermodynamics_An_Engineering_Approach/Chapter5.ipynb')
| -rwxr-xr-x | Thermodynamics_An_Engineering_Approach/Chapter5.ipynb | 614 |
1 files changed, 614 insertions, 0 deletions
diff --git a/Thermodynamics_An_Engineering_Approach/Chapter5.ipynb b/Thermodynamics_An_Engineering_Approach/Chapter5.ipynb new file mode 100755 index 00000000..da9be875 --- /dev/null +++ b/Thermodynamics_An_Engineering_Approach/Chapter5.ipynb @@ -0,0 +1,614 @@ +{
+ "metadata": {
+ "name": ""
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter 5: Mass and Energy Analysis of Control Volumes"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-1 ,Page No.224"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#given data\n",
+ "V=10.0;#volumne of bucket in galon\n",
+ "t=50;#time taken to fill the bucket in sec\n",
+ "p=1;#density of water in kg/L\n",
+ "re=0.8/2/100;#radius of nozzle exit in m\n",
+ "\n",
+ "#calculations\n",
+ "Vd=V/t*3.7854;#factor 0f 3.7854 for gal to L\n",
+ "print'volumne flow rate through hose %f L/s'%round(Vd,3);\n",
+ "m=p*Vd;\n",
+ "print'mass flow rate through hose %f kg/s'%round(m,3);\n",
+ "Ae=math.pi*re**2;\n",
+ "Ve=Vd/Ae/1000;#factor of 1000 for L to m^3\n",
+ "print'average velocity at the nozzle %f m/s'%round(Ve,1);\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "volumne flow rate through hose 0.757000 L/s\n",
+ "mass flow rate through hose 0.757000 kg/s\n",
+ "average velocity at the nozzle 15.100000 m/s\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-2 ,Page No.225"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from scipy.integrate import quad \n",
+ "from pylab import *\n",
+ "\n",
+ "#given data\n",
+ "Dtank=3*12;#diameter of tank in inches\n",
+ "Djet=0.5;#diameter of water jet in inches\n",
+ "h0=2;#bottom reference in ft\n",
+ "h1=4;#height of water tank in ft\n",
+ "\n",
+ "#constants used \n",
+ "g=32.2;#in ft/s^2\n",
+ "\n",
+ "#calculations\n",
+ "#min - mout = dmCV/dt\n",
+ "#mout = p*(2*g*h*Ajet)^2\n",
+ "#mCV = p*Atank*h\n",
+ "#from these we get dt = Dtank^2/Djet^2 * (dh/(2*g*h)^2)\n",
+ "def intgrnd1(h): \n",
+ " return (Dtank**2/Djet**2*(1/sqrt(2*g*h)))\n",
+ "t, err = quad(intgrnd1, h0, h1) \n",
+ "t=(t/60);#in min\n",
+ "print'time taken to drop to 2 ft %f min'%round(t,1)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "time taken to drop to 2 ft 12.600000 min\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-3 ,Page No.229"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "P=150;#operating pressure in kPa\n",
+ "Vliquid=0.6/1000.0;#amount of liquid in the cooker in m^3\n",
+ "t=40.0*60;#period of operation in sec\n",
+ "Ac=8*10**-6;#exit area of opening in m^2\n",
+ "\n",
+ "#from Table A-5\n",
+ "#from P = 150 kPa\n",
+ "h=2693.1;\n",
+ "ug=2519.2;\n",
+ "vf=0.001053;\n",
+ "vg=1.1594;\n",
+ "\n",
+ "#calculations\n",
+ "m=Vliquid/vf;\n",
+ "md=m/t;\n",
+ "print'mass flow rate %f kg/s'%md;\n",
+ "V=md*vg/(Ac);\n",
+ "print'exit velocity %f m/s'%round(V,1);\n",
+ "Eflow=h-ug;\n",
+ "Et=h;\n",
+ "print'flow energy %f kJ/kg'%round(Eflow,1);\n",
+ "print'total energy %f kJ/kg'%round(Et,1);\n",
+ "Emass=md*Et;\n",
+ "print'rate at which energy leaves the cooker %f kW'%round(Emass,3)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "mass flow rate 0.000237 kg/s\n",
+ "exit velocity 34.400000 m/s\n",
+ "flow energy 173.900000 kJ/kg\n",
+ "total energy 2693.100000 kJ/kg\n",
+ "rate at which energy leaves the cooker 0.639000 kW\n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-4 ,Page No.234"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "T1=283;#temperature of air in K\n",
+ "P1=80;#entering pressure in kPa\n",
+ "V1=200;#velocity ar inlet in m/s\n",
+ "A1=0.4;#inlet area in m^2\n",
+ "\n",
+ "#constants used\n",
+ "R=0.287;#in kPa-m^3/kg-K\n",
+ "\n",
+ "#calulations\n",
+ "v1=R*T1/P1;\n",
+ "m=V1*A1/v1;\n",
+ "print'mass flow rate of air %f kg/s'%round(m,1);\n",
+ "# Ein - Eout = dEsystem / dt\n",
+ "#from Table A-17\n",
+ "h1=283.14;\n",
+ "V2=0;\n",
+ "h2=h1-(V2^2 - V1^2)/2/1000;#factor of 1000 to convert to kJ/kg\n",
+ "#from Table A-17 at this value of h2\n",
+ "T2=303;\n",
+ "print'the temperature %i K is'%T2;\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "mass flow rate of air 78.800000 kg/s\n",
+ "the temperature 303 K is\n"
+ ]
+ }
+ ],
+ "prompt_number": 15
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-5 ,Page No.235"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "P1=250;#inlet pressure in psia\n",
+ "T1=700;#inlet temp in F\n",
+ "A1=0.2;#inlet area in ft^2\n",
+ "qout=1.2;#heat losses in Btu/lbm\n",
+ "m=10;#mass flow rate in lbm/sec\n",
+ "P2=200;#nozzle pressure in kPa\n",
+ "V2=900;#nozzle velocity in m/s\n",
+ "\n",
+ "#from Table A-6E\n",
+ "v1=2.6883;\n",
+ "h1=1371.4;\n",
+ "\n",
+ "#calculations\n",
+ "V1=m*v1/A1;\n",
+ "print'the inlet velocit %f f/s'%round(V1,1);\n",
+ "# Ein - Eout = dEsystem / dt\n",
+ "h2=h1-qout-(V2**2 - V1**2)/2/25037;#factor of 25037 to convert to Btu/lbm\n",
+ "#at this value h2, from Tablw A-6E\n",
+ "T2=662;\n",
+ "print'exit temperature %i F'%T2\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the inlet velocit 134.400000 f/s\n",
+ "exit temperature 662 F\n"
+ ]
+ }
+ ],
+ "prompt_number": 19
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-6 ,Page No.236"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "T1=280;#intial pressure in kPa\n",
+ "P1=100;#intial temp in K\n",
+ "m=0.02;#mass flow rate in kg/s\n",
+ "qout=16;#heat losses in kJ/kg\n",
+ "P2=600;#final pressure in kPa\n",
+ "T2=400;#final temp in K\n",
+ "\n",
+ "#from Table A-17\n",
+ "h1=280.13;\n",
+ "h2=400.98;\n",
+ "\n",
+ "#calculations\n",
+ "# Ein - Eout = dEsystem / dt\n",
+ "Win=m*qout+m*(h2-h1);\n",
+ "print'the input power of compressor %f kW'%round(Win,2)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the input power of compressor 2.740000 kW\n"
+ ]
+ }
+ ],
+ "prompt_number": 20
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-7 ,Page No.237"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "P1=2;#inlet pressure in MPa\n",
+ "T1=400;#inlet temp in K\n",
+ "V1=50.0;#inlet velocity in m/s\n",
+ "z1=10;#inlet height in m\n",
+ "P2=15;#final pressure in MPa\n",
+ "x2=0.9;#final dryness fraction\n",
+ "V2=180.0;#final velocity in m/s\n",
+ "z2=6;#final height in m\n",
+ "Wout=5*1000;#power output in kW\n",
+ "\n",
+ "#from Table A-6\n",
+ "h1=3248.4;\n",
+ "#similarly for P2\n",
+ "hf=225.94;\n",
+ "hfg=2372.3;\n",
+ "\n",
+ "#constants used \n",
+ "g=9.8;#in m/s^2\n",
+ "\n",
+ "#calcualtions\n",
+ "h2=hf+x2*hfg;\n",
+ "print'difference in enthalpies %f kJ/kg'%round((h2-h1),2);\n",
+ "print'difference in kinetic energy %f kJ/kg'%round((V2**2-V1**2)/2/1000,2);#factor of 1000 to convert to kJ/kg\n",
+ "print'difference in potential energy %f kJ/kg'%round(g*(z2-z1)/1000,2);#factor of 1000 to convert to kJ/kg\n",
+ "wout=-((h2-h1)+(V2**2-V1**2)/2/1000+g*(z2-z1)/1000);#factor of 1000 to convert to kJ/kg\n",
+ "print'work done per unit of mass %f kJ/kg'%round(wout,2);\n",
+ "m=Wout/wout;\n",
+ "print'mass flow rate %f kg/s'%round(m,2)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "difference in enthalpies -887.390000 kJ/kg\n",
+ "difference in kinetic energy 14.950000 kJ/kg\n",
+ "difference in potential energy -0.040000 kJ/kg\n",
+ "work done per unit of mass 872.480000 kJ/kg\n",
+ "mass flow rate 5.730000 kg/s\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-8 ,Page No.239"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "P1=0.8;#entering pressure in MPa\n",
+ "P2=0.12;#throttled pressure in MPa\n",
+ "\n",
+ "#from Table A-12\n",
+ "#sat. liq at P1\n",
+ "T1=31.31;\n",
+ "h1=95.47;\n",
+ "#since process is insentropic and at P2\n",
+ "h2=h1;\n",
+ "hf=22.49;\n",
+ "hg=236.97;\n",
+ "T2=-22.32;\n",
+ "\n",
+ "#calculations\n",
+ "x2=(h2-hf)/(hg-hf);\n",
+ "print'the final state is %f'%round(x2,3);\n",
+ "dT=T2-T1;\n",
+ "print'temperature drop %f C'%round(dT,2)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the final state is 0.340000\n",
+ "temperature drop -53.630000 C\n"
+ ]
+ }
+ ],
+ "prompt_number": 26
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-9 ,Page No.241"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "T1=140;#inlet temp of hot water in F\n",
+ "T2=50;#inlet temp of cold water in F\n",
+ "T3=110;#steady state output temp in F\n",
+ "P=20;#operating pressure in psia\n",
+ "\n",
+ "#for a compressed liq at given temp\n",
+ "h1=107.99;\n",
+ "h2=18.07;\n",
+ "h3=78.02;\n",
+ "\n",
+ "#calculations\n",
+ "#Mass balance min = mout So, m1+m2 = m3\n",
+ "#Energy balance Ein = Eout So, m1*h1 + m2*h2 = m3*h3\n",
+ "#combining realations\n",
+ "#m1*h1 + m2*h2 = (m1+m2)*h3\n",
+ "#dividing by m2 and y=m1/m2\n",
+ "#we get, yh1 + h2 = (y+1)*h3\n",
+ "y=(h3-h2)/(h1-h3);\n",
+ "print'the ratio of mass flow rates %f'%round(y,1)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the ratio of mass flow rates 2.000000\n"
+ ]
+ }
+ ],
+ "prompt_number": 27
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-10 ,Page No.242"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "T1=15;#water inlet temp in C\n",
+ "P1=300;#water inlet pressure in kPa\n",
+ "T2=25;#water outlet temp in C\n",
+ "T3=70;#R-134a inlet temp in C\n",
+ "P3=1000;#R-134a inlet pressure in kPa\n",
+ "T4=35;#R-134a outlet temp in C\n",
+ "mr=6;#mass flow rate in kg/min\n",
+ "\n",
+ "#from Table A-4, A-13 and A-11\n",
+ "h1=62.982;\n",
+ "h2=104.83;\n",
+ "h3=303.85;\n",
+ "h4=100.87;\n",
+ "\n",
+ "#calculations\n",
+ "#mass balance m1=m2=mw and m3=m4=mr\n",
+ "#energy balance m1*h1 + m3*h3 = m2*h2 + m4*h4\n",
+ "#combining them mw*(h1-h2) = mr*(h4-h3)\n",
+ "mw= mr*(h4-h3)/(h1-h2);\n",
+ "print'mass flow rate of cooling water %f kg/min'%round(mw,1);\n",
+ "Qin=mw*(h2-h1);\n",
+ "print'heat transfer rate %i kJ/min'%round(Qin)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "mass flow rate of cooling water 29.100000 kg/min\n",
+ "heat transfer rate 1218 kJ/min\n"
+ ]
+ }
+ ],
+ "prompt_number": 29
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-11 ,Page No.245"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "T1=17+273;#inlet temp in K\n",
+ "P1=100;#inlet pressure in kPa\n",
+ "V1=150;#inlet volumetric rate in m^3/min\n",
+ "Win=15;#rated power in kW\n",
+ "Qout=200/1000;#heat lost in kJ/s\n",
+ "\n",
+ "#constants used\n",
+ "R=0.287;#in kPa-m^3/kg-K\n",
+ "cp=1.005;#in kJ/kg C\n",
+ "\n",
+ "#calculations\n",
+ "v1=R*T1/P1;\n",
+ "m=V1/v1/60;#factor of 6 to convert to s\n",
+ "# Win - Qout = m*cp*(T2-T1)\n",
+ "T2= T1 + (Win - Qout)/(m*cp);\n",
+ "print'exit temperature %f C'%round((T2-273),2)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "exit temperature 21.970000 C\n"
+ ]
+ }
+ ],
+ "prompt_number": 33
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5-13 ,Page No.250"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#given data\n",
+ "V=6.0/1000;#volumne of cooker in m^3\n",
+ "Pgage=75;#operating pressure in kPa\n",
+ "Patm=100;#atmospheric pressure in kPa\n",
+ "m1=1;#mass of water in kg\n",
+ "Qind=0.5;#heat supplying rate in kJ/sec ; d stands for .\n",
+ "t=30*60;#operating time in s\n",
+ "\n",
+ "#calculation\n",
+ "Pabs=Pgage+Patm;\n",
+ "#from Table A-5, ths saturation temp \n",
+ "T=116.04;\n",
+ "print'the temperature at which cooking takes place %f C'%T;\n",
+ "#mass balance me=(m1-m2)cv\n",
+ "#energy balance Qin - mehe = (m2u2 - m1u1)cv\n",
+ "Qin=Qind*t;\n",
+ "#from Table A-5\n",
+ "he=2700.2;\n",
+ "vf=0.001;\n",
+ "vg=1.004;\n",
+ "uf=486.82;\n",
+ "ufg=2037.7;\n",
+ "v1=V/m1;\n",
+ "x1=(v1-vf)/(vg-vf);\n",
+ "u1=uf+x1*ufg;\n",
+ "U=m1*u1;\n",
+ "#Qin = (m1 - V/v2)*he + (V/v2*u2 - m1*u1)\n",
+ "#v2=vf + x2*(vg-vf)\n",
+ "#u2=uf + x2*ufg\n",
+ "#combining these equations we get\n",
+ "#solved using EES\n",
+ "x2=0.009;\n",
+ "v2=vf + x2*(vg-vf);\n",
+ "m2=V/v2;\n",
+ "print'amount of water left %f kg'%round(m2,1)\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the temperature at which cooking takes place 116.040000 C\n",
+ "amount of water left 0.600000 kg\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |