initial commit / add all books

9 files changed, 291 insertions, 0 deletions

diff --git a/2774/CH8/EX8.1/Ex8_1.sce b/2774/CH8/EX8.1/Ex8_1.sce
new file mode 100755
index 000000000..f5e50f0fe
--- /dev/null
+++ b/2774/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,40 @@
+clc

+// solution

+//initialization of variables

+Ra=0.287 // specific gas constant for air

+P=100 // pressure of room in kPa

+V=150 // volume of room in m^3

+T=25+273 // temperature of air in kelvin

+phi=0.6 // relative humidity

+Pg=3.29 // saturation vapour pressure in kPa at 25 *C from table C.1

+Mv= 18 // molecular mass of water vapor

+Ma=28.97 // molecular mass of air

+

+Pv=Pg*phi // partial pressure of water vapour

+

+Pa=P-Pv // partial pressure of air

+

+w=0.622*(Pv/Pa) // humidity ratio in Kg of water/ Kg of dry air

+Tdp=17.4 // dew point temperature from interpolation in table C.2 corresponding to partial pressure Pv=1.98 kPa

+

+ma=Pa*V/(Ra*T) // mass of air

+mv=w*ma // mass of water vapour in  kg

+

+// now we find volume percentage

+Nv=mv/Mv // moles of vapour

+Na=ma/Ma // moles of air

+

+Vw= Nv/(Na+Nv) // fraction of volume occupied by water vapour

+

+printf(" The humidity ratio is %0.3f Kg water/ kg of dry air \n",w)

+printf("The dew point is %0.1f degree celsius \n ",Tdp)

+printf("The mass of water vapour in the air is %0.2f kg \n",mv)

+printf("The volume percentage of the room that is water vapor is %0.2f %%",Vw*100)

+// The answers are correct within given limits 

+// The variation in answers is due to approximations made by

+// textbook while scilab is precise

+

+

+

+

+

diff --git a/2774/CH8/EX8.10/Ex8_10.sce b/2774/CH8/EX8.10/Ex8_10.sce
new file mode 100755
index 000000000..3cdbedaa6
--- /dev/null
+++ b/2774/CH8/EX8.10/Ex8_10.sce
@@ -0,0 +1,39 @@
+clc

+//solution

+// initialization of variables

+mdotw3=10000 // mass flow rate of water entering in cooling tower in kg/min

+Tw1=40+273 // temperature of water entering cooling tower in kelvin

+Ta1=20+273 // temperature of air entering cooling tower in kelvin

+phi1=0.5// relative humidity of entering air

+Tw2=25+273 // temperature of water leaving cooling tower in kelvin

+Ta2=32+273 // temperature of air leaving cooling tower in kelvin

+phi2=0.98 // relative humidity of leaving air

+// from psychrometric chart

+h1=37// specific enthalpy of air @ 20*C DBT and 50% humidity

+w1=0.0073 // specific humidity of air @ 20*C DBT and 50% humidity

+h2=110// specific enthalpy of air @ 32*C DBT and 98% humidity

+w2=0.030 // specific humidity of air @ 32*C DBT and 98% humidity

+

+h3=167.5 // specific enthalpy of water from steam table at 40 degree celsius

+h4=104.9 // specific enthalpy of water from steam table at 25 degree celsius

+

+mdota=(mdotw3*(h4-h3))/(h1-h2+(w2-w1)*h4) // by energy balance

+

+

+v1=0.84 // specific volume of air entering tower from psychrometric chart

+

+Qf=mdota*v1 // volume flow rate in m^3/min

+printf("The volume flow rate of air into the cooling tower is %i m^3/min \n",Qf)

+

+mdot4=mdotw3-(w2-w1)*mdota // by mass balance

+printf("The mass flow rate of water that leaves the cooling tower is %i kg/min",mdot4) 

+// The answers is slightly different in textbook due to approximations in calculations while in scilab solution is precise

+

+

+

+

+

+

+

+

+

diff --git a/2774/CH8/EX8.2/Ex8_2.sce b/2774/CH8/EX8.2/Ex8_2.sce
new file mode 100755
index 000000000..ee749290d
--- /dev/null
+++ b/2774/CH8/EX8.2/Ex8_2.sce
@@ -0,0 +1,23 @@
+clc

+clear

+// solution

+//initialization of variables

+Ra=0.287 // specific gas constant for air

+P=100 // pressure of room in kPa

+w1=0.0126 // old humidity ratio of example 8.1-

+Pg=3.29 // saturation vapour pressure in kPa at 25 *C from table C.1

+mv=2.17 // initial mass of water vapour in example 8.1

+T=25+273 // temperature after reheat

+V=150 // volume of room in m^3 

+Pv=1.228 // saturation vapour pressure in kPa @ 10 degree celsius from table C.1

+Pa=P-Pv // partial pressure of air

+w2=0.622*(Pv/Pa) // new humidity ratio in Kg of water/ Kg of dry air

+deltaw=w1-w2 // difference in humidity ratio

+ma=Pa*V/(Ra*T) // mass of air

+deltamv=deltaw*ma // mass of water vapour condensed

