summaryrefslogtreecommitdiff
path: root/779/CH17/EX17.5/17_5.sce
diff options
context:
space:
mode:
Diffstat (limited to '779/CH17/EX17.5/17_5.sce')
-rwxr-xr-x779/CH17/EX17.5/17_5.sce22
1 files changed, 22 insertions, 0 deletions
diff --git a/779/CH17/EX17.5/17_5.sce b/779/CH17/EX17.5/17_5.sce
new file mode 100755
index 000000000..e0910babd
--- /dev/null
+++ b/779/CH17/EX17.5/17_5.sce
@@ -0,0 +1,22 @@
+Ax = 18.75; A_ = 12.50; // A_= A*
+AA_ = 1.5; // A/A*
+Mx = 1.86; Pxox = 0.159; R = 0.287;
+Pox = 0.21e03; // in kPa
+Px = Pxox*Pox;
+// from the gas table on normal shock
+Mx = 1.86; My = 0.604; Pyx = 3.87; Poyx = 4.95; Poyox = 0.786;
+Py = Pyx*Px;
+Poy = Poyx*Px;
+My = 0.604;
+Ay_ = 1.183;
+A2 = 25; Ay = 18.75;
+A2_ = (A2/Ay)*Ay_;
+// From isentropic table
+M2 = 0.402;
+P2oy = 0.895;
+P2 = P2oy*Poy;
+syx = -R*log(Poy/Pox); // sy-sx
+disp(M2,"Exit mach number is M2")
+disp("kPa",P2,"Exit pressure is")
+disp("kPa",Pox-Poy,"Exit Stagnation pressure is")
+disp("kJ/kg K",syx,"Entropy increase is") \ No newline at end of file