summaryrefslogtreecommitdiff
path: root/1187/CH11/EX11.2/2.sce
diff options
context:
space:
mode:
Diffstat (limited to '1187/CH11/EX11.2/2.sce')
-rwxr-xr-x1187/CH11/EX11.2/2.sce37
1 files changed, 37 insertions, 0 deletions
diff --git a/1187/CH11/EX11.2/2.sce b/1187/CH11/EX11.2/2.sce
new file mode 100755
index 000000000..e82a3c491
--- /dev/null
+++ b/1187/CH11/EX11.2/2.sce
@@ -0,0 +1,37 @@
+clc
+
+disp("(a) the angle through which the airstream is deflected")
+
+y=1.4;
+R=287; // J/kg.K
+T1=238; // K
+u1=773; // m/s
+beta1=38; // degrees
+cp=1005; // J/kg.K
+
+a1=sqrt(y*R*T1);
+M1=u1/a1;
+
+beta2=atand(tand(beta1)*((2+(y-1)*M1^2*(sind(beta1))^2)/((y+1)*M1^2*(sind(beta1))^2)));
+
+deflection_angle=beta1-beta2;
+disp("Deflection angle =")
+disp(deflection_angle)
+disp("degrees")
+
+disp("(b) the final Mach number")
+
+u2=u1*cosd(beta1)/cosd(beta2);
+
+T2=T1+1/(2*cp)*(u1^2-u2^2);
+a2=sqrt(y*R*T2);
+
+M2=u2/a2;
+
+disp("Final Mach number =")
+disp(M2)
+
+disp("(c) the pressure ratio across the wave.")
+ratio=T2/T1*(tand(beta1)/tand(beta2));
+disp("Pressure ratio =")
+disp(ratio) \ No newline at end of file