+X=deltamv*100/mv // percentage of water vapour condensed

+printf("The percentage that condenses is %0.1f %% \n",X)

+// AFTER REHEATING

+phi=1.608*w2*Pa/Pg

+printf("The relative humidity is %0.1f %%",phi*100)

+

diff --git a/2774/CH8/EX8.3/Ex8_3.sce b/2774/CH8/EX8.3/Ex8_3.sce
new file mode 100755
index 000000000..422105566
--- /dev/null
+++ b/2774/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,24 @@
+clc

+// solution

+//initialization of variables

+T1=40 // dry bulb temperature in degree celsius

+T2=20 // wet bulb temperature in degree celsius

+Cp=1.0 // specific heat

+P=100 // pressure of air stream in kPa

+pg1=7.383 //saturation pressure @ 40 degree celsius

+hfg2=2454 // latent heat for 20 degree celsius

+Pg2=2.338 // saturation pressure @ 20 degree celsius

+w2=0.622*Pg2/(P-Pg2) // specific humidity for wet bulb condition

+hg1=2574 // specific enthalpy of saturated vapour @ 40 degree celsius

+hf2=83.9 //spedific enthalpy of saturated liquid @ 20 degree celsius

+w1=((w2*hfg2)+Cp*(T2-T1))/(hg1-hf2)// specific humidity for 40 degree celsius

+printf("The humidity ratio is %0.5f kg water/ Kg dry air \n",w1)

+pv1=100*w1/(0.622+w1) // partial pressure of vapour

+phi=pv1/pg1 // relative humidity

+printf("The relative humidity is %0.1f %% \n",phi*100)

+

+hv=hg1 // temperature is at DBT=40 degree celsius

+h=Cp*T1+w1*hv // specific enthalpy  of air

+printf("The specific enthalpy is %0.1f kJ/kg dry air",h)

+

+

diff --git a/2774/CH8/EX8.5/Ex8_5.sce b/2774/CH8/EX8.5/Ex8_5.sce
new file mode 100755
index 000000000..caf93c1b9
--- /dev/null
+++ b/2774/CH8/EX8.5/Ex8_5.sce
@@ -0,0 +1,18 @@
+clc

+//solution

+// initialization of variables

+

+T1=40 // inlet temperature in degree celsius

+T2=27 // outlet temperature in degree celsius

+phi1= 10 // relative humidity at inlet

+// as no heat transfer takes place thus isenthalpic process

+//Thus following the enthalpy line at DBT=40 and Relative humidity=10

+phi2=45 // by interpolation of constant enthalpy line

+w1=0.0046// specific humidity @ T=40 and phi1=10

+w2=0.010 // specific humidity at outlet

+W=w2-w1 // amount of water added

+Tmin=18.5 // minimum temperature at 100% relative humidity

+printf("The relative humidity is %i %% \n ",phi2)

+printf("The added water is %0.04f kg water/kg dry air \n",W)

+printf("The lowest possible temperature is %0.1f *C ",Tmin)

+

diff --git a/2774/CH8/EX8.6/Ex8_6.sce b/2774/CH8/EX8.6/Ex8_6.sce
new file mode 100755
index 000000000..e53620f5f
--- /dev/null
+++ b/2774/CH8/EX8.6/Ex8_6.sce
@@ -0,0 +1,32 @@
+clc

+//solution

+// initialization of variables

+T1=5+273 // outside air temperature in kelvin

+P=100 // pressure in kPa

+Ra=0.287 // specific gas constant for air

+phi=0.7 // relative humidity outside

+Qf=50/60 // volume flow rate in m^3/sec

+Pg1=0.872 // saturation pressure at 278 K

+Pv1=phi*Pg1 // partial pressure of water vapour

+Pa1=P-Pv1 // partial pressure of air

+

+rhoa=Pa1/(Ra*T1) // density of dry air

+

+mdota= Qf*rhoa // mass flow rate of dry air

+

+// using psychrometric chart at T1=5*C and phi1=70% 

+h1=14 // inlet enthalpy in kJ/kg

+h2=35 // enthalpy after heating in kJ/kg

+

+Qdot=mdota*(h2-h1) // heat transfer rate

+// from psychrometric chart for T=25 *C and 35 kJ/kg enthalpy

+phi2=19 // realtive humidity

+printf("The heat transfer rate is %0.1f kJ/s \n",Qdot)

+printf("The final relative humidity is %i %% ",phi2)

+

+

+

+

+

+

+

diff --git a/2774/CH8/EX8.7/Ex8_7.sce b/2774/CH8/EX8.7/Ex8_7.sce
new file mode 100755
index 000000000..a718c657a
--- /dev/null
+++ b/2774/CH8/EX8.7/Ex8_7.sce
@@ -0,0 +1,36 @@
+clc

+//solution

+// initialization of variables

+//DATA TAKEN FROM PSYCHROMETRIC CHART

+T1=5+273 // outside temperature in kelvin

+h1=10// enthalpy in kJ/kg @ T=5 *C and 40 % relative humidity

+Pg1=0.872 // saturaion pressure in kPa for 5 degree celsius DBT

+phi1=0.4

