diff options
Diffstat (limited to '689/CH3')
-rw-r--r-- | 689/CH3/EX3.1/1.sce | 16 | ||||
-rw-r--r-- | 689/CH3/EX3.2/2.sce | 18 | ||||
-rw-r--r-- | 689/CH3/EX3.3/3.sce | 18 | ||||
-rw-r--r-- | 689/CH3/EX3.4/4.sce | 14 | ||||
-rw-r--r-- | 689/CH3/EX3.5/5.sce | 13 | ||||
-rw-r--r-- | 689/CH3/EX3.6/6.sce | 16 | ||||
-rw-r--r-- | 689/CH3/EX3.7/7.sce | 14 |
7 files changed, 109 insertions, 0 deletions
diff --git a/689/CH3/EX3.1/1.sce b/689/CH3/EX3.1/1.sce new file mode 100644 index 000000000..331eb6908 --- /dev/null +++ b/689/CH3/EX3.1/1.sce @@ -0,0 +1,16 @@ +clc; funcprot(0);
+//Example 3.1 Equation of states
+
+// Initialisation of variables
+T = 45+459.4;
+P = 25.93;
+P0 = 29.92;
+T0 = 518.4;
+rho_0 = 0.002378
+
+// Calculations
+rho = P*rho_0*T0/(P0*T);
+
+//Results
+disp(rho , "Density of dry air when pressure is 25.93 inch (pound/inch3):");
+
diff --git a/689/CH3/EX3.2/2.sce b/689/CH3/EX3.2/2.sce new file mode 100644 index 000000000..8ca4c518e --- /dev/null +++ b/689/CH3/EX3.2/2.sce @@ -0,0 +1,18 @@ +clc; funcprot(0);
+//Example 3.2 Equation of states
+
+// Initialisation of variables
+T = -10+459.4;
+P = 16.38;
+P0 = 29.92;
+T0 = 518.4;
+rho_0 = 0.002378
+g = 32.1740;
+
+// Calculations
+rho = P*rho_0*T0/(P0*T);
+W = rho*g;
+
+//Results
+disp(W , "Spwcific gravity of dry air (lb/ft3):");
+
diff --git a/689/CH3/EX3.3/3.sce b/689/CH3/EX3.3/3.sce new file mode 100644 index 000000000..593cdd39d --- /dev/null +++ b/689/CH3/EX3.3/3.sce @@ -0,0 +1,18 @@ +clc; funcprot(0);
+//Example 3.3 On Adiabatic Process
+
+// Initialisation of variables
+gma = 1.4;
+rho_0 = 0.002378;
+P1 = 2*2116.2; // Pressure ion lb per sq ft
+P0 = 1*2116.2; // Pressure ion lb per sq ft
+T0 = 59+459.4;
+
+// Calculations
+rho1 = rho_0*(P1/P0)^(1/gma);
+T1 =T0*(rho_0/rho1)*(P1/P0);
+
+//Results
+
+disp(T1-459.4 , "(b)Temperature if air is compressed adiabatically to 2 atm (degree farenheit):",rho1, "(a)Density if air is compressed adiabatically to 2 atm (Slug per cu ft):");
+
diff --git a/689/CH3/EX3.4/4.sce b/689/CH3/EX3.4/4.sce new file mode 100644 index 000000000..f110b6f19 --- /dev/null +++ b/689/CH3/EX3.4/4.sce @@ -0,0 +1,14 @@ +clc; funcprot(0);
+//Example 3.4 On Speed of sound
+
+// Initialisation of variables
+gma = 1.4;
+g = 32.174;
+R = 53.351;
+T = 59+459.4;
+
+// Calculations
+a = sqrt(gma*g*R*T);
+//Results
+disp(a , "Speed of sound in standard temperature 59 degree farenheit (ft/sec):");
+
diff --git a/689/CH3/EX3.5/5.sce b/689/CH3/EX3.5/5.sce new file mode 100644 index 000000000..95c64eeb7 --- /dev/null +++ b/689/CH3/EX3.5/5.sce @@ -0,0 +1,13 @@ +clc; funcprot(0);
+//Example 3.5 Speed of sound
+
+// Initialisation of variables
+gma = 1.4;
+P = 2116.2;
+rho = 0.002378;
+
+// Calculations
+a = sqrt(gma*P/rho);
+//Results
+disp(a , "Speed of sound in standard pressure 2116.2 lbper sq ft, and standard density 0.002378slug per cu ft. (in ft/sec):");
+
diff --git a/689/CH3/EX3.6/6.sce b/689/CH3/EX3.6/6.sce new file mode 100644 index 000000000..477408485 --- /dev/null +++ b/689/CH3/EX3.6/6.sce @@ -0,0 +1,16 @@ +clc; funcprot(0);
+//Example 3.6 Bernoulli Equation for compressible flow
+
+// Initialisation of variables
+gma = 1.4;
+P0 = 14.7;
+rho = 0.002378;
+V0= 500;
+P1 = 13.5;
+
+// Calculations
+a0 = sqrt(gma*P/rho);
+V1 =sqrt(V0^2 + 2*a0^2*(1-(P1/P0)^(1-1/gma))/(gma-1));
+
+//Results
+disp(V1 , "Speed of sound when pressure is 13.5 lb per sq in (in ft/sec):");
diff --git a/689/CH3/EX3.7/7.sce b/689/CH3/EX3.7/7.sce new file mode 100644 index 000000000..0f48b741c --- /dev/null +++ b/689/CH3/EX3.7/7.sce @@ -0,0 +1,14 @@ +clc; funcprot(0);
+//Example 3.7 Similar Flows
+
+// Initialisation of variables
+c = 3/12; //chord length in feet
+V = 100*1.467; // velocity in ft/sec
+R = 0.002378;
+mu = 0.000000373;
+
+// Calculations
+RN = rho*V*c/mu;
+//Results
+disp(RN , "Reynolds no:");
+
|