initial commit / add all books

1 files changed, 31 insertions, 0 deletions

diff --git a/3733/CH18/EX18.3/Ex18_3.sce b/3733/CH18/EX18.3/Ex18_3.sce
new file mode 100644
index 000000000..6b80cccd7
--- /dev/null
+++ b/3733/CH18/EX18.3/Ex18_3.sce
@@ -0,0 +1,31 @@
+// Example 18_3

+clc;funcprot(0);

+//Given data 

+T_1=45;//The temperature of water at inlet in °C

+m_w1=360;// kg/min

+V=10;// The air circulated in the tower in m^3/sec

+Q=4900;// The amount of heat absorbed by the air in watts

+DBT=20;// Dry bulb temperature in °C

+RH=60;// Relative humidity in %

+T_a2=26;// The temperature of air leaves the tower at saturated condition in °C

+C_pw=4.2;// kJ/kg°C

+

+//Calculation

+// The condtions of air entering and leaving the tower are represented onn psychrometric chart as shown in Fig.Prob.18.3

+// From psychrometric chart,

+H_a1=45;// kJ/kg

+H_a2=81;// kJ/kg

+w_1=9.6;// grams/kg

+w_2=21.6;// grams/kg

+v_s1=0.848;// m^3/kg

+m_a=V/v_s1;// kg/sec

+Q=Q/1000;// kW=kJ/sec

+//T_2=y(1)

+function[X]=Temperature(y);

+    X(1)=(((m_w1*C_pw*(T_1-y(1)))/60)+Q)-(m_a*((H_a2-H_a1)-(((w_2-w_1)/1000)*C_pw*y(1))));

+endfunction

+y=[10]

+z=fsolve(y,Temperature);

+T_2=z(1);// The temperature of water coming out of the tower in °C

+m_m=m_a*((w_2-w_1)/1000);// The make up water required per hour in kg/min

+printf('\nThe temperature of water coming out of the tower=%0.0f°C \nThe make up water required per hour=%0.3f kg/min',T_2,m_m);