diff options
Diffstat (limited to '3720/CH12/EX12.3/Ex12_3.sce')
-rw-r--r-- | 3720/CH12/EX12.3/Ex12_3.sce | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/3720/CH12/EX12.3/Ex12_3.sce b/3720/CH12/EX12.3/Ex12_3.sce new file mode 100644 index 000000000..dd635aeb0 --- /dev/null +++ b/3720/CH12/EX12.3/Ex12_3.sce @@ -0,0 +1,24 @@ +//Example 12_3
+clc;clear;funcprot(0);
+// Given values
+m=3;//Mass flow rate in kg/s
+T_0=473;// T_0=T_1 in K
+P_0=1400;// P_0=P_1 in kPa
+P=1200;// kPa
+// Properties
+C_p=0.846;// kJ/(kg.K)
+k=1.289;
+R=0.1889;// kJ/(kg.K)
+
+//Calculation
+T=T_0*(P/P_0)^((k-1)/k);// k
+V=sqrt(2*C_p*(T_0-T)*1000);// m/s
+printf('Velocity ,V=%0.1f m/s\n',V);
+rho=P/(R*T);// kg/m^3
+printf('Density ,rho=%0.1f kg/m^3\n',rho);
+A=(m/(rho*V))*10000;//cm^2
+printf('Area ,A=%0.1f cm^2\n',A);
+c=sqrt(k*R*T*1000);// m/s
+Ma=V/c;
+printf('Mach number ,Ma=%0.3f \n',Ma);
+// The answer vary due to round off error
|