summaryrefslogtreecommitdiff
path: root/629/CH12/EX12.6
diff options
context:
space:
mode:
Diffstat (limited to '629/CH12/EX12.6')
-rw-r--r--629/CH12/EX12.6/ex12_6.txt3
-rw-r--r--629/CH12/EX12.6/example12_6.sce20
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