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 /629/CH12/EX12.9/example12_9.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 '629/CH12/EX12.9/example12_9.sce')
| -rw-r--r-- | 629/CH12/EX12.9/example12_9.sce | 19 |
1 files changed, 19 insertions, 0 deletions
diff --git a/629/CH12/EX12.9/example12_9.sce b/629/CH12/EX12.9/example12_9.sce new file mode 100644 index 000000000..c7d9e299b --- /dev/null +++ b/629/CH12/EX12.9/example12_9.sce @@ -0,0 +1,19 @@ +clear
+clc
+//Example 12.9 MASS FLOW RATE IN SUPERSONIC WIND TUNNEL
+k=1.4;
+R=287; //[J/kg.K]
+M=3; //Mach number
+A=225*10^-4; //[m^2]
+//Throat area
+Ao=A/((1/M)*{(1+[(k-1)/2]*M^2)/((k+1)/2)}^((k+1)/(2*(k-1)))) //[m^2]
+//Static T,p
+T=273+(-20) //[K]
+p=50; //[kPa]
+//Total temperature
+Tt=T*(1+(k-1)*M^2/2) //[K]
+//Total pressure
+pt=p*(1+(k-1)*M^2/2)^(k/(k-1)) //[kPa]
+//Mass flow rate
+m=k^(1/2)*[(2/(k+1))^((k+1)/(2*(k-1)))]*pt*10^3*Ao/((R*Tt)^(1/2)) //[kg]
+printf("\n The mass flow rate = %.1f kg/s.\n",m)
\ No newline at end of file |