diff options
Diffstat (limited to '629/CH12/EX12.5')
| -rw-r--r-- | 629/CH12/EX12.5/ex12_5.txt | 9 | ||||
| -rw-r--r-- | 629/CH12/EX12.5/example12_5.sce | 16 |
2 files changed, 25 insertions, 0 deletions
diff --git a/629/CH12/EX12.5/ex12_5.txt b/629/CH12/EX12.5/ex12_5.txt new file mode 100644 index 000000000..fd0cafe91 --- /dev/null +++ b/629/CH12/EX12.5/ex12_5.txt @@ -0,0 +1,9 @@ +
+The Mach number for downstream of the shock wave, M = 0.668.
+
+
+The pressure for downstream of the shock wave, P = 282 kPa, absolute.
+
+
+The temperature for downstream of the shock wave, T = 400 K or 127°C.
+
\ No newline at end of file diff --git a/629/CH12/EX12.5/example12_5.sce b/629/CH12/EX12.5/example12_5.sce new file mode 100644 index 000000000..1848c77e1 --- /dev/null +++ b/629/CH12/EX12.5/example12_5.sce @@ -0,0 +1,16 @@ +clear
+clc
+//Example 12.5 PROPERTY CHANGES ACROSS NORMAL SHOCK WAVE
+k=1.4;
+M1=1.6; //Upstream Mach number
+p1=100; //pressure[kPa]
+T1=15+273; //[K]
+//Downstream Mach number
+M2=sqrt(((k-1)*M1^2+2)/(2*k*M1^2-(k-1)))
+printf("\nThe Mach number for downstream of the shock wave, M = %.3f.\n\n",M2)
+//Downstream pressure
+p2=p1*[(1+k*M1^2)/(1+k*M2^2)] //[kPa]
+printf("\nThe pressure for downstream of the shock wave, P = %.f kPa, absolute.\n\n",p2)
+//Downstream temperature
+T2=T1*{(1+[(k-1)/2]*M1^2)/(1+[(k-1)/2]*M2^2)} //[K]
+printf("\nThe temperature for downstream of the shock wave, T = %.f K or %.f°C.\n",T2,T2-273)
\ No newline at end of file |