diff options
Diffstat (limited to '629/CH12/EX12.6')
-rw-r--r-- | 629/CH12/EX12.6/ex12_6.txt | 3 | ||||
-rw-r--r-- | 629/CH12/EX12.6/example12_6.sce | 20 |
2 files changed, 23 insertions, 0 deletions
diff --git a/629/CH12/EX12.6/ex12_6.txt b/629/CH12/EX12.6/ex12_6.txt new file mode 100644 index 000000000..51179500f --- /dev/null +++ b/629/CH12/EX12.6/ex12_6.txt @@ -0,0 +1,3 @@ +
+ The change in entropy across the wave = 20.5 J/kg.K.
+
\ No newline at end of file diff --git a/629/CH12/EX12.6/example12_6.sce b/629/CH12/EX12.6/example12_6.sce new file mode 100644 index 000000000..f50a6ce45 --- /dev/null +++ b/629/CH12/EX12.6/example12_6.sce @@ -0,0 +1,20 @@ +clear
+clc
+//Example 12.6 ENTROPY INCREASE ACROSS SHOCK WAVE
+//To find Approx Value
+function [A]= approx (V,n)
+ A= round(V*10^n)/10^n; //V-Value, n-to what place
+ funcprot (0)
+endfunction
+k=1.4;
+M1=1.5;
+//Downstream Mach number
+M2=approx(sqrt(((k-1)*M1^2+2)/(2*k*M1^2-(k-1))),3)
+//Pressure ratio, (p21=p2/p1)
+p21=approx([(1+k*M1^2)/(1+k*M2^2)],2)
+//Temperature ratio, (T21=T2/T1)
+T21=approx((1+[(k-1)/2]*M1^2)/(1+[(k-1)/2]*M2^2),2)
+R=287; //[J/kg.K]
+//Entropy change
+del_s=R*log((T21^(k/(k-1))/p21)) //[J/kg.K]
+printf("\n The change in entropy across the wave = %.1f J/kg.K.\n",del_s)
\ No newline at end of file |