diff options
Diffstat (limited to '617/CH16')
| -rwxr-xr-x | 617/CH16/EX12.3/Example16_3.sci | 17 | ||||
| -rwxr-xr-x | 617/CH16/EX16.1/Example16_1.sci | 13 | ||||
| -rwxr-xr-x | 617/CH16/EX16.2/Example16_2.sci | 13 | ||||
| -rwxr-xr-x | 617/CH16/EX16.4/Example16_4.sci | 12 | ||||
| -rwxr-xr-x | 617/CH16/EX16.5/Example16_5.sci | 26 | ||||
| -rwxr-xr-x | 617/CH16/EX16.6/Example16_6.sci | 32 | ||||
| -rwxr-xr-x | 617/CH16/EX16.7/Example16_7.sci | 34 | ||||
| -rwxr-xr-x | 617/CH16/EX16.8/Example16_8.sci | 18 | ||||
| -rwxr-xr-x | 617/CH16/EX16.9/Example16_9.sci | 18 |
9 files changed, 183 insertions, 0 deletions
diff --git a/617/CH16/EX12.3/Example16_3.sci b/617/CH16/EX12.3/Example16_3.sci new file mode 100755 index 000000000..2e61d6961 --- /dev/null +++ b/617/CH16/EX12.3/Example16_3.sci @@ -0,0 +1,17 @@ +clc();
+clear;
+
+// To compute the ammonia diffusing through the stagnant air
+
+x=0.1/12; // thickness of still air layer in ft
+T=77+460; // temperature in degR
+p=1; // Atmospheric pressure in atm
+pa1=0.3; // Pressure of ammonia in still air in atm
+pb1=p-pa1; // pressure of air in atm
+pa2=0; // pressure of ammonia in the absorption plane
+pb2=p-pa2; // pressure of air in absorption plane
+pbm=(pb2-pb1)/(log(pb2/pb1)); // Logarithmic mean pressure
+D=0.914; // Diffusion coefficient for ammonia
+R=0.729; // Gas constant in ft^3-atm/lb-mole-degR
+N=D*p*(pa1-pa2)/(R*T*x*pbm);
+printf("The amount of ammonia diffusing through the stagnant air is %.1f lb-mol/hr-ft^2",N);
\ No newline at end of file diff --git a/617/CH16/EX16.1/Example16_1.sci b/617/CH16/EX16.1/Example16_1.sci new file mode 100755 index 000000000..bbcc250b8 --- /dev/null +++ b/617/CH16/EX16.1/Example16_1.sci @@ -0,0 +1,13 @@ +clc();
+clear;
+
+// To compute the diffusion coeffiient for water vapour in air
+
+T=25+273; // Temperature in degK
+p=1; // Pressure in atm
+Va=18.9; // Molecular volume of water vapour in cm^3/gm-mol
+Vb=29.9; // Molecular volume of air in cm^3/gm-mol
+Ma=18; // Molecular weight of water vapour in gm/mol
+Mb=29; // Molecular weight of air in gm/mol
+Dab=0.0043*(T^1.5)*sqrt((1/Ma)+(1/Mb))/(p*(Va^(1/3)+Vb^(1/3))^2);
+printf("The diffusion coefficient is %.3f cm^3/sec ",Dab);
\ No newline at end of file diff --git a/617/CH16/EX16.2/Example16_2.sci b/617/CH16/EX16.2/Example16_2.sci new file mode 100755 index 000000000..32375dc82 --- /dev/null +++ b/617/CH16/EX16.2/Example16_2.sci @@ -0,0 +1,13 @@ +clc();
+clear;
+
+// To compute the diffusion coeffiient for benzene in air
+
+T=25+273; // Temperature in degK
+p=1; // Pressure in atm
+Va=96; // Molecular volume of benzene in cm^3/gm-mol
+Vb=29.9; // Molecular volume of air in cm^3/gm-mol
+Ma=78; // Molecular weight of benzene in gm/mol
+Mb=29; // Molecular weight of air in gm/mol
+Dab=0.0043*(T^1.5)*sqrt((1/Ma)+(1/Mb))/(p*(Va^(1/3)+Vb^(1/3))^2);
+printf("The diffusion coefficient is %.3f cm^3/sec ",Dab);
\ No newline at end of file diff --git a/617/CH16/EX16.4/Example16_4.sci b/617/CH16/EX16.4/Example16_4.sci new file mode 100755 index 000000000..4ddfede85 --- /dev/null +++ b/617/CH16/EX16.4/Example16_4.sci @@ -0,0 +1,12 @@ +clc();
+clear;
+
+// To compute the hydrogen loss per unit pipe by diffusion
+
+ri=3/96; // Inner radius of pipe in ft
+ro=1/24; // Outer radius of pipe in ft
+Ca1=0.0003; // Concentration at the inner hose of pipe in lb-mol/ft^2
+Ca2=0; // Concentration at the outer surface
+D=0.7*10^-5; // Diffusion coefficient of hydrogen in rubber in ft^2/hr
+N=2*%pi*D*(Ca1-Ca2)/log(ro/ri); // Rate of diffusion in lb-mol/hr
+printf("The rate of diffusion iof hydrogen in rubber is %.2f*10^-8 lb-mole/hr",N*10^8);
\ No newline at end of file diff --git a/617/CH16/EX16.5/Example16_5.sci b/617/CH16/EX16.5/Example16_5.sci new file mode 100755 index 000000000..62eb97ce7 --- /dev/null +++ b/617/CH16/EX16.5/Example16_5.sci @@ -0,0 +1,26 @@ +clc();
+clear;
+
+// To calculate the amount of water evaporated per hour per square foot of surface area
+
+u=0.0437; // Viscosity in lb/hr-ft
+rho=0.077; // Density in lb-ft^2
+D=0.992; // Diameter of pipe in ft
+v=4*3600; // Velocity in ft/sec
+L=6/12; // Length of pipe parallel to direction of air flow in ft
+p=14.7; // Atmospheric pressure in psi
+T=460+65; // Temperature in degR
+
+// Heat transfer equation for laminar flow of a flat surface
+Nre=L*v*rho/u; // Reynolds number
+Ns=u/(rho*D); // Schimdt mumber
+Nnu=0.662*(Ns)^(1/3)*sqrt(Nre); // Nusselt number
+hmc=Nnu*D/L; // Heat transfer coefficient
+pv1=0.144; // Vapour pressure at 40% humidity
+pv2=0.252; // Vapour pressure at saturation
+pa1=p-pv1; // Absolute pressure of air at 40% rel. humidity in psi
+pa2=p-pv2; // Absolute pressure of saturated air in psi
+pbm=(pa1+pa2)/2; // Log mean pressure in psi
+R=1544; // Universal gas constant in ft^3-psi/lbmol-degR
+N=hmc*p*(pa1-pa2)*144/(R*T*pbm);
+printf("The amount of water evaporated per hour is %.4f lb mol/hr-ft^2",N);
\ No newline at end of file diff --git a/617/CH16/EX16.6/Example16_6.sci b/617/CH16/EX16.6/Example16_6.sci new file mode 100755 index 000000000..d814b3803 --- /dev/null +++ b/617/CH16/EX16.6/Example16_6.sci @@ -0,0 +1,32 @@ +clc();
+clear;
+
+// To estimate the amount of water transferred
+
+u=0.047; // Viscosity in lb/hr-ft
+rho=0.069; // Density in lb-ft^2
+D=0.992; // Diameter of pipe in ft
+v=7.5*3600; // Velocity in ft/sec
+L=2; // Length of pipe parallel to direction of air flow in ft
+M=0.992; // Molecular weight
+p=14.696; // Atmospheric pressure in psi
+T=460+65; // Temperature in degR
+M=29; // molecular weight of air
+M2=18; // Molecular weight of water vapour
+A=4; // Area of water surface in ft^2
+// Heat transfer equation for laminar flow of a flat surface
+Nre=L*v*rho/u; // Reynolds number
+
+// Assuming the case that of a fluid flowing parallel to a flat plate , jm=0.0039
+jm=0.0039;
+Ns=u/(rho*D); // Schimdt mumber
+Gm=v*rho/M; // Mole flow rate
+pv1=0.672; // Vapour pressure at 40% humidity
+pv2=0.600; // Vapour pressure at saturation
+pa1=p-pv1; // Absolute pressure of air at 40% rel. humidity in psi
+pa2=p-pv2; // Absolute pressure of saturated air in psi
+pbm=(pa1+pa2)/2; // Log mean pressure in psi
+hmp=jm*Gm/(pbm*144*Ns^(2/3)); // Heat transfer coefficient in lbmol/ft^2-hr-psi
+N=hmp*(pv1-pv2)*144; // Mass transfer rate in lb mol/hr-ft^2
+W=N*A*M2;
+printf("The amount of water evaporated per hour is %.3f lb mol/hr-ft^2",W);
\ No newline at end of file diff --git a/617/CH16/EX16.7/Example16_7.sci b/617/CH16/EX16.7/Example16_7.sci new file mode 100755 index 000000000..e8b91d617 --- /dev/null +++ b/617/CH16/EX16.7/Example16_7.sci @@ -0,0 +1,34 @@ +clc();
+clear;
+
+// To calculate the amount of water evaporated in per hour for a square foot of water surface
+
+u=3.82*10^-7; // Viscosity in lb-sec/ft^2
+rho=2.3*10^-3; // Density in lbsec^2/ft^4
+A=1; // Area in ft^2
+Cp=0.24; // Specific heat capacity in abtu/lbm-degF
+v=4*3600; // Velocity in ft/sec
+k=0.015; // Thermal conductivity in Btu/hr-ft-degF
+p=14.7; // Atmospheric pressure in psi
+M=29; // Avg. molecular weight of air
+T1=70+460; // Temperature of still air in degF
+T2=90+460; // temperature of surface of water in degF
+L=1; // For characteristic of 1 ft
+D=0.992; // Diffusivity in ft^2/sec
+
+// Heat transfer equation for laminar flow of a flat surface
+Ngr=32.2*L^3*((T2/T1)-1)/(u/rho)^2; // Grasshops number
+Npr=u*3600*Cp*32.2/k; // Prandtls number
+Nnu=0.75*(Ngr*Npr)^.25; // Nusselt number
+h=Nnu*k/L; // Heat transfer coefficient
+Ns=u*3600/(rho*D); // Schimdt mumber
+hmc=h*D*(Ns/Npr)^0.25/k; // Heat transfer coe
+pv1=0.18; // Vapour pressure at 40% humidity
+pv2=0.69; // Vapour pressure at saturation
+pa1=p-pv1; // Absolute pressure of air at 40% rel. humidity in psi
+pa2=p-pv2; // Absolute pressure of saturated air in psi
+pbm=(pa1+pa2)/2; // Log mean pressure in psi
+R=1544; // Universal gas constant in ft^3-psi/lbmol-degR
+T=(T1+T2)/2; // Average temperature in degR
+N=hmc*p*(pv2-pv1)*144/(R*T*pbm)*18; // mass transfer rate in lbmol/hr-ft^2
+printf("The amount of water evaporated per hour is %.4f lb mol/hr-ft^2",N);
diff --git a/617/CH16/EX16.8/Example16_8.sci b/617/CH16/EX16.8/Example16_8.sci new file mode 100755 index 000000000..26928ec34 --- /dev/null +++ b/617/CH16/EX16.8/Example16_8.sci @@ -0,0 +1,18 @@ +clc();
+clear;
+
+// To know the moisture content of air
+
+Td=70+460; // Dry bulb temperature in degR
+Tw=60+460; // Wet bulb temperature in degR
+a=0.26; // Ratio of coefficients ie. h/hmw from table
+L=1059.9; // Latent heat Btu/lbmol
+p=14.7; // Atmospheric pressure in psi
+pa=0.259; // Partial pressure of water in psi
+Ma=18; // Molecular weight of water vapour
+Mb=29; // Molecular weight of air
+
+Wwb=pa*Ma/(Mb*(p-pa)); // Absolte dry bulb humidity of air
+Wdb=Wwb-(a*(Td-Tw)/L); // Absolte dry bulb humidity of air
+printf("The humidity of air at dry conditions is %.5f lbm/lbm of dry air",Wdb);
+
\ No newline at end of file diff --git a/617/CH16/EX16.9/Example16_9.sci b/617/CH16/EX16.9/Example16_9.sci new file mode 100755 index 000000000..f75a49d59 --- /dev/null +++ b/617/CH16/EX16.9/Example16_9.sci @@ -0,0 +1,18 @@ +clc();
+clear;
+
+// To estimate the mass transfer coefficient
+
+v=20; // Velocity of air ammonia mixture in ft/sec
+Npr=0.72; // Prandtls number
+Ns=0.60; // Schimdt number
+pbm=14.7; // log mean pressure in psi
+Mm=29; // Molecular weight of mixture
+Mv=17; // Molecular weight of ammonia
+Ma=29; // Molecular weight of air
+Cp=0.24; // specific heat capacity in Btu/lbm-degF
+h=8; // Heat transfer coefficient
+p=1; // Atospheric pressure in atm
+
+hmp=h*Mv*(Npr/Ns)^(2/3)/(Cp*p*Ma); // Mass transfer coefficient based on pressure
+printf("The mass transfer coefficient based on pressure is %.1f lbm/hr-ft^2-atm",hmp);
\ No newline at end of file |