summaryrefslogtreecommitdiff
path: root/896/CH2
diff options
context:
space:
mode:
Diffstat (limited to '896/CH2')
-rwxr-xr-x896/CH2/EX2.1/1.sce10
-rwxr-xr-x896/CH2/EX2.10/10.sce10
-rwxr-xr-x896/CH2/EX2.11/11.sce18
-rwxr-xr-x896/CH2/EX2.12/12.sce9
-rwxr-xr-x896/CH2/EX2.13/13.sce12
-rwxr-xr-x896/CH2/EX2.14/14.sce14
-rwxr-xr-x896/CH2/EX2.15/15.sce9
-rwxr-xr-x896/CH2/EX2.16/16.sce11
-rwxr-xr-x896/CH2/EX2.17/17.sce21
-rwxr-xr-x896/CH2/EX2.18/18.sce9
-rwxr-xr-x896/CH2/EX2.19/19.sce10
-rwxr-xr-x896/CH2/EX2.2/2.sce11
-rwxr-xr-x896/CH2/EX2.20/20.sce11
-rwxr-xr-x896/CH2/EX2.3/3.sce14
-rwxr-xr-x896/CH2/EX2.4/4.sce25
-rwxr-xr-x896/CH2/EX2.5/5.sce26
-rwxr-xr-x896/CH2/EX2.6/6.sce10
-rwxr-xr-x896/CH2/EX2.7/7.sce13
-rwxr-xr-x896/CH2/EX2.8/8.sce19
-rwxr-xr-x896/CH2/EX2.9/9.sce11
20 files changed, 273 insertions, 0 deletions
diff --git a/896/CH2/EX2.1/1.sce b/896/CH2/EX2.1/1.sce
new file mode 100755
index 000000000..17e908ce9
--- /dev/null
+++ b/896/CH2/EX2.1/1.sce
@@ -0,0 +1,10 @@
+clc
+g=32.2;//ft/s^2
+rho_water=62.3;//lbm/ft^3
+//specific weoight=(density)*(acceleration due to gravity)
+specific_wt=rho_water*g;//lbm.ft/ft^3.s^2
+//1 lbf=32.2 lbm.ft/s^2
+specific_wt=specific_wt/32.2;//lbf/ft^3
+disp("Specific weight of water is")
+disp(specific_wt)
+disp("lbf/ft^3") \ No newline at end of file
diff --git a/896/CH2/EX2.10/10.sce b/896/CH2/EX2.10/10.sce
new file mode 100755
index 000000000..9a9073ab3
--- /dev/null
+++ b/896/CH2/EX2.10/10.sce
@@ -0,0 +1,10 @@
+clc
+//calc thickness of a storage tank
+p_working=250;//lbf/in^2
+//diameter of the cylinder = 10ft = 120in
+d=120;//in
+sigma_tensile=20000;//lbf/in^2
+t=p_working*d/sigma_tensile/2;//in
+disp("Thichness of the storage tank is")
+disp(t)
+disp("in") \ No newline at end of file
diff --git a/896/CH2/EX2.11/11.sce b/896/CH2/EX2.11/11.sce
new file mode 100755
index 000000000..e06943e4b
--- /dev/null
+++ b/896/CH2/EX2.11/11.sce
@@ -0,0 +1,18 @@
+clc
+//calc payload of a helium balloon
+p_atm=1;//atm
+T=293;//K
+d=3;//m (diameter of the balloon)
+//buoyant force=(density of air)*g*(volume of balloon)
+//weight of balloon = (density of helium)*g*(volume of balloon)
+//density for gases = PM/RT
+//payload of balloon = buoyant force - weight
+V_balloon=(%pi)*d^3/6;//m^3
+R=8.2*10^(-2);//m^3.atm/mol/K
+M_air=29;//Kg/Kmol
+M_he=4;//Kg/Kmol
+g=9.81;//m/s^2
+payload=(V_balloon)*g*p_atm*(M_air-M_he)/R/T;//N
+disp("Payload of the balloon is")
+disp(payload)
+disp("N") \ No newline at end of file
diff --git a/896/CH2/EX2.12/12.sce b/896/CH2/EX2.12/12.sce
new file mode 100755
index 000000000..135182d70
--- /dev/null
+++ b/896/CH2/EX2.12/12.sce
@@ -0,0 +1,9 @@
+clc
+//wooden block floating in two phase mix of water and gasoline
+//calc fraction of block in water
+SG_wood=0.96;//Specific gravity
+SG_gasoline=0.72;
+//Let r be the ratio - V_water/V_wood
+r=(SG_wood-SG_gasoline)/(1-SG_gasoline);
+disp("Fraction of wood in water")
+disp(r) \ No newline at end of file
diff --git a/896/CH2/EX2.13/13.sce b/896/CH2/EX2.13/13.sce
new file mode 100755
index 000000000..3304f06e6
--- /dev/null
+++ b/896/CH2/EX2.13/13.sce
@@ -0,0 +1,12 @@
+clc
+//calc gauge pressure of cylinder in a manometer
+//height of water above pt.C = 2.5ft
+rho_water=62.3;//lbm/ft^3;
+h1=2.5;//ft
+rho_gas=0.1;//lbm/ft^3
+h2=0.5;//ft (height of gas)
+g=32.2;//ft/s^2
+gauge_pressure=[(rho_water)*g*h1+(rho_gas)*g*h2]/144/32.2;//lbf/in^2
+disp("Gauge pressure is")
+disp(gauge_pressure)
+disp("lbf/in^2") \ No newline at end of file
diff --git a/896/CH2/EX2.14/14.sce b/896/CH2/EX2.14/14.sce
new file mode 100755
index 000000000..188d2f016
--- /dev/null
+++ b/896/CH2/EX2.14/14.sce
@@ -0,0 +1,14 @@
+clc
+//calc pressure diff between two tanks in a two liquid manometer
+rho_water=62.3;//lbm/ft^3
+SG_oil=1.1;
+rho_oil=SG_oil*(rho_water);
+g=32.2;//ft/s^2
+h1_1=1;//ft
+h1_2=2;//ft
+h2_1=2;//ft
+h2_2=1;//ft
+p_diff=[(rho_water)*g*(h1_1-h1_2)+(rho_oil)*g*(h2_1-h2_2)]/32.2/144;//lbf/in^2
+disp("The pressure difference is")
+disp(p_diff)
+disp("lbf/in^2") \ No newline at end of file
diff --git a/896/CH2/EX2.15/15.sce b/896/CH2/EX2.15/15.sce
new file mode 100755
index 000000000..97376172a
--- /dev/null
+++ b/896/CH2/EX2.15/15.sce
@@ -0,0 +1,9 @@
+clc
+//calc pressure of gauge through a spring piston system
+k=10000;//N/m (spring constant)
+x=0.025;//m (displacement in spring)
+A=0.01;//m^2 (area of piston)
+gauge_pressure=k*x/A/1000;//KPa
+disp("The gauge pressure is")
+disp(gauge_pressure)
+disp("KPa") \ No newline at end of file
diff --git a/896/CH2/EX2.16/16.sce b/896/CH2/EX2.16/16.sce
new file mode 100755
index 000000000..14481525a
--- /dev/null
+++ b/896/CH2/EX2.16/16.sce
@@ -0,0 +1,11 @@
+clc
+//calc pressure diff at the mouth of the fire place
+g=32.2;//ft/s^2
+h=20;//ft (height of fireplace)
+rho_air=0.075;//lbm/ft^3
+T_air=293;//K (surrounding temperature)
+T_fluegas=422;//K
+p_diff=g*h*(rho_air)*[1-(T_air/T_fluegas)]/32.2/144;//lbf/in^2
+disp("The pressure difference is")
+disp(p_diff)
+disp("lbf/in^2") \ No newline at end of file
diff --git a/896/CH2/EX2.17/17.sce b/896/CH2/EX2.17/17.sce
new file mode 100755
index 000000000..1f75b7d02
--- /dev/null
+++ b/896/CH2/EX2.17/17.sce
@@ -0,0 +1,21 @@
+clc
+rho_water=1000;//Kg/m^3
+g=9.81;//m/s^2
+h=5;//m (depth of water)
+//for elevator not accelerated
+p_gauge=(rho_water)*g*h/1000;//KPa
+disp("THe gauge pressure is")
+disp(p_gauge)
+disp("KPa")
+//for elevator accelerated at 5m/s^2 in upward direction
+a=5;//m/s^2
+p_gauge=(rho_water)*(g+a)*h/1000;//KPa
+disp("THe gauge pressure is")
+disp(p_gauge)
+disp("KPa")
+//for elevator accelerated at 5m/s^2 in downward direction
+a=5;//m/s^2
+p_gauge=(rho_water)*(g-a)*h/1000;//KPa
+disp("THe gauge pressure is")
+disp(p_gauge)
+disp("KPa") \ No newline at end of file
diff --git a/896/CH2/EX2.18/18.sce b/896/CH2/EX2.18/18.sce
new file mode 100755
index 000000000..fcc28ed86
--- /dev/null
+++ b/896/CH2/EX2.18/18.sce
@@ -0,0 +1,9 @@
+clc
+//angle free surface makes with the horizontal in an accelerated body
+a=1;//ft/s^2
+g=32.2;//ft/s^2
+theta=atan(a/g);//radians
+theta=theta*180/%pi;//degrees
+disp("The angle made by free surface with the horizontal is")
+disp(theta)
+disp("degrees") \ No newline at end of file
diff --git a/896/CH2/EX2.19/19.sce b/896/CH2/EX2.19/19.sce
new file mode 100755
index 000000000..74c027404
--- /dev/null
+++ b/896/CH2/EX2.19/19.sce
@@ -0,0 +1,10 @@
+clc
+//calc the height to which liq in a cylinder rises when rotated
+f=78/60;//rps
+r=0.15;//m
+g=9.81;//m/s^2
+//omega=2*(%pi)*f
+z=[(2*(%pi)*f)^2]*r^2/2/g;//m
+disp("The liquid in the cylinder rises to a height of")
+disp(z)
+disp("m") \ No newline at end of file
diff --git a/896/CH2/EX2.2/2.sce b/896/CH2/EX2.2/2.sce
new file mode 100755
index 000000000..384adaa61
--- /dev/null
+++ b/896/CH2/EX2.2/2.sce
@@ -0,0 +1,11 @@
+clc
+//calc pressure at depth of 304.9m
+d=304.9;//m
+rho_water=1024;//Kg/m^3
+g=9.81;//m/s^2
+p_atm=101.3;//KPa
+//gauge pressure=(desity)*(acc. due to gravity)*(depth)
+p_depth=p_atm+rho_water*g*d/1000;//KPa
+disp("pressure at the depth is")
+disp(p_depth)
+disp("KPa") \ No newline at end of file
diff --git a/896/CH2/EX2.20/20.sce b/896/CH2/EX2.20/20.sce
new file mode 100755
index 000000000..74a3310f4
--- /dev/null
+++ b/896/CH2/EX2.20/20.sce
@@ -0,0 +1,11 @@
+clc
+//calc thickness of liquid strip at the bottom of the industrial centrifuge
+//Let difference between heights at bottom and top be d
+d=20;//in
+r_a=14;//in
+f=1000/60;//rps
+g=32.2;//ft/s^2
+r_b=[(r_a)^2-2*(d)*g*12/(2*(%pi)*f)^2]^0.5;//in
+disp("The thickness of water strip at bottom of industrial centrifuge")
+disp(r_b)
+disp("in") \ No newline at end of file
diff --git a/896/CH2/EX2.3/3.sce b/896/CH2/EX2.3/3.sce
new file mode 100755
index 000000000..1c289941b
--- /dev/null
+++ b/896/CH2/EX2.3/3.sce
@@ -0,0 +1,14 @@
+clc
+//gauge pressure=(density)*(acc. due to gravity)*(depth)
+rho_oil=55;//lbm/ft^3
+g=32.2;//ft/s^2
+d=60;//ft (depth of oil cylinder)
+gauge_pressure=rho_oil*g*d/32.2;//lbf/ft^2
+disp("Gauge pressure is")
+disp(gauge_pressure)
+disp("lbf/ft^2")
+//1 ft=12 in
+gauge_pressure=gauge_pressure/144;//lbf/in^2
+disp("Gauge pressure is")
+disp(gauge_pressure)
+disp("lbf/in^2") \ No newline at end of file
diff --git a/896/CH2/EX2.4/4.sce b/896/CH2/EX2.4/4.sce
new file mode 100755
index 000000000..e234cba99
--- /dev/null
+++ b/896/CH2/EX2.4/4.sce
@@ -0,0 +1,25 @@
+clc
+//calc of density of air at a certain height
+p_atm=14.7;//psia
+T=289;//K
+//P2=P1*exp^(-(acc. due to gravity)*(mass of air)*(height)/(universal gas const.)/(temp.))
+g=9.81;//m/s^2
+R=8314;//N.m^2/Kmol/K
+//for height of 1000 ft=304.8m
+h=304.8;//m
+p_1000=14.7*exp(-g*29*h/R/289);
+disp("pressure at 1000ft is")
+disp(p_1000)
+disp("psia")
+//for height of 10000 ft=3048m
+h=3048;//m
+p_10000=p_atm*exp(-g*29*h/R/289);
+disp("pressure at 10000ft is")
+disp(p_10000)
+disp("psia")
+//for height of 100000 ft=30480m
+h=30480;//m
+p_100000=14.7*exp(-g*29*h/R/289);
+disp("pressure at 100000ft is")
+disp(p_100000)
+disp("psia") \ No newline at end of file
diff --git a/896/CH2/EX2.5/5.sce b/896/CH2/EX2.5/5.sce
new file mode 100755
index 000000000..03c23ebf3
--- /dev/null
+++ b/896/CH2/EX2.5/5.sce
@@ -0,0 +1,26 @@
+clc
+//calc pressuer at different heights considering on density change in air
+p_atm=14.7;//psia
+g=9.81;//m/s^2
+//P2=P1*[1-(acc. due to gravity)*(mass of air)*(height)/(univ. gas const.)/(temp.)]
+T=289;//K
+R=8314//N.m^2/Kmol/K
+//for height of 1000ft=304.8m
+h=304.8//m
+p_1000=p_atm*[1-g*29*h/R/T];
+disp("pressure at 1000ft is")
+disp(p_1000)
+disp("psia")
+//for height of 10000ft=3048m
+h=3048//m
+p_10000=p_atm*[1-g*29*h/R/T];
+disp("pressure at 10000ft is")
+disp(p_10000)
+disp("psia")
+//for height of 100000ft=30480m
+h=30480//m
+p_100000=p_atm*[1-g*29*h/R/T];
+disp("pressure at 100000ft is")
+disp(p_100000)
+disp("psia")
+//NOTE that the pressure comes out to be negative at 100000ft justifying that density of air changes with altitude \ No newline at end of file
diff --git a/896/CH2/EX2.6/6.sce b/896/CH2/EX2.6/6.sce
new file mode 100755
index 000000000..86af25787
--- /dev/null
+++ b/896/CH2/EX2.6/6.sce
@@ -0,0 +1,10 @@
+clc
+//calc atm pressure on a storage tank roof
+p_atm=14.7;//psia
+//diameter of roof is 120ft
+d_roof=120;//ft
+//force=(pressure)*(area)
+f_roof=p_atm*(%pi)*d_roof^2/4*144;//lbf ;144 because 1ft=12inch
+disp("Force exerted by atmosphere on the roof is")
+disp(f_roof)
+disp("lbf") \ No newline at end of file
diff --git a/896/CH2/EX2.7/7.sce b/896/CH2/EX2.7/7.sce
new file mode 100755
index 000000000..525a8b8e1
--- /dev/null
+++ b/896/CH2/EX2.7/7.sce
@@ -0,0 +1,13 @@
+clc
+//calc atm pressure on a storage tank roof
+p_atm=14.7;//psia
+//diameter of roof is 120ft
+d_roof=120;//ft
+//force=(atm. pressure + gauge pressure)*(area)
+//gauge pressure=(desity)*(acc. due to gravity)*(depth)
+rho_water=62.3//lbm/ft^3
+g=32.2;//ft/s^2
+//depth of water on roof=8 inch=o.667 ft
+h=0.667;//ft
+gauge_pressure=rho_water*g*h/32.2*(%pi)*d_roof^2/4;//lbf
+disp(gauge_pressure)
diff --git a/896/CH2/EX2.8/8.sce b/896/CH2/EX2.8/8.sce
new file mode 100755
index 000000000..b8ba35919
--- /dev/null
+++ b/896/CH2/EX2.8/8.sce
@@ -0,0 +1,19 @@
+clc
+//calc the total force on a lock gate
+//lock gate has water on one side and air on the other at atm. pressure
+w=20;//m (width of the lock gate)
+h=10;//m (height of the lock gate)
+p_atm=1;//atm
+rho_water=1000;//Kg/m^3
+g=9.81//m/s^2
+//for a small strip of dx height at the depth of x on the lock gate
+//net pressure on strip = (p_atm+(rho_water)*g*x) - p_atm
+//thus, net pressure on strip = (rho_water)*g*x
+//force on strip = (rho_water*g*x)*w.dx = (rho_water)*g*w*(x.dx)
+//force on lock gate = integration of force on strip fromm h=0 to h=10
+//integration(x.dx) = x^2/2
+//for h=0 to h=10; integration (x.dx) = h^2/2
+force_lockgate=(rho_water)*g*w*h^2/2;
+disp("The net force on the lock gate is")
+disp(force_lockgate/10^6)
+disp("MN") \ No newline at end of file
diff --git a/896/CH2/EX2.9/9.sce b/896/CH2/EX2.9/9.sce
new file mode 100755
index 000000000..f9f452a9d
--- /dev/null
+++ b/896/CH2/EX2.9/9.sce
@@ -0,0 +1,11 @@
+clc
+//calc thickness of an oil storage
+sigma_tensile=20000;//lbf/in^2 (tensile stress is normally 1/4 rupture stress)
+//max pressure is observed at the bottom of the storage
+p_max=22.9;//lbf/in^2
+//diameter of storaeg tank = 120ft =1440in
+d=1440;//in
+t=(p_max)*d/sigma_tensile/2;//in
+disp("Thichness of the storage tank is")
+disp(t)
+disp("in") \ No newline at end of file
generated by cgit (git 2.34.1) at 2024-12-20 22:37:45 +0000