Introduction to Basic Principles of Fluid Mechanics

2.016 Hydrodynamics Reading #3 2.016 Hydrodynamics Prof. A.H. Techet Introduction to basic principles of fluid mechanics I. Flow Descriptions 1. Lagrangian (following the particle): In rigid body mechanics the motion of a body is described in terms of the body’s position in time. This body can be translating and possibly rotating, but not deforming. This description, following a particle in time, is a Lagrangian description, with velocity vector V = ui +v j + wz . (3.1) Using the Lagrangian approach, we can describe a particle located at point xo = ( xo , yo , zo ) for some time t = to, such that the particle velocity is V= ∂x , ∂t (3.2) a= ∂V . ∂t (3.3) and particle acceleration is We can use Newton’s Law of motion ( F = ma ) on the body to determine the acceleration and thus, the velocity and position. However, in fluid mechanics, it is difficult to track a single fluid particle. But in the lab we can observe many particles passing by one single location. 2. Eulerian (observing at one location): In the lab, we can easily observe many particles passing a single location, and we can make measurements such as drag on a stationary model as fluid flows past. Thus it is useful to use the Eulerian description, or control volume approach, and describe the flow at every fixed point in space (x, y, z) as a function of time, t . version 3.0 updated 8/30/2005 -1- ©2005 A. Techet 2.016 Hydrodynamics Reading #3 z x w u Figure 1: An Eulerian description gives a velocity vector at every point in x,y,z as a function of time. In an Eulerian velocity field, velocity is a function of the position vector and time, V ( x, t ) . For example: V ( x, t ) = 6tx 2 i + 3 zy j + 10 xyt 2 z 3. Reynolds Transport Theorem (the link between the two views): In order to apply Newton’s Laws of motion to a control volume, we need to be able to link the control volume view to the motion of fluid particles. To do this, we use the Reynolds Transport Theorem, which you’ll derive in graduate fluids classes, like 2.25. Suffice it to say that the theorem exists. For this class, we’ll use control volumes to describe fluid motion. 4. Description of Motion: Streamlines: Line everywhere tangent to velocity (Eulerian) (No velocity exists perpendicular to the streamline!) Streaklines: instantaneous loci of all fluid particles that pass through a given point xo. version 3.0 updated 8/30/2005 -2- ©2005 A. Techet 2.016 Hydrodynamics Reading #3 Particle Pathlines: Trajectory of fluid particles (“more” lagrangian) In steady flow stream, streak, and pathlines are identical!! (Steady flow has no time dependence.) II. Governing Laws The governing laws of fluid motion can be derived using a control volume approach. This is equivalent to a “fluidic black box” where all we know about the flow is what is going in and what is coming out of the control volume: mass, momentum, and energy. The control volume (CV) can be fixed or move with the fluid. For simplicity it is often ideal to fix the CV, but this does not always provide the easiest solution in all cases. For most of this class the CV will be fixed. When analyzing a control volume problem there are three laws that are always true: 1. Conservation of Mass 2. Conservation of Momentum 3. Conservation of Energy 1. Conservation of Mass: Basic fluid mechanics laws dictate that mass is conserved within a control volume for constant density fluids. Thus the total mass entering the control volume must equal the total mass exiting the control volume plus the mass accumulating within the control volume. mass in – mass out = mass accumulating min − mout = macc (3.4) Let us consider three cases: version 3.0 updated 8/30/2005 -3- ©2005 A. Techet 2.016 Hydrodynamics Reading #3 Case I: In control volume 1, water is poured into a tub with a drain. The mass flow rate into the tub is min = ρVin A in . Similarly, the mass flow rate out of the tub is mout = ρVout A out . If the mass flow rate in is greater than the mass flow rate out, water will accumulate in the tub. If the mass flow rate out is greater, then the tub will drain. If the two are equal, then no water will accumulate in the tub. Think about it next time you’re in the shower! Case II: Control Volume 2 is a section of a pipe full of water. Since the CV is full, the only way more mass can accumulate is if it becomes denser. Remember, in this class, we treat water as incompressible, so the density cannot change and we never have mass accumulating in a full control volume. Thus, the mass flow rate in equals the mass flow rate out. Furthermore, since CV2 is drawn in a pipe of constant area, then the velocity in must equal the velocity out. For an incompressible fluid, there is no change in velocity through a pipe of constant area! version 3.0 updated 8/30/2005 -4- ©2005 A. Techet 2.016 Hydrodynamics Reading #3 Case III: Let’s now consider a general control volume immersed in a fluid: We can write a 2D mass balance equation for the fluid entering and exiting the control volume. min − mout = macc (3.5) ρudydz + ρvdxdz − ρ (u + −ρ Which simplifies to, ∂u ∂v dx)dydz − ρ (v + dy)dxdz = 0 ∂x ∂y (3.6) ∂u ∂v dxdydz − ρ dydxdz = 0 ∂x ∂y (3.7) ∂u ∂v + =0 ∂x ∂y (3.8) In three dimensions, the derivation is the same, and we have: ∂u ∂v ∂w + + =0 ∂x ∂y ∂z (3.9) or in vector notation, recalling the gradient operator: ∇ = ( ∂x∂ , ∂y∂ , ∂z∂ ) , we have ∇ ⋅V = 0 version 3.0 updated 8/30/2005 -5- (3.10) ©2005 A. Techet 2.016 Hydrodynamics Reading #3 2. Conservation of Momentum: Newton’s second law is simply the law of conservation of momentum. It states that the time rate of change of momentum of a system of particles is equal to the sum of external forces acting on that body. ΣFi = d {MV} dt (3.11) where M = ρδ xδ z is the mass of the fluid parcel (in two dimensions, ie mass per unit length) and MV is the linear momentum of the system (V is the velocity vector). Since the fluid density is constant, the time-rate of change of linear momentum can be written as d dV {MV} = ρδ xδ z . dt dt (3.12) The rate of change of velocity of the fluid parcel can be found, for small δ t , as dV lim 1 = {V( x + δ x p , z + δ z p , t + δ t p ) − V( x, z, t )} dt δ t → 0 δ t (3.13) We can substitute, δ x p = uδ t , and, δ z p = wδ t , into equation (3.13) and cancel terms to arrive at a more familiar form of the momentum equation. The total derivative of the velocity is written as: DV = ∂V + ∂V u + ∂V w ∂z ∂t ∂x Dt (3.14) which can be simplified using the vector identity, V ⋅∇ = u ∂ + v ∂ + w ∂ ∂x ∂y ∂z version 3.0 updated 8/30/2005 -6- (3.15) ©2005 A. Techet 2.016 Hydrodynamics Reading #3 The total (material) derivative of the velocity is the sum of the conventional acceleration, ∂∂tV , and the advection term, (V ⋅∇)V : DV ∂V = + (V ⋅∇)V. ∂t Dt (3.16) Finally, the momentum equation, from equation (3.11), can be rewritten in two dimensions as ΣFi = ρ DV δ xδ z . Dt (3.17) 3. Forces The LHS of equation (3.17) is the sum of the forces acting on the control volume. Contributions from gravity and pressure both play a role in this term as well as any applied external forces. These forces are found as follows: 1. Force on a fluid volume due to gravity: Fg = −( ρ gδ xδ z )kˆ (3.18) 2. Pressure Forces due: Force due to pressure is simply pressure times the surface area it acts on FP = P ⋅ A . (3.19) 1 ∂p ⎞ 1 ∂p ∂p ⎞ ∂p ⎛ ⎛ FPx = ⎜ p + δ z ⎟δ z − ⎜ p + δ z + δ x ⎟ δ z = − δ xδ z ∂x 2 ∂z ⎠ 2 ∂z ∂x ⎠ ⎝ ⎝ (3.20) The pressure force in x-direction is and the pressure force in z-direction is FPz = − version 3.0 updated 8/30/2005 ∂p δ xδ z . ∂z -7- (3.21) ©2005 A. Techet 2.016 Hydrodynamics Reading #3 Thus the total pressure force in two dimensions is FP = −( ∂p ∂p , )δ xδ z = −∇pδ xδ z. ∂x ∂z (3.22) 4. Euler Equation Substituting relations (3.18) and (3.22) for the gravity and pressure forces acting on the body, into the momentum equation (3.17) we arrive at ⎧ ∂V ⎫ + (V ⋅∇)V ⎬ δ xδ z = (− ρ gδ xδ z )kˆ − ∇pδ xδ z ⎩ ∂t ⎭ ρ⎨ (3.23) for any δ x, δ z . The final result is the Euler equation in vector form: ⎧ ∂V ⎫ + (V ⋅∇)V ⎬ = − ρ gk̂ − ∇p. ⎩ ∂t ⎭ ρ⎨ (3.24) We can further manipulate this equation with the vector identity (V ⋅∇)V = 12 ∇(V ⋅ V ) , (3.25) such that the Euler equation becomes ⎧ ∂V 1 ⎫ + ∇(V ⋅ V) ⎬ = − ρ gk̂ − ∇p. ⎩ ∂t 2 ⎭ ρ⎨ version 3.0 updated 8/30/2005 -8- (3.26) ©2005 A. Techet 2.016 Hydrodynamics Reading #3 5. Bernoulli’s Equation: Application of Newton’s Second Law along a streamline: p1 + 1 2 ρV12 + ρ gz1 = p2 + 1 2 ρV22 + ρ gz2 = C (3.27) Assuming the following conditions: 1) Points 1 and 2 are on the same streamline! 2) Fluid density is constant dV = 0 (no time dependence or turbulence) 3) Flow is steady: dt 4) Fluid is “inviscid” or can be approximated as inviscid. No frictional effects 5) No Work Added! We can derive this through a Lagrangian derivation: Looking at a small elemental volume along a streamline d∀ = dn ds dx (dx is the depth into the paper). Fluid weight in the (-z) direction.: ρ g dn ds dx (3.28) Component of weight acting in the s-direction: − ρ g sin β dn ds dx Where sin β = version 3.0 (3.29) dz so that the weight in the s-direction is: ds updated 8/30/2005 -9- ©2005 A. Techet 2.016 Hydrodynamics Reading #3 −ρ g dz dn ds dx . ds (3.30) The force due to pressure in the s-direction is found similarly: ∂p ⎞ ∂p ds ⎞ ∂z ⎞ ⎛ ⎛ ∂p ⎛ Fs = ⎜ p − ⎟ dn dx − ⎜ p + ⎟ dn dx = ⎜ − − ρ g ⎟ dn ds dx ∂s ⎠ ∂s 2 ⎠ ∂s ⎠ ⎝ ⎝ ⎝ ∂s (3.31) The force accelerates the fluid along the streamline such that the rate of change in momentum, per unit volume, is ∂V ⎛ ⎜ V + ∂s ds − V ρ⎜ dt ⎜ ⎝ where V = ⎞ ⎟ ∂V ⎟ = ρV ∂s ⎟ ⎠ (3.32) ∂s . ∂t So Euler’s equation in one dimension along a streamline becomes: ρV ∂V ∂p ∂z + + ρg = 0 . ∂s ∂s ∂s (3.33) ∂p ds ∂s ∂V ds Change in Velocity along a streamline: dV = ∂s ∂z Change in height along a streamline: dz = ds ∂s Change in Pressure along a streamline: dp = Multiplying equation 23 through by ds gives us ρV dV + dp + ρ g dz = 0 dp ρ + V dV + g dz = 0 (3.34) (3.35) If density is constant along the streamline then we can integrate along the streamline to get: version 3.0 updated 8/30/2005 -10- ©2005 A. Techet 2.016 Hydrodynamics Reading #3 p ρ + 12 V 2 + g z = C Along a streamline Bernoulli’s equation relates pressure, height and velocity at two points: p1 ρ + 1 2 V12 + g z1 = p2 ρ + 1 2 V2 2 + g z2 = C (3.36) This equation also assumes that NO additional heat or work is added to the system along the streamline. 6. Irrotational flow An irrotational flow is defined as a flow for which each and every fluid particle is not rotating. Mathematically speaking, the curl of the velocity is identically zero. iˆ ∂ ω = ∇×V = ∂x u ω = iˆ( ĵ ∂ ∂y v k̂ ∂ =0 ∂z w ∂w ∂v ∂u ∂w ∂v ∂u − ) + ĵ( − ) + k̂( − ) = 0. ∂y ∂z ∂z ∂x ∂x ∂y (3.37) (3.38) For 2D flow this reduces to ∂u = ∂w . ∂z ∂x version 3.0 updated 8/30/2005 -11- ©2005 A. Techet