* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download Review Notes 2
Survey
Document related concepts
Flow conditioning wikipedia , lookup
Magnetorotational instability wikipedia , lookup
Wind-turbine aerodynamics wikipedia , lookup
Airy wave theory wikipedia , lookup
Aerodynamics wikipedia , lookup
Cnoidal wave wikipedia , lookup
Lattice Boltzmann methods wikipedia , lookup
Reynolds number wikipedia , lookup
Euler equations (fluid dynamics) wikipedia , lookup
Computational fluid dynamics wikipedia , lookup
Bernoulli's principle wikipedia , lookup
Navier–Stokes equations wikipedia , lookup
Fluid dynamics wikipedia , lookup
Derivation of the Navier–Stokes equations wikipedia , lookup
Transcript
Equations of Change for Isothermal Systems • In the previous lecture, we showed how to derive the velocity distribution for simple flows by the application of the shell momentum balance or the force balance. • It is however more reliable to start with general equations for – the conservation of mass (continuity equation) – the conservation of momentum (equation of motion, N2L) to describe any flow problem and then simplify these equations for the case at hand. • For non-isothermal fluids (heat transfer + flow problems), the same technique can be applied combined with the use of the equation for the conservation of energy Time Derivatives Any quantity c which depends on time as well as position can be written as c=f(t, x, y, z) then it differential dc is: c c c c dt + dx + dz dc = dy + t x z y • Partial Time Derivative c t The partial time derivative is the derivative of the function c with time holding x, y, z constant. (fixed observer) dc • Total Time Derivative dt The Total Time Derivative accounts for the fact that the observer is moving (how c varies with t because of changing location). dc c c dx c dy c dz = + + + dt t x dt y dt z dt • Substantial Time Derivative Dc Dt The substantial Time Derivative is a particular case of the Total Time Derivative for which the velocity v of the observer is the same as the velocity of the flow. It is also called the Derivative Following the Motion c Dc c c c = + v x + vy + vz Dt t x z y vx v = v where the local fluyid velocity is defined by: y v z • Analogy of the river, the fisherman and the boat... Equation of Continuity • Write the mass balance over a stationary elementary volume ∆x∆y∆z through which the fluid is flowing z position (x+∆x, y+∆y, z+∆z) ρvx |x y ∆y position (x, y, z) ∆z ρvx |x+∆x ∆x x • Rate of Mass accumulation = (Rate of Mass In) - (Rate of Mass Out) Rate of Mass In through face at x is ρvx |x ∆y ∆z Rate of Mass Out through face at x+∆x is ρvx |x+∆x ∆y ∆z Same Thing as Above for other faces Rate of Mass Accumulation is ∆x ∆y ∆z t ∆x∆y∆z = ∆y ∆z (ρvx |x - ρvx |x+∆x ) + ∆x ∆z (ρvy |y - ρvy |y+∆y) t + ∆y ∆x (ρv | - ρv | z z z z+∆z ) • Dividing on both sides by ∆x∆y∆z and taking the limit as ∆x∆y∆z goes to zero yields: = − t ( vx ) + x ( v y ) + y ( vz ) = − ∇ • ( v) z ( ) which is known as the Continuity Equation (mass balance for fixed observer). The above equation can be rewritten as: v x vy vz + v x + vy + vz = − + + t y z x y z x which is equivalent for an observer moving along the flow to: D Note that for a fluid of constant density =− ∇• v Dρ/Dt = 0 and Div(v) = 0 (incompressibility) Dt ( ) The Equation of Motion • We will again consider a small volume element ∆x∆y∆z at position x, y, z and write the momentum balnce on the fluid element (similar approach to mass balance) • Rate of Momentum Accumulation = (Rate of Momentum In) (Rate of Momentum Out) + (Sum of Forces Acting on System) (Note that this is a generalization of what we did in the case of the Shell Momentum Balance for unsteady state conditions) for transport of xτyx |y +∆y momentum through τxz |z +∆z z each of the 6 faces (the transport of yτxx |x τxx |x +∆x and z-momentum ∆z ∆y through each of the y ∆x τyx |y faces is handled similarly) x τzx |z • Note that momentum flows into and out of the volume element by two different mechanisms: 1- Convection (bulk fluid flow) and 2Molecular Transfer (velocity gradient in Newton’s or related laws) • Molecular Transfer: Remembering that τyx is the flux of x-momentum through a face perpendicular to the y-axis, the rate at which x-momentum enters the face at y (perpendicular to the y-axis) is ∆x∆z τyx |y and the rate at which x-momentum leaves the face at y+∆y is ∆x∆z τyx |y+∆y . Similarly, the rate at which x-momentum enters the face at x by molecular transfer is ∆y∆z τxx |x and the rate at which x-momentum leaves the face at x+∆x is ∆y∆z τxx |x+∆x . Similarly, there is another term for the x-momentum entering by the face at z and leaving by face at z+ ∆z ( ) ∆y∆z( xx|x − xx|x+ ∆x ) + ∆x∆z yx|y − yx|y+∆y + ∆y∆x( zx|z − zx|z+∆z ) The same should be done with y-momentum and with z-momentum. • Note that τxx is the normal stress on the x-face Note that τyx is the x-directed tangential (shear) stress on the yface resulting from viscous forces Note that τzx is the x-directed tangential (shear) stress on the zface resulting from viscous forces • Convection: Rate at which x-momentum enters face at x by convection Rate at which x-momentum leaves face at x+∆x by convection Rate at which x-momentum enters face at y by convection Rate at which x-momentum leaves face at y+∆y by convection Rate at which x-momentum enters face at z by convection Rate at whic x-momentum leaves face at z+∆z by convection ∆y ∆z (ρvxvx|x) ∆y ∆z (ρvxvx|x+∆x) ∆x ∆z (ρvyvx|y) ∆x ∆z (ρvyvx|y+∆y) ∆y ∆x (ρvzvx|z) ∆y ∆x (ρvzvx|z+∆z) The same can be done for yand z- momenta • The other important contributions arise from forces acting on the fluid (fluid pressure and gravity) in the x-direction these forces contribute: ∆y∆z( p x|x − p x|x+∆x ) + ∆x∆y∆z ( g x ) p is a scalar quantity, which is a function of ρ and T g is a vectorial quantity with components (gx, gy, gz) • Summing all contributions and dividing by ∆x∆y∆z and taking the limit of ∆x∆y∆z to zero yields ( v x ) = − ( v x v x ) + v y v x + ( v z v x ) t x y z ( − xx + x yx y + ) zx − p + g x z x Similarly for the contributions in the yu- and z-directions: • Note that ρvx, ρvy, ρvz are the components of the vector ρv Note that gx, gy, gz are the components of g p p p are the components of grad(p) , , x y z Note that Note that ρvxvx , ρvyvx , ρvzvx , ρvxvy , .... are the nine components of the convective momentum flux ( a dyadic product of ρv and v (not the dot product). Note that τxx , τxy , τxz , τyx , τyy , ....are the nine components of the stress tensor τ. Gravitational ( v) = − ∇ • ( v v ) − ∇p − ∇ • + g force per unit t volume [ ] [ ] Rate of increase Convection Pressure force of momentum contribution per unit volume per unit volume Viscous Transfer Contribution • The previous equation is Newton’s second law (equation of motion) expressed for a stationary volume element • This equation can be rewritten for a small volume element of fluid moving along the flow. [ ]+ Dv = −∇p − ∇ • Dt g • To determine the velocity distribution, one now needs to insert expressions for the various stresses in terms of velocity gradients and fluid properties. For Newtonian Fluids one has: xx = −2 yy = −2 zz = −2 vx 2 + x 3 vy 2 + y 3 (∇ • v ) vy vx + yx = xy = − y x (∇ • v ) vy vz + yz = zy = − y z vz 2 + z 3 ( vx vz = = − + zx xz x z ∇•v ) • Combining Newton’s second law with the equations describing the relationships between stresses and viscosity allows to derive the velocity profile for the flow system. One may also need the fluid equation of state (P = f(ρ,T)) and the density dependence of viscosity (µ = f(ρ))along Boundary and Initial conditions. • To solve a flow problem, write the Continuity equation and the Equation of Motion in the appropriate coordinate system and for the appropriate symmetry (cartesian, cylindrical, spherical), then discard all terms that are zero. Use your intuition, while keeping track of the terms you are ignoring (check your assumptions at the end). Use the Newtonian or Non-Newtonian relationship between velocity gradient and shear stresses. Integrate differential equation using appropriate boundary conditions General Equation of Motion in Cartesian Coordinate System Equation of Motion In Cartesian Coordinate System For Newtonian Incompressible Fluids General Equation of Motion in the Cylindrical Coordinate System Equation of Motion In the Cylindrical Coordinate System For Newtonian Incompressible Fluids General Equation of Motion in the Spherical Coordinate System Equation of Motion In the Cylindrical Coordinate System For Newtonian Incompressible Fluids Stress Tensor for Newtonian Fluids in the Cylindrical Coordinate System Stress Tensor For Newtonian Fluids In Spherical Coordinates The function : ∇v = Φ v for Newtonian Fluids Axial Flow of an Incompressible Fluid in a Circular Tube of Length L and Radius R • • • • Use Cylindrical Coordinates Set vθ = vr = 0 (flow along the z-axis) vz is not a function of θ because of cylindrical symmetry Worry only about the z-component of the equation of motion 2 1 v vz vz p z vz =− + gz + r + 2 z z r z r r 2v vz z Continuity equation: =0⇒ =0 2 z z vz p 0 =− + gz + r Integrate twice w/respect to r using z r r r v =0 at r=R and v finite at r=0 z z [ P − P ] R2 r2 0 L 1 − v z = 2 4L R