DOCSLIB.ORG

BEM for Compressible Fluid Dynamics

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
BEM for Compressible Fluid Dynamics Transactions on Modelling and Simulation vol 35, © 2003 WIT Press, www.witpress.com, ISSN 1743-355X BEM for compressible fluid dynamics L. gkerget & N. Samec Faculty of Mechanical Engineering, Institute of Power, Process and Environmental Engineering, University of Maribor, Slovenia Abstract The fully developed boundary element method (BEM) numerical model of compressible fluid dynamics is presented. In particular, the singular boundary domain integral approach, which has been established for the viscous incompressible flow problem, is modified and extended to capture the compressible fluid state. As test cases a natural convection of compressible fluid in a closed cavity and an L-shape cavity are studied. 1 Introduction This paper deals with a BEM numerical scheme developed for simulation of motion of compressible viscous fluid. The method is based on the approximate solution of the set of Navier-Stokes equations in the velocity-vorticity formulation. Particular attention is given to proper transformation of the governing differential equations into corresponding integral representations, which satisfy the continuity equation exactly, e.g. velocity and mass density field functions. The pressure field is computed by using the Poisson equation for pressure for known velocity, vorticity and mass density functions. 2 Conservation laws The analytical description of the continuous homogenous medium motion is based on conservation of mass, momentum and energy with associated rheological models and equation of state. The present development will be focused on laminar flow of compressible viscous isotropic fluid in solution domain $2 bounded by boundary T. The field functions of interest are velocity Transactions on Modelling and Simulation vol 35, © 2003 WIT Press, www.witpress.com, ISSN 1743-355X 124 Boutrdary Elcmatr~~XXV vector field vi(rj,t), pressure field p(rj,t), mass density field p(rj,t) and temperature field T(rjt) such that the mass, momentum and energy equations are satisfied, written in the Cartesian tensor notation xi , where p and c,, are fluid mass density and isobaric specific heat per unit mass, t is time, gi is gravitational acceleration vector, qi denotes the components of the total stress tensor, g; stands for heat flux vector. Using the Stokes mass time derivative in general form --- ~t at the following alternative form of the conservation laws can be stated where c is the isobaric specific heat per unit volume, c cpP . -- For an incompressible fluid, the rate of change of mass density following the motion is zero, that is and the mass conservation equation takes simple form V4=0, expressing the solenoidality constraint for the velocity vector. The set of field eqns (4)-(6) has to be closed and solved in conjunction with appropriate rheological models of the fluid and boundary and initial conditions Transactions on Modelling and Simulation vol 35, © 2003 WIT Press, www.witpress.com, ISSN 1743-355X of the flow problem. Boundary conditions in general depend on the dependent variables applied, i.e. primitive or velocity-vorticity variables formulation. For a compressible fluid, the Cauchy total stress aucan be decomposed into a pressure contribution plus an extra deviatoric stress tensor field function 0.. = -p6.. +r.. 11 'I II ' (9) where fiij is the Kronecker delta function. 3 Rheological models In general, real fluid in motion sustains shear stresses. The most general relationship between the extra stress tensor qj and the strain rate tensor E, is given by the Reiner-Rivlin model r0 = a6,i+ P&ii + y2ik&kj, (10) where the coefficients a, p and y are functions of three scalar invariants of strain rate tensor Eij. For a simple viscous shear compressible fluid in motion one can consider the relations a = -2qD / 3 and P = 2q , such that the following constitutive model can be stated n where quantity D = div6 = Cii, and r] is dynamic viscosity. For a most heat transfer problems of practical importance Fouier model of heat diffusion is accurate enough where k is heat conductivity. 4 Summary of governing equations Combining constitutive models for stress tensor and heat diffusion flux, eqns (1 1) and (12) in conservation eqns (5) and (6) the following system of nonlinear equations is developed Transactions on Modelling and Simulation vol 35, © 2003 WIT Press, www.witpress.com, ISSN 1743-355X 126 Bnurzdary Elcmatzt~XXV Using an extended form of the operator divz , i.e. the momentum eqn (14) can be written in the form appropriate for development of the velocity-vorticity formulation DU 4 p- = -grad p -I-pi - rot(r]G)+-grad(q~)+ 2grad6 gradv Dt 3 (17) +2 grad 17 X c3 - 2213 grad 7. 5 Velocity-vorticity formulation The divergence and the curl of a vector field function are fundamental differential operators in vector analysis. Applied to the velocity vector field they give local rate of expansion D and local vorficity vectoor representing a solenoidal vector by definition, the fluid motion computation procedure is partitioned into its kinetics and kinematics [l]. The vorticity transport in fluid domain is governed by non-linear parabolic difhsive-convective equation obtained as a curl ot the momentum eqn (17), i.e. written in general vector form D ac5 1 - -=-+(~V)CZ=V,AW+(U~~V)~~-WD+-VXP, (20) ~t at PO or in Cartesian notation form The pseudo body force vector Fincludes the effects of variable material properties F=~"(~-Z)-~~VXG-GXV~+~VU.V~-~DV~,(22) the following tensor notation form is also valid Transactions on Modelling and Simulation vol 35, © 2003 WIT Press, www.witpress.com, ISSN 1743-355X where a = DV l Dt represents acceleration vector. For the two dimensional plane notation the vorticity vector i7i has just one component perpendicular to the plane of the flow, and it can be treated as a scalar field function. The stretching-twisting term is identically zero, reducing the vector vorticity eqn (20) to a scalar one for the vorticity w where pseudo body force term is am a17 avi aq F;. = p(gi- -a,)-r7e..-+e..m-+2---2D-, 317 (25) rJ lJ axj axj ax, axi or by putting together last three terms represented also as In eqns above, the material properties are considered as a sum of constant and variable part, i.e. 17 =% +ff and p=po+P. (27) Applying the curl operator to the vorticity deffinition PX~~~=TX(V~V)=V(V.~)-A~~, (28) and by using the continuity eqn (13), the following elliptic Poisson eqn is obtained AU+VX~~~-VD=O, (29) or in tensor notation form a *V, am, m +eij, -0 a.,' axj axi The eqn (30) represents the kinematics of a compressible fluid motion expressing the compatibility and restriction conditions between velocity and solenoidal vorticity vector field functions at a given point in space and time. TO accelerate the convergence of the coupled velocity-vorticity iterative scheme the false transient approach is applied. Thus, in the solution scheme the eqn (30) is rewritten as parabolic diffusion eqn for velocity vector with a as a relaxation parameter. It is obvious that the governing velocity eqn (30) is exactly satisfied only at the steady state (t+-) when the false time derivative or false accumulation term vanishes. Transactions on Modelling and Simulation vol 35, © 2003 WIT Press, www.witpress.com, ISSN 1743-355X 6 Pressure equation Let us rewrite momentum cqn (17) for pressure gradient gradp vp=fp=P(X-ri)-~x(rlui)+~~(rl~)+2~~~~0+2~~xoi-2~~~1, (32) 3 whcrc in vector function fr, inertia, gravitational, diffusion and non-linear malcrial effects are incorporated. To derive pressure equation dependent on known field function values the divergence of eqn (32) should bc considered =O. A~-V.&, (33) Considering the normal component of the eqn (32) thc Neuman pressure boundary conditions ar specified. 7 Integral representations The unique advantage of BEM originates from the application of Green fundamental solutions as particular weighting functions. Sincc they only consider the linear transport phenomena, an appropriate selection of a linear differential operator is of main importance in establishing stable and accurate singular integral representations of thc original differential conservation equations. 7.1 Kinematics Consider an integral representation of false transient velocity eqn (31), which can be recognized as a non-homogenous parabolic PDE of the form the following corresponding boundary-domain integral eqn can be obtained by applying weighted residual statement, e.g. written in a time incremental form with a timc step At = tF - tF-, where U* is the parabolic diffusion fundamental solution Transactions on Modelling and Simulation vol 35, © 2003 WIT Press, www.witpress.com, ISSN 1743-355X Equating pseudo body force term b with rotational and compressible part of fluid motion and assuming constant variation of all field functions within individual time increment, one can derive the following integral statement c(<@({,r,)+ J(vu* = ~(vu*xii)xi#+ JGXVU*~Q r I- R (37) DVU*~~+ J"F-lu;-,d~, R R t~ where U* =a ju*dt t~-l The kinematics of planar fluid motion is given by Eqn (37) is equivalent to continuity equation also recognized as compatibility and restriction conditions between velocity and mass density field functions, and vorticity definition expressing the kinematics of general compressible fluid motion in the integral form. Velocity boundary conditions are incorporated in the boundary integrals, while the first two domain integrals express the influence of the local vorticity and expansion rate on the velocity field. The last domain integral is taking into account the influence of initial velocity conditions of false transient phenomena. In eqn (37) normal and tangential derivative of fundamental solution aU * /an and aU * /at are employed. To compute boundary values of field fuctions normal or tangential form of vector eqn (37) [2] is demanded. The boundary vorticity values are expressed in integral form within the domain integral.
Load more

