From 92cca121f959c6616e3da431c1e2d23c4fa5e886 Mon Sep 17 00:00:00 2001 From: hardythe1 Date: Tue, 7 Apr 2015 15:58:05 +0530 Subject: added books --- .../Chapter_4.ipynb | 466 +++++++++++++++++++++ 1 file changed, 466 insertions(+) create mode 100755 Thermodynamics:_From_concepts_to_applications/Chapter_4.ipynb (limited to 'Thermodynamics:_From_concepts_to_applications/Chapter_4.ipynb') diff --git a/Thermodynamics:_From_concepts_to_applications/Chapter_4.ipynb b/Thermodynamics:_From_concepts_to_applications/Chapter_4.ipynb new file mode 100755 index 00000000..b4cbb2bf --- /dev/null +++ b/Thermodynamics:_From_concepts_to_applications/Chapter_4.ipynb @@ -0,0 +1,466 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:f794a4d18f36c13ddab8d6024721da6bc9aca8b98b0260eaf3e1fb73a8bd345d" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Chapter 4-simple systems" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example1-pg43" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#caluclate the heat interaction and temperature rise\n", + "##initialisation of variables\n", + "m= 5000. ##kg\n", + "cp= 1.4 ##kJ/kg K\n", + "T2= 27.6 ##K\n", + "T1= 22. ##K\n", + "t= 40. ##min\n", + "P= 20. ##kW\n", + "##CALCULATIONS\n", + "H= m*cp*(T2-T1)\n", + "W= -P*t*60\n", + "Q= H+W\n", + "dT= -W/(m*cp)\n", + "##RESULTS\n", + "print'%s %.f %s'%('heat interaction =',Q,'kJ')\n", + "print'%s %.2f %s'%('Temperature rise = ',dT,'C')\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "heat interaction = -8800 kJ\n", + "Temperature rise = 6.86 C\n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 2-pg 49\n" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#caluclate the volumes and internal energy and enthaply\n", + "##initialisation of variables\n", + "T= 300. ##C\n", + "p= 2. ##Mpa\n", + "T1= 300. ##C\n", + "p1= 20.##Mpa\n", + "T2= 300. ##C\n", + "p2= 8.501 ##Mpa\n", + "##CALCULATIONS\n", + "v= 0.12547\n", + "v1= 0.00136\n", + "u= 2772.6\n", + "u1= 1306.1\n", + "h= 3023.5\n", + "h1= 1333.3\n", + "##RESULTS\n", + "print'%s %.5f %s'%('volume =',v,'m^3/kg')\n", + "print'%s %.5f %s'%('volume =',v1,'m^3/kg')\n", + "print'%s %.1f %s'%('internal energy = ',u,'kJ/kg')\n", + "print'%s %.1f %s'%('internal energy =',u1,'kJ/kg')\n", + "print'%s %.1f %s'%('enthalpy =',h,'kJ/kg')\n", + "print'%s %.f %s'%('enthalpy = ',h1,'kJ/kg')\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "volume = 0.12547 m^3/kg\n", + "volume = 0.00136 m^3/kg\n", + "internal energy = 2772.6 kJ/kg\n", + "internal energy = 1306.1 kJ/kg\n", + "enthalpy = 3023.5 kJ/kg\n", + "enthalpy = 1333 kJ/kg\n" + ] + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example3-pg49\n" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#calculate the volume and internal energy and enthalpy\n", + "##initialisation of variables\n", + "vf= 0.001404 ##m^3/kg\n", + "x= 0.8\n", + "vg= 0.02167 ##m^3/kg\n", + "uf= 1332. ##kJ/kg\n", + "ug= 1231. ##kJ/kg\n", + "hf= 1344. ##kJ/kg\n", + "hg= 1404.9 ##kJ/kg\n", + "##CALCULATIONS\n", + "v= vf+x*(vg-vf)\n", + "u= uf+x*ug\n", + "h= hf+x*hg\n", + "##RESULTS\n", + "print'%s %.5f %s'%('volume =',v,'m^3/kg')\n", + "print'%s %.1f %s'%('internal energy =',u,'kJ/kg')\n", + "print'%s %.1f %s'%('enthalpy =',h, 'kJ/kg')\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "volume = 0.01762 m^3/kg\n", + "internal energy = 2316.8 kJ/kg\n", + "enthalpy = 2467.9 kJ/kg\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example4-pg 52" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#calculate the specific volume and internal energy and enthalpy and entropy\n", + "##initialisation of variables\n", + "T= 296. ##K\n", + "T1= 250. ##K\n", + "T2= 300. ##K\n", + "v= 0.1257 ##m^3/kg\n", + "v1= 0.11144 ##m^3/kg\n", + "u1= 27772.6 ##kJ/kg\n", + "u2= 2679.6 ##kJ/kg\n", + "h1= 3023.5 ##kJ/kg\n", + "h2= 2902.5 ##kJ/kg\n", + "s1= 6.7664 ##kJ/kg K\n", + "s2= 6.5433 ##kJ/kg K\n", + "##CALCULATIONS\n", + "a1= (T-T1)/(T2-T1)\n", + "a2= 1-a1\n", + "V= a1*v+a2*v1\n", + "U= a1*u1+a2*u2\n", + "H= a1*h1+a2*h2\n", + "S= a1*s1+a2*s2\n", + "##RESULTS\n", + "print'%s %.3f %s'%('a2 = ',a2,'')\n", + "print'%s %.5f %s'%('specific volume =',V,'m^3/kg')\n", + "print'%s %.1f %s'%('internal energy =',U,'kJ/kg')\n", + "print'%s %.1f %s'%('enthalpy =',H,'kJ/kg')\n", + "print'%s %.1f %s'%('Entropy =',S,'kJ/kg')\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "a2 = 0.080 \n", + "specific volume = 0.12456 m^3/kg\n", + "internal energy = 25765.2 kJ/kg\n", + "enthalpy = 3013.8 kJ/kg\n", + "Entropy = 6.7 kJ/kg\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Exapmle5-pg54" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#calculate a2 and internal enerjy\n", + "##initialisation of variables\n", + "v= 0.15 ##m^3/kg\n", + "v1= 0.13857 ##m^3/kg\n", + "v2= 0.1512 ##m^3/kg\n", + "v3= 0.050 ##m^3##kg\n", + "vf= 0.001177 ##m^3/kg\n", + "vg= 0.09963 ##m^3/kg\n", + "uf= 906.44 ##kJ/kg\n", + "ufg= 1693.8 ##kJ/kg\n", + "##CALCULATIONS\n", + "a1= (v-v1)/(v2-v1)\n", + "a2= 1-a1\n", + "x= (v3-vf)/(vg-vf)\n", + "u= uf+x*ufg\n", + "##RESULTS\n", + "print'%s %.3f %s'%('a2 =',a2,'')\n", + "print'%s %.1f %s'%('internal energy =',u,'kJ/kg')\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "a2 = 0.095 \n", + "internal energy = 1746.4 kJ/kg\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example6-pg52" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#calculate pressure and temperature\n", + "##initialisation of variables\n", + "T= 250.##C\n", + "T1= 300. ##C\n", + "v300= 0.6548 ##m^3/kg\n", + "v250= 0.591 ##m^3/kg\n", + "v= 0.6 ##m^3/kg\n", + "u= 3000. ##kJ/kg\n", + "u250= 2726.1 ##kJ/kg\n", + "u300= 2804.8 ##kJ/kg\n", + "T2= 510.30 \n", + "u2= 3145.26 ##kJ/kg\n", + "p= 0.4 ##Mpa\n", + "p2= 0.2 ##Mpa\n", + "##CALCULATIONS\n", + "T1= T1+((v-v250)/(v300-v250))*(T2-T)\n", + "u1= u250+((v-v250)/(v300-v250))*(u300-u250)\n", + "du= u1-u\n", + "p1= p+((u-u1)/(u2-u1))*p2\n", + "##RESULTS\n", + "print'%s %.3f %s'%('pressure =',p1,'Mpa')\n", + "print'%s %.f %s'%('temperature =',T2,'C')\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "pressure = 0.529 Mpa\n", + "temperature = 510 C\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 7-pg54" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#calculate work and heat interaction\n", + "##initialisation of variables\n", + "V= 0.2 ##m^3\n", + "v1= 0.02995 ##m^3/kg\n", + "u2= 2826.7 ##kJ/kg\n", + "u1= 2747.7 ##kJ/kg\n", + "h2= 3092.5 ##kJ/kg\n", + "h1= 2987.3 ##kJ/kg\n", + "p= 4 ##Mpa\n", + "v2= 0.06645 ##m^3/kg\n", + "v1= 0.02995 ##m^3/kg\n", + "##CALCULATIONS\n", + "m= V/v1\n", + "U= m*(u2-u1)\n", + "H= m*(h2-h1)\n", + "W= m*p*10*10*10*(v2-v1)\n", + "Q= U+W\n", + "##RESULTS\n", + "print'%s %.1f %s'%('work =',W,'kJ') \n", + "print'%s %.1f %s'%('heat interaction =',Q,'kJ')\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "work = 975.0 kJ\n", + "heat interaction = 1502.5 kJ\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example8-pg 55" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#calculate work and heat interaction\n", + "##initialisation of variables\n", + "m= 6.678 ##kg\n", + "u2= 2826.7 ##kJ/kg\n", + "u1= 2747.7 ##kJ/kg\n", + "p1= 8. ##Mpa\n", + "p2= 7. ##Mpa\n", + "p3= 6. ##Mpa\n", + "p4= 5. ##Mpa\n", + "p5= 4. ##Mpa\n", + "v1= 29.95 ##L/kg\n", + "v2= 35.24 ##L/kg\n", + "v3= 42.23 ##L/kg\n", + "v4= 51.94 ##L/kg\n", + "v5= 66.45 ##L/kg\n", + "##CALCULATIONS\n", + "U= m*(u2-u1)\n", + "W= m*0.5*((p1+p2)*(v2-v1)+(p2+p3)*(v3-v2)+(p3+p4)*(v4-v3)+(p4+p5)*(v5-v4))\n", + "Q=U+W\n", + "##RESULTS\n", + "print'%s %.1f %s'%('work =',W,'kJ') \n", + "print'%s %.f %s'%('heat interaction =',Q,'kJ') \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "work = 1361.0 kJ\n", + "heat interaction = 1889 kJ\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 9-pg56\n" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#calculate final pressure and enthalpy and etropy\n", + "##initialisation of variables\n", + "p0= 100. ##kpa\n", + "A= 0.1 ##m^2\n", + "F= 20. ##kN\n", + "m3= 0.8873 ##kg\n", + "m1= 1.1384 ##kg\n", + "m2= 0.2511 ##kg\n", + "u1= 3116.2 ##kJ/kg\n", + "u2= 2728.7 ##kJ/kg\n", + "v3= 0.9942 ##m^3/kg\n", + "##CALCULATIONS\n", + "pe= (p0+(F/A))/1000\n", + "h3= (m1*u1-m2*u2)/m3\n", + "z3= m3*v3/A\n", + "##RESULTS\n", + "print'%s %.1f %s'%('final pressure =',pe,'Mpa') \n", + "print'%s %.1f %s'%('enthalpy = ',h3,'kJ/kg') \n", + "print'%s %.2f %s'%('piston rise =',z3,'m') \n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "final pressure = 0.3 Mpa\n", + "enthalpy = 3225.9 kJ/kg\n", + "piston rise = 8.82 m\n" + ] + } + ], + "prompt_number": 14 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file -- cgit