Conservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
|
|
- Elmer Francis
- 5 years ago
- Views:
Transcription
1 Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation of the equations governing fluid flow in integral form Conservation of Momentum Conservation of Energy Summary Incompressible flows Inviscid compressible flows Conservation of Mass In general, mass can be added or removed: Rate of Net influx increase of = of mass mass time t The conservation law must be stated as: Final Mass = Original Mass + time t+ t Mass Mass - Added Removed We now apply this statement to an arbitrary control volume in an arbitrary flow field
2 volume volume The rate of change of total mass in the control volume is given by: Rate of change ρdv Total mass surface S To find the mass flux through the control surface, let s examine a small part of the surface where the velocity is normal to the surface surface S Density ρ elocity U x Select a small rectangle outside the boundary such that during time t it flows into the C: A U t time t time t+ t Mass flow through the boundary A during t is: ρu x ΔtΔA where U x = -u n The sign of the normal component of the velocity determines whether the fluid flows in or out of the control volume. We will take the outflow to be positive: The net in-flow through the boundary of the control volume is therefore: Negative since this is net inflow Integral over the boundary u n < 0 S ρu n ds Inflow Mass flow normal to boundary Equation in Integral Form ρdv = ρu nds S Rate of change of mass Net inflow of mass Conservation of momentum
3 Conservation of momentum Rate of change of momentum Momentum per volume: ρu Momentum in control volume: ρu dv Rate of increase of = Net influx of Body momentum + momentum forces + Surface forces The rate of change of momentum in the control volume is given by: Rate of change ρu dv Total momentum volume surface S Density ρ elocity u In/out flow of momentum Select a small rectangle outside the boundary such that during time t it flows into the C: A U x U x t time t time t+ t U x = -u n so the net influx of momentum per unit time (divided by t ) is: Momentum flow through the boundary A during t is: ( ρu)u x ΔtΔA ρuu nds Body forces, such as gravity act on the fluid in the control volume. Body force dv f ρ f dv volume Total body force iscous force iscous flux of momentum The viscous force is given by the dot product of the normal with the stress tensor T ds T n T nds S volume surface S Stress Tensor T = ( p + λ u)i + 2µD Deformation Tensor ( ) D = 1 u + ut 2 Whose components are: D ij = 1 u i + u j 2 x j x i For incompressible flows u = 0 And the stress tensor is T = pi + 2µD Total stress force
4 The pressure produces a normal force on the control surface. The total force is: Pressure force ds pn pnds S volume surface S Total pressure force Conservation of Momentum Gathering the terms ρudv = ρuu nds Tnds + ρ f dv Substitute for the stress tensor ρudv = ρuu nds pnds + 2µD nds + ρ f dv Rate of change of momentum Net inflow of momentum Total pressure T = pi + 2µD Total viscous force Total body force Conservation of energy Conservation of energy The energy equation in integral form Rate of change of kinetic+internal energy ρ e u2 dv = ρ e u2 u nds + u f dv + Work done by body forces Net inflow of kinetic+internal energy n (ut)ds n qds Net work done by the stress tensor Net heat flow Differential Form of the Governing Equations The Divergence or Gauss Theorem can be used to convert surface integrals to volume integrals a dv = a n ds S
5 Start with the integral form of the mass conservation equation ρdv = ρu nds S Using Gauss s theorem ρu nds = (ρu)dv S The mass conservation equations becomes ρdv = (ρu)dv Rearranging ρdv + (ρu)dv = 0 Since the control volume is fixed, the derivative can be moved under the integral sign ρ dv + (ρu)dv = 0 Rearranging ρ + (ρu) dv = 0 ρ + (ρu) dv = 0 This equation must hold for ANY control volume, no matte what shape and size. Therefore the integrand must be equal to zero ρ + (ρu) = 0 Expanding the divergence term: (ρu) = u ρ + ρ u The mass conservation equation becomes: or where ρ ρ + (ρu) = + u ρ + ρ u Dρ Dt + ρ u = 0 D() Dt = () + u () Convective derivative The differential form of the momentum equation is derived in the same way: Start with ρu dv = ρ f dv + (nt ρu(u n))ds Rewrite as: ρu dv = ρ f dv + ( T ρuu )dv To get: ρu = ρ f + (T ρuu) Using the mass conservation equation the advection part can be rewritten: ρu + ρuu = u ρ u + u uρ + ρ + ρu u = u ρ + uρ + ρ Du Dt = ρ Du Dt =0, by mass conservation
6 The momentum equation equation ρu = ρ f + (T ρuu) Can therefore be rewritten as ρ Du Dt = ρ f + T The energy equation equation can be converted to a differential form in the same way. It is usually simplified by subtracting the mechanical energy where ρ Du Dt = ρ u + u u The mechanical energy equation is obtained by taking the dot product of the momentum equation and the velocity: u ρ u + ρu u = ρ f + T The result is: ρ u 2 = ρu u 2 + ρu f 2 2 b + u T The mechanical energy equation ( ) The final result is Subtract the mechanical energy equation from the energy equation ρ De T u + q = 0 Dt Using the constitutive relation presented earlier for the stress tensor and Fourie s law for the heat conduction: T = ( p + λ u)i + 2µD q = k T The final result is where ρ De Dt + p u = Φ + k T Φ = λ( u) 2 + 2µD D is the dissipation function and is the rate at which work is converted into heat. Generally we also need: p = p(e, ρ); T = T(e, ρ); and equations for µ, λ, k Summary of governing equations
7 Integral Form ρ dv = ρu nds S ρudv = ρ f dv + (nt ρu(u n))ds ρ(e u2 )dv = u ρ f dv + n (ut ρ(e u2 ) q)ds Conservative Form ρ + (ρu) = 0 ρu = ρ f + (T ρuu) ρ(e u2 ) = ρ(e u2 )u ut +q Convective Form Dρ Dt + ρ u = 0 ρ Du Dt = ρ f + T ρ De = T u q Dt Special Cases Compressible inviscid flows Incompressible flows Stokes flow Not in this course! Potential flows Inviscid compressible flows U = U + E x + F y + G z = 0 ρ ρu ρv ρw ρe E = ρu ρu 2 + p ρuv ρuw ( ρe + p)u F = Where we have taken ρv ρuv ρv 2 + p ρvw ( ρe + p)v T = pi ρu ρuw G = ρvw ρw 2 + p ( ρe + p)u Incompressible flow p = ρrt; e = c v T; h = c p T; γ = c p /c v c v = R (γ 1)e ; p = (γ 1)ρe; T = ; γ 1 R Assuming Ideal Gas
8 Incompressible flows: Dρ Dt = 0 Dρ Dt + ρ u = 0 u = 0 Navier-Stokes equations (conservation of momentum) ρu = ρuu- P + ρ f + µ u + b T u For constant viscosity ρ u + ρu u = - P + ρ f b + µ 2 u ( ) Objectives: The 2D Navier-Stokes Equations for incompressible, homogeneous flow: u + u u x + v u y = 1 P ρ x + u x + v y = 0 µ 2 u ρ x + 2 u 2 y 2 v + u v x + v v y = 1 P ρ y + µ 2 v ρ x + 2 v 2 y 2 u nds = 0 S ud vd Incompressibility (conservation of mass) Navier-Stokes equations (conservation of momentum) = uu nds +ν u nds 1 ρ pn x ds S S = vu nds +ν v n S ds 1 S ρ pn S y ds S Advection: The role of pressure: Pressure and Advection Newton s law of motion: In the absence of forces, a fluid particle will move in a straight line Needed to accelerate/decelerate a fluid particle: Easy to use if viscous forces are small Needed to prevent accumulation/depletion of fluid particles: Use if there are strong viscous forces Pressure Pressure The pressure opposes local accumulation of fluid. For compressible flow, the pressure increases if the density increases. For incompressible flows, the pressure takes on whatever value necessary to prevent local accumulation: Increasing pressure slows the fluid down and decreasing pressure accelerates it Pressure increases Pressure decreases High Pressure Low Pressure u < 0 u > 0 Slowing down Speeding up
9 iscosity Diffusion of fluid momentum is the result of friction between fluid particles moving at uneven speed. The velocity of fluid particles initially moving with different velocities will gradually become the same. Due to friction, more and more of the fluid next to a solid wall will move with the wall velocity. Instability and unsteadyness For all but the lowest Reynolds numbers, the fluid flow is unstable and unsteady, forming transient whorls and vortices that greatly increase the transfer of momentum over what viscosity alone can do. The Kelvin Helmholtz instability of a slip line between flows in different directions is one of the fundamental ways in which steady flow becomes unstable. Time Nondimensional Numbers the Reynolds number U 2 Nondimensional Numbers u + u u x + v u y = 1 P ρ x + u, v = u,v U ; x, y = x, y L ; µ 2 u + 2 u ρ x 2 y 2 t = Ut L u L t + u u x + v u y = 1 P Lρ x + U µ 2 u L 2 ρ x + 2 u 2 y 2 u t + u u x + v u y = L 1 U 2 P Lρ x + L U µ 2 u U 2 L 2 ρ x + 2 u 2 y 2 u t + Nondimensional Numbers u u x + v u y = L 1 U 2 P Lρ x + L 1 P U 2 Lρ x = 1 P ρu 2 x = P x L U µ 2 u U 2 L 2 ρ x + 2 u 2 y 2 where P = P ρu 2 Nondimensional Numbers The momentum equations are: u t + u u x + v u y = P x u Re x + 2 u 2 y 2 v t + u v x + v v y = P y v Re x + 2 v 2 y 2 L U µ U 2 L 2 ρ = µ ULρ = 1 Re where Re = ρul µ The continuity equation is unchanged u x + v y = 0
10 Outline Derivation of the equations governing fluid flow in integral form Conservation of Momentum Conservation of Energy Summary Incompressible flows Inviscid compressible flows Integral versus differential form of the governing equations When using FINITE OLUME approximations, we work directly with the integral form of the conservation principles. The average values of f in a small volume are stored f 1 +1 x j-1 x j x j+1 x j-1/2 x j+1/2 x To derive an equation that governs the evolution of the average value in each cell, consider: d dt f + F x = 0 Integrating the equation over a small control volume of size h= Δx=x j+1 -x j yields: f dx + F dx = 0 x Δx Δx f dx + F j +1/ 2 F j 1/ 2 = 0 Δx The volume is fixed in space or: d dt (h ) = F j 1/2 F j +1/2 Notation: Average value of f in cell j where the average in each cell is defined by: = 1 h Δx f dx F j 1/ 2 Flow of f into cell j x x j j-1/2 x j+1/2 F j +1/ 2 Flow of f out of cell j F j-1/2 is the flow of f into cell j F j+1/2 is the flow of f out off cell j F is usually called the flux function Left boundary of cell j Right boundary of cell j Center of cell j
11 The exact form of F depends on the problem: F = Uf F = D f x F = Uf D f x Advection Diffusion Advection/Diffusion The statement d dt (h) = F j 1/2 F j +1/2 Is exactly the original conservation law The rate of increase of the average value of f is equal to the inflow into the control volume minus the outflow from the control volume Rate of increase in f = inflow of f -outflow of f Example: advection/diffusion equation Rewrite: were f + U f 2 x = D f x 2 f + F x = 0 F = Uf D f x Assume here that U is a constant For each cell Average value of f in cell j x x j j-1/2 d dt (h) = F j 1/2 F j +1/2 x j+1/2 F j +1/ 2 = 1 h Δx f dx F j +1/2 = U +1/2 D f x j +1/2 x x j j-1/2 x j+1/2 Average value of f in cell j Approximate: F j +1/ 2 +1 / ( +1 + ) f at a cell boundary Need to find F j +1/2 = U +1/2 D f x j +1/2 and f +1 x j +1 / 2 h gradient of f at a cell boundary Approximate: F j +1/2 = U +1/2 D f 1 x j +1/2 2 U( f + f ) D f +1 j j +1 j h F j 1/2 =U 1/2 D f 1 x j 1/2 2 U( f + f ) D 1 j j 1 h and d dt 1 Δt ( f n+1 j f n j ) time derivative of f
12 Substituting: h Δt ( f n+1 j f n j ) = U( 1 ( f + f ) 2 j +1 j 1 ( f 2 j + 1 )) + D f +1 j h Or, rearranging the terms: n +1 n 1 h + U 2h ( f f ) = D j +1 j 1 ( f h 2 j ) Δt Which is exactly the same as the finite difference equation if we take the average value to be the same as the value in the center of the cell Example: U=1; D=0.05; k=1 N=21 Δt=0.05 Exact Numerical The finite volume point of view has two main advantages over the finite difference point of view 1. Working directly with the conservation principle (no implicit assumption of smoothness and differentiability) Most physical laws are based on CONSERATION principles: In the absence of explicit sources or sinks, f is neither created nor destroyed. For a control volume fixed in space, we can state the conservation of f as: 2. Easier visualization of how the solution is updated (the fluxes have a physical meaning) Analysis of accuracy and stability is, however, usually easier using the finite difference point of view Amount of f in a control volume after a given time interval Δt Amount of f in the control volume at the beginning of the time interval Δt Amount of f that flows into the control volume during Δt = + - Amount of f that flows out of the control volume during Δt Rate of increase of f in a control volume during a given time interval Δt This can be restated in rate form as: = Net rate of flow of f into the control volume during the time interval Δt The mathematical equivalent of this statement is: d dt C f dv =! CS fu nds While we can derive a partial differential expression from this equation, often it is better to work directly with the conservation principle in integral form
Getting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationAE/ME 339. K. M. Isaac. 9/22/2005 Topic 6 FluidFlowEquations_Introduction. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 1...in the phrase computational fluid dynamics the word computational is simply an adjective to fluid dynamics.... -John D. Anderson 2 1 Equations of Fluid
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationSIMPLE Algorithm for Two-Dimensional Channel Flow. Fluid Flow and Heat Transfer
SIMPLE Algorithm for Two-Dimensional Channel Flow Fluid Flow and Heat Transfer by Professor Jung-Yang San Mechanical Engineering Department National Chung Hsing University Two-dimensional, transient, incompressible
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationLecture 3: 1. Lecture 3.
Lecture 3: 1 Lecture 3. Lecture 3: 2 Plan for today Summary of the key points of the last lecture. Review of vector and tensor products : the dot product (or inner product ) and the cross product (or vector
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationIntroduction to Fluid Dynamics
Introduction to Fluid Dynamics Roger K. Smith Skript - auf englisch! Umsonst im Internet http://www.meteo.physik.uni-muenchen.de Wählen: Lehre Manuskripte Download User Name: meteo Password: download Aim
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationIntroduction to Physical Acoustics
Introduction to Physical Acoustics Class webpage CMSC 828D: Algorithms and systems for capture and playback of spatial audio. www.umiacs.umd.edu/~ramani/cmsc828d_audio Send me a test email message with
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationWhere does Bernoulli's Equation come from?
Where does Bernoulli's Equation come from? Introduction By now, you have seen the following equation many times, using it to solve simple fluid problems. P ρ + v + gz = constant (along a streamline) This
More information1 Equations of motion
Part A Fluid Dynamics & Waves Draft date: 21 January 2014 1 1 1 Equations of motion 1.1 Introduction In this section we will derive the equations of motion for an inviscid fluid, that is a fluid with zero
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
More informationIntroduction to Heat and Mass Transfer. Week 10
Introduction to Heat and Mass Transfer Week 10 Concentration Boundary Layer No concentration jump condition requires species adjacent to surface to have same concentration as at the surface Owing to concentration
More informationUNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence
UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Examination in: Day of examination: Friday 15. June 212 Examination hours: 9. 13. This problem set consists of 5 pages. Appendices: Permitted
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
More informationChapter 3: Newtonian Fluids
Chapter 3: Newtonian Fluids CM4650 Michigan Tech Navier-Stokes Equation v vv p t 2 v g 1 Chapter 3: Newtonian Fluid TWO GOALS Derive governing equations (mass and momentum balances Solve governing equations
More informationFundamentals of compressible and viscous flow analysis - Part II
Fundamentals of compressible and viscous flow analysis - Part II Lectures 3, 4, 5 Instantaneous and averaged temperature contours in a shock-boundary layer interaction. Taken from (Pasquariello et al.,
More informationMathematical Theory of Non-Newtonian Fluid
Mathematical Theory of Non-Newtonian Fluid 1. Derivation of the Incompressible Fluid Dynamics 2. Existence of Non-Newtonian Flow and its Dynamics 3. Existence in the Domain with Boundary Hyeong Ohk Bae
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationBasic concepts in viscous flow
Élisabeth Guazzelli and Jeffrey F. Morris with illustrations by Sylvie Pic Adapted from Chapter 1 of Cambridge Texts in Applied Mathematics 1 The fluid dynamic equations Navier-Stokes equations Dimensionless
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationThe Shallow Water Equations
The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationLecture: Wave-induced Momentum Fluxes: Radiation Stresses
Chapter 4 Lecture: Wave-induced Momentum Fluxes: Radiation Stresses Here we derive the wave-induced depth-integrated momentum fluxes, otherwise known as the radiation stress tensor S. These are the 2nd-order
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationMAE 3130: Fluid Mechanics Lecture 7: Differential Analysis/Part 1 Spring Dr. Jason Roney Mechanical and Aerospace Engineering
MAE 3130: Fluid Mechanics Lecture 7: Differential Analysis/Part 1 Spring 2003 Dr. Jason Roney Mechanical and Aerospace Engineering Outline Introduction Kinematics Review Conservation of Mass Stream Function
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More information2. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations
. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations We need to review the governing equations of fluid mechanics before examining the methods of computational aerodynamics in detail.
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationModule 2: Governing Equations and Hypersonic Relations
Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More information7 EQUATIONS OF MOTION FOR AN INVISCID FLUID
7 EQUATIONS OF MOTION FOR AN INISCID FLUID iscosity is a measure of the thickness of a fluid, and its resistance to shearing motions. Honey is difficult to stir because of its high viscosity, whereas water
More informationFluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows
Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 17 Laminar and Turbulent flows Welcome back to the video course on fluid mechanics. In
More informationNumerical Solution of Partial Differential Equations governing compressible flows
Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1
AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationGoverning Equations of Fluid Dynamics
Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present
More information2 Law of conservation of energy
1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationFLUID MECHANICS EQUATIONS
FLUID MECHANIC EQUATION M. Ragheb 11/2/2017 INTRODUCTION The early part of the 18 th -century saw the burgeoning of the field of theoretical fluid mechanics pioneered by Leonhard Euler and the father and
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationComputer Fluid Dynamics E181107
Computer Fluid Dynamics E181107 2181106 Transport equations, Navier Stokes equations Remark: foils with black background could be skipped, they are aimed to the more advanced courses Rudolf Žitný, Ústav
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationChapter 2: Basic Governing Equations
-1 Reynolds Transport Theorem (RTT) - Continuity Equation -3 The Linear Momentum Equation -4 The First Law of Thermodynamics -5 General Equation in Conservative Form -6 General Equation in Non-Conservative
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 5
.9 Numerical Fluid Mechanics Fall 011 Lecture 5 REVIEW Lecture 4 Roots of nonlinear equations: Open Methods Fixed-point Iteration (General method or Picard Iteration), with examples Iteration rule: x g(
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationThe vorticity field. A dust devil
The vorticity field The vector ω = u curl u is twice the local angular velocity in the flow, and is called the vorticity of the flow (from Latin for a whirlpool). Vortex lines are everywhere in the direction
More informationResearch Article Analysis of Mechanical Energy Transport on Free-Falling Wedge during Water-Entry Phase
Journal of Applied Mathematics Volume 212, Article ID 73882, 21 pages doi:1.1155/212/73882 Research Article Analysis of Mechanical Energy Transport on Free-Falling Wedge during Water-Entry Phase Wen-Hua
More information3.5 Vorticity Equation
.0 - Marine Hydrodynamics, Spring 005 Lecture 9.0 - Marine Hydrodynamics Lecture 9 Lecture 9 is structured as follows: In paragraph 3.5 we return to the full Navier-Stokes equations (unsteady, viscous
More informationCONSERVATION OF MASS AND BALANCE OF LINEAR MOMENTUM
CONSERVATION OF MASS AND BALANCE OF LINEAR MOMENTUM Summary of integral theorems Material time derivative Reynolds transport theorem Principle of conservation of mass Principle of balance of linear momentum
More informationTotal energy in volume
General Heat Transfer Equations (Set #3) ChE 1B Fundamental Energy Postulate time rate of change of internal +kinetic energy = rate of heat transfer + surface work transfer (viscous & other deformations)
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
More informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More information6. Laminar and turbulent boundary layers
6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary
More informationThe Orchestra of Partial Differential Equations. Adam Larios
The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline
More informationVortex stretching in incompressible and compressible fluids
Vortex stretching in incompressible and compressible fluids Esteban G. Tabak, Fluid Dynamics II, Spring 00 1 Introduction The primitive form of the incompressible Euler equations is given by ( ) du P =
More informationBasic hydrodynamics. David Gurarie. 1 Newtonian fluids: Euler and Navier-Stokes equations
Basic hydrodynamics David Gurarie 1 Newtonian fluids: Euler and Navier-Stokes equations The basic hydrodynamic equations in the Eulerian form consist of conservation of mass, momentum and energy. We denote
More informationThermodynamics ENGR360-MEP112 LECTURE 7
Thermodynamics ENGR360-MEP11 LECTURE 7 Thermodynamics ENGR360/MEP11 Objectives: 1. Conservation of mass principle.. Conservation of energy principle applied to control volumes (first law of thermodynamics).
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
More informationIntroduction to Turbulence AEEM Why study turbulent flows?
Introduction to Turbulence AEEM 7063-003 Dr. Peter J. Disimile UC-FEST Department of Aerospace Engineering Peter.disimile@uc.edu Intro to Turbulence: C1A Why 1 Most flows encountered in engineering and
More information2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ v = 0;
Accepted for publication in J. Fluid Mech. 1 Viscous Potential Flow By D.D. Joseph Department of Aerospace Engineering and Mechanics, University of Minnesota, MN 55455 USA Email: joseph@aem.umn.edu (Received
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More information