1 files changed, 22 insertions, 0 deletions
diff --git a/896/CH8/EX8.6/6.sce b/896/CH8/EX8.6/6.sce
new file mode 100755
index 000000000..85bd4184a
--- /dev/null
+++ b/896/CH8/EX8.6/6.sce
@@ -0,0 +1,22 @@
+clc
+//Example 8.6
+//Calculate the cross sectional area, pressure, temperature and mach number at a pt in duct where air velocity is 1400ft/s
+P1=30//psia
+T1=660//R (Rankine temperature scale)
+m=10//lbm/s mass flow rate
+v1=1400//ft/s
+R=4.98*10^4//(ft^2/s^2)*(lbm/lbmol.R)^0.5
+k=1.4//dimentionless
+M=29//lbm/lbmol
+T2=T1-v1^2*((k-1)/k)*M/2/R//R (Rankine temperature scale)
+printf("The temperature at the pt in the duct where air velocity is 1400 ft/s is %f R\n",T2);
+c=223*(k*T2/M)^0.5//ft/s
+Ma=v1/c//dimentionless (Mach number)
+printf("The mach number at the pt in the duct where air velocity is 1400 ft/s is %f\n",Ma);
+P2=P1/(T1/T2)^(k/(k-1))//psia
+printf("The pressure at the pt in the duct where air velocity is 1400 ft/s is %f psia\n",P2);
+//1 lbf.s^2 = 32.2 lbm.ft
+A0=m/(P1*(M*k)^0.5*32.2/223/(T1)^0.5/((k-1)/2+1)^((k+1)/2/(k-1)))//in^2
+ratio_A=((Ma^2*(k-1)/2+1)/((k-1)/2+1))^((k+1)/2/(k-1))/Ma//dimentionless
+A=ratio_A*A0//in^2
+printf("The cross sectional at the pt in the duct where air velocity is 1400 ft/s is %f in^2",A);
+\ No newline at end of file