summaryrefslogtreecommitdiff
path: root/Fluid_Mechanics_/Chapter8.ipynb
diff options
context:
space:
mode:
Diffstat (limited to 'Fluid_Mechanics_/Chapter8.ipynb')
-rw-r--r--Fluid_Mechanics_/Chapter8.ipynb33
1 files changed, 0 insertions, 33 deletions
diff --git a/Fluid_Mechanics_/Chapter8.ipynb b/Fluid_Mechanics_/Chapter8.ipynb
index 51f3481e..6972a7f7 100644
--- a/Fluid_Mechanics_/Chapter8.ipynb
+++ b/Fluid_Mechanics_/Chapter8.ipynb
@@ -28,9 +28,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "# Determine maximum velocity and shear stress\n",
"\n",
- "# Given\n",
"\n",
"from math import *\n",
"\n",
@@ -46,7 +44,6 @@
"\n",
"g = 9.81 # acceleration due to gravity in m/s**2\n",
"\n",
- "# properties of kerosene\n",
"\n",
"mu = 19.1*10**-4 # viscosity of kerosene at 20 deg C\n",
"\n",
@@ -54,17 +51,13 @@
"\n",
"rho = 1000 # density in kg/m**3\n",
"\n",
- "# Solution\n",
"\n",
- "# calculating direction of flow\n",
"\n",
"p1 = (P1+g*z*S)*1000 # point 1\n",
"\n",
"p2 = (P2)*1000 # point 2\n",
"\n",
- "# direction of flow is from point 1-2\n",
"\n",
- "# shear stress\n",
"\n",
"Sp = -((p1-p2)/sqrt(L**2+z**2))\n",
"\n",
@@ -74,19 +67,16 @@
"\n",
"print \"(a) Maximum shear stress =\",round(Tau_max,3),\"N/m**2\"\n",
"\n",
- "# maximum velocity\n",
"\n",
"Vmax = r**2*Sp/(4*mu)\n",
"\n",
"print \"(b) Maximum velocity =\",round(Vmax,3),\"m/s\"\n",
"\n",
- "# discharge\n",
"\n",
"Q = pi*r**4*Sp/(8*mu)\n",
"\n",
"print \"(c) Discharge = \",round(Q,7),\"m**3/s\"\n",
"\n",
- "# calculate reynolds number\n",
"\n",
"V = Vmax/2\n",
"\n",
@@ -122,9 +112,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "# Determine the head loss\n",
"\n",
- "# Given\n",
"\n",
"d = 0.02 # diameter of the pipe in m\n",
"\n",
@@ -134,11 +122,9 @@
"\n",
"g = 9.81 # acceleration due to gravity in m/s**2\n",
"\n",
- "# for water at 5 deg C\n",
"\n",
"nu = 1.54*10**-6 # kinematic viscosity of water in m**2/s\n",
"\n",
- "# Solution\n",
"\n",
"R = v*d/nu\n",
"\n",
@@ -178,13 +164,10 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "# Horsepower required to pump 50 tons of oil\n",
"\n",
- "# Given\n",
"\n",
"from math import *\n",
"\n",
- "# oil properties\n",
"\n",
"S = 0.92 # specific gravity\n",
"\n",
@@ -198,7 +181,6 @@
"\n",
"g = 32.2 # acceleration due to gravity in ft/s**2\n",
"\n",
- "# Solution\n",
"\n",
"Q = W*2000/(gma*3600) # discharge in ft**3/s\n",
"\n",
@@ -244,9 +226,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "# Viscosity of the liquid in poise\n",
"\n",
- "# Given\n",
"\n",
"from math import *\n",
"\n",
@@ -268,9 +248,7 @@
"\n",
"a = pi*d1**2/4\n",
"\n",
- "# Solution\n",
"\n",
- "# From derivation we obtain a equation for T\n",
"\n",
"H2 = H1-(V/A)\n",
"\n",
@@ -303,13 +281,10 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "# Velocity distribution; Discharge ; shear on the upper plate\n",
"\n",
- "# Given\n",
"\n",
"from math import *\n",
"\n",
- "# properties of kerosene oil at 20 deg C\n",
"\n",
"S = 0.81 # specific gravity of oil\n",
"\n",
@@ -327,35 +302,27 @@
"\n",
"theta = pi/6 # angle of inclination\n",
"\n",
- "# Solution\n",
"\n",
- "# point 1\n",
"\n",
"P1 = p1*144 + gma*l*sin(theta)\n",
"\n",
- "# point 2\n",
"\n",
"P2 = p2*144\n",
"\n",
- "# flow is taking from poont 2-1\n",
"\n",
"Sp = (P2-P1)/4\n",
"\n",
- "# equation for u = 2154.75*y-359125*y**2\n",
"\n",
"y = h\n",
"\n",
- "# discharge per ft width\n",
"\n",
"q = (2154.75*y**2/2) - (359125*y**3/3)\n",
"\n",
"print \"Discharge q = \",round(q,3),\"per unit ft of the plate\"\n",
"\n",
- "# to find shear at the top of the plate take du/dy = 0\n",
"\n",
"dV = 2154.75 - 718250*h\n",
"\n",
- "# shear stress\n",
"\n",
"T = -mu*dV\n",
"\n",