diff options
Diffstat (limited to '2126/CH5/EX5.14/14.sce')
| -rwxr-xr-x | 2126/CH5/EX5.14/14.sce | 35 |
1 files changed, 35 insertions, 0 deletions
diff --git a/2126/CH5/EX5.14/14.sce b/2126/CH5/EX5.14/14.sce new file mode 100755 index 000000000..b3bace8ac --- /dev/null +++ b/2126/CH5/EX5.14/14.sce @@ -0,0 +1,35 @@ +clc
+clear
+
+//Input data
+T=300 //Temperature in K
+P=1.01325*10^5 //Absolute pressure in Pa
+k=1.4 //Adiabatic constant
+R=287 //Specific gas constant in J/kg-K
+C_1=60 //Velocity of air in m/s
+C_2=200 //Velocity of air in m/s
+C_3=500 //Velocity of air in m/s
+d_hg=13600 //Density of mercury in kg/m^3
+g=9.81 //Acceleration due to gravity in m/s^2
+
+//Calculation
+a=sqrt(k*R*T) //Sound velocity in m/s
+M_1=C_1/a //Mach number
+dP1=(P*C_1^2)/(2*R*T) //Difference in mercury levels in Pa
+dP1_hg=(dP1/(d_hg*g))*1000 //Difference in mercury levels in mm of Hg
+M_2=C_2/a //Mach number
+p1=(1+((k-1)/2)*M_2^2)^(k/(k-1)) //Stagnation to static pressure ratio
+Po=p1*P //Stagnation pressure in Pa
+dP2=abs(Po-P) //Difference in mercury levels in Pa
+dP2_hg=(dP2/(d_hg*g))*1000 //Difference in mercury levels in mm of Hg
+M_3=C_3/a //Mach number & M_3=Mach number just before shock
+My=0.723 //Mach number just after shock
+p1=2.2530 //Ratio of pressure after shock to before shock from gas tables @My
+Py=p1*P //Pressure after shock in Pa
+p2=0.706 //Ratio of pressure after shock to Stagnation pressure from gas tables @My
+Po=Py/p2 //Stagnation pressure in Pa
+dP3=Po-Py //Difference in mercury levels in Pa
+dP3_hg=(dP3/(d_hg*g))*1000 //Difference in mercury levels in mm of Hg
+
+//Output
+printf('Difference in mercury levels at velocity equal to:\n (A)%2i m/s is %3.3f mm of Hg\n (B)%3i m/s is %3.1f mm of Hg\n (C)%3i m/s is %3i mm of Hg',C_1,dP1_hg,C_2,dP2_hg,C_3,dP3_hg)
|