blob: d09bc6010803f366f41ec674e053372190fa1b33 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
clc;
p1=3; // Pressure of fluid at inlet in bar
T1=150; // Temperature of fluid at inlet in degree celcius
V1=90; // Velocity of fluid at inlet in m/s
eff_nozzle=0.85; // Nozzle efficiency
k=1.4; // Index of reversible adiabatic process
p2=1/3*p1;
// (a).Steam
// Following are taken from steam table
h1=2761; // specific enthalpy in kJ/kg
s1=7.0778;// specific entropy in kJ/kg K
s2s=s1; // Isentropic process
sf2s=1.3026; sfg2s=6.0568;// specific entropy in kJ/kg K
hf2=417.46; hfg2=2258; // specific enthalpy in kJ/kg
x2s=(s2s-sf2s)/sfg2s; // Quality of steam
h2s=hf2+x2s*hfg2;
V2s=sqrt (2000*(h1-h2s)+V1^2); // Isentropic Velocity
V2=sqrt (eff_nozzle) *V2s; // Actual nozzle exit velocity
disp ("m/s (round off error)",V2," Actual nozzle exit velocity = ","(a).Steam");
// (b).Air
Cpo=1.0035; // Specific heat at constant pressure in kJ/kg K
T2s=(T1+273)*(p2/p1)^((k-1)/k); // Isentropic temperature
V2s=sqrt ((2000*Cpo*((T1+273)-T2s))+V1^2); // Isentropic Velocity and (answer mentioned in the textbook is wrong)
V2=sqrt (eff_nozzle) *V2s; // Actual nozzle exit velocity
disp ("m/s (answer mentioned in the textbook is wrong)",V2," Actual nozzle exit velocity = ","(b).Air");
|
