summaryrefslogtreecommitdiff
path: root/812/CH4
diff options
context:
space:
mode:
Diffstat (limited to '812/CH4')
-rwxr-xr-x812/CH4/EX4.01/4_01.sce19
-rwxr-xr-x812/CH4/EX4.02/4_02.sce20
-rwxr-xr-x812/CH4/EX4.03/4_03.sce9
-rwxr-xr-x812/CH4/EX4.04/4_04.sce24
-rwxr-xr-x812/CH4/EX4.05/4_05.sce16
-rwxr-xr-x812/CH4/EX4.06/4_06.sce10
-rwxr-xr-x812/CH4/EX4.07/4_07.sce16
-rwxr-xr-x812/CH4/EX4.08/4_08.sce7
-rwxr-xr-x812/CH4/EX4.09/4_09.sce8
-rwxr-xr-x812/CH4/EX4.10/4_10.sce19
-rwxr-xr-x812/CH4/EX4.11/4_11.sce16
-rwxr-xr-x812/CH4/EX4.12/4_12.sce11
-rwxr-xr-x812/CH4/EX4.14/4_14.sce15
-rwxr-xr-x812/CH4/EX4.16/4_16.sce11
-rwxr-xr-x812/CH4/EX4.17/4_17.sce10
15 files changed, 211 insertions, 0 deletions
diff --git a/812/CH4/EX4.01/4_01.sce b/812/CH4/EX4.01/4_01.sce
new file mode 100755
index 000000000..e50c6df62
--- /dev/null
+++ b/812/CH4/EX4.01/4_01.sce
@@ -0,0 +1,19 @@
+//Velocity//
+pathname=get_absolute_file_path('4.01.sce')
+filename=pathname+filesep()+'4.01-data.sci'
+exec(filename)
+//If I=integral of(pV.dA):
+//For system: Ics=IA1+IA2+IA3+IA4.
+//For area 1
+IA1=-d*V1*A1
+//For area 3: IA2=d*V3*A3=m3
+IA3=m3
+//For area 4: IA4=-d*V4*A4=-d*Q4
+IA4=-d*Q4
+//For area 2:
+IA2=-IA1-IA3-IA4
+//Velocity at section 2(in ft/sec):
+V2=IA2/d/A2
+//V2 is in the negative y direction
+printf("\n\nRESULTS\n\n")
+printf("\n\nVelocity at section 2: -%.0fj ft/sec\n\n",V2)
diff --git a/812/CH4/EX4.02/4_02.sce b/812/CH4/EX4.02/4_02.sce
new file mode 100755
index 000000000..a209b74e0
--- /dev/null
+++ b/812/CH4/EX4.02/4_02.sce
@@ -0,0 +1,20 @@
+//Mass flow//
+pathname=get_absolute_file_path('4.02.sce')
+filename=pathname+filesep()+'4.02-data.sci'
+exec(filename)
+//If I=integral of(pV.dA):
+//For system: ICS=Iab+Ibc+Icd+Ida
+//But ICS=0
+
+//For Aab:
+function p=f(y),p=-d*U*w*y^0,endfunction
+IAab=intg(0,t,f)
+
+//For Acd:
+function q=g(y),q=d*U*w*(2*y/t-(y/t)^2),endfunction
+IAcd=intg(0,t,g)
+
+//Mass flow rate across surface bc(in kg/sec):
+mbc=(-IAab-IAcd)/1000
+printf("\n\nRESULTS\n\n")
+printf("\n\nMass flow rate across surface bc: %.4f kg/sec\n\n",mbc)
diff --git a/812/CH4/EX4.03/4_03.sce b/812/CH4/EX4.03/4_03.sce
new file mode 100755
index 000000000..d9ad2ace3
--- /dev/null
+++ b/812/CH4/EX4.03/4_03.sce
@@ -0,0 +1,9 @@
+//density//
+pathname=get_absolute_file_path('4.03.sce')
+filename=pathname+filesep()+'4.03-data.sci'
+exec(filename)
+//Rate of change of air density in tank(in (kg/m^3)/s):
+r=-d*v*A/V/10^6
+printf("\n\nRESULTS\n\n")
+printf("\n\nRate of change of air density in tank: %.3f kg/m^3\n\n",r)
+printf("\n\nThe density decreases as is indicated by the negative sign\n\n")
diff --git a/812/CH4/EX4.04/4_04.sce b/812/CH4/EX4.04/4_04.sce
new file mode 100755
index 000000000..40e6d7721
--- /dev/null
+++ b/812/CH4/EX4.04/4_04.sce
@@ -0,0 +1,24 @@
+//Horizontal force//
+pathname=get_absolute_file_path('4.04.sce')
+filename=pathname+filesep()+'4.04-data.sci'
+exec(filename)
+//1) Control Volume selected so that area of left surface is equal to the area of the right surface
+u1=15;
+//Force of support on control volum(in kN):
+function y=f(A),y=-u1*d*V,endfunction
+Rx1=intg(0,0.01,f)
+//Horizontal force on support(in kN):
+Kx=-Rx1
+//2) Control volumes are selected do that the area of the left and right surfaces are equial to the area of the plate
+
+function z=g(A),z=-u1*d*V,endfunction
+Fsx=intg(0,0.01,g)
+//Net force on plate:Fx=0=-Bx-pa*Ap+Rx
+// Rx=pa*Ap+Bx
+//From the above, it is obtained that:
+Rx2=-2.25
+//Horizontal force on support(in kN):
+Kx2=-Rx2
+printf("\n\nRESULTS\n\n")
+printf("\n\nHorizontal force on support: %.3f kN\n\n",Kx/1000)
+
diff --git a/812/CH4/EX4.05/4_05.sce b/812/CH4/EX4.05/4_05.sce
new file mode 100755
index 000000000..27c799ef0
--- /dev/null
+++ b/812/CH4/EX4.05/4_05.sce
@@ -0,0 +1,16 @@
+//Scale//
+pathname=get_absolute_file_path('4.05.sce')
+filename=pathname+filesep()+'4.05-data.sci'
+exec(filename)
+//Weight of water in the tank(in lbf):
+d1=62.4;
+WH2O=d1*A*h
+v=-5;
+//Total body force in negative y direction(lbf):
+function y=f(A),y=-v*d2*V1,endfunction
+F=intg(0,A1,f)
+//Force of scale on control volume(in kN):
+Ry=W+WH2O-F
+printf("\n\nRESULTS\n\n")
+printf("\n\nScale Reading: %.3f lbf\n\n",Ry)
+
diff --git a/812/CH4/EX4.06/4_06.sce b/812/CH4/EX4.06/4_06.sce
new file mode 100755
index 000000000..8220dc73e
--- /dev/null
+++ b/812/CH4/EX4.06/4_06.sce
@@ -0,0 +1,10 @@
+//force exerted per unt//
+pathname=get_absolute_file_path('4.06.sce')
+filename=pathname+filesep()+'4.06-data.sci'
+exec(filename)
+//X-component of reaction force per unit width of the gate(in N/m):
+Rxw=(d*(V2^2*D2-V1^2*D1))-(d*g/2*(D1^2-D2^2))
+//Horizontal force exerted per unt width on the gate(in N/m):
+Kxw=-Rxw
+printf("\n\nRESULTS\n\n")
+printf("\n\nHorizontal force exerted per unt width on the gate: %.3f kN/m\n\n",Kxw/1000)
diff --git a/812/CH4/EX4.07/4_07.sce b/812/CH4/EX4.07/4_07.sce
new file mode 100755
index 000000000..6f0000427
--- /dev/null
+++ b/812/CH4/EX4.07/4_07.sce
@@ -0,0 +1,16 @@
+//Force to hold//
+pathname=get_absolute_file_path('4.07.sce')
+filename=pathname+filesep()+'4.07-data.sci'
+exec(filename)
+//Velocity at section 1(in m/sec):
+V1=V2*A2/A1
+//Gauge pressure(in kPa):
+p1g=p1-patm
+u1=V1;u2=-V2;
+//Reaction force component in the x direction(in N):
+Rx=-p1g*A1-u1*d*V1*A1
+//Reaction force component in the y direction(in N):
+Ry=u2*d*V2*A2
+printf("\n\nRESULTS\n\n")
+printf("\n\nForce to hold elbow acting to the left: %.3f kN\n\n",Rx/1000)
+printf("\n\nForce to hold elbow acting downwards: %.3f N\n\n",Ry)
diff --git a/812/CH4/EX4.08/4_08.sce b/812/CH4/EX4.08/4_08.sce
new file mode 100755
index 000000000..fd8e5eddf
--- /dev/null
+++ b/812/CH4/EX4.08/4_08.sce
@@ -0,0 +1,7 @@
+//Tension//
+filename=pathname+filesep()+'4.08-data.sci'
+exec(filename)
+//Tension required to pull the belt(in lbf):
+T=Vbelt*m/32.2
+printf("\n\nRESULTS\n\n")
+printf("\n\nTension required to pull the belt: %.3f lbf\n\n",T)
diff --git a/812/CH4/EX4.09/4_09.sce b/812/CH4/EX4.09/4_09.sce
new file mode 100755
index 000000000..3415f9e67
--- /dev/null
+++ b/812/CH4/EX4.09/4_09.sce
@@ -0,0 +1,8 @@
+//pressure required//
+pathname=get_absolute_file_path('4.09.sce')
+filename=pathname+filesep()+'4.09-data.sci'
+exec(filename)
+//Minimum gauge pressure required(in lbf/in^2):
+pg=8/%pi^2*d/D1^4*Q^2*((D1/D2)^4-1)*144
+printf("\n\nRESULTS\n\n")
+printf("Minimum gauge pressure required: %.3f lbf/in^2",pg)
diff --git a/812/CH4/EX4.10/4_10.sce b/812/CH4/EX4.10/4_10.sce
new file mode 100755
index 000000000..39073d8e6
--- /dev/null
+++ b/812/CH4/EX4.10/4_10.sce
@@ -0,0 +1,19 @@
+//Net force//
+pathname=get_absolute_file_path('4.10.sce')
+filename=pathname+filesep()+'4.10-data.sci'
+exec(filename)
+u1=V-U
+u2=(V-U)*cosd(theta)
+v2=(V-U)*sind(theta)
+V1=V-U
+V2=V1
+//X component of moment equation(in N):
+function y=f(A),y=u1*-(d*V1),endfunction
+function z=g(A),z=u2*d*V2,endfunction
+Rx=intg(0,A,f)+intg(0,A,g)
+
+//Y component of moment equation(in N):
+function a=h(A),a=v2*d*V1,endfunction
+Ry=intg(0,A,h) //This is after neglecting weight of vane and the water.
+printf("\n\nRESULTS\n\n")
+printf("\n\nNet force on the vane: %.3f i+%.2f j kN\n\n",Rx/1000,Ry/1000)
diff --git a/812/CH4/EX4.11/4_11.sce b/812/CH4/EX4.11/4_11.sce
new file mode 100755
index 000000000..718c78372
--- /dev/null
+++ b/812/CH4/EX4.11/4_11.sce
@@ -0,0 +1,16 @@
+//PLOTTING//
+pathname=get_absolute_file_path('4.11.sce')
+filename=pathname+filesep()+'4.11-data.sci'
+exec(filename)
+//Evaluating the value of Vb:
+Vb=V*(1-cosd(theta))*d*A/M
+//Value of U/V for various values of t
+t=0:20;
+[m n]=size(t)
+for i=1:n
+ U_V(i)=Vb*t(i)/(1+Vb*t(i));
+end
+
+//Plotting U/V vs t:
+plot(t,U_V)
+xtitle('U/V vs t','t (in sec)','U/V')
diff --git a/812/CH4/EX4.12/4_12.sce b/812/CH4/EX4.12/4_12.sce
new file mode 100755
index 000000000..86464c549
--- /dev/null
+++ b/812/CH4/EX4.12/4_12.sce
@@ -0,0 +1,11 @@
+//Velocity of rocket//
+pathname=get_absolute_file_path('4.12.sce')
+filename=pathname+filesep()+'4.12-data.sci'
+exec(filename)
+//Acceleration of rocket at t=0(in m/sec^2):
+Ve*me/M0-g
+//Velocity of rocket at t=10 (in m/sec):
+function y=f(t),y=Ve*me/(M0-me*t)-g,endfunction
+Vcv=intg(0,t,f)
+printf("\n\nRESULTS\n\n")
+printf("\n\nVelocity of rocket at t=10: %.1f m/sec\n\n",Vcv)
diff --git a/812/CH4/EX4.14/4_14.sce b/812/CH4/EX4.14/4_14.sce
new file mode 100755
index 000000000..755aa6047
--- /dev/null
+++ b/812/CH4/EX4.14/4_14.sce
@@ -0,0 +1,15 @@
+//Relative speed and friction//
+pathname=get_absolute_file_path('4.14.sce')
+filename=pathname+filesep()+'4.14-data.sci'
+exec(filename)
+//Area of jet(in mm^2):
+Ajet=%pi/4*D^2
+//Jet speed relative to the nozzle(in m/sec):
+Vrel=Q/2/Ajet*10^6/60/1000
+//Value of w*R in m/sec:
+wR=w*R*2*%pi/60/1000
+//Friction torque at pivot(in N-m):
+Tf=R*(Vrel*cosd(alpha)-wR)*d*Q/1000/60/1000
+printf("\n\nRESULTS\n\n")
+printf("\n\nJet speed relative to each nozzle: %.2f m/sec\n\n",Vrel)
+printf("\n\nFriction torque at pivot: %.5f N-m\n\n",Tf)
diff --git a/812/CH4/EX4.16/4_16.sce b/812/CH4/EX4.16/4_16.sce
new file mode 100755
index 000000000..4feb7092b
--- /dev/null
+++ b/812/CH4/EX4.16/4_16.sce
@@ -0,0 +1,11 @@
+//Rate of heat//
+pathname=get_absolute_file_path('4.16.sce')
+filename=pathname+filesep()+'4.16-data.sci'
+exec(filename)
+//Velocity at exit(in ft/sec):
+V2=m*R*(T2+460)/A2/p2/144
+//As power input is to CV, Ws=-600
+//Rate of heat transfer(in Btu/sec):
+Q=Ws*550/778+m*cp*(T2-T1)+m*V2^2/2/32.2/778
+printf("\n\nRESULTS\n\n")
+printf("\n\nRate of heat transfer: %.3f Btu/sec\n\n",Q)
diff --git a/812/CH4/EX4.17/4_17.sce b/812/CH4/EX4.17/4_17.sce
new file mode 100755
index 000000000..38b8f55c8
--- /dev/null
+++ b/812/CH4/EX4.17/4_17.sce
@@ -0,0 +1,10 @@
+//Mass flow rate//
+pathname=get_absolute_file_path('4.17.sce')
+filename=pathname+filesep()+'4.17-data.sci'
+exec(filename)
+//Density of tank(in kg/m^3):
+d=(p1+patm)/R/T
+//Mass flow rate of air into the tank(in kg/sec):
+m=d*V*cv*r/R/T*1000
+printf("\n\nRESULTS\n\n")
+printf("\n\nMass flow rate of air into the tank: %.3f g/sec\n\n",m)
generated by cgit (git 2.34.1) at 2025-01-06 12:51:41 +0000