+h2=33 // specific enthalpy at 25 *C and 40 % relatuve humidity

+h3=45 // specific enthalpy at state 3

+P=100 // atmospheric pressure in kPa

+Ra=0.287 // specific gas constant for air

+Qf=60/60 // volume flow rate in m^3/s

+Pv1=phi1*Pg1 // partial presure of water vapour 

+Pa1=P-Pv1 // partial pressure of air

+w2=0.0021 // specific humidity @ 40 % relative humidity and 25*C temperature

+w3=0.008 // final specific humidity

+rhoa1=Pa1/(Ra*T1) // air density

+mdota=Qf*rhoa1 // mass flow rate of dry air

+

+Qdot=mdota*(h2-h1) // heat transfer rate

+

+// as the process is isothermal thus

+mdots=mdota*(w3-w2)// mass flow rate of steam by conservation of mass

+printf("the rate of steam supplied is %0.4f kg/s \n",mdots)

+// also using energy balance

+hs=(mdota*(h3-h2))/mdots // enthalpy of steam

+hf=604.7 // enthalpy of saturated liquid @ 400 kPa

+hg=2738.5 // enthalpy of saturated vapour @ 400 kPa

+xs=(hs-hf)/(hg-hf)

+printf("The quality of steam is %0.2f ",xs)

+

+

+

+

+

diff --git a/2774/CH8/EX8.8/Ex8_8.sce b/2774/CH8/EX8.8/Ex8_8.sce
new file mode 100755
index 000000000..a672683ba
--- /dev/null
+++ b/2774/CH8/EX8.8/Ex8_8.sce
@@ -0,0 +1,33 @@
+clc

+//solution

+// initialization of variables

+// REFER TO FIG. 8.4

+T1=30 // outside temperature in degree celsius

+phi1=0.9 // outside relative humidity

+T2=23 // room temperature in degree celsius

+phi2=0.4 // relative humidity in room

+

+// using psychrometric chart

+w1=0.0245 // specific humidity @ 30 *C and relative humidity 0.9

+h1=93 // specific enthalpy @ 30 *C and relative humidity 0.9

+w2=w1 // during cooling humidity remains constant 

+w3=0.007 // specific humidity @ 23 *C and relative humidity 0.4

+h4=41 // final specific enthalpy

+h3=26 // specific enthalpy @ 23 *C and relative humidity 0.4

+deltaw=w3-w2 // moisture removed

+printf("the amount of moisture removed is %0.4f kg \n",deltaw)

+

+

+qout=h3-h1 // heat removed F-G-H process

+

+printf(" the heat removed is %i kJ/kg \n ",qout)

+

+qin=h4-h3 // heat added to bring to desired state

+

+printf(" the heat added is %i kJ/kg ",qin)

+

+

+

+

+

+

diff --git a/2774/CH8/EX8.9/Ex8_9.sce b/2774/CH8/EX8.9/Ex8_9.sce
new file mode 100755
index 000000000..ea67971dd
--- /dev/null
+++ b/2774/CH8/EX8.9/Ex8_9.sce
@@ -0,0 +1,46 @@
+clc

+//solution

+// initialization of variables

+P=100 // atospheric pressure in kPa

+R=0.287 // specific gas constant for air

+T1=15+273 // outside temperature in kelvin

+phi1=0.4// outside air relative humidity

+Qf1=40 // outside air flow rate in m^3/min

+T2=32+273 // inside temperature in kelvin

+phi2=0.7 // inside air relative humidity

+Qf2=20 // outside air flow rate in m^3/min

+Ps1=1.7 // saturation pressure @ 15 degree celsius and 40% humidity

+Ps2=4.9 // saturation pressure @ 32 degree celsius and 70% humidity

+

+Pv1=Ps1*phi1 // partial pressure of water vapour outside

+

+Pv2=Ps2*phi2 // partial pressure of water vapour inside

+

+Pa1=P-Pv1 //partial pressure of dry air outside

+Pa2=P-Pv2 //partial pressure of dry air inside

+

+rhoa1=Pa1/(R*T1) // density of outside air

+mdota1=Qf1*rhoa1 // mass flow rate of air outside

+

+rhoa2=Pa2/(R*T2) // density of inside air

+mdota2=Qf2*rhoa2 // mass flow rate of inside air

+ // using psychrometric chart locating state 1 and 2

+h1=37 // specific enthalpy @ DBT 15*C and 40 % humidity

+w1=0.0073 // specific humidity @ DBT 15*C and 40 % humidity

+h2=110 // specific enthalpy @ DBT 32*C and 70 % humidity

+w2=0.0302 // specific humidity @ DBT 32*C and 70 % humidity

+ratio=mdota1/mdota2 // ratio of distance between states 

+//  using this ratio state 3 is located on psychrometric chart

+T3=(mdota1*T1+mdota2*T2)/(mdota1+mdota2)-273 // final temparature in celsius

+

+phi3=65// final relative humidity at T3 from psychrometric chart

+

+printf("The relative humidity is %i %% \n",phi3)

+printf(" The resultant temperature is %i degree celsius",T3)

+

+

+

+

+

+

+