1 files changed, 27 insertions, 0 deletions
diff --git a/746/CH12/EX12.02/12_02.sce b/746/CH12/EX12.02/12_02.sce
new file mode 100755
index 000000000..9410fab06
--- /dev/null
+++ b/746/CH12/EX12.02/12_02.sce
@@ -0,0 +1,27 @@
+//Mass flow//
+pathname=get_absolute_file_path('12.02.sce')
+filename=pathname+filesep()+'12.02-data.sci'
+exec(filename)
+//Checking for chocking:
+c=pb/p0;
+if(c<=0.528)
+  //choked
+else
+  //Not choked
+  //Therefore pressure at exit = back pressure
+  pe=pb;
+  //Mach number at exit:
+  Me=(((p0/pe)^((k-1)/k)-1)*(2/(k-1)))^0.5
+  //Temperature at exit(in K):
+  Te=T0/(1+(k-1)/2*Me^2)
+  //Velocity at exit(in m/sec):
+  Ve=Me*sqrt(k*R*Te)
+  //Density at exit(in kg/m^3):
+  de=pe*10^3/R/Te
+  //Mass flow rate of air(kg/sec):
+  m=de*Ve*Ae
+end;
+printf("\n\nRESULTS\n\n")
+printf("\n\nMach number at exit: %.3f\n\n",Me)
+printf("\n\nMass flow rate of air: %.3f kg/sec\n\n",m)
+