diff options
Diffstat (limited to '3720/CH12/EX12.12/Ex12_12.sce')
| -rw-r--r-- | 3720/CH12/EX12.12/Ex12_12.sce | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/3720/CH12/EX12.12/Ex12_12.sce b/3720/CH12/EX12.12/Ex12_12.sce new file mode 100644 index 000000000..1f2ce61eb --- /dev/null +++ b/3720/CH12/EX12.12/Ex12_12.sce @@ -0,0 +1,23 @@ +// Example 12_12
+clc;clear;funcprot(0);
+//Given values
+Ma_1=2.0;// Mach number
+P_1=230;// kPa
+delta=10;// degree
+//Properties
+k=1.4//The specific heat ratio
+
+//Calculation
+theta=delta;
+v_1=(sqrt((k+1)/(k-1))*atand(sqrt(((k-1)*(Ma_1^2-1))/(k+1))))-atand(sqrt(Ma_1^2-1));// degree
+v_2=theta+v_1;// degree
+// Ma_2=y(1);
+function[X]=Machnumber(y);
+ X(1)=((sqrt((k+1)/(k-1))*atand(sqrt(((k-1)*(y(1)^2-1))/(k+1))))-atand(sqrt(y(1)^2-1))-v_2);
+endfunction
+y=[1];
+z=fsolve(y,Machnumber);
+printf('The downstream Mach number Ma_2=%0.3f\n',z(1));
+Ma_2=z(1);
+P_2=((((1+(((k-1)/2)*Ma_2^2)))^(-k/(k-1)))/(((1+(((k-1)/2)*Ma_1^2)))^(-k/(k-1))))*(P_1);
+printf('The downstream pressure,P_2=%0.0f kPa\n',P_2);
|