diff options
Diffstat (limited to 'Thermodynamics:_From_concepts_to_applications/Chapter14_2.ipynb')
| -rwxr-xr-x | Thermodynamics:_From_concepts_to_applications/Chapter14_2.ipynb | 313 |
1 files changed, 313 insertions, 0 deletions
diff --git a/Thermodynamics:_From_concepts_to_applications/Chapter14_2.ipynb b/Thermodynamics:_From_concepts_to_applications/Chapter14_2.ipynb new file mode 100755 index 00000000..5bcffb75 --- /dev/null +++ b/Thermodynamics:_From_concepts_to_applications/Chapter14_2.ipynb @@ -0,0 +1,313 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:dceed5d123e7adb1d73db762a18a52615a5b487342c94704535d968cf751d6c0"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "chapter14-equations of state and generlized charts"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example2-pg"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "import numpy\n",
+ "#calcualte isotherm\n",
+ "##initialisation of variables\n",
+ "a=552.6 ##kPa m^6/kmol^2\n",
+ "b= 0.03402 ##m^3/kmol\n",
+ "p= 100. ##kPa\n",
+ "R= 8.314 ##J/mol K\n",
+ "##CALCULATIONS\n",
+ "vec=([p,-a,2*a*b])\n",
+ "vector= numpy.roots(vec)\n",
+ "x=vec[0]\n",
+ "T= 2*a*(x-b)*(x-b)/(R*x*x*x)\n",
+ "##RESULTS\n",
+ "print'%s %.1f %s'% (' isotherm=',T,'K')\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ " isotherm= 1.3 K\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example1-pg307"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "#calculate volumes of given temperatures\n",
+ "##initialisation of variables\n",
+ "R= 8.314 ##J/mol K\n",
+ "T= 400. ##C\n",
+ "T1= 500. ##C\n",
+ "M= 18.015 ##kg/k mol\n",
+ "p1= 30. ##Mpa\n",
+ "##CALCULATIONS\n",
+ "v1= R*(273.15+T)/(M*p1*10*10*10)\n",
+ "v2= R*(273.15+T1)/(M*p1*10*10*10)\n",
+ "##RESULTS\n",
+ "print'%s %.5f %s'% ('volume=',v1,'m^3/kg')\n",
+ "print'%s %.5f %s'% ('volume=',v2,'m^3/kg')\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "volume= 0.01036 m^3/kg\n",
+ "volume= 0.01189 m^3/kg\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example2-pg310"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "#calculate specific enthalpy at each and every point their differences\n",
+ "##initialisation of variables\n",
+ "h1= 3892.2 ##kJ/kg\n",
+ "h2= 4102.2 ##kJ/kg\n",
+ "dh= 1015.4 ##kJ/kg\n",
+ "dh1= 448. ##kJ/kg\n",
+ "h3= 2151.1 ##kJ/kg\n",
+ "h4= 3081.1 ##kJ/kg\n",
+ "##RESULTS\n",
+ "print'%s %.1f %s'% ('Specific Enthalpy=',h1,'kJ/kg')\n",
+ "print'%s %.1f %s'% ('Specific Enthalpy=',h2,'kJ/kg')\n",
+ "print'%s %.1f %s'% ('Specific Enthalpy=',h3,'kJ/kg')\n",
+ "print'%s %.1f %s'% ('Specific Enthalpy=',h4,'kJ/kg')\n",
+ "print'%s %.f %s'% ('Enthalpy difference=',dh,'kJ/kg')\n",
+ "print'%s %.f %s'% ('Enthalpy difference=',dh1,'kJ/kg')\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Specific Enthalpy= 3892.2 kJ/kg\n",
+ "Specific Enthalpy= 4102.2 kJ/kg\n",
+ "Specific Enthalpy= 2151.1 kJ/kg\n",
+ "Specific Enthalpy= 3081.1 kJ/kg\n",
+ "Enthalpy difference= 1015 kJ/kg\n",
+ "Enthalpy difference= 448 kJ/kg\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example3-pg313"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "#calculate entropy at every point\n",
+ "##initialisation of variables\n",
+ "s2= 5.7905 ##kJ/kg K\n",
+ "s1= 4.4728 ##kJ/kg K\n",
+ "s3= 4.64437 ##kJ/kg K\n",
+ "s4= 5.7883 ##kJ/kg K\n",
+ "s5= 6.2036 ##kJ/kg K\n",
+ "s6= 5.9128 ##kJ/kg K\n",
+ "##CALCULATIONS\n",
+ "S1= s2-s1\n",
+ "S2= s4-s3\n",
+ "S3= s5-s6\n",
+ "##RESULTS\n",
+ "print'%s %.4f %s'% ('Entropy=',S1,'kJ/kg K')\n",
+ "print'%s %.4f %s'% ('Entropy=',S2,'kJ/kg K')\n",
+ "print'%s %.4f %s'% ('Entropy=',S3,'kJ/kg K')\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Entropy= 1.3177 kJ/kg K\n",
+ "Entropy= 1.1439 kJ/kg K\n",
+ "Entropy= 0.2908 kJ/kg K\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example4-pg315"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "#calculate entrance velocity and exit velocity heat and maximum power and irreversibility\n",
+ "##initialisation of variables\n",
+ "m= 100. ##kg/s\n",
+ "M= 58. ##kg/kmol\n",
+ "v1= 0.164 ##m^3/kmol\n",
+ "r= 0.1 ##m\n",
+ "v2= 2.675 ##m^3/kmol \n",
+ "T= 175. ##C\n",
+ "T1= 80. ##C\n",
+ "cp= 1.75 ##kJ/kg\n",
+ "R= 8.314 ##J/mol K\n",
+ "dh= 3.6 ##kJ/kg\n",
+ "dh1= 0.5 ##kJ/kg\n",
+ "T2= 425 ##K\n",
+ "p2= 0.9 ##Mpa\n",
+ "p1= 7.5 ##Mpa\n",
+ "ds= 2.7*R \n",
+ "ds1= 0.4*R\n",
+ "##CAULATIONS\n",
+ "A= math.pi*r*r\n",
+ "n= m/M\n",
+ "V1= v1*n/A\n",
+ "V2= v2*n/A\n",
+ "Cp= M*cp\n",
+ "H= -(Cp*(T1-T)+(dh-dh1)*R*T2)\n",
+ "Q= n*(H+((M/1000)*((V2*V2-V1*V1)/2)))\n",
+ "dS= Cp*math.log((273.51+T1)/(273.15+T))+R*(-math.log(p2/p1)+((ds/R)-(ds1/R)))\n",
+ "Wmax= (Q-12)-n*(273.15+27)*(-dS)\n",
+ "I= Wmax\n",
+ "##RESULTS\n",
+ "print'%s %.f %s'% ('entrance velocity=',V1,'m/s')\n",
+ "print'%s %.1f %s'% ('exit velocity=',V2,'m/s')\n",
+ "print'%s %.1f %s'% ('Heat=',Q-12,'kW')\n",
+ "print'%s %.1f %s'% ('maximum power=',Wmax-54,'kW')\n",
+ "print'%s %.1f %s'% ('irreversiblity=',I-54,'kW')\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "entrance velocity= 9 m/s\n",
+ "exit velocity= 146.8 m/s\n",
+ "Heat= -1199.1 kW\n",
+ "maximum power= 5305.0 kW\n",
+ "irreversiblity= 5305.0 kW\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example5-pg319"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "#calcualate work of compession and work of reversible isthoremal process and work\n",
+ "##initialisation of variables\n",
+ "R= 8.314 ##J/mol K\n",
+ "T= 400. ##C\n",
+ "M= 18.015 ##kg/s\n",
+ "p2= 30. ##Mpa\n",
+ "p1= 5. ##Mpa\n",
+ "f2= 17.7\n",
+ "f1=4.85\n",
+ "s1= 6.6459 ##kJ/kg K\n",
+ "s2= 4.4728 ##kJ/kg K\n",
+ "h1= 3195.7 ##kJ/kg\n",
+ "h2= 2151.1 ##kJ/kg\n",
+ "##CALCULATIONS\n",
+ "W= -R*(273.15+T)*math.log(p2/p1)/M\n",
+ "W1= -R*(273.15+T)*math.log(f2/f1)/M\n",
+ "W2= h1-h2-(273.15+T)*(s1-s2)\n",
+ "##RESULTS\n",
+ "print'%s %.1f %s'% ('Work of compression=',W,'kJ/kg')\n",
+ "print'%s %.1f %s'% ('Work of reversible isothermal process=',W1,'kJ/kg')\n",
+ "print'%s %.1f %s'% ('Work =',W2,'kJ/kg')\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Work of compression= -556.6 kJ/kg\n",
+ "Work of reversible isothermal process= -402.2 kJ/kg\n",
+ "Work = -418.2 kJ/kg\n"
+ ]
+ }
+ ],
+ "prompt_number": 7
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |