Topics in Shear Flow Chapter 1 - Introduction
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Chapter 4: Immersed Body Flow [Pp
MECH 3492 Fluid Mechanics and Applications Univ. of Manitoba Fall Term, 2017 Chapter 4: Immersed Body Flow [pp. 445-459 (8e), or 374-386 (9e)] Dr. Bing-Chen Wang Dept. of Mechanical Engineering Univ. of Manitoba, Winnipeg, MB, R3T 5V6 When a viscous fluid flow passes a solid body (fully-immersed in the fluid), the body experiences a net force, F, which can be decomposed into two components: a drag force F , which is parallel to the flow direction, and • D a lift force F , which is perpendicular to the flow direction. • L The drag coefficient CD and lift coefficient CL are defined as follows: FD FL CD = 1 2 and CL = 1 2 , (112) 2 ρU A 2 ρU Ap respectively. Here, U is the free-stream velocity, A is the “wetted area” (total surface area in contact with fluid), and Ap is the “planform area” (maximum projected area of an object such as a wing). In the remainder of this section, we focus our attention on the drag forces. As discussed previously, there are two types of drag forces acting on a solid body immersed in a viscous flow: friction drag (also called “viscous drag”), due to the wall friction shear stress exerted on the • surface of a solid body; pressure drag (also called “form drag”), due to the difference in the pressure exerted on the front • and rear surfaces of a solid body. The friction drag and pressure drag on a finite immersed body are defined as FD,vis = τwdA and FD, pres = pdA , (113) ZA ZA Streamwise component respectively. -
Derivation of Fluid Flow Equations
TPG4150 Reservoir Recovery Techniques 2017 1 Fluid Flow Equations DERIVATION OF FLUID FLOW EQUATIONS Review of basic steps Generally speaking, flow equations for flow in porous materials are based on a set of mass, momentum and energy conservation equations, and constitutive equations for the fluids and the porous material involved. For simplicity, we will in the following assume isothermal conditions, so that we not have to involve an energy conservation equation. However, in cases of changing reservoir temperature, such as in the case of cold water injection into a warmer reservoir, this may be of importance. Below, equations are initially described for single phase flow in linear, one- dimensional, horizontal systems, but are later on extended to multi-phase flow in two and three dimensions, and to other coordinate systems. Conservation of mass Consider the following one dimensional rod of porous material: Mass conservation may be formulated across a control element of the slab, with one fluid of density ρ is flowing through it at a velocity u: u ρ Δx The mass balance for the control element is then written as: ⎧Mass into the⎫ ⎧Mass out of the ⎫ ⎧ Rate of change of mass⎫ ⎨ ⎬ − ⎨ ⎬ = ⎨ ⎬ , ⎩element at x ⎭ ⎩element at x + Δx⎭ ⎩ inside the element ⎭ or ∂ {uρA} − {uρA} = {φAΔxρ}. x x+ Δx ∂t Dividing by Δx, and taking the limit as Δx approaches zero, we get the conservation of mass, or continuity equation: ∂ ∂ − (Aρu) = (Aφρ). ∂x ∂t For constant cross sectional area, the continuity equation simplifies to: ∂ ∂ − (ρu) = (φρ) . ∂x ∂t Next, we need to replace the velocity term by an equation relating it to pressure gradient and fluid and rock properties, and the density and porosity terms by appropriate pressure dependent functions. -
Literature Review
SECTION VII GLOSSARY OF TERMS SECTION VII: TABLE OF CONTENTS 7. GLOSSARY OF TERMS..................................................................... VII - 1 7.1 CFD Glossary .......................................................................................... VII - 1 7.2. Particle Tracking Glossary.................................................................... VII - 3 Section VII – Glossary of Terms Page VII - 1 7. GLOSSARY OF TERMS 7.1 CFD Glossary Advection: The process by which a quantity of fluid is transferred from one point to another due to the movement of the fluid. Boundary condition(s): either: A set of conditions that define the physical problem. or: A plane at which a known solution is applied to the governing equations. Boundary layer: A very narrow region next to a solid object in a moving fluid, and containing high gradients in velocity. CFD: Computational Fluid Dynamics. The study of the behavior of fluids using computers to solve the equations that govern fluid flow. Clustering: Increasing the number of grid points in a region to better resolve a geometric or flow feature. Increasing the local grid resolution. Continuum: Having properties that vary continuously with position. The air in a room can be thought of as a continuum because any cube of air will behave much like any other chosen cube of air. Convection: A similar term to Advection but is a more generic description of the Advection process. Convergence: Convergence is achieved when the imbalances in the governing equations fall below an acceptably low level during the solution process. Diffusion: The process by which a quantity spreads from one point to another due to the existence of a gradient in that variable. Diffusion, molecular: The spreading of a quantity due to molecular interactions within the fluid. -
Chapter 4: Immersed Body Flow [Pp
MECH 3492 Fluid Mechanics and Applications Univ. of Manitoba Fall Term, 2017 Chapter 4: Immersed Body Flow [pp. 445-459 (8e), or 374-386 (9e)] Dr. Bing-Chen Wang Dept. of Mechanical Engineering Univ. of Manitoba, Winnipeg, MB, R3T 5V6 When a viscous fluid flow passes a solid body (fully-immersed in the fluid), the body experiences a net force, F, which can be decomposed into two components: a drag force F , which is parallel to the flow direction, and • D a lift force F , which is perpendicular to the flow direction. • L The drag coefficient CD and lift coefficient CL are defined as follows: FD FL CD = 1 2 and CL = 1 2 , (112) 2 ρU A 2 ρU Ap respectively. Here, U is the free-stream velocity, A is the “wetted area” (total surface area in contact with fluid), and Ap is the “planform area” (maximum projected area of an object such as a wing). In the remainder of this section, we focus our attention on the drag forces. As discussed previously, there are two types of drag forces acting on a solid body immersed in a viscous flow: friction drag (also called “viscous drag”), due to the wall friction shear stress exerted on the • surface of a solid body; pressure drag (also called “form drag”), due to the difference in the pressure exerted on the front • and rear surfaces of a solid body. The friction drag and pressure drag on a finite immersed body are defined as FD,vis = τwdA and FD, pres = pdA , (113) ZA ZA Streamwise component respectively. -
Part III: the Viscous Flow
Why does an airfoil drag: the viscous problem – André Deperrois – March 2019 Rev. 1.1 © Navier-Stokes equations Inviscid fluid Time averaged turbulence CFD « RANS » Reynolds Averaged Euler’s equations Reynolds equations Navier-stokes solvers irrotational flow Viscosity models, uniform 3d Boundary Layer eq. pressure in BL thickness, Prandlt Potential flow mixing length hypothesis. 2d BL equations Time independent, incompressible flow Laplace’s equation 1d BL Integral 2d BL differential equations equations 2d mixed empirical + theoretical 2d, 3d turbulence and transition models 2d viscous results interpolation The inviscid flow around an airfoil Favourable pressure gradient, the flo! a""elerates ro# $ero at the leading edge%s stagnation point& Adverse pressure gradient, the low decelerates way from the surface, the flow free tends asymptotically towards the stream air freestream uniform flow flow inviscid ◀—▶ “laminar”, The boundary layer way from the surface, the fluid’s velocity tends !ue to viscosity, the asymptotically towards the tangential velocity at the velocity field of an ideal inviscid contact of the foil is " free flow around an airfoil$ stream air flow (magnified scale) The boundary layer is defined as the flow between the foil’s surface and the thic%ness where the fluid#s velocity reaches &&' or &&$(' of the inviscid flow’s velocity. The viscous flow around an airfoil at low Reynolds number Favourable pressure gradient, the low a""elerates ro# $ero at the leading edge%s stagnation point& Adverse pressure gradient, the low decelerates +n adverse pressure gradients, the laminar separation bubble forms. The flow goes flow separates. The velocity close to the progressively turbulent inside the bubble$ surface goes negative. -
Incompressible Irrotational Flow
Incompressible irrotational flow Enrique Ortega [email protected] Rotation of a fluid element As seen in M1_2, the arbitrary motion of a fluid element can be decomposed into • A translation or displacement due to the velocity. • A deformation (due to extensional and shear strains) mainly related to viscous and compressibility effects. • A rotation (of solid body type) measured through the midpoint of the diagonal of the fluid element. The rate of rotation is defined as the angular velocity. The latter is related to the vorticity of the flow through: d 2 V (1) dt Important: for an incompressible, inviscid flow, the momentum equations show that the vorticity for each fluid element remains constant (see pp. 17 of M1_3). Note that w is positive in the 2 – Irrotational flow counterclockwise sense Irrotational and rotational flow ij 0 According to the Prandtl’s boundary layer concept, thedomaininatypical(high-Re) aerodynamic problem at low can be divided into outer and inner flow regions under the following considerations: Extracted from [1]. • In the outer region (away from the body) the flow is considered inviscid and irrotational (viscous contributions vanish in the momentum equations and =0 due to farfield vorticity conservation). • In the inner region the viscous effects are confined to a very thin layer close to the body (vorticity is created at the boundary layer by viscous stresses) and a thin wake extending downstream (vorticity must be convected with the flow). Under these hypotheses, it is assumed that the disturbance of the outer flow, caused by the body and the thin boundary layer around it, is about the same caused by the body alone. -
Continuity Equation in Pressure Coordinates
Continuity Equation in Pressure Coordinates Here we will derive the continuity equation from the principle that mass is conserved for a parcel following the fluid motion (i.e., there is no flow across the boundaries of the parcel). This implies that δxδyδp δM = ρ δV = ρ δxδyδz = − g is conserved following the fluid motion: 1 d(δM ) = 0 δM dt 1 d()δM = 0 δM dt g d ⎛ δxδyδp ⎞ ⎜ ⎟ = 0 δxδyδp dt ⎝ g ⎠ 1 ⎛ d(δp) d(δy) d(δx)⎞ ⎜δxδy +δxδp +δyδp ⎟ = 0 δxδyδp ⎝ dt dt dt ⎠ 1 ⎛ dp ⎞ 1 ⎛ dy ⎞ 1 ⎛ dx ⎞ δ ⎜ ⎟ + δ ⎜ ⎟ + δ ⎜ ⎟ = 0 δp ⎝ dt ⎠ δy ⎝ dt ⎠ δx ⎝ dt ⎠ Taking the limit as δx, δy, δp → 0, ∂u ∂v ∂ω Continuity equation + + = 0 in pressure ∂x ∂y ∂p coordinates 1 Determining Vertical Velocities • Typical large-scale vertical motions in the atmosphere are of the order of 0. 01-01m/s0.1 m/s. • Such motions are very difficult, if not impossible, to measure directly. Typical observational errors for wind measurements are ~1 m/s. • Quantitative estimates of vertical velocity must be inferred from quantities that can be directly measured with sufficient accuracy. Vertical Velocity in P-Coordinates The equivalent of the vertical velocity in p-coordinates is: dp ∂p r ∂p ω = = +V ⋅∇p + w dt ∂t ∂z Based on a scaling of the three terms on the r.h.s., the last term is at least an order of magnitude larger than the other two. Making the hydrostatic approximation yields ∂p ω ≈ w = −ρgw ∂z Typical large-scale values: for w, 0.01 m/s = 1 cm/s for ω, 0.1 Pa/s = 1 μbar/s 2 The Kinematic Method By integrating the continuity equation in (x,y,p) coordinates, ω can be obtained from the mean divergence in a layer: ⎛ ∂u ∂v ⎞ ∂ω ⎜ + ⎟ + = 0 continuity equation in (x,y,p) coordinates ⎝ ∂x ∂y ⎠ p ∂p p2 p2 ⎛ ∂u ∂v ⎞ ∂ω = − ⎜ + ⎟ ∂p rearrange and integrate over the layer ∫p ∫ ⎜ ⎟ 1 ∂x ∂y p1⎝ ⎠ p ⎛ ∂u ∂v ⎞ ω(p )−ω(p ) = (p − p )⎜ + ⎟ overbar denotes pressure- 2 1 1 2 ⎜ ⎟ weighted vertical average ⎝ ∂x ∂y ⎠ p To determine vertical motion at a pressure level p2, assume that p1 = surface pressure and there is no vertical motion at the surface. -
Chapter 3 Newtonian Fluids
CM4650 Chapter 3 Newtonian Fluid 2/5/2018 Mechanics Chapter 3: Newtonian Fluids CM4650 Polymer Rheology Michigan Tech Navier-Stokes Equation v vv p 2 v g t 1 © Faith A. Morrison, Michigan Tech U. Chapter 3: Newtonian Fluid Mechanics TWO GOALS •Derive governing equations (mass and momentum balances •Solve governing equations for velocity and stress fields QUICK START V W x First, before we get deep into 2 v (x ) H derivation, let’s do a Navier-Stokes 1 2 x1 problem to get you started in the x3 mechanics of this type of problem solving. 2 © Faith A. Morrison, Michigan Tech U. 1 CM4650 Chapter 3 Newtonian Fluid 2/5/2018 Mechanics EXAMPLE: Drag flow between infinite parallel plates •Newtonian •steady state •incompressible fluid •very wide, long V •uniform pressure W x2 v1(x2) H x1 x3 3 EXAMPLE: Poiseuille flow between infinite parallel plates •Newtonian •steady state •Incompressible fluid •infinitely wide, long W x2 2H x1 x3 v (x ) x1=0 1 2 x1=L p=Po p=PL 4 2 CM4650 Chapter 3 Newtonian Fluid 2/5/2018 Mechanics Engineering Quantities of In more complex flows, we can use Interest general expressions that work in all cases. (any flow) volumetric ⋅ flow rate ∬ ⋅ | average 〈 〉 velocity ∬ Using the general formulas will Here, is the outwardly pointing unit normal help prevent errors. of ; it points in the direction “through” 5 © Faith A. Morrison, Michigan Tech U. The stress tensor was Total stress tensor, Π: invented to make the calculation of fluid stress easier. Π ≡ b (any flow, small surface) dS nˆ Force on the S ⋅ Π surface V (using the stress convention of Understanding Rheology) Here, is the outwardly pointing unit normal of ; it points in the direction “through” 6 © Faith A. -
Geosc 548 Notes R. Dibiase 9/2/2016
Geosc 548 Notes R. DiBiase 9/2/2016 1. Fluid properties and concept of continuum • Fluid: a substance that deforms continuously under an applied shear stress (e.g., air, water, upper mantle...) • Not practical/possible to treat fluid mechanics at the molecular level! • Instead, need to define a representative elementary volume (REV) to average quantities like velocity, density, temperature, etc. within a continuum • Continuum: smoothly varying and continuously distributed body of matter – no holes or discontinuities 1.1 What sets the scale of analysis? • Too small: bad averaging • Too big: smooth over relevant scales of variability… An obvious length scale L for ideal gases is the mean free path (average distance traveled by before hitting another molecule): = ( 1 ) 2 2 where kb is the Boltzman constant, πr is the effective4√2 cross sectional area of a molecule, T is temperature, and P is pressure. Geosc 548 Notes R. DiBiase 9/2/2016 Mean free path of atmosphere Sea level L ~ 0.1 μm z = 50 km L ~ 0.1 mm z = 150 km L ~ 1 m For liquids, not as straightforward to estimate L, but typically much smaller than for gas. 1.2 Consequences of continuum approach Consider a fluid particle in a flow with a gradient in the velocity field : �⃑ For real fluids, some “slow” molecules get caught in faster flow, and some “fast” molecules get caught in slower flow. This cannot be reconciled in continuum approach, so must be modeled. This is the origin of fluid shear stress and molecular viscosity. For gases, we can estimate viscosity from first principles using ideal gas law, calculating rate of momentum exchange directly. -
Flow Over Immerced Bodies
Department of Mechanical Engineering AMEE401 / AUTO400 Instructor: Marios M. Fyrillas Aerodynamics Email: [email protected] HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible for the drag and the lift, the pressure and the shear stress: The drag force and the lift force on an object can be obtained by: Dpcos d Aw sin d A Lpsin d A w cos d A LIFT: Most common lift-generating devices (i.e., airfoils, fans, spoilers on cars, etc.) operate in the large Reynolds number range in which the flow has a boundary layer character, with viscous effects confined to the boundary layers and wake regions. For such cases the wall shear stress, w , contributes little to the lift. Most of the lift comes from the surface pressure distribution, as justified through Bernoulli’s equation. For objects operating in very low Reynolds number regimes (i.e. Re 1) viscous effects are important, and the contribution of the shear stress to the lift may be as important as that of the pressure. Such situations include the flight of minute insects and the swimming of microscopic organisms. DRAG: Similarly, when flow separation occurs, i.e. a blunt body or a streamlined body at a large angle of attack, the major component of the drag force is pressure differential due to the low pressure in the flow separation region. On the ultimately streamlined body (a zero thickness flat plate parallel to the flow) the drag is entirely due to the shear stress at the surface (boundary layers) and, is relatively small. -
CONTINUITY EQUATION Another Principle on Which We Can Derive a New Equation Is the Conservation of Mass
ESCI 342 – Atmospheric Dynamics I Lesson 7 – The Continuity and Additional Equations Suggested Reading: Martin, Chapter 3 THE SYSTEM OF EQUATIONS IS INCOMPLETE The momentum equations in component form comprise a system of three equations with 4 unknown quantities (u, v, p, and ). Du1 p fv (1) Dt x Dv1 p fu (2) Dt y p g (3) z They are not a closed set, because there are four dependent variables (u, v, p, and ), but only three equations. We need to come up with some more equations in order to close the set. DERIVATION OF THE CONTINUITY EQUATION Another principle on which we can derive a new equation is the conservation of mass. The equation derived from this principle is called the mass continuity equation, or simply the continuity equation. Imagine a cube at a fixed point in space. The net change in mass contained within the cube is found by adding up the mass fluxes entering and leaving through each face of the cube.1 The mass flux across a face of the cube normal to the x-axis is given by u. Referring to the picture below, these fluxes will lead to a rate of change in mass within the cube given by m u yz u yz (4) t x xx 1 A flux is a quantity per unit area per unit time. Mass flux is therefore the rate at which mass moves across a unit area, and would have units of kg s1 m2. The mass in the cube can be written in terms of the density as m xyz so that m x y z . -
1 the Continuity Equation 2 the Heat Equation
1 The Continuity Equation Imagine a fluid flowing in a region R of the plane in a time dependent fashion. At each point (x; y) R2 it has a velocity v = v (x; y; t) at time t. Let ρ = ρ(x; y; t) be the density of the fluid 2 −! −! at (x; y) at time t. Let P be any point in the interior of R and let Dr be the closed disk of radius r > 0 and center P . The mass of fluid inside Dr at any time t is ρ dxdy: ZZDr If matter is neither created nor destroyed inside Dr, the rate of decrease of this quantity is equal to the flux of the vector field −!J = ρ−!v across Cr, the positively oriented boundary of Dr. We therefore have d ρ dxdy = −!J −!N ds; dt − · ZZDr ZCr where N~ is the outer normal and ds is the element of arc length. Notice that the minus sign is needed since positive flux at time t represents loss of total mass at that time. Also observe that the amount of fluid transported across a small piece ds of the boundary of Dr at time t is ρ−!v −!N ds. Differentiating under the integral sign on the left-hand side and using the flux form of Green’s· Theorem on the right-hand side, we get @ρ dxdy = −!J dxdy: @t − r · ZZDr ZZDr Gathering terms on the left-hand side, we get @ρ ( + −!J ) dxdy = 0: @t r · ZZDr If the integrand was not zero at P it would be different from zero on Dr for some sufficiently small r and hence the integral would not be zero which is not the case.