initial commit / add all books

26 files changed, 232 insertions, 0 deletions

diff --git a/797/CH3/EX3.1.e/3_01_example.sci b/797/CH3/EX3.1.e/3_01_example.sci
new file mode 100644
index 000000000..afc8905df
--- /dev/null
+++ b/797/CH3/EX3.1.e/3_01_example.sci
@@ -0,0 +1,3 @@
+//Example 3-1 Absolute Pressure of Vacuum Chamber
+P_atm = 100 //Atmospheric pressure [kPa]
+P_vac = 40 //Vacuum Gauge pressure [kPa]
diff --git a/797/CH3/EX3.1.s/3_01_solution.sce b/797/CH3/EX3.1.s/3_01_solution.sce
new file mode 100644
index 000000000..06d546448
--- /dev/null
+++ b/797/CH3/EX3.1.s/3_01_solution.sce
@@ -0,0 +1,7 @@
+//Soultion 3-01
+WD=get_absolute_file_path('3_01_solution.sce');
+datafile=WD+filesep()+'3_01_example.sci';
+clc;
+exec(datafile)
+P_abs = P_atm - P_vac; //Pressure relationship
+printf("Absolute pressure of Vacuum chamber is %1.2f kPa", P_abs);
diff --git a/797/CH3/EX3.10.e/3_10_example.sci b/797/CH3/EX3.10.e/3_10_example.sci
new file mode 100644
index 000000000..1c394d6b9
--- /dev/null
+++ b/797/CH3/EX3.10.e/3_10_example.sci
@@ -0,0 +1,4 @@
+//Example 3-10 Measuring specific gravity by hydrometer
+D = 1 //diameter of hydrometer [cm]
+h = 10 //height of water surface from the bottom of hydrometer [cm]
+rho = 1000 //density of water [kg/m^3]
diff --git a/797/CH3/EX3.10.s/3_10_solution.sce b/797/CH3/EX3.10.s/3_10_solution.sce
new file mode 100644
index 000000000..996f78e9f
--- /dev/null
+++ b/797/CH3/EX3.10.s/3_10_solution.sce
@@ -0,0 +1,9 @@
+//Solution 3-10
+WD=get_absolute_file_path('3_10_solution.sce');
+datafile=WD+filesep()+'3_10_example.sci';
+clc;
+exec(datafile)
+R = D / 2 / 100; //radius of cylinder [m]
+h = h / 100;
+m = rho * (%pi * R**2 * h);
+printf("The mass of lead required for hydrostat to attain depth of 10cm is %1.5f kg", m);
diff --git a/797/CH3/EX3.11.e/3_11_example.sci b/797/CH3/EX3.11.e/3_11_example.sci
new file mode 100644
index 000000000..b9735f893
--- /dev/null
+++ b/797/CH3/EX3.11.e/3_11_example.sci
@@ -0,0 +1,5 @@
+//Example 3-11 Weight loss of an Object in Seawater
+rho_seawater = 1025 //density of seawater [kg/m^3]
+V = 0.4 * 0.4 * 3 //volume of block [m^3]
+rho_block = 2300 //density of block [kg/m^3]
+g = 9.81 //gravitational acceleration [m/s^2]
diff --git a/797/CH3/EX3.11.s/3_11_solution.sce b/797/CH3/EX3.11.s/3_11_solution.sce
new file mode 100644
index 000000000..d1eff4477
--- /dev/null
+++ b/797/CH3/EX3.11.s/3_11_solution.sce
@@ -0,0 +1,12 @@
+//Solution 3-11
+WD=get_absolute_file_path('3_11_solution.sce');
+datafile=WD+filesep()+'3_11_example.sci';
+clc;
+exec(datafile)
+//(a)
+F_Tair = rho_block * V * g;
+printf("Tension in the rope of crane \n1.block suspended in air= %1.2f N", F_Tair/1000);
+//(b)
+F_B = rho_seawater * g * V; //bouyancy force
+F_Twater = F_Tair - F_B; //net force in downward diretion
+printf("\n2.block suspended in seawater= %1.2f kN", F_Twater / 1000);
diff --git a/797/CH3/EX3.12.e/3_12_example.sci b/797/CH3/EX3.12.e/3_12_example.sci
new file mode 100644
index 000000000..1168af4cf
--- /dev/null
+++ b/797/CH3/EX3.12.e/3_12_example.sci
@@ -0,0 +1,9 @@
+//Example 3-12 Overflow From Water Tank During Acceleration
+h = 80 //height of fish tank [cm]
+b_1 = 2 //one of the cross-sectional dimension of fish tank [m]
+b_2 = 0.6 //other cross-sectional dimension of fish tank [m]
+V_0 = 0 //initial velocity of truck [km/h]
+V_1 = 90 //velocity of truck after 10 s [km/h]
+t = 10 //duration of acceleration of truck [s]
+g = 9.81 //gravitaional acceleration [m/s^2]
+a_z = 0 //other acceleration in Z direction [m/s^2]
diff --git a/797/CH3/EX3.12.s/3_12_solution.sce b/797/CH3/EX3.12.s/3_12_solution.sce
new file mode 100644
index 000000000..d1b1745b0
--- /dev/null
+++ b/797/CH3/EX3.12.s/3_12_solution.sce
@@ -0,0 +1,20 @@
+//Solution 3-12
+WD=get_absolute_file_path('3_12_solution.sce');
+datafile=WD+filesep()+'3_12_example.sci';
+clc;
+exec(datafile)
+a_x = (V_1 - V_0) / t; //acceleration = rate of change of velocity  (horizontal)
+a_x = a_x / 3.6 //converting acceleration to [m/s^2]
+theta = atan(a_x / (g + a_z)) //angle made by free surface of water with horizontal [radians]
+printf("Vertical rise at the back of the tank relative to the midplande is")
+//Case 1:
+deltaz_1 = b_1 / 2 * tan(theta);
+printf("\n1.For long side parallel to direction of motion =%1.2f cm", deltaz_1 * 100);
+//Case 2:
+deltaz_2 = b_2 / 2 * tan(theta);
+printf("\n2.For short side parallel to direction of motion =%1.2f cm", deltaz_2 * 100);
+if(deltaz_2 < deltaz_1)
+    printf("\n Hence short side must be parallel to the direction of motion.");
+else
+    printf("\n Hence long side must be parallel to the direction of motion ");
+end
diff --git a/797/CH3/EX3.13.e/3_13_example.sci b/797/CH3/EX3.13.e/3_13_example.sci
new file mode 100644
index 000000000..5598e99a6
--- /dev/null
+++ b/797/CH3/EX3.13.e/3_13_example.sci
@@ -0,0 +1,6 @@
+//Example 3-13 Rising of Liquid During Rotation
+D = 20 //diameter of cylinder [cm]
+H = 60 //height of cylinder [cm]
+h_0 = 50 //height of liquid in the container [cm]
+rho = 850 //density of liquid in the container [kg/m^3]
+g = 9.81 //gravitaional acceleration [m/s^2]
diff --git a/797/CH3/EX3.13.s/3_13_solution.sce b/797/CH3/EX3.13.s/3_13_solution.sce
new file mode 100644
index 000000000..ca2840a60
--- /dev/null
+++ b/797/CH3/EX3.13.s/3_13_solution.sce
@@ -0,0 +1,13 @@
+//Solution 3-13
+WD=get_absolute_file_path('3_13_solution.sce');
+datafile=WD+filesep()+'3_13_example.sci';
+clc;
+exec(datafile)
+H = H / 100; //converting height from [cm] to [m]
+R = D / 2 / 100; //radius of cylinder [m]
+h_0 = h_0 / 100;
+omega = sqrt(4 * g * (H - h_0) / R**2); //from equation for the free surface of liquid
+ndot = omega / (2 * %pi) * 60; //rotational speed [rpm]
+printf("Rotational speed of the container must be restricted to %1.2f rpm to avoid any spill of liquid as a result of centrifugal effect", ndot);
+z_0 = h_0 - omega^2 * R^2 / (4 * g); //height of liquid at the center [m]
+printf("\nHeight of liquid at the center is %1.2f m>0 hence our assumption is valid", z_0);
diff --git a/797/CH3/EX3.2.e/3_02_example.sci b/797/CH3/EX3.2.e/3_02_example.sci
new file mode 100644
index 000000000..55a7b2dca
--- /dev/null
+++ b/797/CH3/EX3.2.e/3_02_example.sci
@@ -0,0 +1,4 @@
+//Example 3-02 Measuring Atmospheric Pressure with Barometer
+h = 740 //height of mercury column [m]
+g = 9.805 //gravitational acceleration [m^2/s]
+rho = 13570 //density of mercury [kg/m^3]
diff --git a/797/CH3/EX3.2.s/3_02_solution.sce b/797/CH3/EX3.2.s/3_02_solution.sce
new file mode 100644
index 000000000..c487614c3
--- /dev/null
+++ b/797/CH3/EX3.2.s/3_02_solution.sce
@@ -0,0 +1,10 @@
+//Soultion 3-02
+WD=get_absolute_file_path('3_02_solution.sce');
+datafile=WD+filesep()+'3_02_example.sci';
+clc;
+exec(datafile)
+h = h / 1000; //converting height of Hg column from [mm] to [m]
+P = rho * g * h; //Basic pressure eqaution [Pa]
+P = P / 1000;
+//result
+printf("Atmospheric pressure is %1.1f kPa", P);
diff --git a/797/CH3/EX3.3.e/3_03_example.sci b/797/CH3/EX3.3.e/3_03_example.sci
new file mode 100644
index 000000000..e33077f2d
--- /dev/null
+++ b/797/CH3/EX3.3.e/3_03_example.sci
@@ -0,0 +1,5 @@
+//Example 3-03 Gravity Driven flow in IV bottle
+rho = 1020 //density of IV fluid [kg/m^3]
+h_bottle1 = 1.2 //height of bottle for blood pressure balance
+P_gauge2 = 20 //gauge pressure required for sufficient flow rate
+g = 9.81 //gravitational acceleration [m^2/s]
diff --git a/797/CH3/EX3.3.s/3_03_solution.sce b/797/CH3/EX3.3.s/3_03_solution.sce
new file mode 100644
index 000000000..9e76c0cdb
--- /dev/null
+++ b/797/CH3/EX3.3.s/3_03_solution.sce
@@ -0,0 +1,13 @@
+//Soultion 3-03
+WD=get_absolute_file_path('3_03_solution.sce');
+datafile=WD+filesep()+'3_03_example.sci';
+clc;
+exec(datafile)
+//(a)
+P_gauge1 = rho * g * h_bottle1; //Basic Pressure formula
+P_gauge1 = P_gauge1 / 1000; //conversion from [Pa] to [kPa]
+printf("Gauge pressure of blood is %1.2f kPa", P_gauge1);
+//(b)
+P_gauge2 = P_gauge2 * 1000;
+h_bottle2 = P_gauge2 / (rho * g);
+printf("\nHeight required for maintaining 20kPa pressure is %1.2f m", h_bottle2);
diff --git a/797/CH3/EX3.4.e/3_04_example.sci b/797/CH3/EX3.4.e/3_04_example.sci
new file mode 100644
index 000000000..915ea4c99
--- /dev/null
+++ b/797/CH3/EX3.4.e/3_04_example.sci
@@ -0,0 +1,5 @@
+//Example 3-04 Hydrostatic pressure in solar pond with variable density
+rho_0 = 1040 //density on water surface [kg/m^3]
+H = 4 //thickness of the gradient zone [m]
+h_1 = 0.8 //thickness of surface zone [m]
+g = 9.81 //gravitational acceleration [m^2/s]
diff --git a/797/CH3/EX3.4.s/3_04_solution.sce b/797/CH3/EX3.4.s/3_04_solution.sce
new file mode 100644
index 000000000..3aefd7c9f
--- /dev/null
+++ b/797/CH3/EX3.4.s/3_04_solution.sce
@@ -0,0 +1,10 @@
+//Soultion 3-04
+WD=get_absolute_file_path('3_04_solution.sce');
+datafile=WD+filesep()+'3_04_example.sci';
+clc;
+exec(datafile)
+P_1 = rho_0 * g * h_1; //Gauge pressure at the bottom of surface zone [Pa]
+P_2 = P_1 + rho_0 * g * 4 *H / %pi * asinh(tan( %pi * H / (4 * H))); //After integrating w.r.t depth s
+P_2 = P_2 / 1000; //conversion from [Pa] to [kPa]
+//result
+printf("Pressure at the bottom of Gradient layer is %1.1f kPa",P_2);
diff --git a/797/CH3/EX3.5.e/3_05_example.sci b/797/CH3/EX3.5.e/3_05_example.sci
new file mode 100644
index 000000000..379012a62
--- /dev/null
+++ b/797/CH3/EX3.5.e/3_05_example.sci
@@ -0,0 +1,6 @@
+//Example 3-5 Measuring Pressure with Manometer
+SG = 0.85 //specific gravity of manometric fluid
+h = 55 //manometer column height [cm]
+P_atm = 96 //Local atmospheric pressure [kPa]
+g = 9.81 //gravitational acceleration [m^2/s]
+rho_water = 1000 //density of water [kg/m^3]
diff --git a/797/CH3/EX3.5.s/3_05_solution.sce b/797/CH3/EX3.5.s/3_05_solution.sce
new file mode 100644
index 000000000..9e05db619
--- /dev/null
+++ b/797/CH3/EX3.5.s/3_05_solution.sce
@@ -0,0 +1,12 @@
+//Soultion 3-05
+WD=get_absolute_file_path('3_05_solution.sce');
+datafile=WD+filesep()+'3_05_example.sci';
+clc;
+exec(datafile)
+rho = SG * rho_water; //definition of specific gravity
+P_atm = P_atm * 1000; //converting from [kPa] to [Pa]
+h = h / 100; //converting from [cm] to [m]
+P = P_atm + rho * g * h; //Pressure in manometer
+P = P / 1000; //converting from [Pa] to [kPa]
+//result
+printf("Absolute pressure in the tank is %1.1f kPa", P);
diff --git a/797/CH3/EX3.6.e/3_06_example.sci b/797/CH3/EX3.6.e/3_06_example.sci
new file mode 100644
index 000000000..bb9bb86e0
--- /dev/null
+++ b/797/CH3/EX3.6.e/3_06_example.sci
@@ -0,0 +1,9 @@
+//Example 3-06 Measuring Fluid with Multifluid Manometer
+P_atm = 85.6 //Atmospheric pressure at 1400m altitude [kPa]
+h_1 = 0.1 //differnce of water and oil level in manometer [m]
+h_2 = 0.2 //difference between water and mercury level in manometer [m]
+h_3 = 0.35 //difference between oil and mercury level in manometer [m]
+rho_water = 1000 //density of water [kg/m^3]
+rho_oil = 850 //density of oil [kg/m^3]
+rho_mercury = 13600 //density of mercury [kg/m^3]
+g = 9.81 //gravitational acceleration [m^2/s]
diff --git a/797/CH3/EX3.6.s/3_06_solution.sce b/797/CH3/EX3.6.s/3_06_solution.sce
new file mode 100644
index 000000000..9869f186e
--- /dev/null
+++ b/797/CH3/EX3.6.s/3_06_solution.sce
@@ -0,0 +1,10 @@
+//Soultion 3-06
+WD=get_absolute_file_path('3_06_solution.sce');
+datafile=WD+filesep()+'3_06_example.sci';
+clc;
+exec(datafile)
+P_atm = P_atm * 1000;
+P_1 = P_atm - rho_water * g * h_1 - rho_oil * g * h_2 + rho_mercury * g * h_3; //pressure equilibrium
+P_1 = P_1 / 1000; //converting from [Pa] to [kPa]
+//result
+printf("Air pressure in the tank is %1.0f kPa", P_1);
diff --git a/797/CH3/EX3.7.e/3_07_example.sci b/797/CH3/EX3.7.e/3_07_example.sci
new file mode 100644
index 000000000..a1c7254a6
--- /dev/null
+++ b/797/CH3/EX3.7.e/3_07_example.sci
@@ -0,0 +1,9 @@
+//Example 3-07 Analyzing Multifluid manometer
+P_atm = 85.6 //Atmospheric pressure at 1400m altitude [kPa]
+P_1 = 130 // Air pressure in tank [kPa]
+h_1 = 0.1 //differnce of water and oil level in manometer [m]
+h_2 = 0.2 //difference between water and mercury level in manometer [m]
+rho_water = 1000 //density of water [kg/m^3]
+rho_oil = 850 //density of oil [kg/m^3]
+rho_seawater = 1030//density of seawater [kg/m^3]
+g = 9.81 //gravitational acceleration [m^2/s]
diff --git a/797/CH3/EX3.7.s/3_07_solution.sce b/797/CH3/EX3.7.s/3_07_solution.sce
new file mode 100644
index 000000000..5d4a26d12
--- /dev/null
+++ b/797/CH3/EX3.7.s/3_07_solution.sce
@@ -0,0 +1,10 @@
+//Soultion 3-07
+WD=get_absolute_file_path('3_07_solution.sce');
+datafile=WD+filesep()+'3_07_example.sci';
+clc;
+exec(datafile)
+//converting pressures into [Pa]
+P_atm = P_atm * 1000;
+P_1 = P_1 * 1000;
+h_3 = (P_1 - P_atm + rho_water * g *h_1 + rho_oil * g * h_2) / (rho_seawater * g); //pressure eqquilibrium
+printf("If mercury is changed to seawater the height of seawater will be %1.2f m", h_3);
diff --git a/797/CH3/EX3.8.e/3_08_example.sci b/797/CH3/EX3.8.e/3_08_example.sci
new file mode 100644
index 000000000..3e8a6a146
--- /dev/null
+++ b/797/CH3/EX3.8.e/3_08_example.sci
@@ -0,0 +1,6 @@
+//Example 3-08 Hydrostatic Force Acting on the Door of Submerged car
+s = 8 //depth of car door top from the water surface [m]
+b = 1.2 //height of car door [m]
+h = 1 //breath of car door [m]
+rho = 1000 //density of water [kg/m^3]
+g = 9.81 //gravitational acceleration [m/s^2]
diff --git a/797/CH3/EX3.8.s/3_08_solution.sce b/797/CH3/EX3.8.s/3_08_solution.sce
new file mode 100644
index 000000000..77c073a66
--- /dev/null
+++ b/797/CH3/EX3.8.s/3_08_solution.sce
@@ -0,0 +1,12 @@
+//Soultion 3-08
+WD=get_absolute_file_path('3_08_solution.sce');
+datafile=WD+filesep()+'3_08_example.sci';
+clc;
+exec(datafile)
+A = b * h; //area of door
+P_avg = rho * g * (s + b / 2); //hydrostatic pressure formula
+F_g = P_avg * A;
+F_g= F_g / 1000; //conversion from [N] to [kN]
+printf("Hydrostatic force on the door is %1.2f kN", F_g);
+y_p = s + b / 2 + b**2 / (12 * (s + b / 2)) //formula for centre of pressure for P_o=0
+printf("\nThe center of pressure is %1.2f m", y_p)
diff --git a/797/CH3/EX3.9.e/3_09_example.sci b/797/CH3/EX3.9.e/3_09_example.sci
new file mode 100644
index 000000000..ab4501b17
--- /dev/null
+++ b/797/CH3/EX3.9.e/3_09_example.sci
@@ -0,0 +1,5 @@
+//Example 3-09 A Gravity Controlled Cylinderical Gate
+R = 0.8 //radius of cylinder [m]
+h_bottom = 5 //maximum level of water in tank [m]
+g = 9.81 //gravitational acceleration
+rho = 1000 //density of water [kg/m^3]
diff --git a/797/CH3/EX3.9.s/3_09_solution.sce b/797/CH3/EX3.9.s/3_09_solution.sce
new file mode 100644
index 000000000..256158edd
--- /dev/null
+++ b/797/CH3/EX3.9.s/3_09_solution.sce
@@ -0,0 +1,18 @@
+//Solution 3-09
+WD=get_absolute_file_path('3_09_solution.sce');
+datafile=WD+filesep()+'3_09_example.sci';
+clc;
+exec(datafile)
+//(a)
+s = h_bottom - R; //distance of cylinder top from water surface
+F_h = rho * g * (s + R / 2) * R; //horizontal force acting on vericle surface of cylinder
+F_y= rho * g * h_bottom * R; //vericle for acting on cylinder
+W = rho * g * R**2 * (1 - %pi / 4); //weight of fluid block per m width
+F_v = F_y - W; //net upward force
+F_R = sqrt(F_v**2 + F_h**2); //magnitude of resultant force
+theta = atan(F_v / F_h) * 180 / %pi; //angle made by resultant with horizontal
+printf("Resultant hydrostatic force acting on cylinder is %1.2f kN", F_R/1000);
+printf("\nAngle made by hydrostatic force with horizontal is %1.2f degrees", theta);
+//(b)
+W_cyl = F_R * sin(theta * %pi / 180); //equating moment at hinge to zero
+printf("\nWeight of cylinder per m length is %1.2f kN", W_cyl/1000);
\ No newline at end of file