summaryrefslogtreecommitdiff
path: root/3831/CH16/EX16.2/Ex16_2.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3831/CH16/EX16.2/Ex16_2.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3831/CH16/EX16.2/Ex16_2.sce')
-rw-r--r--3831/CH16/EX16.2/Ex16_2.sce16
1 files changed, 16 insertions, 0 deletions
diff --git a/3831/CH16/EX16.2/Ex16_2.sce b/3831/CH16/EX16.2/Ex16_2.sce
new file mode 100644
index 000000000..7f40f256f
--- /dev/null
+++ b/3831/CH16/EX16.2/Ex16_2.sce
@@ -0,0 +1,16 @@
+// Example 16_2
+clc;funcprot(0);
+// Given data
+T=20+273.15;// K
+V=25.0;// m/s
+k=1.40;// The specific heat ratio
+p=0.101;// MPa
+g_c=1;// The gravitational constant
+c_p=1.004;// kJ/kg.K
+R=0.286;// kJ/kg.K
+
+// Solution
+p_os=p*(1+((V^2/1000)/(2*g_c*c_p*T)))^(k/(k-1));// The isentropic stagnation pressure in MPa
+rho=(p*10^3)/(R*T);// kg/m^3
+rho_os=rho*(1+((V^2/1000)/(2*g_c*c_p*T)))^(1/(k-1));// The isentropic stagnation density in kg/m^3
+printf("\nThe isentropic stagnation pressure,p_os=%0.4f MPa \nThe isentropic stagnation density,rho_os=%1.4f kg/m^3",p_os,rho_os);