Recommended publications
  • Modeling Conservation of Mass
    Modeling Conservation of Mass How is mass conserved (protected from loss)? Imagine an evening campfire. As the wood burns, you notice that the logs have become a small pile of ashes. What happened? Was the wood destroyed by the fire? A scientific principle called the law of conservation of mass states that matter is neither created nor destroyed. So, what happened to the wood? Think back. Did you observe smoke rising from the fire? When wood burns, atoms in the wood combine with oxygen atoms in the air in a chemical reaction called combustion. The products of this burning reaction are ashes as well as the carbon dioxide and water vapor in smoke. The gases escape into the air. We also know from the law of conservation of mass that the mass of the reactants must equal the mass of all the products. How does that work with the campfire? mass – a measure of how much matter is present in a substance law of conservation of mass – states that the mass of all reactants must equal the mass of all products and that matter is neither created nor destroyed If you could measure the mass of the wood and oxygen before you started the fire and then measure the mass of the smoke and ashes after it burned, what would you find? The total mass of matter after the fire would be the same as the total mass of matter before the fire. Therefore, matter was neither created nor destroyed in the campfire; it just changed form. The same atoms that made up the materials before the reaction were simply rearranged to form the materials left after the reaction.
    [Show full text]
  • Law of Conversation of Energy
    Law of Conservation of Mass: "In any kind of physical or chemical process, mass is neither created nor destroyed - the mass before the process equals the mass after the process." - the total mass of the system does not change, the total mass of the products of a chemical reaction is always the same as the total mass of the original materials. "Physics for scientists and engineers," 4th edition, Vol.1, Raymond A. Serway, Saunders College Publishing, 1996. Ex. 1) When wood burns, mass seems to disappear because some of the products of reaction are gases; if the mass of the original wood is added to the mass of the oxygen that combined with it and if the mass of the resulting ash is added to the mass o the gaseous products, the two sums will turn out exactly equal. 2) Iron increases in weight on rusting because it combines with gases from the air, and the increase in weight is exactly equal to the weight of gas consumed. Out of thousands of reactions that have been tested with accurate chemical balances, no deviation from the law has ever been found. Law of Conversation of Energy: The total energy of a closed system is constant. Matter is neither created nor destroyed – total mass of reactants equals total mass of products You can calculate the change of temp by simply understanding that energy and the mass is conserved - it means that we added the two heat quantities together we can calculate the change of temperature by using the law or measure change of temp and show the conservation of energy E1 + E2 = E3 -> E(universe) = E(System) + E(Surroundings) M1 + M2 = M3 Is T1 + T2 = unknown (No, no law of conservation of temperature, so we have to use the concept of conservation of energy) Total amount of thermal energy in beaker of water in absolute terms as opposed to differential terms (reference point is 0 degrees Kelvin) Knowns: M1, M2, T1, T2 (Kelvin) When add the two together, want to know what T3 and M3 are going to be.
    [Show full text]
  • Key Concepts: Conservation of Mass, Momentum, Energy Fluid: a Material
    Key concepts: Conservation of mass, momentum, energy Fluid: a material that deforms continuously and permanently under the application of a shearing stress, no matter how small. Fluids are either gases or liquids. (Under very specialized conditions, a phase of intermediate properties can be stable, but we won’t consider that possibility.) In liquids, the molecules are relatively closely spaced, allowing the magnitude of their (attractive, electrically-based) interaction energy to be of the same magnitude as their kinetic energy. As a result, they exist as a loose collection of clusters. In gases, the molecules are much more widely separated, so the kinetic energy (at a given temperature, identical to that in the liquid) is far greater than the interaction energy (much less than in the liquid), and molecule do not form clusters. In a liquid, the molecules themselves typically occupy a few percent of the total space available; in a gas, they occupy a few thousandths of a percent. Nevertheless, for our purposes, all fluids are considered to be continua (no voids or holes).The absence of significant intermolecular attraction allows gases to fill whatever volume is available to them, whereas the presence of such attraction in liquids prevents them from doing so. The attractive forces in liquid water are unusually strong, compared to other liquids. Properties of Fluids: m Density is mass/volume: ρ = . The density of liquid water is V 3 o −3 3 ~1.0 kg/m ; that of air at 20 C is ~1.2x10 kg/m . mg W Specific weight is weight/volume: γ = = = ρg V V C:\Adata\CLASNOTE\342\Class Notes\Key concepts_Topic 1.doc 1 γ Specific gravity is density normalized to the density of water: sg..= i γ w V 1 Specific volume is volume/mass: V = = m ρ Bulk modulus or modulus of elasticity is the pressure change per dp dp fractional change in volume or density: E =− = .
    [Show full text]
  • Chapter 15 - Fluid Mechanics Thursday, March 24Th
    Chapter 15 - Fluid Mechanics Thursday, March 24th •Fluids – Static properties • Density and pressure • Hydrostatic equilibrium • Archimedes principle and buoyancy •Fluid Motion • The continuity equation • Bernoulli’s effect •Demonstration, iClicker and example problems Reading: pages 243 to 255 in text book (Chapter 15) Definitions: Density Pressure, ρ , is defined as force per unit area: Mass M ρ = = [Units – kg.m-3] Volume V Definition of mass – 1 kg is the mass of 1 liter (10-3 m3) of pure water. Therefore, density of water given by: Mass 1 kg 3 −3 ρH O = = 3 3 = 10 kg ⋅m 2 Volume 10− m Definitions: Pressure (p ) Pressure, p, is defined as force per unit area: Force F p = = [Units – N.m-2, or Pascal (Pa)] Area A Atmospheric pressure (1 atm.) is equal to 101325 N.m-2. 1 pound per square inch (1 psi) is equal to: 1 psi = 6944 Pa = 0.068 atm 1atm = 14.7 psi Definitions: Pressure (p ) Pressure, p, is defined as force per unit area: Force F p = = [Units – N.m-2, or Pascal (Pa)] Area A Pressure in Fluids Pressure, " p, is defined as force per unit area: # Force F p = = [Units – N.m-2, or Pascal (Pa)] " A8" rea A + $ In the presence of gravity, pressure in a static+ 8" fluid increases with depth. " – This allows an upward pressure force " to balance the downward gravitational force. + " $ – This condition is hydrostatic equilibrium. – Incompressible fluids like liquids have constant density; for them, pressure as a function of depth h is p p gh = 0+ρ p0 = pressure at surface " + Pressure in Fluids Pressure, p, is defined as force per unit area: Force F p = = [Units – N.m-2, or Pascal (Pa)] Area A In the presence of gravity, pressure in a static fluid increases with depth.
    [Show full text]
  • Conservation of Mass
    Conservation of mass Henryk Kudela Contents 1 Principle of Conservation of mass. Transport theorem. 1 1 Principle of Conservation of mass. Transport theorem. The conservation of mass principle is one of the most fundamental principles in nature. We are all familiar with this principle, and it is not difficult to understand. As the saying goes: you cannot have your cake and eat it too!. For closed system Ω(t), which we mean a collection of unchanging contents, so all particles move together inside the region Ω(t) with boundary S(t). The mass of the system can be express by the density of fluid: Msys = ρ(x,t) (1) ZΩ(t) The conservation of mass can be express as follows: Definition 1. Let the density ρ(t,x),x ⊂ Ω, be a positive, smooth function and let v be a vector field with with map motion Φ(t,x). We say ρ,v satisfy the principle of conservation of mass if d ρ(t,x) dυ = 0 (2) dt ZΩ(t) The (2) can be rewrite as: d M = 0 (3) dt sys It is worth to emphasize that Ω(t),ρ and v may change with time, but they must do so in a way that leaves Msys unchanged if the conservation of mass is fulfil. Theorem 1. The principle of conservation of mass is satisfied by ρ,v if and only if any of the following equivalent condition hold: d ρ ρ 1. dt + div v = 0 ∂ρ ρ 2. ∂t + div ( v)= 0 1 ∂ ρ υ ρ 3.
    [Show full text]
  • Fluid Properties and the Conservation of Mass Lab Lecture the Week of Sep 21 Lab Held in Marston 10 the Week of Sep 28
    9/21/2015 Fluid Properties and the Conservation of Mass Lab Lecture the week of Sep 21 Lab held in Marston 10 the week of Sep 28 In-Lab Rules Must wear : Personal Protective Equipment: lab coat, gloves, goggles Closed-toed shoes and pants Label all vessels (beakers, test tubes, flasks): Contents, date, initials, class When in doubt, label! No food or drink in the lab Safety first Handle chemicals with care, always clean broken glass, don’t put flammables in close proximity to flame 1 9/21/2015 Fluid Properties 1. Density of a substance: the quantity of matter contained in a unit volume of the substance Mass density ρ(kg/m3)=m/V Specific weight ω(N/m3)=ρg Relative density σ=ρs/ρH20 2. Viscosity: property of fluid, due to cohesion and interaction between molecules, which offers resistance to deformation. Dynamic viscosity μ Kinematic viscosity ν Reynold’s Number Ratio of the inertial forces (ρv2/L ) to viscous forces (μv/L2) . Re= ρvL/μ=vL/υ=vD/υ ν=μ/ρ μ is the dynamic viscosity of the fluid (kg/(m.s)), v is the maximum velocity of an object relative to a fluid (m/s) or mean fluid velocity, L is the traveled length of the fluid (m) (the symbol D is used sometimes instead of L as the hydraulic diameter), and υ is kinematic viscosity (m2/s). 2 9/21/2015 Conservation of Mass Antoine Lavoisier’s Law (1789): min m mass is neither created nor Control out volume destroyed: min=mout Qin Qout Objectives Measure the density of water Check fluid velocity using Reynold’s number criterion Estimate terms in conservation of water
    [Show full text]
  • 3 Stress and the Balance Principles
    3 Stress and the Balance Principles Three basic laws of physics are discussed in this Chapter: (1) The Law of Conservation of Mass (2) The Balance of Linear Momentum (3) The Balance of Angular Momentum together with the conservation of mechanical energy and the principle of virtual work, which are different versions of (2). (2) and (3) involve the concept of stress, which allows one to describe the action of forces in materials. 315 316 Section 3.1 3.1 Conservation of Mass 3.1.1 Mass and Density Mass is a non-negative scalar measure of a body’s tendency to resist a change in motion. Consider a small volume element Δv whose mass is Δm . Define the average density of this volume element by the ratio Δm ρ = (3.1.1) AVE Δv If p is some point within the volume element, then define the spatial mass density at p to be the limiting value of this ratio as the volume shrinks down to the point, Δm ρ(x,t) = lim Spatial Density (3.1.2) Δv→0 Δv In a real material, the incremental volume element Δv must not actually get too small since then the limit ρ would depend on the atomistic structure of the material; the volume is only allowed to decrease to some minimum value which contains a large number of molecules. The spatial mass density is a representative average obtained by having Δv large compared to the atomic scale, but small compared to a typical length scale of the problem under consideration. The density, as with displacement, velocity, and other quantities, is defined for specific particles of a continuum, and is a continuous function of coordinates and time, ρ = ρ(x,t) .
    [Show full text]
  • 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 Vu= i+ v j + wz . (3.1) Using the Lagrangian approach, we can describe a particle located at point xo = (,xyo o, z o ) for some time t = to, such that the particle velocity is ∂x V = , (3.2) ∂t and particle acceleration is ∂V a = . (3.3) ∂t 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,,yz) 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.
    [Show full text]
  • Introduction to Fluid Flow
    Stanford University Department of Aeronautics and Astronautics AA200A Applied Aerodynamics Chapter 1 - Introduction to fluid flow Stanford University Department of Aeronautics and Astronautics 1.1 Introduction Compressible flows play a crucial role in a vast variety of man-made and natural phenomena. Propulsion and power systems High speed flight Star formation, evolution and death Geysers and geothermal vents Earth meteor and comet impacts Gas processing and pipeline transfer Sound formation and propagation Stanford University Department of Aeronautics and Astronautics 1.2 Conservation of mass Stanford University Department of Aeronautics and Astronautics Divide through by the volume of the control volume. 1.2.1 Conservation of mass - Incompressible flow If the density is constant the continuity equation reduces to Note that this equation applies to both steady and unsteady incompressible flow Stanford University Department of Aeronautics and Astronautics 1.2.2 Index notation and the Einstein convention Make the following replacements Using index notation the continuity equation is Einstein recognized that such sums from vector calculus always involve a repeated index. For convenience he dropped the summation symbol. Coordinate independent form Stanford University Department of Aeronautics and Astronautics 1.3 Particle paths and streamlines in 2-D steady flow The figure below shows the streamlines over a 2-D airfoil. The flow is irrotational and incompressible Stanford University Department of Aeronautics and Astronautics Streamlines Streaklines Stanford University Department of Aeronautics and Astronautics A vector field that satisfies can always be represented as the gradient of a scalar potential or If the vector potential is substituted into the continuity equation the result is Laplaces equation.
    [Show full text]
  • Fluids – Lecture 6 Notes 1
    Fluids – Lecture 6 Notes 1. Control Volumes 2. Mass Conservation 3. Control Volume Applications Reading: Anderson 2.3, 2.4 Control Volumes In developing the equations of aerodynamics we will invoke the firmly established and time- tested physical laws: Conservation of Mass Conservation of Momentum Conservation of Energy Because we are not dealing with isolated point masses, but rather a continuous deformable medium, we will require new conceptual and mathematical techniques to apply these laws correctly. One concept is the control volume, which can be either finite or infinitesimal. Two types of control volumes can be employed: 1) Volume is fixed in space (Eulerian type). Fluid can freely pass through the volume’s boundary. 2) Volume is attached to the fluid (Lagrangian type). Volume is freely carried along with the fluid, and no fluid passes through its boundary. This is essentially the same as the free-body concept employed in solid mechanics. V V Eulerian Control Volume Lagrangian Control Volume fixed in space moving with fluid Both approaches are valid. Here we will focus on the fixed control volume. Mass Conservation Mass flow Consider a small patch of the surface of the fixed, permeable control volume. The patch has 1 area A, and normal unit vectorn ˆ. The plane of fluid particles which are on the surface at time t will move off the surface at time t + ∆t, sweeping out a volume given by ∆V = Vn A ∆t where Vn = V~ · nˆ is the component of the velocity vector normal to the area. V n^ Vn V ^ swept volume n ∆ V ∆ t V t n V ∆ t A A The mass of fluid in this swept volume, which evidently passed through the area during the ∆t interval, is ∆m = ρ ∆V = ρVn A ∆t The mass flow is defined as the time rate of this mass passing though the area.
    [Show full text]
  • Chapter 6 Conservation of Momentum: Fluids and Elastic Solids
    Chapter 6 Conservation of Momentum: Fluids and Elastic Solids The description of the motion of fluids plays a fundamental role in applied mathematics. The basic model is a system of partial di↵erential equations of evolution type. 6.1 Equations of Fluid Motion Fluids consists of molecules; thus, on a microscopic level, a fluid is a discrete material. To obtain a useful approximation, we will describe fluid motion on the macroscopic level by taking into account forces that act on a parcel of fluid, which we assume to be a collection of sufficiently many molecules of the fluid so that the continuity assumption is valid. More precisely, we will assume the fluid to be a continuous medium contained in three-dimensional space R3 such that every parcel of the fluid, no matter how small in compar- ison with the whole body of fluid, can be viewed as a continuous material. Mathematically, a parcel of fluid (at every moment of time) is a bounded open subset A of the fluid that is assumed to be in an open set whose closure is the region containing the fluid. To ensure correctness of theR math- ematical operations that follow, we assume that each parcel A has a C1 boundary. At every given moment of time, a particle of fluid is identified with a point in . R Let ⇢ = ⇢(x, t)denotethedensityandu = u(x, t)thevelocityofthefluid at position x and time t R. The motion of the fluid is modeled by 2R 2 257 258 Fluids di↵erential equations for ⇢ and u determined by the fluid’s internal material properties, its container, and the external forces acting on the fluid.
    [Show full text]
  • Lecture 1 Conservation Equations
    Summer 2013 Lecture 1 Conservation Equations Moshe Matalon Conservation equations for a pure substance We focus on a material volume Vm(t) (made of the same set of material n n points). d ⇢dV =0 dt Vm(t1) Vm(t2) VZ(t) Mass Conservation The time rate of a change of a quantity φ(x,t) over a given volume V,moving with velocity VI is d @ φdV = dV + φ(V n) dS dt @t I · VZ(t) VZ(t) SZ(t) This is the multidimensional analogue of Leibnitz’s theorem from calculus; it is a kinematic relation, not a law of fluid mechanics. Moshe Matalon Moshe Matalon University of Illinois at Urbana-Champaign 1 Summer 2013 When VI = v,i.e.,theCVmoveswiththefluidvelocity,theresultisthe Reynolds’ transport theorem D @ φdV = dV + φ v n dS Dt @t · VZ(t) VZ(t) S(Zt) D @ φdV = + (φv) dV Dt @t r· ZV(t) ZV(t) ⇢ D Dφ φdV = + φ v dV Dt Dt r· ZV(t) ZV(t) ⇢ and the common notation for the “time rate of change” is D/Dt,theconvective or material derivative. For example: time rate of change of the D⇢ @⇢ = + v ⇢ density of a fluid particle Dt @t ·r acceleration of a fluid Dv @v = + v v particle Dt @t ·r Moshe Matalon Conservation of mass D @⇢ ⇢dV =0 + ⇢v dV =0 Dt @t r· ZV(t) ) ZV ⇢ @⇢ + ⇢v =0 @t r· @⇢ + v ⇢ +⇢ v =0 D⇢ @t ·r r· + ⇢ v =0 ) Dt r· D⇢/Dt | {z } 1 D⇢ v = r· −⇢ Dt rate of compression or rate of dilation D⇢ =0 v =0 Dt )r· Moshe Matalon Moshe Matalon University of Illinois at Urbana-Champaign 2 Summer 2013 Conservation of momentum Newton’s second law can be applied directly to a material C.V.
    [Show full text]