8 files changed, 261 insertions, 0 deletions
diff --git a/3564/CH1/EX1.1/Ex1_1.sce b/3564/CH1/EX1.1/Ex1_1.sce
new file mode 100644
index 000000000..ec0eee9f8
--- /dev/null
+++ b/3564/CH1/EX1.1/Ex1_1.sce
@@ -0,0 +1,30 @@
+
+// Display mode
+mode(0);
+
+// Display warning for floating point exception
+ieee(1);
+
+clc;
+disp("Introduction to Fluid Mechanics, 3rd Ed. William S. Janna Chapter - 1 Example # 1.1 ")
+
+//Weight in lbf
+F = 150;
+
+//Solving part a
+disp("Part a)")
+//Acceleration due to gravity in ft/s2
+a = 32.2;
+
+//Calculating mass of the person in slug
+disp("Mass of person in lbf.s2/ft (or slug) is")
+m = F/a
+//Answer varies slightly because of round-off error
+
+//Solving part b
+disp("Part b)")
+//New aceleration due to gravity in ft/s2
+a = a/6;
+//Calculating new weight of the person in lbf
+disp("Weight of person on the moon in lbf is")
+F = m*a
diff --git a/3564/CH1/EX1.2/Ex1_2.sce b/3564/CH1/EX1.2/Ex1_2.sce
new file mode 100644
index 000000000..cbe2ce887
--- /dev/null
+++ b/3564/CH1/EX1.2/Ex1_2.sce
@@ -0,0 +1,31 @@
+
+// Display mode
+mode(0);
+
+// Display warning for floating point exception
+ieee(1);
+
+clc;
+disp("Introduction to Fluid Mechanics, 3rd Ed. William S. Janna Chapter - 1 Example # 1.2 ")
+
+//Solving part a
+disp("Part a)")
+//Volume of water in m3
+V = 1;
+//Mass of water in kg
+m = 1000;
+//Acceleration due to gravity in m/s2
+a = 9.81;
+
+//Weight of 1 m3 water in N
+disp("Weight of water in N is")
+F = m*a
+
+//Solving part b
+disp("Part b)")
+//New aceleration due to gravity in m/s2
+a = (2*a)/5;
+//Calculating new weight of the water in N
+disp("Weight of water on Mars in N is")
+F = m*a
+//Answer varies slightly because of round-off error
diff --git a/3564/CH1/EX1.3/Ex1_3.sce b/3564/CH1/EX1.3/Ex1_3.sce
new file mode 100644
index 000000000..1d6ab8d39
--- /dev/null
+++ b/3564/CH1/EX1.3/Ex1_3.sce
@@ -0,0 +1,60 @@
+
+// Display mode
+mode(0);
+
+// Display warning for floating point exception
+ieee(1);
+
+clc;
+disp("Introduction to Fluid Mechanics, 3rd Ed. William S. Janna Chapter - 1 Example # 1.3 ")
+
+//Solving part a
+disp("Part a)")
+//Mass in kg
+m = 0.001;
+//Deltay in mm
+deltay = 5;
+//Acceleration due to gravity in m/s2
+g = 9.81;
+//Area of contact in m2
+A = 0.5;
+
+//Using Appendix table A.5 for properties of linseed oil
+//Viscosity myu in N.s/m2
+myu = 0.0331;
+
+//Force therefore in N is
+F = m*g;
+//Shear stress in N/m2 is
+tau = F/A;
+
+//Since shear stress is myu*(velocity gradient) i.e. myu*(deltaV/deltay)
+//deltaV = V - 0 = V
+
+//Velocity of the plate in mm/s is
+disp("Velocity of the plate in mm/s is")
+V = (tau*deltay)/myu
+//Answer varies slightly because of round-off error
+
+//Solving part b
+disp("Part b)")
+//Using Appendix table A.5 for properties of water
+//Viscosity myu in N.s/m2
+myu = 0.89/1000;
+
+//Velocity of the plate in mm/s is
+V = (tau*deltay)/myu;
+disp("Velocity of the plate in m/s is")
+V = V/1000
+//Answer varies slightly because of round-off error
+
+//Solving part c
+disp("Part c)")
+//Initial shear stress in N/m2 is
+tau0 = 4;
+//Inital viscosity myu0 in N.s/m2
+myu0 = 0.004;
+if tau<tau0 then
+  disp("Applied shear stress is less than initial shear stress")
+  disp("Therefore velocity of plate is 0 m/s")
+end;
diff --git a/3564/CH1/EX1.4/Ex1_4.sce b/3564/CH1/EX1.4/Ex1_4.sce
new file mode 100644
index 000000000..9399f2639
--- /dev/null
+++ b/3564/CH1/EX1.4/Ex1_4.sce
@@ -0,0 +1,31 @@
+
+// Display mode
+mode(0);
+
+// Display warning for floating point exception
+ieee(1);
+
+clc;
+disp("Introduction to Fluid Mechanics, 3rd Ed. William S. Janna Chapter - 1 Example # 1.4 ")
+
+//Solving part a
+disp("Part a)")
+disp("Part a is theoretical and does not require computation")
+disp("Final result is pi - p0 = 2*sigma/R")
+
+
+//Solving part b
+disp("Part b)")
+//Diameter of droplet in cm
+d = 0.01;//Authors have used 0.01 diameter for calculation while the diameter quoted in question is 0.1
+//Using Appendix table A.5 for properties of water
+//Surface tension sigma in N/m
+sigma = 71.97/1000;
+//Atmospheric pressure for droplet in N/m2 is
+p0 = 101300;
+//Radius of droplet in  m
+R = 0.01*(d/2);
+//Calculating pressure inside the droplet in N/m2
+disp("Pressure inside the droplet in N/m2 is")
+pi = p0+(2*sigma)/R
+//Answer varies slightly because of round-off error
diff --git a/3564/CH1/EX1.5/Ex1_5.sce b/3564/CH1/EX1.5/Ex1_5.sce
new file mode 100644
index 000000000..b2b7b6c3e
--- /dev/null
+++ b/3564/CH1/EX1.5/Ex1_5.sce
@@ -0,0 +1,12 @@
+
+// Display mode
+mode(0);
+
+// Display warning for floating point exception
+ieee(1);
+
+clc;
+disp("Introduction to Fluid Mechanics, 3rd Ed. William S. Janna Chapter - 1 Example # 1.5 ")
+
+disp("This example is theoretical and does not require computation")
+disp("Final result is h = (2*sigma*cos(theta))/(rho*R*g)")
diff --git a/3564/CH1/EX1.6/Ex1_6.sce b/3564/CH1/EX1.6/Ex1_6.sce
new file mode 100644
index 000000000..984185d3e
--- /dev/null
+++ b/3564/CH1/EX1.6/Ex1_6.sce
@@ -0,0 +1,37 @@
+
+// Display mode
+mode(0);
+
+// Display warning for floating point exception
+ieee(1);
+
+clc;
+disp("Introduction to Fluid Mechanics, 3rd Ed. William S. Janna Chapter - 1 Example # 1.6 ")
+
+//Solving part a
+disp("Part a)")
+//Mass of air in kg
+m = 4;
+//Using Appendix table A.6 for properties of air
+//Specific heat of air in J/Kg.K
+cp = 1005;
+//Representing temperature difference T2 - T1 as deltaT for computation
+//detaT in degree K
+deltaT = 10;
+
+//Heat required in kJ
+disp("Heat required in kJ is")
+Q = ((m*cp)*deltaT)/1000
+
+//Solving part b
+disp("Part b)")
+//Using Appendix table A.6 for properties of air
+funcprot(0);
+//Specific heat ratio gamma for air is
+gamma = 1.4;
+//Specific heat at constant volume of air in J/Kg.K
+cv = cp/gamma;
+//Heat required at constant volume in kJ
+disp("Heat required at constant volume  in kJ is")
+Q = ((m*cv)*deltaT)/1000
+//Answer varies slightly because of round-off error
diff --git a/3564/CH1/EX1.7/Ex1_7.sce b/3564/CH1/EX1.7/Ex1_7.sce
new file mode 100644
index 000000000..3794091ed
--- /dev/null
+++ b/3564/CH1/EX1.7/Ex1_7.sce
@@ -0,0 +1,29 @@
+
+// Display mode
+mode(0);
+
+// Display warning for floating point exception
+ieee(1);
+
+clc;
+disp("Introduction to Fluid Mechanics, 3rd Ed. William S. Janna Chapter - 1 Example # 1.7 ")
+
+//Using Appendix table A.6 for properties of air
+//Specific heat of air in Btu/slag
+cp = 7.72;
+//Specific heat ratio gamma for air is
+gamma = 1.4;
+//Specific heat at constant volume of air in Btu/slag.R
+cv = cp/gamma;
+
+//Inital temperature in degree F
+T1 = 70;//Symbol is assumed for computation
+//Final temperature in degree F
+T2 = 110;//Symbol is assumed for computation
+//Temperature difference deltaT in degree K
+deltaT = T2-T1;
+
+//Change in internal energy in Btu/slag
+disp("Change in internal energy in Btu/slag is")
+deltau = cv*deltaT
+//Answer varies slightly because of round-off error
diff --git a/3564/CH1/EX1.8/Ex1_8.sce b/3564/CH1/EX1.8/Ex1_8.sce
new file mode 100644
index 000000000..6dde947db
--- /dev/null
+++ b/3564/CH1/EX1.8/Ex1_8.sce
@@ -0,0 +1,31 @@
+
+// Display mode
+mode(0);
+
+// Display warning for floating point exception
+ieee(1);
+
+clc;
+disp("Introduction to Fluid Mechanics, 3rd Ed. William S. Janna Chapter - 1 Example # 1.8 ")
+
+//Using Appendix table A.5 for properties of acetone
+//Density in kg/m3
+rho = 0.787*1000;
+//Speed of sound in m/s
+a = 1174;
+//Isothermal bulk modulus in N/m2
+k = (rho*a)*a;
+
+//Change in volume is given to be -1 percent
+//Assuming V = 1 cubic m so that deltaV/V = -0.01
+
+//Initial volume in m3
+V = 1;
+//Change in volume in m3
+deltaV = -0.01;
+
+
+//Change in pressure in N/m2
+disp("Change in pressure in N/m2 is")
+deltap = -(k*deltaV)/V
+//Answer varies slightly because of round-off error