diff options
Diffstat (limited to 'Fluid_Mechanics_/Chapter1.ipynb')
| -rw-r--r-- | Fluid_Mechanics_/Chapter1.ipynb | 49 |
1 files changed, 0 insertions, 49 deletions
diff --git a/Fluid_Mechanics_/Chapter1.ipynb b/Fluid_Mechanics_/Chapter1.ipynb index e770487c..2869f384 100644 --- a/Fluid_Mechanics_/Chapter1.ipynb +++ b/Fluid_Mechanics_/Chapter1.ipynb @@ -28,9 +28,7 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Density of air\n", "\n", - "#Given \n", "\n", "Mw = 29.0 # Molecular weight of air\n", "\n", @@ -40,7 +38,6 @@ "\n", "p = 50*144*47.88 # Pressure in N/m**2\n", "\n", - "# Solution\n", "\n", "rho = p/(R*T) # from the state law\n", "\n", @@ -71,9 +68,7 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Reduction in volume\n", "\n", - "# Given\n", "\n", "dP = 10**6 # Pressure drop in N/m**2\n", "\n", @@ -81,7 +76,6 @@ "\n", "bta = 2.2*10**9 # Bulk modulus of elasticity in N/m**2\n", "\n", - "# Solution\n", "\n", "dV = -dP*V/bta # Change in volume in m**3\n", "\n", @@ -114,9 +108,7 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Volume reduction\n", "\n", - "# Given\n", "\n", "bta1 = 2.28*10**9 # Bulk modulus of elasticity at 20 deg C and 103.4 N/m**2\n", "\n", @@ -126,7 +118,6 @@ "\n", "p2 = 1034 # Pressure in N/m**2\n", "\n", - "# Solution \n", "\n", "bavg = (bta1+bta2)/2 # bulk modulus average in N/m**2\n", "\n", @@ -168,9 +159,7 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Volume reduction \n", "\n", - "# Given\n", "\n", "patm = 14.6 # Atmospheric pressure in psia\n", "\n", @@ -179,7 +168,6 @@ "p2 = 102 # gauge pressure at point 2 in psia\n", "\n", "V = 1 # volume in m**3\n", - "# solution\n", "\n", "p = p1+patm # absolute pressure in psia\n", "\n", @@ -224,11 +212,9 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Sonic velocity of air\n", "\n", "from math import *\n", "\n", - "# Given\n", "\n", "k = 1.4 # gas constant\n", "\n", @@ -236,7 +222,6 @@ "\n", "T = 68+460 # temperature in *oR\n", "\n", - "# solution\n", "\n", "c = sqrt(k*R*T)\n", "\n", @@ -267,11 +252,9 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Force required to move the piston\n", "\n", "from math import *\n", "\n", - "# Given\n", "\n", "d = 0.05 # diameter of cylinder 1 in m\n", "\n", @@ -287,7 +270,6 @@ "\n", "A = pi*l*d # area in m**2\n", "\n", - "# Solution \n", "\n", "tau = mu*U/Y # Shear stress in N/m**2\n", "\n", @@ -320,27 +302,19 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Distance between the walls; Shear Stress; Location of maximum velocity\n", "\n", "from math import *\n", "\n", "from sympy import *\n", "\n", - "# Given\n", "\n", "mu = 1.005*10**-3 # Viscosity of water in Ns/m**2\n", "\n", - "# Solution \n", "\n", - "# Part a\n", "\n", - "# Velocity is given by the formula u = 10*(0.01*y-y**2)\n", "\n", - "# two boundary conditions must be satisfied\n", "\n", - "# at y=0;u=0 at the bottom of the plate\n", "\n", - "# at y=Y ; u = 0 at top of the plate\n", "\n", "Y = 0.01 # Distance between the walls\n", "\n", @@ -348,11 +322,8 @@ "\n", "print \" (a) Distance between the walls = \",round(Y1,1),\"cm\"\n", "\n", - "# Part b\n", "\n", - "# tau = mu*du/dy # Newtons law of viscosity\n", "\n", - "# differentiate u wrt y\n", "\n", "y = Symbol('y')\n", "\n", @@ -364,8 +335,6 @@ "\n", "U = uprime\n", "\n", - "#print U\n", - "# for y =0 at the bottom plate we get \n", "\n", "U1 = 0.01 # from U\n", "\n", @@ -373,19 +342,14 @@ "\n", "print \" (b) Shear stress = \",round(tau,9),\"N/m**2\"\n", "\n", - "# Part c\n", "\n", - "# Shaer stress at 20um from the wall\n", "\n", "tau1 = mu*10*(0.01-2*20*10**-6) # using the equation of U and y = 20*10**-6 calc shear stress in N/m**2\n", "\n", "print \" (c) Shear Stress at 20 um from the plate = \",round(tau1,9),\"N/m**2\"\n", "\n", - "# Part D\n", "\n", - "# Distance at which shaer stress is zero can be found from the location of maximum velocity \n", "\n", - "# equating uprime = 0\n", "\n", "y1 = 0.01/2 # shear stress location\n", "\n", @@ -423,9 +387,7 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Shaft torque for linear and non linear distribution of velocity \n", "from math import *\n", - "# Given\n", "\n", "d = 0.1 # diameter of shaft in m\n", "\n", @@ -443,9 +405,7 @@ "\n", "mu = 0.44 # Viscosity of SAE-30 oil in Ns/m**2\n", "\n", - "# Solution\n", "\n", - "# part a\n", "\n", "F = 2*pi*r1*l*mu*U/t\n", "\n", @@ -453,7 +413,6 @@ "\n", "print \" (a) For linear distribution of velocity , shaft torque = \",round(T,2),\"m.N\"\n", "\n", - "# part b\n", "\n", "F1 = 2*pi*l*mu*U/log(r2/r1)\n", "\n", @@ -488,9 +447,7 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Watts of energy lost to overcome friction\n", "from math import *\n", - "# Given\n", "\n", "mu = 0.44 # viscosity of the oil in Ns/m**2\n", "\n", @@ -533,15 +490,12 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Excessive pressure inside droplet\n", "\n", - "# Given\n", "\n", "d = 0.01 # diameter in m\n", "\n", "sigma = 0.073 # surface tension in N/m\n", "\n", - "# Solution \n", "\n", "dP = 4*sigma/d # pressure excessive\n", "\n", @@ -572,9 +526,7 @@ "cell_type": "code", "collapsed": false, "input": [ - "# Height to which alcohol will rise\n", "\n", - "# Given\n", "\n", "sigma = 0.022 # surface tension in N/m\n", "\n", @@ -584,7 +536,6 @@ "\n", "d = 0.002 # diameter in m\n", "\n", - "# Solution \n", "\n", "\n", "h =4*sigma*1000/(gma*S*d) # capillary height in m\n", |