From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 629/CH12/EX12.9/example12_9.sce | 19 +++++++++++++++++++
 1 file changed, 19 insertions(+)
 create mode 100644 629/CH12/EX12.9/example12_9.sce

(limited to '629/CH12/EX12.9/example12_9.sce')

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
-- 
cgit