0.2. CONSERVATION LAW FOR FLUID 9
|
|
- Erika Bailey
- 5 years ago
- Views:
Transcription
1 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 the x-component force exerted by the pressure of the neighboring fluid. Applying Green s theorem, we have D(ρu) + ρu( v) dv t = For y and z components, we have D(ρv) + ρv( v) dv t = D(ρw) + ρw( v) dv t = Combining the three equations, we have Since from mass conservation we have D(ρv) In the differential form or more explicitly ( ρg x p ) dv t. (28) x ( ρg y p ) dv t, (29) y ( ρg z p ) dv t. (30) z D(ρv) + ρv( v) dv t = (ρg p) dv t. (31) + ρv( v) = ρ Dv + v ( Dρ + ρ( v) ) = ρ Dv, ρ ( ρ Dv ) + p dv t = ρgdv t. (32) ( v ρ Dv + (v )v + p = ρg, (33) ) + p = ρg (34)
2 10 CONTENTS Conservation of Energy In a fluid, the energy consists of two parts, kinetic energy and internal energy (thermal energy). Assuming that e is the specific internal energy (internal energy per unit mass), the total energy density E (total energy per unit volume) is E = 1 2 ρ v 2 +ρe. The energy flux is, therefore, ( 1 2 ρ v 2 +ρe)v. For a fixed volume, the conservation of energy requires that the rate of the energy change equals the net influx of energy plus the work done at the boundary by neighboring fluid which equals -pv ds per unit area and the work done by bodily force. Thus the energy conservation gives us ( 1 ( 1 W 2 ρ v 2 + ρe) dv = S 2 ρ v 2 + ρe + p) v ds+ ρg vdv. (35) W From Green s theorem, we have W ( E ) + ((E + p)v) dv = W ρg vdv. (36) Since W is aritary, we have the differential form of the energy conservation as E + ((E + p)v) = ρg v. (37) Incompressible Fluid For an incompressible fluid, if one follows the fluid particle in the moving volume, this volume should not be changed dv t = JdV = W W J dv = 0. (38) This requires J = 0. (39) Another way to look at the incompressible fluid is that if one follows the particle motion, the fluid density should not be changed, that is Dρ/ = 0. From conservation of mass (Eq. (25)), this is equivalent to v = 0, (40)
3 0.2. CONSERVATION LAW FOR FLUID 11 therefore the complete set of equations for incompressible fluid is if density is uniform, it is simply ρ + (v )ρ = 0 ( ) v ρ + (v )v = p + ρg (41) v = 0 ρ ( v Isentropic Flow + (v )v ) v = 0 = p + ρg (42) In thermodynamics, the enthalpy of the fluid can be written as dw = T ds + V dp, (43) where T is the absolute temperature, S is entropy, and V is the specific volume V = 1. For isentropic flow, the entropy S is a constant, therefore we ρ have dw = V dp = dp p. It follows that = w, and the momentum equation ρ ρ becomes v + (v )v = w + g. (44) Since normally g can be expressed as the gradient of a potential function as g = φ, and (v )v = 1 2 v 2 v v, we have v v 2 v v = (w + φ). (45) Taking the curl of both sides, we have or v ω = (v v), (46) = (v ω), (47)
4 12 CONTENTS where ω is the vorticity of the fluid. If the fluid is initially irrotational, it will remain irrotational because of Eq. (47). In such case, we simply have For steady flow, / = 0, we have Integrating over space, we have Eq. (50) is called Bernoulli s law. v v 2 = (w + φ). (48) 1 2 v 2 + w + φ = 0. (49) 1 2 v 2 + w + φ = Constant, (50) 0.3 Conservation Law and Wave Equations Linear Wave Equation The general form of first order scalar wave equation has the following form u t + c(u)u x = 0, (51) Here the function c(u) serve as the propagation speed of hyperbolic wave. If c is independent of u, the equation is linear, otherwise it is nonlinear. When compared with conservation law Eq. (8), we have c(u) = F (u). For linear equation F (u) = au and c(u) = a, where a is a constant. We start from the linear wave equation: Du = u t + au x = 0 or u t + (au) x = 0 (52) and discuss the following aspects of the equation (1). The traveling wave solution in (, ): u(x, 0) = φ(x), u(x, t) = φ(x at), (53) where a is the speed of the traveling wave.
5 0.3. CONSERVATION LAW AND WAVE EQUATIONS 13 (2). The characteristics of the equation: dx dt = a, along which du dt = u t + au x = 0. (54) This tells us that along each characteristics, u is a constant. (3). Boundary condition: if a > 0, boundary condition can only be applied from the left side, if a < 0 boundary condition can only be applied from the right side. (4). Discontinuity: with smooth initial condition, linear equation will not produce discontinuity, but the traveling wave solution does allow discontinuous initial condition. The discontinuous solution does NOT satisfy the original equation at the discontinuity. Most of the physical waves, including sound wave, electromagnetic wave, etc, are described by the second order wave equation u tt a 2 u xx = 0 (55) We can factorize the differential operator 2 a 2 2 into ( a )( +a 2 x 2 x In the one dimensional infinite domain, the d Alembert solution gives: x ). u(x, t) = φ(x + at) + ψ(x at) (56) corresponding to the left-traveling and the right-traveling waves with characteristics dx = ±a (57) dt A typical initial value problem in the infinite domain is The d Alembert solution must satisfy u tt a 2 u xx = 0, (58) u(x, 0) = g(x), u t (x, 0) = h(x). (59) φ(x) + ψ(x) = g(x), aφ (x) aψ (x) = h(x). (60) After integration, we have u(x, t) = 1 1 x+at [g(x + at) + g(x at)] + h(ξ)dξ. (61) 2 2a x at
6 14 CONTENTS For solutions in finite domain, boundary condition must be considered. Two of the most frequently used boundary conditions are: (1) fixed boundary condition u(x B, t) = 0, and (2) free boundary condition u x (x B, t) = 0. Consider solution in a semi-infinite domain with a boundary at x B = 0, for the fix boundary condition (1), the solution should be an odd function about the point x = 0, while for free boundary condition (2), the solution should be an even function about the point x = Sound Wave Consider one dimensional Euler equations with no external force ρ t + (ρu) x = 0 (62) ρu t + ρuu x + p x = 0, (63) for small oscillation, we may assume p = p 0 + δp, ρ = ρ 0 + δρ and u = δu. Substitute this into Eqs. (62-63) and neglect all second order terms, we have δρ t + ρ 0 δu x = 0 (64) ρ 0 δu t + δp x = 0. (65) Now assume that p = p(ρ), we have δp x = p (ρ)ρ x. Substituting this into Eqs. (64-65), we have δρ t + ρ 0 δu x = 0 (66) ρ 0 δu t + δp x = 0. (67) Differentiating Eq. (66) with respect to t and Eq. (67) with respect to x and substituting the second one into the first one, we have δρ tt p (ρ)δρ xx = 0. (68) This gives the sound speed square a 2 = p (ρ). During the small oscillation of sound, the pressure responds to the change of density adiabatically, therefore, we have γp0 a = p (ρ) =, (69) ρ 0 where γ is the gas constant, p 0 is the ambient air pressure, and ρ 0 is the ambient air density.
7 0.3. CONSERVATION LAW AND WAVE EQUATIONS Numerical Diffusion Linear diffusion Consider viscous linear equations u t + au x = νu xx (70) with discontinuous initial condition u(x, 0) = { ul if x < 0 u r if x > 0 (71) A substitution of X = x at will transform Eq. (70) into a diffusion equation u t = νu XX (72) Diffusion equation For the parabolic diffusion equation u t = νu XX, and initial condition u(x, 0) = u 0 (X), (73) using Fourier transform here û(k, t) = 1 2π u(x, t)e ikx dx (74) û t = νk 2 û, û(k, t) = û 0 (k)e νk2t, (75) û 0 (k) = 1 2π is the Fourier transform of initial condition. u(x, t) = 1 2π u(x, 0)e ikx dx (76) û(k, t)e ikx dk = 1 u(y, 0)e νk2 t+ik(x Y ) dkdy. 2π (77) Integrating over k (see Appendix), we obtain u(x, t) = 1 4πνt u 0 (Y )e (X Y )2 /4νt dy (78) Solution to singular initial condition u 0 (X) = δ(x), u(x, t) = 1 4πνt e X2 /4νt (79)
8 46 CONTENTS Birkoff and Carter used conformal mapping to predict a 2D bubble has C 1 = Garabedian estimated C Numerical experiments have been in support to such terminal bubble velocities. Youngs and the α Paradox Read [12] and Youngs [15, 16, 17] first presented their expreiment of Rayleigh- Taylor instability with randomly perturbed interface and claimed that the perturbed interface will grow under constant acceleration. h h 0 = αagt 2 (192) His experimental value of α was in the range of , but numerical simulation failed to support this result. Front tracking simulation has a much closer agreement with the experiment. Sharp-Wheeler Merger Model Researchers have tried to explain why the random surface RT instability has constant acceleration while single bubble has constant velocity. The first premitive model was given by Sharp and Wheeler [8] as the so-called bubble merger model. With the merger criterion h 2 h 1 m(r 1 ) = C 2 r 1 (193) with appropriate choice (some how arbitary) of C 2, numerical simulation has suggested Youngs s constant acceleration. More sophisticated models include the renormalization model, buoyancydrag mode, average equation model etc. It is still an open area for imagination and creative thinking Richtmyer-Meshkov Instability The Impulsive Model The impulsive model was proposed by Richtmyer in In this model, the growth rate of the instability ȧ 0 (t) is given by a simple formula ȧ 0 (t) = kū ρ 1 ρ 2 ρ 1 + ρ 2 a 0 (0), (194)
0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
More informationGetting 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 informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
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 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 information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
More information4.1 LAWS OF MECHANICS - Review
4.1 LAWS OF MECHANICS - Review Ch4 9 SYSTEM System: Moving Fluid Definitions: System is defined as an arbitrary quantity of mass of fixed identity. Surrounding is everything external to this system. Boundary
More informationDiffusion Equation and Mean Free Path
Diffusion Equation and Mean Free Path Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Content General Introduction Analytic
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 informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
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 informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
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 informationAST242 LECTURE NOTES PART 5
AST242 LECTURE NOTES PART 5 Contents 1. Waves and instabilities 1 1.1. Sound waves compressive waves in 1D 1 2. Jeans Instability 5 3. Stratified Fluid Flows Waves or Instabilities on a Fluid Boundary
More information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
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 informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
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 informationMAE 101A. Homework 7 - Solutions 3/12/2018
MAE 101A Homework 7 - Solutions 3/12/2018 Munson 6.31: The stream function for a two-dimensional, nonviscous, incompressible flow field is given by the expression ψ = 2(x y) where the stream function has
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 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 informationContinuum Mechanics Lecture 5 Ideal fluids
Continuum Mechanics Lecture 5 Ideal fluids Prof. http://www.itp.uzh.ch/~teyssier Outline - Helmholtz decomposition - Divergence and curl theorem - Kelvin s circulation theorem - The vorticity equation
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationSection 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow
Section 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow Anderson: Chapter 2 pp. 41-54 1 Equation of State: Section 1 Review p = R g T " > R g = R u M w - R u = 8314.4126
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 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 informationVorticity Equation Marine Hydrodynamics Lecture 9. Return to viscous incompressible flow. N-S equation: v. Now: v = v + = 0 incompressible
13.01 Marine Hydrodynamics, Fall 004 Lecture 9 Copyright c 004 MIT - Department of Ocean Engineering, All rights reserved. Vorticity Equation 13.01 - Marine Hydrodynamics Lecture 9 Return to viscous incompressible
More information5. Coupling of Chemical Kinetics & Thermodynamics
5. Coupling of Chemical Kinetics & Thermodynamics Objectives of this section: Thermodynamics: Initial and final states are considered: - Adiabatic flame temperature - Equilibrium composition of products
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 018 Properties of Sound Sound Waves Requires medium for propagation Mainly
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
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 informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
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 informationGasdynamics 1-D compressible, inviscid, stationary, adiabatic flows
Gasdynamics 1-D compressible, inviscid, stationary, adiabatic flows 1st law of thermodynamics ρ const Kontrollfläche 1 2 m u 2 u 1 z Q 12 +P 12 = ṁ } h 2 h {{} 1 Enthalpy Q 12 + 1 2 (u2 2 u2 1 }{{} ) +
More informationSimple waves and a characteristic decomposition of the two dimensional compressible Euler equations
Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Jiequan Li 1 Department of Mathematics, Capital Normal University, Beijing, 100037 Tong Zhang Institute
More information(Jim You have a note for yourself here that reads Fill in full derivation, this is a sloppy treatment ).
Lecture. dministration Collect problem set. Distribute problem set due October 3, 004.. nd law of thermodynamics (Jim You have a note for yourself here that reads Fill in full derivation, this is a sloppy
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
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 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 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 informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationSOLAR MHD Lecture 2 Plan
SOLAR MHD Lecture Plan Magnetostatic Equilibrium ü Structure of Magnetic Flux Tubes ü Force-free fields Waves in a homogenous magnetized medium ü Linearized wave equation ü Alfvén wave ü Magnetoacoustic
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 informationOn Fluid Maxwell Equations
On Fluid Maxwell Equations Tsutomu Kambe, (Former Professor, Physics), University of Tokyo, Abstract Fluid mechanics is a field theory of Newtonian mechanics of Galilean symmetry, concerned with fluid
More informationFundamentals of Atmospheric Modelling
M.Sc. in Computational Science Fundamentals of Atmospheric Modelling Peter Lynch, Met Éireann Mathematical Computation Laboratory (Opp. Room 30) Dept. of Maths. Physics, UCD, Belfield. January April, 2004.
More informationPhysics Dec Time Independent Solutions of the Diffusion Equation
Physics 301 10-Dec-2004 33-1 Time Independent Solutions of the Diffusion Equation In some cases we ll be interested in the time independent solution of the diffusion equation Why would be interested in
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
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 informationBorn of the Sound Wave Equation
Corso di Laurea in Fisica - UNITS ISTITUZIONI DI FISICA PER IL SISTEMA TERRA Born of the Sound Wave Equation FABIO ROMANELLI Department of Mathematics & Geosciences University of Trieste romanel@units.it
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
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 informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
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 informationConservation and dissipation principles for PDEs
Conservation and dissipation principles for PDEs Modeling through conservation laws The notion of conservation - of number, energy, mass, momentum - is a fundamental principle that can be used to derive
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 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 informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationBERNOULLI EQUATION. The motion of a fluid is usually extremely complex.
BERNOULLI EQUATION The motion of a fluid is usually extremely complex. The study of a fluid at rest, or in relative equilibrium, was simplified by the absence of shear stress, but when a fluid flows over
More informationBoundary conditions. Diffusion 2: Boundary conditions, long time behavior
Boundary conditions In a domain Ω one has to add boundary conditions to the heat (or diffusion) equation: 1. u(x, t) = φ for x Ω. Temperature given at the boundary. Also density given at the boundary.
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
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 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 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 information1 Energy dissipation in astrophysical plasmas
1 1 Energy dissipation in astrophysical plasmas The following presentation should give a summary of possible mechanisms, that can give rise to temperatures in astrophysical plasmas. It will be classified
More informationModels of ocean circulation are all based on the equations of motion.
Equations of motion Models of ocean circulation are all based on the equations of motion. Only in simple cases the equations of motion can be solved analytically, usually they must be solved numerically.
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationOffshore Hydromechanics Module 1
Offshore Hydromechanics Module 1 Dr. ir. Pepijn de Jong 4. Potential Flows part 2 Introduction Topics of Module 1 Problems of interest Chapter 1 Hydrostatics Chapter 2 Floating stability Chapter 2 Constant
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 information1 POTENTIAL FLOW THEORY Formulation of the seakeeping problem
1 POTENTIAL FLOW THEORY Formulation of the seakeeping problem Objective of the Chapter: Formulation of the potential flow around the hull of a ship advancing and oscillationg in waves Results of the Chapter:
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationFundamentals of Acoustics
Fundamentals of Acoustics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences
More informationexamples of equations: what and why intrinsic view, physical origin, probability, geometry
Lecture 1 Introduction examples of equations: what and why intrinsic view, physical origin, probability, geometry Intrinsic/abstract F ( x, Du, D u, D 3 u, = 0 Recall algebraic equations such as linear
More informationThe Euler equations in fluid mechanics
The Euler equations in fluid mechanics Jordan Bell jordan.bell@gmail.com Department of Mathematics, niversity of Toronto April 14, 014 1 Continuity equation Let Ω be a domain in R n and let C (Ω R); perhaps
More informationIntroduction to Gas Dynamics All Lecture Slides
Introduction to Gas Dynamics All Lecture Slides Teknillinen Korkeakoulu / Helsinki University of Technology Autumn 009 1 Compressible flow Zeroth law of thermodynamics 3 First law of thermodynamics 4 Equation
More information5.1 Fluid momentum equation Hydrostatics Archimedes theorem The vorticity equation... 42
Chapter 5 Euler s equation Contents 5.1 Fluid momentum equation........................ 39 5. Hydrostatics................................ 40 5.3 Archimedes theorem........................... 41 5.4 The
More informationRocket Thermodynamics
Rocket Thermodynamics PROFESSOR CHRIS CHATWIN LECTURE FOR SATELLITE AND SPACE SYSTEMS MSC UNIVERSITY OF SUSSEX SCHOOL OF ENGINEERING & INFORMATICS 25 TH APRIL 2017 Thermodynamics of Chemical Rockets ΣForce
More informationSummary of the Equations of Fluid Dynamics
Reference: Summary of the Equations of Fluid Dynamics Fluid Mechanics, L.D. Landau & E.M. Lifshitz 1 Introduction Emission processes give us diagnostics with which to estimate important parameters, such
More informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationTurbomachinery & Turbulence. Lecture 2: One dimensional thermodynamics.
Turbomachinery & Turbulence. Lecture 2: One dimensional thermodynamics. F. Ravelet Laboratoire DynFluid, Arts et Metiers-ParisTech February 3, 2016 Control volume Global balance equations in open systems
More informationTime Evolution of Speed
ENERGY & MOMENTUM Newton's 2 nd Law: Time Evolution of Speed Describes how force affects velocity of the CM How do forces affect speed of the CM? Depends on relative direction of force and velocity: F
More informationSAMPLE CHAPTERS UNESCO EOLSS WAVES IN THE OCEANS. Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany
WAVES IN THE OCEANS Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany Keywords: Wind waves, dispersion, internal waves, inertial oscillations, inertial waves,
More informationPart IB Fluid Dynamics
Part IB Fluid Dynamics Based on lectures by P. F. Linden Notes taken by Dexter Chua Lent 2016 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after lectures.
More information(The) Three Linear Partial Differential Equations
(The) Three Linear Partial Differential Equations 1 Introduction A partial differential equation (PDE) is an equation of a function of 2 or more variables, involving 2 or more partial derivatives in different
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More information2.25 Advanced Fluid Mechanics
MIT Department of Mechanical Engineering.5 Advanced Fluid Mechanics Problem 4.05 This problem is from Advanced Fluid Mechanics Problems by A.H. Shapiro and A.A. Sonin Consider the frictionless, steady
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 informationIntroduction to Aerospace Engineering
Introduction to Aerospace Engineering Lecture slides Challenge the future 3-0-0 Introduction to Aerospace Engineering Aerodynamics 5 & 6 Prof. H. Bijl ir. N. Timmer Delft University of Technology 5. Compressibility
More informationRational derivation of the Boussinesq approximation
Rational derivation of the Boussinesq approximation Kiyoshi Maruyama Department of Earth and Ocean Sciences, National Defense Academy, Yokosuka, Kanagawa 239-8686, Japan February 22, 2019 Abstract This
More informationWaves in Uniform Flow on Half Plane
May 1, 6 HB Ju Waves in Uniform Flow on Half Plane Hongbin Ju Department of Mathematics Florida State University, Tallahassee, FL.336 www.aeroacoustics.info Please send comments to: hju@math.fsu.edu The
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
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 informationF1.9AB2 1. r 2 θ2 + sin 2 α. and. p θ = mr 2 θ. p2 θ. (d) In light of the information in part (c) above, we can express the Hamiltonian in the form
F1.9AB2 1 Question 1 (20 Marks) A cone of semi-angle α has its axis vertical and vertex downwards, as in Figure 1 (overleaf). A point mass m slides without friction on the inside of the cone under the
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 27 May, 2004 1.30 to 3.30 PAPER 64 ASTROPHYSICAL FLUID DYNAMICS Attempt THREE questions. There are four questions in total. The questions carry equal weight. Candidates
More informationF11AE1 1. C = ρν r r. r u z r
F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions
More informationTopics in Relativistic Astrophysics
Topics in Relativistic Astrophysics John Friedman ICTP/SAIFR Advanced School in General Relativity Parker Center for Gravitation, Cosmology, and Astrophysics Part I: General relativistic perfect fluids
More information