diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3785/CH12/EX12.5/Ex12_5.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '3785/CH12/EX12.5/Ex12_5.sce')
| -rw-r--r-- | 3785/CH12/EX12.5/Ex12_5.sce | 16 |
1 files changed, 16 insertions, 0 deletions
diff --git a/3785/CH12/EX12.5/Ex12_5.sce b/3785/CH12/EX12.5/Ex12_5.sce new file mode 100644 index 000000000..1563e54cc --- /dev/null +++ b/3785/CH12/EX12.5/Ex12_5.sce @@ -0,0 +1,16 @@ +// Example 12_5
+clc;funcprot(0);
+// Given data
+r=1.3;// The specific heat ratio
+V_e=0.90;// Exit velocity in % of maximum possible velocity
+
+// Calculation
+V_mebya_s=sqrt(2/(r-1));// The maximum possible exit velocity
+function[X]=machnumber(y)
+ X(1)=(V_e*V_mebya_s)-(y(1)*(1+((0.3*y(1)^2)/2))^(-1/2));
+endfunction
+y=[1];
+M_e=fsolve(y,machnumber);// The Mach number of the exit flow
+A_ebyA_c=M_e*((2/(r+1))*((1+(((r-1)/2)*(M_e)^2))))^((r+1)/(2*(r-1)));// The area ratio
+D_r=sqrt(A_ebyA_c);// Corresponding diameter ratio
+printf("\nThe exit Mach number,M_e=%1.3f \nThe area ratio=%1.3e \nThe diameter ratio=%2.2f",M_e,A_ebyA_c,D_r)
|