diff options
Diffstat (limited to 'Examples_in_Thermodynamics_Problems')
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/Chapter1.ipynb | 562 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/Chapter2.ipynb | 150 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/Chapter3.ipynb | 469 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/Chapter4.ipynb | 514 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/Chapter5.ipynb | 542 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/Chapter6.ipynb | 792 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/Chapter7.ipynb | 539 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/Chapter8.ipynb | 502 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/Chapter9.ipynb | 227 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/README.txt | 10 | ||||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/screenshots/Heat_supplied.png | bin | 0 -> 74751 bytes | |||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/screenshots/Ideal_efficiency.png | bin | 0 -> 67751 bytes | |||
| -rwxr-xr-x | Examples_in_Thermodynamics_Problems/screenshots/Thermal_capacity.png | bin | 0 -> 64520 bytes |
13 files changed, 4307 insertions, 0 deletions
diff --git a/Examples_in_Thermodynamics_Problems/Chapter1.ipynb b/Examples_in_Thermodynamics_Problems/Chapter1.ipynb new file mode 100755 index 00000000..0d02f117 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/Chapter1.ipynb @@ -0,0 +1,562 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:48ea781a71227cd2ed588b90ba4f6ceb57773ac789c06adb226b03cd21058126"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 1: HEATING AND EXPANSION OF GASES ENTROPY"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.1 Page 1"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p1=280#lb/in^2\n",
+ "v=2#ft^3\n",
+ "p2=20#lb/in^2\n",
+ "v2=18.03#ft^3\n",
+ "\n",
+ "#calculation\n",
+ "W=144*(p1*v-p2*v2)/(1.2-1)#ft/lb\n",
+ "\n",
+ "#result\n",
+ "print\"The volume and work done during the expansion is\",round(W,2),\"ft/lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The volume and work done during the expansion is 143568.0 ft/lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.2 Page 2"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "v=2#ft^3\n",
+ "v2=20#ft^3\n",
+ "p=100000#ft lb\n",
+ "v2=10.41#lb/in^2\n",
+ "v3=10#lb/in^2\n",
+ "p1=1.3#lb\n",
+ "p2=(v2*199.5)/9.95#lb/in^2\n",
+ "\n",
+ "#calculation\n",
+ "P=(p2/v3-v2)#lb/in^2\n",
+ "\n",
+ "#result\n",
+ "print\"The initial andfinal pressure is\",round(P,4),\"lb/in^2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The initial andfinal pressure is 10.4623 lb/in^2\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.4 Page 3"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "Cp=0.24#lb/in^2\n",
+ "Cv=0.18#ft^3\n",
+ "p1=5#lb/in^2\n",
+ "T1=20#Degree C\n",
+ "T2=150#Degree C\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=p1*Cp*(T2-T1)#C.H.U\n",
+ "H=p1*Cv*(T2-T1)#C.H.U\n",
+ "Gamma=Cp/Cv#lb/in^2\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the constant pressure is\",round(W,2),\"C.H.U\"\n",
+ "print\"the constant volume the value of gas is\",round(Gamma,2),\"lb/in^2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the constant pressure is 156.0 C.H.U\n",
+ "the constant volume the value of gas is 1.33 lb/in^2\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.5 Page 4"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "Gama=1.33#ft/lb\n",
+ "p=100#lb/in^2\n",
+ "p1=20#lb/in^2\n",
+ "v2=10.05#ft^3\n",
+ "v=3#ft/lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=144*(p*v-p1*v2)/0.33#ft lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The work done is\",round(W,2),\"lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The work done is 43200.0 lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.8 Page 5"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=3.74#ft/lb\n",
+ "p1=2.48#ft/lb\n",
+ "v=5.70#ft lb\n",
+ "Cv=0.21#ft/lb\n",
+ "P=440.00#lb/in^2\n",
+ "P1=160.00#lb/in^2\n",
+ "P2=14.00#lb/in^2\n",
+ "T=25.00#degree C\n",
+ "T1=100.00#F\n",
+ "vs=(pi*(p1)**2/4)*(p/1728.00)#ft^3\n",
+ "vc=5.70#ft^3\n",
+ "v1=4.70#ft^3\n",
+ "v2=vs/v1#ft^3\n",
+ "v3=0.01273#ft^3\n",
+ "T2=298.00#F\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=(P2*144.00*v3)/(T2*T1)#lb\n",
+ "T3=((P1*144.0*1.0)/(P2*144.0*7.0)*T2)#Degree C\n",
+ "T4=(P*T3)/P1#Degree C\n",
+ "H=W*Cv*(T4-T3)#C.H.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The heat supplied during explosion is\",round(H,5),\"C.H.U\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat supplied during explosion is 0.15398 C.H.U\n"
+ ]
+ }
+ ],
+ "prompt_number": 39
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.9 Page 6"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "v=10.0#ft^3\n",
+ "p=100.0#lb/in^2\n",
+ "p1=18.0#lb/in^2\n",
+ "v1=50.0#ft^3\n",
+ "n=log(p/p1)/log(5)\n",
+ "gama=1.40#air\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=(144.0*(p*v-p1*v1))/(n-1)#ft lb\n",
+ "H=(gama-n)/(gama-1)*W#ft lb\n",
+ "E=W-H#ft lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The heat supplied and the change of internal energy\",round(E,2),\"ft lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat supplied and the change of internal energy 36000.0 ft lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 37
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.11 Page 7"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "v=2.0#ft^3\n",
+ "p=1100.0#lb/in^2\n",
+ "t1=44.0#Degree C\n",
+ "t2=15.0#Degree C\n",
+ "p1=300.0#lb/in^2\n",
+ "t3=3.0#Degree c\n",
+ "Cv=0.17#ft/lb\n",
+ "T=273.0#F\n",
+ "R=96.0#ft lb\n",
+ "p3=300.0#lb/in^2\n",
+ "n=1.12#lb\n",
+ "gama=1.404#lb\n",
+ "W=((144*p*v)/(T+t1))/R#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "Wc=W*Cv*(t1-t2)#C.H.U\n",
+ "p2=p*(T+t2)/(T+t1)#lb /in^2\n",
+ "A=(144*p3*v)/(R*276)#lb\n",
+ "W1=(A/W)*v#ft^3\n",
+ "H=((gama-n)/(gama-1))*(144*(p*0.65-p1*v)/(n-1))#ft lb\n",
+ "H1=H/1400#C.H.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the heat was lost by all the air in the vessel before leakage began\",round(Wc,4),\"C.H.U\"\n",
+ "print\"the heat was lost or gainned leakage by the air\",round(H1,4),\"C.H.U\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the heat was lost by all the air in the vessel before leakage began 51.3218 C.H.U\n",
+ "the heat was lost or gainned leakage by the air 69.2928 C.H.U\n"
+ ]
+ }
+ ],
+ "prompt_number": 42
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.13 Page 9"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "h=0.218#ft^3\n",
+ "h1=0.156#ft^3\n",
+ "n=0.249#lb\n",
+ "h2=0.178#lb\n",
+ "c=0.208#lb\n",
+ "c1=0.162#lb\n",
+ "w1=1.0#ft^3\n",
+ "p=150.0#lb/in^2\n",
+ "T=100.0#Degree C\n",
+ "T1=373.0#F\n",
+ "Cp=(h*0.2312)+(n*0.3237)+(c*0.4451)#C.H.U/lb\n",
+ "Cv=(h1*0.2312)+(h2*0.3237)+(c1*0.4451)#C.H.U//lb\n",
+ "R=1400*(Cp-Cv)#ft lb units\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=(144*p*w1)/(R*T1)#lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The characteristic constant of the gas is\",round(W,4),\"lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The characteristic constant of the gas is 0.7157 lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 43
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.20 Page 12"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "T=200.0#Degree C\n",
+ "p=150.0#lb/in^2\n",
+ "v=12.0#ft^3\n",
+ "R=96.0#Lb\n",
+ "T1=473.0#F\n",
+ "T2=273.0#F\n",
+ "j=1400.0#lb\n",
+ "Cv=0.169#lb/in^2\n",
+ "v1=(R*T1)/(p*144)#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "Fhi=(R/j)*log(v/v1)+Cv*log(T2/T1)#rank\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The change of entropy is\",round(Fhi,4),\"rank\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The change of entropy is 0.0266 rank\n"
+ ]
+ }
+ ],
+ "prompt_number": 44
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.22 Page 13"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "v=10.0#ft^3\n",
+ "T=20.0#Degree C\n",
+ "p=15.0#lb in^2\n",
+ "p1=200.0#lb//in^2\n",
+ "gama=1.41 #lb\n",
+ "Cv=0.169#lb\n",
+ "v2=1.153#ft^3\n",
+ "j=1400.0#lb\n",
+ "T1=293.0#F\n",
+ "T2=451.0#F\n",
+ "T1=((p1*v2)/(p*v))*T1#Degree C\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "R=Cv*j*(gama-1)\n",
+ "W=0.816#lb\n",
+ "Fhi=Cv*((gama-1.2)/(1.2-1))*log(T1/T2)*W#rnak\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The change of entropy is\",round(Fhi,5),\"rank\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The change of entropy is -0.00018 rank\n"
+ ]
+ }
+ ],
+ "prompt_number": 46
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.23 Page 15"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "p=1.0#lb\n",
+ "T=200.0#Degree C\n",
+ "p1=15.0#lb/in^2\n",
+ "v1=4.0#ft^3\n",
+ "gama=1.41#lb\n",
+ "Cv=0.169#lb\n",
+ "J=1400.0#lb\n",
+ "n=1.2\n",
+ "T=473.0#F\n",
+ "v2=16.1#ft^3\n",
+ "T1=473.0#F\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T2=(p1*v2)/(p*v1)*T1#Degree C\n",
+ "R=Cv*J*(gama-p)#lb/in^2\n",
+ "Fhi=0.1772*log(1.317)#rank\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the change of entropy from intial condition is\",round(Fhi,5),\"rank\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the change of entropy from intial condition is 0.04879 rank\n"
+ ]
+ }
+ ],
+ "prompt_number": 47
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 1.26 Page 16"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "w=0.066#ft^3\n",
+ "p=14.7#lb/in^2\n",
+ "w1=14.2#lb/in^2\n",
+ "w2=2780.0#lb/in^2\n",
+ "g=0.038#lb\n",
+ "a=28.9#lb\n",
+ "R=w2/w1#for gas\n",
+ "R1=93.0#for air\n",
+ "T=273.0#F\n",
+ "V=0.4245#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=(p*144*w)/(T*R)#lb\n",
+ "m=(g-W)#lb gas\n",
+ "T2=(V+w)#ft^3\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The volume of mixture is\",round(T2,4),\"ft^3\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The volume of mixture is 0.4905 ft^3\n"
+ ]
+ }
+ ],
+ "prompt_number": 48
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/Chapter2.ipynb b/Examples_in_Thermodynamics_Problems/Chapter2.ipynb new file mode 100755 index 00000000..d7128b61 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/Chapter2.ipynb @@ -0,0 +1,150 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:8ff0faa4ee8798b8afb099ec4225e3d1dba96fda9be0b2724a96827e87a8c2b0"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 2: AIR CYCLE EFFICIENCIES"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 2.2 Page 19"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "T1=(100+273)#Degree C\n",
+ "T2=(300+273)#degree C\n",
+ "T=(1-T1/T2)*100#F\n",
+ "lam=0.41#in\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "R=log(T2)-log(T1)#lb/in^2\n",
+ "r=2.849#ratio of compression\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The ideal efficiency and the compression ratio is\",round(r,3)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The ideal efficiency and the compression ratio is 2.849\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 2.4 Page 20"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "r=0.60#in\n",
+ "v=3#in\n",
+ "p=15.4#lb\n",
+ "r=5#in\n",
+ "P=2000#r p m\n",
+ "V=19000#B.Th.U Per lb\n",
+ "lam=1.41 #lb\n",
+ "n=0.4831#percent\n",
+ "P=15.4/4#lb\n",
+ "H=P*V#B.Th.U\n",
+ "l=4.5#lb\n",
+ "A=9#lb\n",
+ "S=1000#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "R=0.60*n*100#percent\n",
+ "C=H*R#B.Th.U\n",
+ "I=(C*778)/(60*33000)#lb\n",
+ "P1=(I*12*4*33)/(l*A*pi)#lb/in^2\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The mean efficity pressure is\",round(P1,3)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The mean efficity pressure is 10372.115\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 2.5 Page 21"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "v=15#in\n",
+ "S=(5*14/100)#ln\n",
+ "lam=1.4#in\n",
+ "v1=1.7#in\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "N=(1-0.38)*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The ideal effiecncy for an engine is\",round(N,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The ideal effiecncy for an engine is 62.0 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 7
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/Chapter3.ipynb b/Examples_in_Thermodynamics_Problems/Chapter3.ipynb new file mode 100755 index 00000000..cdd25fa7 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/Chapter3.ipynb @@ -0,0 +1,469 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:a9b498c87a8bc557e26090cded88daaf4b9544d15a7aef988c258f78451cb8ae"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 3: PROPERTIES OF STEAM"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3.1 Page 30"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "p=100#lb/in^2\n",
+ "x=0.8#lb\n",
+ "t1=164#degree C\n",
+ "t2=4.45 #ft^3\n",
+ "p1=0.016#ft^3\n",
+ "h1=493.4#C.H.U/lb\n",
+ "h2=165.9#C.H.U/lb\n",
+ "S=h2+h1#C.H.U/lb\n",
+ "w=(144*p)/1400*(t2-p1)#C.H.U/lb\n",
+ "H=h2+(x*h1)#C.H.U//lb\n",
+ "w1=(x*144*p)/1400*(t2-p1)#C.H.U\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "E=S-w#C.H.U/lb\n",
+ "IE=H-w1#C.H.U/lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The steam is total heat dry and satured is\",round(E,3),\"C.H.U/lb\"\n",
+ "print\"Total heat of wet steam\",round(IE,3),\"C.H.U/lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The steam is total heat dry and satured is 614.96 C.H.U/lb\n",
+ "Total heat of wet steam 524.135 C.H.U/lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3.2 Page 31"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "t1=35#degree C\n",
+ "p=100#lb/in^2\n",
+ "L=435#C.H.U\n",
+ "L2=539.3#C.H.U\n",
+ "h1=165.9#H.C.U/lb\n",
+ "h2=493.4#C.H.U/lb\n",
+ "S=(h1-t1)#C.H.U\n",
+ "h3=304.1#C.H.U\n",
+ "h4=335#C.H.U/lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "X1=h3/h2#C.H.U/lb\n",
+ "X2=h4/L2#C.H.U/lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The heat giving to the water and steam is\",round(X2,3),\"C.H.U/lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat giving to the water and steam is 0.621 C.H.U/lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3.3 Page 32"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "p=35#lb/in^2\n",
+ "w=1425#lb\n",
+ "q=1474#lb\n",
+ "s1=126.7#C.H.U/lb\n",
+ "s2=28#C.H.U/lb\n",
+ "t1=5#degree C\n",
+ "t2=28#degree C\n",
+ "L1=521.4#C.H.U/lb\n",
+ "w1=245#lb\n",
+ "w2=0.2#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=(s1-s2)+L1#C.H.U/lb\n",
+ "H=q*(t2-t1)#C.H.U/lb\n",
+ "T=H/W#lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The total equivalent is\",round(T,3),\"lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The total equivalent is 54.672 lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3.4 Page 33"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "p=100#lb/in^2\n",
+ "w=2400#lb\n",
+ "t1=15#degree C\n",
+ "s1=165.9#C.H.U/lb\n",
+ "x=0.9#lb\n",
+ "L2=493.4#C.H.U/lb\n",
+ "t2=65#degree\n",
+ "x4=0.8#lb\n",
+ "s3=64.8#C.H.U/lb\n",
+ "w1=2000#lb\n",
+ "w2=2400#lb\n",
+ "b1=12400#lb\n",
+ "b2=22000#lb\n",
+ "p1=4400#lb\n",
+ "n=421.65#lb\n",
+ "h1=w2*[s1+(x*L2)]#C.H.U/hr\n",
+ "h2=w1*[s1+(x4+L2)]#C.H.U/hr\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T=w*[(s1-t1)+(x*L2)]#C.H.U/hr\n",
+ "T1=w1*[(s1-s3)+(x4*L2)]#C.H.U/hr\n",
+ "H=T+T1#C.H.U/hr\n",
+ "X=n/L2#C.H.U/lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The thermal capacity of the pipe is\",round(X,3),\"C.H.u/hr\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The thermal capacity of the pipe is 0.855 C.H.u/hr\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3.5 Page 35"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "w1=4.5#lb\n",
+ "s1=45.5#lb\n",
+ "p1=165#lb/in^2\n",
+ "T=140#Degree C\n",
+ "h1=30#in\n",
+ "h2=4#in\n",
+ "p2=0.49#ln/in^2\n",
+ "T1=(w1+s1)#lb\n",
+ "T2=103.5#Degree C\n",
+ "T3=140#Degree\n",
+ "h3=0.48#in\n",
+ "x=0.988#berfore throttling\n",
+ "T=[103.12+537.1+h3*(T3-T2)]#C.H.U/lb\n",
+ "x1=0.012#lb of water\n",
+ "X=s1*x1#lb water\n",
+ "w=50#lb of steam\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "P=h2+h1#in of mercury\n",
+ "P1=s1*x1#lb/in^2\n",
+ "T4=w1+P1#lb\n",
+ "D=(w-T4)/w#lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The dryness of steam with a combined is\",round(D,3),\"lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The dryness of steam with a combined is 0.899 lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3.6 Page 36"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "w=40.0#lb\n",
+ "w1=380.0#lb\n",
+ "t1=80.0#Degree\n",
+ "p=85.0#lb/in^2\n",
+ "p1=15.0#lb/in^2\n",
+ "W=w+w1#lb/hr\n",
+ "P=p+p1#lb/in^2\n",
+ "T=659.3#C.H.U/lb\n",
+ "d=10.0#h.p\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=W*T-w1*t1#C.H.U/hr\n",
+ "H1=(d*33000.0*60.0)/1400.0#C.H.U/hr\n",
+ "T1=H1/H*100.0#percent\n",
+ "D=w1/(w1+w)#C.H.U/hr\n",
+ "H2=(W*(99.6+D*539.3)-w1*t1)#C.H.U/hr\n",
+ "T2=H-H2#C.H.U/hr\n",
+ "H3=T2-H1#C.H.U/hr\n",
+ "E=(1400.0*H3)/(60.0*33000.0)#h.p\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The amount of radiations from the engine is\",round(E,3),\"h.p\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The amount of radiations from the engine is 11.311 h.p\n"
+ ]
+ }
+ ],
+ "prompt_number": 9
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3.10 Page 37"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "w=40.0#lb\n",
+ "w1=380.0#lb\n",
+ "t1=80.0#Degree\n",
+ "p=85.0#lb/in^2\n",
+ "p1=15.0#lb/in^2\n",
+ "W=w+w1#lb/hr\n",
+ "P=p+p1#lb/in^2\n",
+ "T=659.3#C.H.U/lb\n",
+ "d=10.0#h.p\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=W*T-w1*t1#C.H.U/hr\n",
+ "H1=(d*33000*60)/1400#C.H.U/hr\n",
+ "T1=H1/H*100#percent\n",
+ "D=w1/(w1+w)#C.H.U/hr\n",
+ "H2=(W*(99.6+D*539.3)-w1*t1)#C.H.U/hr\n",
+ "T2=H-H2#C.H.U/hr\n",
+ "H3=T2-H1#C.H.U/hr\n",
+ "E=(1400*H3)/(60*33000)#h.p\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The amount of radiations from the engine is\",round(E,3),\"h.p\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The amount of radiations from the engine is 11.311 h.p\n"
+ ]
+ }
+ ],
+ "prompt_number": 11
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3.12 Page 39"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "p=120#lb/in^2\n",
+ "ts=264#degree C\n",
+ "T1=(273+130.6)#F\n",
+ "v=0.0171#ft^3/lb\n",
+ "L1=518.4#lb\n",
+ "T2=(273+171.9)#F\n",
+ "L2=487.4#lb\n",
+ "Cp=0.48#lb\n",
+ "L=0.0894/Cp#lb\n",
+ "Ts=T2*1.205#degree\n",
+ "ta=536-273#Degree C\n",
+ "T=649.9#C.H.U\n",
+ "S=131.2#C.H.U\n",
+ "w=(144*40)/1400*(10.49-v)#C.H.U\n",
+ "C=T-S#C.H.U\n",
+ "I=C-w#C.H.U\n",
+ "E=(704.7-57.8)#C.H.U\n",
+ "E1=E-606.5#C.H.U\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "E1=E-606.5#C.H.U\n",
+ "H=(704.7-T)#C.H.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"Heat and internal energy after each operation is\",round(H,3),\"h.p\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Heat and internal energy after each operation is 54.8 h.p\n"
+ ]
+ }
+ ],
+ "prompt_number": 12
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 3.17 Page 42"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "A=28.1#in Hg vacuum\n",
+ "a=0.93#lb/in^2\n",
+ "T=33#Degree\n",
+ "p=0.729#lb/in^2\n",
+ "P=-p+a#lb/in^2\n",
+ "p1=120000#lb\n",
+ "p2=28.1#in\n",
+ "a1=0.9#ln\n",
+ "p3=1000#lb\n",
+ "t=15#degree C\n",
+ "A1=(a1*(p1/(60*p3)))#lb/mim\n",
+ "v=(A1*96*306)/(144*P)#ft^3 of air per min\n",
+ "V=37.3+a1*610#C.H.U/lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=(V-T)#C.H.U\n",
+ "W=(H/t)*(p1/60)#gal/min\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The water per minute in cubic feet per minute passing to air extractor\",round(W,3),\"gal/min\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The water per minute in cubic feet per minute passing to air extractor 73773.333 gal/min\n"
+ ]
+ }
+ ],
+ "prompt_number": 13
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/Chapter4.ipynb b/Examples_in_Thermodynamics_Problems/Chapter4.ipynb new file mode 100755 index 00000000..a68b7f49 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/Chapter4.ipynb @@ -0,0 +1,514 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:cafa634c5fbcfd5447496e6358027f61ba038aded335f02d8eee875e161267c7"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 4: THE STEAM ENGINE"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.1 Page 46"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=90#lb/in^2\n",
+ "x1=0.9#lb\n",
+ "p1=10#lb/in^2\n",
+ "x2=0.81#lb\n",
+ "s1=161.5#lb.in^2\n",
+ "s2=89.1#lb.in^2\n",
+ "L1=496.8#lb.in^2\n",
+ "L2=545.5#lb.in^2\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "bc=(s1-s2)+(x1*L1)#C.H.U/lb\n",
+ "da=x2*L2#C.H.U/lb\n",
+ "W=bc-da#C.H.U/lb\n",
+ "R=W/bc*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the work done per\",round(R,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the work done per 14.949 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.2 Page 47"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "h=1600.0#i.h.p\n",
+ "h1=20000.0#lb\n",
+ "h2=230.0#lb/in^2\n",
+ "T1=293.3#Degree C\n",
+ "x=25.91#in\n",
+ "v=30.0#in\n",
+ "T2=201#Degree C\n",
+ "T=T1-T2#degree C\n",
+ "x2=0.845#lb\n",
+ "L2=566.51#lb\n",
+ "s1=724#lb\n",
+ "h3=1400#C.H.U/hr\n",
+ "x=33000#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=671.48#C.H.U/lb\n",
+ "ea=x2*L2#C.H.U/lb\n",
+ "W=H-ea#C.H.U/lb\n",
+ "R=W/H*100#percent\n",
+ "S=h2*s1#C.H.U\n",
+ "I=((h*x*60)/(h3*h1*s1))*100#percent\n",
+ "R1=I/R*100#pecent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The indicated thermal efficiency ratio is\",round(R1,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The indicated thermal efficiency ratio is 54.433 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 8
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.3 Page 49"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "h1=180#lb/in^2\n",
+ "h2=3#lb/in^2\n",
+ "r1=60#percent\n",
+ "r2=90#percent\n",
+ "p3=100#lb/in^2\n",
+ "p4=10#lb/in^2\n",
+ "v1=4.4#ft^3/lb\n",
+ "v2=2*v1#ft^3\n",
+ "p=44#lb/in^2\n",
+ "x2=0.95#ft^3\n",
+ "s1=165.9#lb\n",
+ "s2=89.1#lb\n",
+ "L1=493.4#lb\n",
+ "H=(s1-s2)+L1#C.H.U/lb\n",
+ "W=65.8#C.H.U/lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "R=W/H*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The rankine efficiency of the engine is\",round(R,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The rankine efficiency of the engine is 11.54 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 9
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.10 Page 52"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=85.0#lb/in^2\n",
+ "h=210.0#i.p.m\n",
+ "h1=8.0#in\n",
+ "h2=2.5#in\n",
+ "h3=20.0#in\n",
+ "x=0.75#in\n",
+ "p1=100.0#ln/in^2\n",
+ "x1=33000.0#in\n",
+ "p2=15.0#lb/in^2\n",
+ "v2=pi/4*(h1/12)**2*(h3/12)#ft^3\n",
+ "A=144*(29.08*1.6931-8.724)#ft/lb\n",
+ "d=x*A#ft/lb\n",
+ "v3=0.5816#ft^3\n",
+ "P=d/(144*v3)#lb/in^2\n",
+ "P1=pi/4*64#in^2\n",
+ "r=25*pi/16#in^2\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=P*(h3/12)*P1*h/(x1)\n",
+ "I=(P*(h3/12)*(P1-r)*h)/(x1)#I.h.P\n",
+ "T=H+I#I.h.p\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the h.p diameter of the piston and piston rod is\",round(T,3),\"I.h.p\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the h.p diameter of the piston and piston rod is 52.982 I.h.p\n"
+ ]
+ }
+ ],
+ "prompt_number": 11
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.14 Page 55"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "a=1.025#in^2\n",
+ "h=18.0#in\n",
+ "h1=24.0#in\n",
+ "x=8.2#percent\n",
+ "v=15.0#in\n",
+ "v2=6.9#ft^3\n",
+ "p=0.74#lb/in^2\n",
+ "p1=50.0#lb/in^2\n",
+ "p2=83.0#lb/in^2\n",
+ "P3=48.0#lb/in\n",
+ "P1=29.8#lb/in^2\n",
+ "P2=14.6#lb/in^2\n",
+ "h2=29.8#in\n",
+ "D=(pi/4)*(3.0/2.0)**2*2#ft^3\n",
+ "v1=23400.0#ft.lb\n",
+ "W=a*v1#ft.lb\n",
+ "V=0.082*D#ft^3\n",
+ "Q=1.530#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "I=V+Q#ft^3\n",
+ "P=P3+P2#lb/in^2\n",
+ "V1=p*v2#ft^3\n",
+ "W1=I/V1#lb\n",
+ "S=p2+P2#l/in^2\n",
+ "H=659.06#C.H.U/lb\n",
+ "T=W/(H*W1*1400)*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The thermal efficiency of the engine is\",round(T,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The thermal efficiency of the engine is 7.294 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.16 Page 57"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "v=4.6#ft^3\n",
+ "h=5.0#percent\n",
+ "p=60.0#lb/in^2\n",
+ "p1=0.8#ft^3\n",
+ "p2=19.0#lb/in^2\n",
+ "a=100.0#r.p.m\n",
+ "h1=5920.0#lb\n",
+ "W=h1/(2*a*p)#lb\n",
+ "V=(0.25*v)#ft^3\n",
+ "v1=21.07#ft^3\n",
+ "w=V/v1#lb\n",
+ "H=W+w#lb\n",
+ "v2=H*7.17#ft^3\n",
+ "P=w*v2#ft^3\n",
+ "P1=0.675*v#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "DP=P1/v2#ft^2\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The assumptions do you make in working out the dryness of the steam is\",round(DP,3),\"ft^3\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The assumptions do you make in working out the dryness of the steam is 0.79 ft^3\n"
+ ]
+ }
+ ],
+ "prompt_number": 16
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.17 Page 58"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "h=0.08#lb\n",
+ "p=60#lb/in^2\n",
+ "p1=0.50#lb/in^2\n",
+ "v=0.5#ft^3\n",
+ "v1=7.17#ft^3\n",
+ "V=h*v1#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=p1/v1#lb\n",
+ "I=v/v1#lb\n",
+ "M=h-I#lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The dryness of the steam at this pressure and missing quantity is\",round(M,6),\"lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The dryness of the steam at this pressure and missing quantity is 0.010265 lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 18
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.19 Page 59"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p1=120#lb/in^2\n",
+ "p2=15#lb/in^2\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "v=1.65#lb\n",
+ "D=sqrt(v)#lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The above pressure are by gauge is\",round(D,3),\"lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The above pressure are by gauge is 1.285 lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 19
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.21 Page 61"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=150.0#lb/in^2\n",
+ "x=198.0#r.p.m\n",
+ "x1=33000.0#lb\n",
+ "h=2700.0#lb\n",
+ "h1=1400.0#lb\n",
+ "h2=51600.0#lb\n",
+ "r=165.0#C.H.U/lb\n",
+ "s=60.0#lb\n",
+ "t=48.0#Degree C\n",
+ "t1=11.0#degree C\n",
+ "t2=36.0#Degree C\n",
+ "P1=(40*75*t2*x)/(12*x1)#lb\n",
+ "P2=(38*70*t2*x)/(12*x1)#lb\n",
+ "L1=(t1*300*t2*x)/(12*x1)#lb\n",
+ "L2=(12*295*t2*x)/(12*x1)#lb\n",
+ "T=P1*P2*L1*L2#lb\n",
+ "H=5294#C.H.U/min\n",
+ "T1=h/s#lb/min\n",
+ "H1=T1*663#lb/min\n",
+ "H2=(h2/s*(36-11)+(h/s)*(t))#C.H.U\n",
+ "H3=(h/60)*t#C.H.U\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "TE=H/H1-H3*100#percent\n",
+ "R=r/(663-t)*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The rankine efficiency is\",round(R,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The rankine efficiency is 26.829 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 21
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.23 Page 64"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p1=100#ln/in^2\n",
+ "p2=2.5#lb/in^2\n",
+ "p3=20#lb/in^2\n",
+ "d=0.75#lb\n",
+ "p=0.5#lb\n",
+ "r=16#in\n",
+ "p4=p1/r#lb/in^2\n",
+ "P5=50#lb/in^2\n",
+ "W1=13960#ft/lb\n",
+ "W2=19040#ft/lb\n",
+ "T=33000#ft/lb\n",
+ "v=4.43#ft^3\n",
+ "v1=v*d#ft^3\n",
+ "W3=T*v1#ft/lb\n",
+ "Hp=3416#ft/lb\n",
+ "Lp=3416#ft/lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=Lp*v1#ft lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The thermal efficiency of a compound steam and work done is\",round(W,3),\"ft lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The thermal efficiency of a compound steam and work done is 11349.66 ft lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 22
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/Chapter5.ipynb b/Examples_in_Thermodynamics_Problems/Chapter5.ipynb new file mode 100755 index 00000000..e3603cbc --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/Chapter5.ipynb @@ -0,0 +1,542 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:983f6479fd716219923fc039ee4060b375eda2b85fbdeb1a9e63c25f153f7ec9"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 5: AIR COMPRESSORS AND MOTORS REFRIGERATIONS"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.1 Page 66"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "a=7.0#in\n",
+ "b=10.0#in\n",
+ "c=12.0#in\n",
+ "r=96.0#in\n",
+ "p1=15.0#lb/in^2\n",
+ "p2=100.0#lb/in^2\n",
+ "T=16.0#Degree C\n",
+ "gama=1.4#in\n",
+ "h=120.0#r.p.m\n",
+ "T1=T+273#C absolute\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "v1=(pi/4)*(a/c)**2*(b/c)#ft^3\n",
+ "w=(p1*144*v1)/(r*T1)#lb\n",
+ "w1=h*w#lb\n",
+ "W=1680*(1.72-1)#ft lb\n",
+ "I=144*p1*v1*log(p2/p1)#ft lb\n",
+ "E=I/W*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The ideal efficiency is defined as the ratio of tthis work is\",round(E,4),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The ideal efficiency is defined as the ratio of tthis work is 75.4482 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.2 Page 68"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "h1=16#i.h.p\n",
+ "p1=100#lb/in^2\n",
+ "p2=15#lb/in^2\n",
+ "R=275#R.p.m\n",
+ "h=550#ft/min\n",
+ "q=33000#in^2\n",
+ "v1=4.85#lb\n",
+ "B=8.53#in\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "M=(p1/v1)-p2+(p1/v1-p2)*1/0.2\n",
+ "S=h/(2*R)#ft\n",
+ "I=(q*h1)/(M*S*R)#in^2\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The effect of the clearance volume is\",round(I,3),\"in^2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The effect of the clearance volume is 56.954 in^2\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.3 Page 69"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "h=100#ft^3\n",
+ "t=15#degree C\n",
+ "p=120#lb/in^2\n",
+ "gama=1.3#in\n",
+ "t1=15#Degree C\n",
+ "M=((144*t*h*2.6)/(0.3)*(1.271-1))#ft lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "V=sqrt(p/t)#ft lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"Compare the values of the two cylinders is\",round(V,3),\"ft lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Compare the values of the two cylinders is 2.828 ft lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.5 Page 71"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "h=0.2#ft^3\n",
+ "v=10#percent\n",
+ "T=15#degree c\n",
+ "p=30#lb/in62\n",
+ "t1=15#Degree C\n",
+ "p1=60#lb/in^2\n",
+ "v1=2.2#ft^3\n",
+ "v3=0.328#ft^3\n",
+ "A=(v1-v3)#ft^3\n",
+ "v2=1.341#ft^3\n",
+ "V=v2-h#ft^\n",
+ "t2=288#Degree C\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T2=(t2*p*v2)/(t1*v1)#Degree C absolute\n",
+ "v5=(t2/T2)*V#ft^3\n",
+ "v7=0.164#ft^3\n",
+ "v8=v5-(v7/11)*v5\n",
+ "v6=v8/(1-v7/11)#ft^3\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The required volume of the H.P cylinder including clearance is\",round(v6,3),\"ft^3\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The required volume of the H.P cylinder including clearance is 0.936 ft^3\n"
+ ]
+ }
+ ],
+ "prompt_number": 7
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.6 Page 73"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "p1=80#lb/in^2\n",
+ "p2=20#lb/in^2\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "P=sqrt(p1*p2)#lb/in62\n",
+ "V=P/p1#stroke\n",
+ "W=p2/P#stroke\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the ratio of cut off to length of stroke is\",round(W,3),\"stroke\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the ratio of cut off to length of stroke is 0.5 stroke\n"
+ ]
+ }
+ ],
+ "prompt_number": 8
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.9 Page 75"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "p1=25.0#lb/in^2\n",
+ "p2=50.0#lb/in^2\n",
+ "p3=75.0#lb/in^2\n",
+ "p4=100.0#lb/in^2\n",
+ "v1=29.2#ft^3\n",
+ "v2=28.8#ft^3\n",
+ "v3=28.1#ft^3\n",
+ "v4=27.2#ft^3\n",
+ "h=14.7#lb/in^2\n",
+ "v=3.0#percent\n",
+ "s=5.0#stroke\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "V=(pi*p1)/(4)*4#in^3\n",
+ "V1=v/p4*V#in^3\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The volume of efficiency of pressure is\",round(V1,3),\"in^3\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The volume of efficiency of pressure is 2.356 in^3\n"
+ ]
+ }
+ ],
+ "prompt_number": 10
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.12 Page 78"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "p1=15.0#lb/in^2\n",
+ "p2=60.0#lb/in^2\n",
+ "t=16.0#Degree C\n",
+ "Ta=273.0+t#Degree C absolute\n",
+ "T=1.486#lb/in^2\n",
+ "Td=Ta/T#Degree C absolute\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "P=Td/(Ta-Td)#Degree C absolute\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The lowest temperature and coefficient of per formance is\",round(P,3),\"Degree C absolute\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The lowest temperature and coefficient of per formance is 2.058 Degree C absolute\n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.14 Page 79"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "T1=30#Degree c\n",
+ "T2=-10#degree C\n",
+ "t1=263#F\n",
+ "t2=303#F\n",
+ "h1=20#Units\n",
+ "h2=79#C.H.U/lb\n",
+ "h=24#hours\n",
+ "T3=1#Degree C\n",
+ "p=2.2046#C.H.U/sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "P=h1*p#C.H.U/sec\n",
+ "T=t1/(t2-t1)#F\n",
+ "H=P*60#C.H.U\n",
+ "W=(H*1400)/T#ft/lb\n",
+ "hp=W/33000#h.p\n",
+ "W1=(H*60*h)/(80*2240)#tons\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the cycle is a perfect one\",round(W1,3),\"tons\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the cycle is a perfect one 21.259 tons\n"
+ ]
+ }
+ ],
+ "prompt_number": 15
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.15 Page 80"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "p1=930#lb/in^2\n",
+ "p2=440#lb/in^2\n",
+ "T=268#F\n",
+ "t1=25#F\n",
+ "t2=5#F\n",
+ "h1=19.4#C.H.U\n",
+ "h2=-1.8#C.H.U\n",
+ "h3=29#C.H.U\n",
+ "h4=58.6#C.H.U\n",
+ "d=0.6#C.H.U\n",
+ "d1=0.06#lb\n",
+ "d2=-0.01#lb\n",
+ "c=40#percent\n",
+ "h=24#hour\n",
+ "t3=10#C\n",
+ "d3=15#lb\n",
+ "h5=80#C.H.U\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "A=(h1-(h2))-(d1-(d2))*T#C.H.U\n",
+ "FD=A/T#units of entropy\n",
+ "AD=(d*h4/T-0.07-A/T)*T#C.H.U\n",
+ "W=4.28#C.H.U\n",
+ "T=AD/W#C.H.U\n",
+ "P=0.4*T#C.H.U\n",
+ "H=P*W*d3#C.H.U\n",
+ "H1=P*W*d3*60*h#C.H.U\n",
+ "H2=t3+h5#C.H.U\n",
+ "W1=H1/(H2*2240)#tond\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The many tons of ice would a machine working between the same limit and having a relative coefficient is\",round(W1,3),\"tons\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The many tons of ice would a machine working between the same limit and having a relative coefficient is 0.598 tons\n"
+ ]
+ }
+ ],
+ "prompt_number": 16
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.16 Page 82"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "t1=20.0#Degeree C\n",
+ "t2=-10.0#degree C\n",
+ "h=0.95#dry\n",
+ "t3=35.0#Degree C\n",
+ "h1=0.066#lb\n",
+ "h2=1.089#lb\n",
+ "v1=-0.033#lb\n",
+ "v2=1.193#lb\n",
+ "v3=0.508#lb\n",
+ "T1=263.0#F\n",
+ "T2=293.0#F\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T=T1/(T2-T1)#F\n",
+ "E=h1-(v1)#lb\n",
+ "C=0.1079#lb\n",
+ "CP=E/C#lb\n",
+ "A=CP*(T2-T1)-E*T1#C.H.U\n",
+ "F=A/T1#units of entropy\n",
+ "H=254.212#C.H.U\n",
+ "H2=274.447#C.H.U\n",
+ "W=(CP*(T2-T1)+h*1.023*(T2-T1)-E*T1)#C.H.U\n",
+ "P=H/W#C.H.U\n",
+ "V=A+v3*15-T1*v3*0.0507#C.H.U\n",
+ "H1=T1*(v3*0.0507+0.05*1.023)#C.H.U\n",
+ "N=H2/(W+V)#C.H.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The upper and lower temperature limits respectively are\",round(T,4),\"F\"\n",
+ "print\"The vapour compression cycle work done is\",round(H,3),\"C.H.U\"\n",
+ "print\"The vapour is now additional work done is\",round(N,3),\"C.H.U\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The upper and lower temperature limits respectively are 8.7667 F\n",
+ "The vapour compression cycle work done is 254.212 C.H.U\n",
+ "The vapour is now additional work done is 8.322 C.H.U\n"
+ ]
+ }
+ ],
+ "prompt_number": 20
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 5.18 Page 85"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan,log\n",
+ "h=0.8#dry\n",
+ "p=120#lb/in^2\n",
+ "p1=1#lb/in^2\n",
+ "t=100#Degree C\n",
+ "A=99.6-38.6-0.178*311.8#C.H.U\n",
+ "G=311.8#units of entropy\n",
+ "AF=440.52#C.H.U\n",
+ "H=399.82#lb/in^2\n",
+ "p=307#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T=H/p#C.H.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"theoretical coefficient pf performance as a refrigeratior is\",round(T,3),\"C.H.U\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "theoretical coefficient pf performance as a refrigeratior is 1.302 C.H.U\n"
+ ]
+ }
+ ],
+ "prompt_number": 21
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/Chapter6.ipynb b/Examples_in_Thermodynamics_Problems/Chapter6.ipynb new file mode 100755 index 00000000..63b34ea8 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/Chapter6.ipynb @@ -0,0 +1,792 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:d4d40e7a0bce77d845a19fb5891f8dca6717a9164c9a3127715de6d72e594b43"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 6: FLOW THROUGH NOZZLES-STEAM TURBINES"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.1 Page 87"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p1=150#lb/in^2\n",
+ "p2=10#lb/in^2\n",
+ "n=10#percent\n",
+ "T=183.6+479.4#C.H.U\n",
+ "x2=0.852#C.H.U\n",
+ "H=553.9#C.H.U/lb\n",
+ "h1=T-H#C.H.U/lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "V=sqrt(2*32.2*1400*h1)#ft/sec\n",
+ "V1=sqrt(2*32.2*1400*0.9*h1)#ft/sec\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the neglecting friction is\",round(V,3),\"ft/sec\"\n",
+ "print\"the frictional drop in the nozzle is 10 recent of the total heat drop is\",round(V1,3),\"ft/sec\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the neglecting friction is 3136.312 ft/sec\n",
+ "the frictional drop in the nozzle is 10 recent of the total heat drop is 2975.367 ft/sec\n"
+ ]
+ }
+ ],
+ "prompt_number": 23
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.2 Page 88"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "v=((3140*pi*60*60)/(4*4*144))#ft/sec\n",
+ "v1=0.852*38.37#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=(v/v1)#lb\n",
+ "V=(2970*pi*60*60)/(4*4*144)#ft^3\n",
+ "W1=(V/v1)#lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the weight of steam per hour is\",round(W,3),\"lb\"\n",
+ "print\"the weight of steam per hour is\",round(W1,3),\"lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the weight of steam per hour is 471.485 lb\n",
+ "the weight of steam per hour is 445.959 lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 24
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.4 Page 89"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p1=300#lb\n",
+ "p=75#lb/in^2\n",
+ "p2=8#lb/in^2\n",
+ "h=90#C.H.U/lb\n",
+ "Pt=0.58*p#lb/in^2 absolute\n",
+ "h1=24#lb/C.H.U\n",
+ "D=0.968#C.H.U\n",
+ "D1=0.886#C.H.U\n",
+ "v=9.7#ft^3\n",
+ "v1=47.24#ft^3\n",
+ "V=sqrt(2*32.2*1400*24)#ft/sec\n",
+ "V1=sqrt(2*32.2*1400*90)#ft/sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=(p1*v*D/3600)#ft^3\n",
+ "V2=(p1*v1*D1/3600)#ft^3\n",
+ "A=0.768#in^2\n",
+ "A1=1.72#in^2\n",
+ "d=sqrt(4*0.768/pi)#in\n",
+ "d1=sqrt((4*A1)/(pi))#in\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the diameters at the throat and exit of the nozzle is\",round(d1,3),\"in\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the diameters at the throat and exit of the nozzle is 1.48 in\n"
+ ]
+ }
+ ],
+ "prompt_number": 25
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.5 Page 90"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "d=2.15#in^2\n",
+ "a=0.98#dry\n",
+ "p=100#lb/in^2\n",
+ "p1=11000#lb\n",
+ "P=0.58*p#lb/in^2\n",
+ "H=24#C.H.U/lb\n",
+ "D=0.947#lb\n",
+ "s=7.407#ft^3\n",
+ "\n",
+ "#CALCULATION\n",
+ "V=sqrt(2*32.2*1400*H)#ft/sec\n",
+ "V1=V*(d/144)#ft^3\n",
+ "T=V1/(s*D)#lb\n",
+ "A=(p1/3600)#lb\n",
+ "C=A/T#lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the coefficient of discharge for the nozzles is\",round(C,3),\"lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the coefficient of discharge for the nozzles is 0.958 lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 26
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.6 Page 91"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=9.5#lb\n",
+ "p1=120.0#lb\n",
+ "e=0.88#in\n",
+ "p2=80.0#lb/in^2\n",
+ "d=25.0#in\n",
+ "d1=0.125#in\n",
+ "t=14.0#degree C\n",
+ "T=e*19.0#C.H.U/lb\n",
+ "D=0.975#in\n",
+ "V=sqrt(2.0*32.2*1400.0*T)#ft/sec\n",
+ "S=5.467#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "V1=p*S*D#ft^3\n",
+ "T1=(V1*144.0/V)#in ^2\n",
+ "C=25.0*pi#in\n",
+ "N=C/2.5#in\n",
+ "P=C/31.0#in\n",
+ "W=d1/sin(t*pi/180)#in\n",
+ "L=P-W#in\n",
+ "W1=L*sin(t*pi/180)#in\n",
+ "T2=(T1)/(31*W1)#in\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The height of nozzle is\",round(T2,3),\"in\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The height of nozzle is 0.393 in\n"
+ ]
+ }
+ ],
+ "prompt_number": 44
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.8 Page 94"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p1=100#lb/in^2\n",
+ "p2=15#lb/in^2\n",
+ "d1=95#percent\n",
+ "d2=30#percent\n",
+ "P=0.58*p1#lb/in^2\n",
+ "H=0.95*25#C.H.U/lb\n",
+ "H1=0.95*76.5#C.H.U/lb\n",
+ "D=0.97#in\n",
+ "D1=0.905#in\n",
+ "V=7.407#ft^3\n",
+ "V1=sqrt(2*32.2*1400*H)#ft/sec\n",
+ "V2=sqrt(2*32.2*1400*H1)#ft/sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "V3=(2*pi*1*V1)/(64*4*144)#ft^3\n",
+ "W=(V3*3600)/(V*D)#lb\n",
+ "K=V2/(2*32.2)#ft lb sec\n",
+ "E=(((V2)**2*W)/(2*32.2*3600))#ft.lb\n",
+ "W1=(E*d2)/(p1*550)#ft.lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the quantity of steam used per hour and horse power developed is\",round(W1,3),\"ft lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the quantity of steam used per hour and horse power developed is 1.927 ft lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 45
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.10 Page 95"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "d=0.15#lb\n",
+ "p=20#lb/in^2\n",
+ "p1=100#lb/iN62\n",
+ "t=200#degree C\n",
+ "f=10#percent\n",
+ "Pt=0.5457*p1#lb/in^2\n",
+ "x1=0.996#in\n",
+ "x2=0.952#in\n",
+ "h=29#C.H.U/lb\n",
+ "h1=65#C.H.U/lb\n",
+ "v=7.73#ft^3\n",
+ "v1=20.12#ft^3\n",
+ "T=0.364#in\n",
+ "T1=0.465#in\n",
+ "v2=sqrt(2*32.2*1400*h)#ft/sec\n",
+ "v3=sqrt(2*32.2*1400*h1)#ft/sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "V1=d*v*x1#ft^3\n",
+ "V2=d*v1*x2#ft^3\n",
+ "A1=(V1/v2)*144#in^2\n",
+ "A2=(V2/v3)*144#in^2\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the throat area of the nozzle is\",round(A1,3),\"and exit area is\",round(A2,3),\"in^2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the throat area of the nozzle is 0.103 and exit area is 0.171 in^2\n"
+ ]
+ }
+ ],
+ "prompt_number": 51
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.11 Page 96"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "h=0.5#lb\n",
+ "p1=2.5#lb/in^2\n",
+ "p2=100#lb/in^2\n",
+ "t=250#Degree C\n",
+ "pv=1.3#constant\n",
+ "pt=0.5457*p2#lb/in^2\n",
+ "t1=18#degree C\n",
+ "h1=32#C.H.U/lb\n",
+ "h2=151#C.H.U/lb\n",
+ "D=0.887#in\n",
+ "V1=sqrt(2*32.2*1400*h1)#ft/sec\n",
+ "V2=sqrt(2*32.2*1400*h2)#ft.sec\n",
+ "s1=8.74#ft^3\n",
+ "s2=140.8#ft^3\n",
+ "T1=0.687#in\n",
+ "T1=1.77#in\n",
+ "V3=h*s1#ft^3/sec\n",
+ "V4=h*s2#ft^3/sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "A1=(V3/V1)*144#in^2\n",
+ "A2=(V4/V2)*144#in^2\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the size os nozzle is\",round(A2,3),\"in^2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the size os nozzle is 2.748 in^2\n"
+ ]
+ }
+ ],
+ "prompt_number": 52
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.13 Page 97"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "h=500.0#gallons\n",
+ "p1=150.0#lb/in^2\n",
+ "p2=0.6#lb/in^2\n",
+ "P=p2*p1#lb/in^2\n",
+ "h=25.0#C.H.U/lb\n",
+ "p=62.4#lb/ft^2\n",
+ "V=sqrt(2.0*32.2*1400.0*h)#ft/sec\n",
+ "D=0.996#in^2\n",
+ "d=4.898#in^2\n",
+ "v1=1.2#in\n",
+ "vi=163.2#ft/sec\n",
+ "m=V/32.2#ft.lb.sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=V/vi-1#lb\n",
+ "W1=(5000.0)/(3600.0*W)#ft/sec\n",
+ "V1=W1*d*D#ft^3\n",
+ "A=V1/V*144.0#in^2\n",
+ "I=(50.0/36.0+W1)#lb/sec\n",
+ "A1=(I*144.0)/(62.4*vi)#in^2\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the aera of the stream and water orifices is\",round(A1,3),\"in^2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the aera of the stream and water orifices is 0.022 in^2\n"
+ ]
+ }
+ ],
+ "prompt_number": 54
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.15 Page 100"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "a=50#degree c\n",
+ "v=2000#ft/sec\n",
+ "p=800#ft/sec\n",
+ "b=20#Degree C\n",
+ "v1=0.9#in^2\n",
+ "v2=513#ft/sec\n",
+ "W=(1/32.2)*(1810-(-313))*p#ft/lb lb stream /sec\n",
+ "K=(v**2)/(2*32.2)#ft/lb sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "D=(W/K)*100#percent/lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the work done per lb is\",round(D,3),\"%/lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the work done per lb is 84.92 %/lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 56
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.16 Page 102"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "t=65#B.Th.U per lb\n",
+ "n=0.98#dry\n",
+ "p=105#lb/in^2\n",
+ "a=14#Degree C\n",
+ "b=20#Degree C\n",
+ "p1=800#ft/sec\n",
+ "v=0.80#ft/lb\n",
+ "p2=3.5#lb/sec\n",
+ "q=1400#in\n",
+ "V=sqrt(2*32.2*778*t)#ft/sec\n",
+ "W=(p1)*(1750-b)/32.2#ft lb/lb stream/sec\n",
+ "H=(W*p2/550)#ft/lb\n",
+ "E=1/64.4*((1053)**2-(825)**2)#ft.lb steam /sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "Hd=(E/q)#C.H.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the steam as it leaves the blades and horse power is\",round(Hd,3),\"C.H.U\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the steam as it leaves the blades and horse power is 4.749 C.H.U\n"
+ ]
+ }
+ ],
+ "prompt_number": 57
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.18 Page 103"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=300#ft/sec\n",
+ "W=880#ft/sec\n",
+ "a=18#degree C\n",
+ "g=32.2#ft\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "Wd=(p*W)/g#ft lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the work done /lb steam sec is\",round(Wd,3),\"ft lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the work done /lb steam sec is 8198.758 ft lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 58
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.19 Page 104"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "a=35.0#Degree C\n",
+ "b=20.0#degree C\n",
+ "f=2.0#ft\n",
+ "w=422.0#ft\n",
+ "w1=222.0#ft\n",
+ "g=32.2#ft\n",
+ "s=1500.0#r p m\n",
+ "j=0.8#ft\n",
+ "p=3.0#lb/sec\n",
+ "h=80.0#percent\n",
+ "i=1400.0#ft\n",
+ "P=(pi*(31.0/12.0)*(s/60))#ft/sec\n",
+ "W=P/g*(w-(-w1))#ft lb\n",
+ "H=(p*W)/550#ft lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "E=W/(j*i)#C.H.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the house -power developed per pair of rings if the steam is\",round(E,3),\"ft lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the house -power developed per pair of rings if the steam is 3.623 ft lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 61
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.23 Page 106"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "d=73#ft\n",
+ "h=2#in\n",
+ "s=750#r p m\n",
+ "s1=31.3#lb/sec\n",
+ "h1=1.5#in\n",
+ "a=25#Degree c\n",
+ "p=5.7#lb/in^2\n",
+ "d1=0.97#in\n",
+ "h2=370#ft/sec\n",
+ "j=32.2#in\n",
+ "k=1400#in\n",
+ "e=0.75#percent\n",
+ "w=326#in\n",
+ "p=290#in\n",
+ "vi=155#ft/sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "P=(pi*7.69*s)/(60)#ft/sec\n",
+ "H=(P*h2*s1)/(550*j)#ft/sec\n",
+ "E=(P*h2)/(j*e*k)#C.H.U/lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the drop in pressure while the steam is passing through the turbine is\",round(E,3),\"C.H.U/lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the drop in pressure while the steam is passing through the turbine is 3.305 C.H.U/lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 63
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.25 Page 108"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=300#lb/in^2\n",
+ "ab=100#degree C\n",
+ "w=26.4#C\n",
+ "t=40#lb/in^2\n",
+ "t1=180#Degree C\n",
+ "p1=0.5#lb/in^2\n",
+ "T=732.38#C.H.U\n",
+ "W=26.2#C.H.U/lb\n",
+ "W1=102#C.H.U/lb\n",
+ "x=0.963#in\n",
+ "d=335#C.H.U/lb\n",
+ "E=743.85#C.H.U/lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=T-w#C.H.U/lb\n",
+ "h=T-W1#C.H.U/lb\n",
+ "H1=E-h#C.H.U/lb\n",
+ "T1=H+H1#C.H.U/lb\n",
+ "Wd=W1+d#C.H.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The total work done per lb steam is\",round(Wd,3),\"C.H.U\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The total work done per lb steam is 437.0 C.H.U\n"
+ ]
+ }
+ ],
+ "prompt_number": 64
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.28 Page 110"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=100#lb/in62\n",
+ "p1=0.5#lb/in^2\n",
+ "T1=659.3#C.H.U/lb\n",
+ "T2=26.2#C H U/lb\n",
+ "W=181#C H U/lb\n",
+ "H1=66#C H U/lb\n",
+ "H2=115#C H U /lb\n",
+ "D=0.912#C H U/lb\n",
+ "H3=533.4#C H U/lb\n",
+ "T3=108.5 #Degree C\n",
+ "T4=26.4#Degree C\n",
+ "W1=82.1/(D*H3)#lb\n",
+ "s=1-W1#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T=W/(T1-T2)*100#percent\n",
+ "Wd=H1+(H2*s)#C H U/lb\n",
+ "H=T1-T3#C H U//lb\n",
+ "TE=Wd/H*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the without bleeding is\",round(T,3),\"%\"\n",
+ "print\"the proper weight of steam is bled is\",round(TE,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the without bleeding is 28.589 %\n",
+ "the proper weight of steam is bled is 29.338 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 66
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/Chapter7.ipynb b/Examples_in_Thermodynamics_Problems/Chapter7.ipynb new file mode 100755 index 00000000..c07f5f58 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/Chapter7.ipynb @@ -0,0 +1,539 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:3a5f70d3bf4cb3a553b5fa0666639707f168933e92d2f8256330898382eb0c39"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 7: COMBUSTION BOILER TRIALS"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.1 Page112"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "C=86#percent\n",
+ "h=4.2#percent\n",
+ "w=20#lb\n",
+ "a=w+0.902#lb\n",
+ "C2=44/12#lb\n",
+ "N=0.77#lb\n",
+ "CO2=3.15\n",
+ "H2O=0.042*9#lb\n",
+ "N2=w*N#lb\n",
+ "Ox=a-CO2-H2O-N2#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "Co2=CO2/a*100#percent\n",
+ "H2o=H2O/a*100#percent\n",
+ "n2=N2/a#percent\n",
+ "o2=Ox/a*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The composition of the products of combutions by weight is\",round(o2,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The composition of the products of combutions by weight is 9.444 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.2 Page 113"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "g=0.05#percent\n",
+ "n=0.35#percent\n",
+ "c=0.5#percent\n",
+ "h=10#percent\n",
+ "m=167#C H U\n",
+ "h1=162#C H U\n",
+ "v=1#ft^3\n",
+ "H2=0.5#ft^3\n",
+ "Co=0.05#ft^3\n",
+ "v2=3#ft\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "G=(g*c)+(n*H2)#ft^3\n",
+ "Tv=(g*h1)+(n*m)#C H U\n",
+ "M=Tv/v2#C H U/ft^3\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The gas with twice its volume of air is\",round(M,3),\"C.H.U/ft^3\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The gas with twice its volume of air is 22.183 C.H.U/ft^3\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.4 Page 114"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "g=8#percent\n",
+ "f=88#percent\n",
+ "C=12#percent\n",
+ "w=20#lb\n",
+ "C1=11/3#lb\n",
+ "CO2=3/11#lb\n",
+ "e=0.08#lb\n",
+ "D=0.0218#lb C\n",
+ "w1=0.88#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W1=w1/D#lb lb fuel\n",
+ "T=w1/D*w#lb/hr\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the total weight of exaust gas leaving the engine per hour is\",round(T,3),\"lb/hr\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the total weight of exaust gas leaving the engine per hour is 807.339 lb/hr\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.6 Page 116"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "a=30#percent\n",
+ "b=20#percent\n",
+ "c=8#percent\n",
+ "h=42#percent\n",
+ "t1=20#degree C\n",
+ "g=0.24#in\n",
+ "t2=320#degree c\n",
+ "M=7.654#lb/lb fuel\n",
+ "A=3*M#lb/lb fuel\n",
+ "W=0.08+0.04#lb\n",
+ "T=A+0.8#lb\n",
+ "w1=0.72+0.3#lb\n",
+ "w=T-w1#lb\n",
+ "d=w*0.24*(t2-b)#C H U/lb fuel\n",
+ "H=1.02*(639+0.49*220-t1)#C H U/lb fuel\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T1=d+H#C H U/lb fuel\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"total heat carried away by flue gases is\",round(T1,3),\"C H U/lb fuel\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "total heat carried away by flue gases is 2378.76 C H U/lb fuel\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.7 Page 117"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "h=40#percent\n",
+ "g=30#percent\n",
+ "c=8#percent\n",
+ "n=10#percent\n",
+ "w=6#percent\n",
+ "g1=10#percent\n",
+ "g2=4.14#ft^3\n",
+ "Ch4=4.562#ft^3 of air\n",
+ "Co2=0.44#ft\n",
+ "H2o=1.18#ft^3\n",
+ "N2=3.7#ft63\n",
+ "x=41.4/11#ft63\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T=Ch4+x#ft^3\n",
+ "v=1+T#ft^3\n",
+ "V=x+g2#ft^3\n",
+ "D=v-V#ft^3\n",
+ "P=D/v*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the volueme of air suplied per is\",round(P,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the volueme of air suplied per is 15.248 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.9 Page 119"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "Ox=2.679#lb\n",
+ "O2=Ox-0.03#lb O2/lb fuel\n",
+ "o2=O2*100.0/23.0#lb air lb fuel\n",
+ "E=o2/2.0#lb\n",
+ "a=17.325#lb /lb fuel\n",
+ "Co2=3.294#lb\n",
+ "H2o=0.315#lb\n",
+ "N2o=13.34#lb\n",
+ "O2=23.0/100.0*E#lb\n",
+ "So2=0.005*2#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "W=Co2+N2o+O2+So2#lb /lb fuel\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the total weight of dry products is\",round(W,3),\"lb /lb fuel\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the total weight of dry products is 17.969 lb /lb fuel\n"
+ ]
+ }
+ ],
+ "prompt_number": 9
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.11 Page 120"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "l=8.7#percent\n",
+ "Co2=42.0#percent\n",
+ "N=28.0#percent\n",
+ "O2=32.0#percent\n",
+ "x=27.65#lb air\n",
+ "W=(O2/12.0)*(100.0/23.0)#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "A=x-W#lb\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the air to flues /lb carbon is\",round(A,3),\"lb\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the air to flues /lb carbon is 16.056 lb\n"
+ ]
+ }
+ ],
+ "prompt_number": 11
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.13 Page 121"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "Co=2420#C H U\n",
+ "a=3400/6#C H U\n",
+ "R=Co/3246#C H U\n",
+ "T=1+0.745#lb\n",
+ "n=1.12 #lb\n",
+ "O2=1.33/1.745#lb\n",
+ "C=O2*100/23#lb\n",
+ "CB=n/T#lb\n",
+ "m=1.74#lb\n",
+ "k=2.33#lb\n",
+ "l=1.33#lb\n",
+ "c=77#lb\n",
+ "d=23#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "Y=l*c/d#N2\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the weight of air and steam is\",round(Y,3),\"N2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the weight of air and steam is 4.453 N2\n"
+ ]
+ }
+ ],
+ "prompt_number": 12
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.15 Page 123"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "w=20#lb\n",
+ "t=320#degree C\n",
+ "t1=22#Degree C\n",
+ "w1=0.0807#lb\n",
+ "A=0.03901#AH\n",
+ "W=0.07469#AH\n",
+ "g=5.2#A\n",
+ "Q=W-A#A\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=(g*0.625)/(Q)#ft\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"weight of equal column of external air is\",round(H,3),\"ft\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "weight of equal column of external air is 91.087 ft\n"
+ ]
+ }
+ ],
+ "prompt_number": 13
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.16 Page 124"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=120#lb/in^2\n",
+ "h=30#in\n",
+ "t=48#degree C\n",
+ "C=1000#lb\n",
+ "t1=26#degree C\n",
+ "m=2.2#percent\n",
+ "g=18#lb\n",
+ "f=127#lb\n",
+ "j=33000#in\n",
+ "q=1400#in\n",
+ "L=0.978*8000#C.H.U\n",
+ "b=50#in\n",
+ "t2=320#degree C\n",
+ "g1=0.24#in\n",
+ "d=0.90#in\n",
+ "a=0.4912*30#lb/in^2\n",
+ "P=p+a#lb/in^2 abs\n",
+ "T=178.62+d*483.45#C.H.U/lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "Wt=C/f#lb\n",
+ "H=Wt*(T-t)#C.H.U\n",
+ "F=0.022*(638.9+0.48*220-t1)#C.H.U\n",
+ "G=g*0.24*(t2-t1)#C.H.U\n",
+ "E=H/L*100#percent\n",
+ "E1=b*j*60/(L*f*q)*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the heat balance for the boiler and find its efficiency and the overall efficiency of the plant is\",round(E1,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the heat balance for the boiler and find its efficiency and the overall efficiency of the plant is 7.117 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 7.17 Page 126"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "v=7950#lb C.H.U /lb\n",
+ "w=15#percent\n",
+ "c=0.85#lb\n",
+ "w1=14#percent\n",
+ "w2=9#percent\n",
+ "t1=15#degree C\n",
+ "t2=325#degree C\n",
+ "g=0.25#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=c*v#C.H.U\n",
+ "H1=0.15*(638.9+0.48*225-15)#C.H.U\n",
+ "C=c*c#lb\n",
+ "A=19.2#lb\n",
+ "Wt=A+C#lb\n",
+ "P=Wt*g*(t2-t1)#C.H.U/lb coal\n",
+ "R=0.14*H#C.H.U\n",
+ "R1=H-H1-P-R#C.H.U\n",
+ "B=R1/H*100#percent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the efficiency of a boiler is\",round(B,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the efficiency of a boiler is 61.527 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 15
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/Chapter8.ipynb b/Examples_in_Thermodynamics_Problems/Chapter8.ipynb new file mode 100755 index 00000000..417f17af --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/Chapter8.ipynb @@ -0,0 +1,502 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:7a4f335b4806f25fef0bf94c3a05e4c422a4dc1cb839f62bebcb4060f38f97a2"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 8: INTERNAL COMBUSTION ENGINES: VARIABLE SPECIFIC HEATS"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.1 Page 128"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "b=6#in\n",
+ "b1=9.0#in\n",
+ "r1=4.0#ratio\n",
+ "r2=1.0#ratio\n",
+ "p=50#lb/in^2\n",
+ "s=300#r p m\n",
+ "e=30#per cent\n",
+ "v=260#C.H.U\n",
+ "a=1.41\n",
+ "h=0.30#in\n",
+ "g=33000#in\n",
+ "g1=1400#in\n",
+ "A=1-(r2/r1)**0.41#lb/in^2\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "I=(p*pi*36/4*9/12*s/2)*1/g#ft^3\n",
+ "X=(I*g)/(g1*v*h)#ft^3\n",
+ "C=X*60/I#ft^3\n",
+ "R=h/A*100#per cent\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The fuel consumption in ft^3/h p hr and the efficiency relative to the air standard cycle is\",round(R,3),\"%\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The fuel consumption in ft^3/h p hr and the efficiency relative to the air standard cycle is 69.195 %\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.3 Page 129"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "h=200.0#r p m\n",
+ "h1=50.0#i h p\n",
+ "P4=33.4#lb/in^2\n",
+ "W=9000.0#ft lb\n",
+ "x=33000.0#ft.lb\n",
+ "p=1728.0#ft/lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "w=h1*x/100.0#ft lb\n",
+ "T=w/W#ft^3\n",
+ "V =13.0/14.0*T#ft^3\n",
+ "D=((V*p*8)/(3*pi))**(1.0/3.0)#in\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The diameter of the cylinder of a single acting and swept volume is\",round(D,3),\"in\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The diameter of the cylinder of a single acting and swept volume is 13.567 in\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.6 Page 132"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "h=12#in\n",
+ "h1=18#in\n",
+ "v=19000#B.Th.U/lb\n",
+ "T=12600#lb/in^2\n",
+ "m=90#lb/in^2\n",
+ "w=120#gal\n",
+ "t1=140#F\n",
+ "t2=60#F\n",
+ "t3=570#F\n",
+ "Cv=0.24#ft/lb\n",
+ "q=810#ft/lb\n",
+ "n=16.9#lb\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "H=(n/t2*v)#B.Th.U\n",
+ "H1=(m*pi*(144/4)*(h1/h)*(T/t2))/(778*2)#B.TH.U/min\n",
+ "H2=1750#B.Th.U\n",
+ "H3=(H1-H2)#B.Th.U\n",
+ "W=(w*10/t2)*(t1-t2)#B,Th.U\n",
+ "G=((q+n)/(t2))*(t3-t2)*Cv#B.TH.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The heat balance showing heat quantities received and the discharged per min is\",round(G,3),\"B.TH.u\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The heat balance showing heat quantities received and the discharged per min is 1686.876 B.TH.u\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.8 Page 133"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "v=12.5 #i.p.h\n",
+ "p1=8.25#in\n",
+ "p2=12.0#in\n",
+ "t=110.0#per min\n",
+ "g1=280.0#C.H.U/ft^3\n",
+ "g2=215.0#ft^3\n",
+ "V=25.0#percent\n",
+ "e=0.875#in\n",
+ "T=33000.0#in\n",
+ "v1=0.4170#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "M=(T*v)/((pi*(p1)**2.0)/(4.0)*(p2/p2)*(t))#lb.in^2\n",
+ "V1=pi*(p1)**2.0/4.0*p2/1728.0*e#ft^3\n",
+ "V2=(pi*(p1)**2.0*p2)/(4.0*4.0*1728.0)#ft^3\n",
+ "G=(g2/60.0*1.0/t)#ft^3\n",
+ "T1=G*g1#C.H.U\n",
+ "T2=(T1/v1)#C.H.U\n",
+ "F=(M/T2)#C.H.U\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The value of the Tookey factor for gas engine is\",round(F,3),\"C.H.U\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The value of the Tookey factor for gas engine is 3.207 C.H.U\n"
+ ]
+ }
+ ],
+ "prompt_number": 8
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.10 Page 135 "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p1=140#lb/in^2\n",
+ "p2=6.6#lb/in^2\n",
+ "v1=122#r.p.m\n",
+ "v2=1250#b.h.p\n",
+ "t=1425#i.h.p\n",
+ "p3=77.8#lb/in^2\n",
+ "h=0.356#lb\n",
+ "v=10000#C.H.U/lb\n",
+ "h2=2400#lb\n",
+ "q=33000#in\n",
+ "j=1400#in\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "t=(v2*q*60)/(j*h*v2*v)*100#percent\n",
+ "V=(p3*144*v1)/(q*2)#V\n",
+ "V1=(p2*144*v1)/q#V\n",
+ "T=24.16#V\n",
+ "V2=t/T#ft^3\n",
+ "I=V*V2#ft^3\n",
+ "I1=V1*V2#ft^3\n",
+ "H=24904#C/.H.U//mim\n",
+ "T=(I*q*60)/(j*h*v2*v)*100#percent\n",
+ "T1=(I1*q)/(j*H)*100#percent\n",
+ "T2=(h*v2*v)/(60)#C.H.U\n",
+ "H1=(v2*q)/(j)#C.H.U/mim\n",
+ "H2=H-(I1*q*v2)/(j*t)#C.H.U/mim\n",
+ "T3=H1+H2#C.H.U/mim\n",
+ "Tn=T2-T3#C.H.U/mim\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the overall thermal effciency is\",round(t,3),\"%\"\n",
+ "print\"the cylinder volume is\",round(V,3),\"ft^3\"\n",
+ "print\"the thermal efficiency of steam engine is\",round(T1,3),\"%\"\n",
+ "print\"total heat in oil.min is\",round(Tn,3),\"C.H.U/min\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the overall thermal effciency is 39.727 %\n",
+ "the cylinder volume is 20.709 ft^3\n",
+ "the thermal efficiency of steam engine is 0.547 %\n",
+ "total heat in oil.min is 24083.671 C.H.U/min\n"
+ ]
+ }
+ ],
+ "prompt_number": 10
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.12 Page 138"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "r=14#in\n",
+ "r1=1.8#in\n",
+ "t=30.4#lb\n",
+ "e=0.6#lb\n",
+ "lam=1.4\n",
+ "d=12#in\n",
+ "d1=18#in\n",
+ "v=10000#C.H.U/lb\n",
+ "P=200#r m p\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "A=1-(1/(lam*(r)**0.4))*((r1)**lam-1)/(r1-1)#percent\n",
+ "T=e*A#percent\n",
+ "H=t/60*v#C.H.U\n",
+ "H1=H*T#C.H.U\n",
+ "I=(H1*1400)/(33000)#ln/in^2\n",
+ "M=(I*33000)/(2*pi*144/4*d1/12*P/2)#lb/in^2\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the indicated hourse-power and the mean effiective pressure of the engine is\",round(M,3),\"lb/in^2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the indicated hourse-power and the mean effiective pressure of the engine is 75.666 lb/in^2\n"
+ ]
+ }
+ ],
+ "prompt_number": 11
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.19 Page 140"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "cv=0.1714#C.H.U\n",
+ "R=100.3#ft.lb\n",
+ "T=500#degree c\n",
+ "J=1400#in\n",
+ "Lam=R/J#C.H.U percent C\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "Cp=Lam+cv#C.H.U percent C\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The specific heat at constant volume of a gaseous mixture is\",round(Cp,3),\"C.H.U %\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The specific heat at constant volume of a gaseous mixture is 0.243 C.H.U %\n"
+ ]
+ }
+ ],
+ "prompt_number": 12
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.20 Page 141"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "a=0.124#in\n",
+ "b=0.000025#in\n",
+ "R=0.0671#heat units\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "Cp=(R+a+b)+b#T\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the specific heat of a gas at constant volume is\",round(Cp,3),\"T\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the specific heat of a gas at constant volume is 0.191 T\n"
+ ]
+ }
+ ],
+ "prompt_number": 13
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.21 Page 142"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "v=18#ft^3\n",
+ "p=14#lb/in^2\n",
+ "p1=150#lb/in^2\n",
+ "Cp=0.242#T\n",
+ "Cv=0.171#T\n",
+ "j=1400#ft\n",
+ "R=j*(Cp-Cv)#ft.lb\n",
+ "p2=144#ft\n",
+ "I1=137500#ft/lb\n",
+ "I2=6.37#ft/lb\n",
+ "v2=3.282#ft^3\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T=(p2*p*v)/R#Degree C\n",
+ "T2=(p2*p1*v2)/(R)#Degree c\n",
+ "W=Cp*(T2-T)+0.00002*((T2)**2-(T)**2)#C.H.U/lb\n",
+ "C=v/v2#ratio\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The work done the temperatures at the beginning and end of compression ratio is\",round(C,3),\"ratio\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The work done the temperatures at the beginning and end of compression ratio is 5.484 ratio\n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 8.22 Page 144"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "r=12.5#rario\n",
+ "p=0.39*10**6#ft.lb\n",
+ "p1=14#lb/in^2\n",
+ "t=373#Degree C\n",
+ "g=18#ft^3\n",
+ "t1=100#Degree C\n",
+ "V=g/r#ft^3\n",
+ "I=0.2*10**6#ft lb/lb\n",
+ "T=0.59*10**6#ft.lb/lb\n",
+ "D=0.221*10**6#ft.lb/lb\n",
+ "A=0.095*10**6#ft.lb/lb\n",
+ "E=0.264*10**6#ft.lb/lb\n",
+ "E1=0.390*10**6#t.lb/lb\n",
+ "\n",
+ "#ALCULATIONS\n",
+ "W=(E/E1)*100#ercent\n",
+ "M=(E)/(144*(g-V))#b.in^2\n",
+ "\n",
+ "#RSULTS\n",
+ "print\"the efficiency of the engine and the m e p on the assumption that the specific heats is\",round(M,3),\"lb in^2\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the efficiency of the engine and the m e p on the assumption that the specific heats is 110.709 lb in^2\n"
+ ]
+ }
+ ],
+ "prompt_number": 16
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/Chapter9.ipynb b/Examples_in_Thermodynamics_Problems/Chapter9.ipynb new file mode 100755 index 00000000..320477a8 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/Chapter9.ipynb @@ -0,0 +1,227 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:023bb2358c5a7a6740b1b9a08290e11f954e790555a353817bbe84615e4679e6"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 9: VALVE DIAGRAMS AND VALVE GEARS"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 9.5 Page 150"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=20#in\n",
+ "l=100#in\n",
+ "r=120#r.p.m\n",
+ "v=3.5#in\n",
+ "l2=1#in\n",
+ "l3=1/8#in\n",
+ "v1=1.44#omega in/sec\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "V=p*(1.06/1.166)#omega in./sec\n",
+ "R=(V/v1)#omega in/sec\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"The ratio of velocity of the piston to the velocity is\",round(R,3),\"Omega in/sec\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The ratio of velocity of the piston to the velocity is 12.626 Omega in/sec\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 9.7 Page 151"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "v=0.6#in\n",
+ "m=1.0#in\n",
+ "t=0.75#in\n",
+ "p=4#in\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "D=t/m#in\n",
+ "A=(p*m/D)#in\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the travel and laps of the value is\",round(A,3),\"in\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the travel and laps of the value is 5.333 in\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 9.10 Page 154"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "l=1.5#in\n",
+ "p=4.0#in\n",
+ "v=0.98#in\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "T=(l*p/v)#in\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the particulars of a value and it eccentric is\",round(T,3),\"in\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the particulars of a value and it eccentric is 6.122 in\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 9.12 Page 156"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "p=1/10#in\n",
+ "v1=3/4#in\n",
+ "v2=3/5#in\n",
+ "m=1*1/2#in\n",
+ "l=4#cranks\n",
+ "a1=1.25#in\n",
+ "a2=0.7#in\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "C=a1/a2#in\n",
+ "A=l*a1/a2#in\n",
+ "S=(A/2-a1)#in\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the travel of the value is\",round(S,3),\"in\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the travel of the value is 2.321 in\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 9.17 Page 161"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#initialisation of variable\n",
+ "from math import pi,sqrt,acos,asin,atan,cos,sin,tan\n",
+ "v=3*1/2#in\n",
+ "a=30#degree\n",
+ "l=0.8#in\n",
+ "v1=0.2#in\n",
+ "L=0.13#in\n",
+ "m=1.075#in\n",
+ "d=0.58#in\n",
+ "p=1.875#in\n",
+ "\n",
+ "#CALCULATIONS\n",
+ "V=(p-d)#in\n",
+ "P=V+1.25#in\n",
+ "\n",
+ "#RESULTS\n",
+ "print\"the main value and the maximum opening to steam is\",round(P,3),\"in\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "the main value and the maximum opening to steam is 2.545 in\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/README.txt b/Examples_in_Thermodynamics_Problems/README.txt new file mode 100755 index 00000000..b203a218 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/README.txt @@ -0,0 +1,10 @@ +Contributed By: Tarun Kumar Das +Course: btech +College/Institute/Organization: College of Engineering +Department/Designation: Industrial Engineering +Book Title: Examples in Thermodynamics Problems +Author: W. R. Crawford +Publisher: Sir Isaac Pitman, London +Year of publication: 1937 +Isbn: 978-0471862567 +Edition: 3
\ No newline at end of file diff --git a/Examples_in_Thermodynamics_Problems/screenshots/Heat_supplied.png b/Examples_in_Thermodynamics_Problems/screenshots/Heat_supplied.png Binary files differnew file mode 100755 index 00000000..da1352c8 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/screenshots/Heat_supplied.png diff --git a/Examples_in_Thermodynamics_Problems/screenshots/Ideal_efficiency.png b/Examples_in_Thermodynamics_Problems/screenshots/Ideal_efficiency.png Binary files differnew file mode 100755 index 00000000..09dc192c --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/screenshots/Ideal_efficiency.png diff --git a/Examples_in_Thermodynamics_Problems/screenshots/Thermal_capacity.png b/Examples_in_Thermodynamics_Problems/screenshots/Thermal_capacity.png Binary files differnew file mode 100755 index 00000000..75f29503 --- /dev/null +++ b/Examples_in_Thermodynamics_Problems/screenshots/Thermal_capacity.png |
