Webster s horn model on Bernoulli flow

Size: px
Start display at page:

Download "Webster s horn model on Bernoulli flow"

Transcription

1 Webster s horn model on Bernoulli flow Aalto University, Dept. Mathematics and Systems Analysis January 5th, 2018

2 Incompressible, steady Bernoulli principle Consider a straight tube Ω R 3 havin circular intersections. By A(x) denote the areas of the intersection along the length variable x R. Let p 0 and v 0 denote the pressure and flow velocity as a function of x. Bernoulli principle: p 0 (x) + ρ 0 2 v 0(x) 2 = p stag, x R where p stag denotes the stagnation pressure. Steadyness everything is independent of time t. Incompressibility density ρ 0 is independent of x, too. Then the volume flow U 0 > 0 is constant and v 0 (x) = U 0 A(x), x R.

3 Lossless Webster s horn model For small pressure and velocity variations (p, v ) in a tubular domain Ω containing stationary fluid, the acoustic approximation for the longitudinal dynamics can be used: 2 φ t 2 c2 A(x) ( A(x) φ ) = 0. The sound pressure and the acoustic perturbation velocity are obtained from the velocity potential φ = φ(t, x) by the partial derivatives v = φ and p φ = ρ 0 t. This model can be derived for the Wave Equation in a tubular domain by averaging the 3D velocity potential ψ = ψ(t, x, y, z) over the transversal intersections having the area A(x).

4 Purpose of this work If the fluid column is moving in the tube Ω with a constant volume velocity U 0 in the longitudinal direction (i.e., x-axis), the fluid movement couples with the propagation of longitudinal acoustic waves. At narrow parts of the tube, the fluid velocity is higher, and the acoustic field gets transported faster. Can we write a PDE that would encompass both the incompressible Bernoulli flow and the Webster horn model at once? Doing so would help us studying the coupling between the transport and acoustic subsystems if the PDE can be kept simple enough for analysis. Let us see about that.

5 Equation of continuity (1) We assume that ρ = ρ 0 + ρ, v = v 0 i + v, p = p 0 + p where ρ 0 > 0 is constant, and v 0 = v 0 (x), p 0 = p 0 (x) are as above. The unit vector i is the direction of the steady flow. The perturbations ρ, v, p depend on both t and (x, y, z) Ω R 3. They define the acoustic field. Since the mass (i.e., the fluid volume) is conserved, we get by omitting the second order perturbation terms the Equation of Continuity 1 ρ ρ 0 t + v + v 0 = 0.

6 Euler s equation (2) Omitting the second order advection term and perturbations in the density, the Conservation of the Momentum takes the form v (v t + v ) 0 + v v 0 x i + 1 p = 0 ρ 0 where v x = v i.

7 Thermodynamics and the speed of sound (1) Assuming the locally isentropic thermodynamics, the perturbation pressure and density are connected by p = c 2 ρ for the speed of sound 0 < c <. Of course, having c < requires compressibility. Neglecting the density and temperature variation (i.e., compressibility) but taking into consideration the pressure variation in the underlying flow, we get c 2 = γ p 0 ρ 0 = c 2 0 γu2 0 2A(x) 2 0. Here c0 2 = γ pstag ρ 0 is the speed of sound at stagnation, and γ is the adiabatic index of the gas. For diatomic ideal gas, γ = 7/5.

8 Thermodynamics and the speed of sound (2) We can now eliminate ρ from the Equation of Continuity, and hence 1 p ρ 0 c 2 t + v + v 0 = 0. It does not complicate the final model much to have c = c(x). It is of some interest to compute the lower limit for constriction A(x) 0 where the underlying flow remains subsonic... but not today. You cannot see Mach = 1 condition from this argument since the density, the absolute temperature, and the speed of sound decrease significantly when A(x) 0 then you must consider the underlying steady flow as compressible to make sense.

9 Wave equation for pressure and velocity Proceeding as usual, the combination of the Continuity Equation, Euler s Equation and Isentropic Equation of State produce a variant of the Wave Equation for the sound pressure p and acoustic perturbation velocity component v x in the direction of the flow: 2 p t 2 c2 2 p + v 0 2 p t = ρ 0c 2 ( 2 v x ) v 0 + (v x v 0 ) 2 v 0 2. Bad news: It is not possible to directly eliminate v x = v i from this equation since the Continuity Equation only gives radially symmetric information about the divergence of v. Good news: We can do tricks by defining a velocity potential candidate and writing the wave equation in terms of it.

10 Velocity potential (1) Taking the i-component from Euler s Equation gives v x t + ) (v 0 v x + p = 0. (1) Note that both v x and v 0 v x + p ρ 0 are functions of t, x, y, z even though the differentiations concern only the variables x, t. For any twice differentiable ψ = ψ(t, x, y, z) satisfying v x = ψ we get Eq. (1) since vx t + ρ 0 and v 0 v x + p = ψ ρ 0 t ( ) v 0 v x + p ρ 0 = 2 ψ t + 2 ψ t = 0. Thus, any solution pair (v x, p) of Eq. (1) can be represented in terms of such ψ, but it is not a true velocity potential in R 3 since we have instead of there.

11 Averaged velocity potential For any f : Ω R, the planar averages are defined by the integral f (x 0 ) = 1 f (x 0, y, z) dx dy. A(x 0 ) Γ(x 0 ) By averaging the i-component of Euler s Equation, v x t + ) (v 0 v x + p = 0 (2) where v x = v x (t, x) and p = p (t, x) come from (v x, p ). Now Eq. (2) allows the 1D velocity potential φ = φ(t, x) v x = φ ( ) ( ) φ φ and p = ρ 0 t v 0 v x = ρ 0 t + v φ 0. ρ 0

12 Averaging the Wave Equation Similarly, averaging the Wave Equation for (v x, p ) introduced above, we get ) (A(x) p 2 p + v 0 t 2 p t 2 c2 A(x) ( = ρ 0 c 2 2 v x v 0 + ( v x v 0 ) 2 v 0 2 Observe that Webster s operator on the LHS W A := 1 ( A(x) ) A(x) ). is what you get by averaging the Laplacian on a tube Ω R 3 with the area function A(x) in the longitudinal direction.

13 Bernoulli Webster model (1) If (v x, p ) satisfy the Wave Equation, what is the PDE satisfied by the 1D velocity potential φ = φ(t, x)? There are several variants. First variant: ( 2 t 2 2 ) ( ) φ c2 W A + v 0 t t + v φ 0 c 2 v 0 B(x) ( W 1/A φ + W 1/B φ ) = c 2 v0 2 ( B (x) B(x) 2) where B(x) := A (x) A(x) is the logarithmic derivative of the area function.

14 Bernoulli Webster model (2) Second variant: ( ) ( t + v 2 φ 0 t 2 2 ) φ c2 W A φ + v 0 t c 2 v 0 B(x) ( 3W B φ + W 1/A φ ) = c 2 v0 2 ( B (x) B(x) 2). And then there are others, but none of them looks any nicer. Note that the equations become homogenous if B (x) = B(x) 2 which takes place if and only if A(x) = a b x for a, b R, a > 0. This defines a hyperbolic horn having its flare at x = b.

15 The combined flow and acoustic field If φ = φ(t, x) is a solution of the Bernoulli Webster model, then the acoustically perturbed steady Bernoulli flow is described by ( ρ(t, x) = ρ ( φ c 2 t (t, x) + U )) 0 φ (t, x), A(x) v(t, x) = U 0 A(x) φ (t, x), and p(t, x) = p stag + ρ 0 ( φ t (t, x) + U 0 A(x) ( φ (t, x) U )) 0. 2A(x)

16 Further developments: 1 The Lagrangian image: i.e., the acoustics in the coordinate system that moves along the flow? (Leads to complications at least if the model is not simplified further.) 2 The underlying steady flow modelled by the compressible Bernoulli principle? (Taking the variations in density and temperature of the underlying flow into account seems technically complicated however, necessary if the constrictions are severe; i.e., A(x) 0 for some x.) 3 Magnetohydrodynamical variants? Acoustical compression waves in moving plasma bounded by a static magnetic field? So little time, so much to do.

17 That s all, folks. Questions?

The Euler Equation of Gas-Dynamics

The 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 information

PHYS 643 Week 4: Compressible fluids Sound waves and shocks

PHYS 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 information

Introduction to Fluid Mechanics. Chapter 13 Compressible Flow. Fox, Pritchard, & McDonald

Introduction to Fluid Mechanics. Chapter 13 Compressible Flow. Fox, Pritchard, & McDonald Introduction to Fluid Mechanics Chapter 13 Compressible Flow Main Topics Basic Equations for One-Dimensional Compressible Flow Isentropic Flow of an Ideal Gas Area Variation Flow in a Constant Area Duct

More information

Continuum Mechanics Lecture 5 Ideal fluids

Continuum 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 information

Notes #4a MAE 533, Fluid Mechanics

Notes #4a MAE 533, Fluid Mechanics Notes #4a MAE 533, Fluid Mechanics S. H. Lam lam@princeton.edu http://www.princeton.edu/ lam October 23, 1998 1 The One-dimensional Continuity Equation The one-dimensional steady flow continuity equation

More information

Introduction to Turbomachinery

Introduction to Turbomachinery 1. Coordinate System Introduction to Turbomachinery Since there are stationary and rotating blades in turbomachines, they tend to form a cylindrical form, represented in three directions; 1. Axial 2. Radial

More information

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION

CHAPTER 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 information

AA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations

AA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations AA210A Fundamentals of Compressible Flow Chapter 5 -The conservation equations 1 5.1 Leibniz rule for differentiation of integrals Differentiation under the integral sign. According to the fundamental

More information

Chapter Two. Basic Thermodynamics, Fluid Mechanics: Definitions of Efficiency. Laith Batarseh

Chapter Two. Basic Thermodynamics, Fluid Mechanics: Definitions of Efficiency. Laith Batarseh Chapter Two Basic Thermodynamics, Fluid Mechanics: Definitions of Efficiency Laith Batarseh The equation of continuity Most analyses in this book are limited to one-dimensional steady flows where the velocity

More information

Macroscopic plasma description

Macroscopic plasma description Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion

More information

Compressible Duct Flow with Friction

Compressible Duct Flow with Friction Compressible Duct Flow with Friction We treat only the effect of friction, neglecting area change and heat transfer. The basic assumptions are 1. Steady one-dimensional adiabatic flow 2. Perfect gas with

More information

In which of the following scenarios is applying the following form of Bernoulli s equation: steady, inviscid, uniform stream of water. Ma = 0.

In which of the following scenarios is applying the following form of Bernoulli s equation: steady, inviscid, uniform stream of water. Ma = 0. bernoulli_11 In which of the following scenarios is applying the following form of Bernoulli s equation: p V z constant! g + g + = from point 1 to point valid? a. 1 stagnant column of water steady, inviscid,

More information

To study the motion of a perfect gas, the conservation equations of continuity

To 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 information

Stream Tube. When density do not depend explicitly on time then from continuity equation, we have V 2 V 1. δa 2. δa 1 PH6L24 1

Stream Tube. When density do not depend explicitly on time then from continuity equation, we have V 2 V 1. δa 2. δa 1 PH6L24 1 Stream Tube A region of the moving fluid bounded on the all sides by streamlines is called a tube of flow or stream tube. As streamline does not intersect each other, no fluid enters or leaves across the

More information

Fundamentals of Fluid Dynamics: Waves in Fluids

Fundamentals of Fluid Dynamics: Waves in Fluids Fundamentals of Fluid Dynamics: Waves in Fluids Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/ tzielins/ Institute

More information

Final: Solutions Math 118A, Fall 2013

Final: Solutions Math 118A, Fall 2013 Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or

More information

Introduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)

Introduction 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 information

1. (20 pts total 2pts each) - Circle the most correct answer for the following questions.

1. (20 pts total 2pts each) - Circle the most correct answer for the following questions. ME 50 Gas Dynamics Spring 009 Final Exam NME:. (0 pts total pts each) - Circle the most correct answer for the following questions. i. normal shock propagated into still air travels with a speed (a) equal

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (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 information

In this lecture... Centrifugal compressors Thermodynamics of centrifugal compressors Components of a centrifugal compressor

In this lecture... Centrifugal compressors Thermodynamics of centrifugal compressors Components of a centrifugal compressor Lect- 3 In this lecture... Centrifugal compressors Thermodynamics of centrifugal compressors Components of a centrifugal compressor Centrifugal compressors Centrifugal compressors were used in the first

More information

Sound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur

Sound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-17 HARMONIC PLANE WAVES Introduction In this lecture, we discuss propagation of 1-D planar waves as they travel

More information

20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma:

20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma: Phys780: Plasma Physics Lecture 20. Alfven Waves. 1 20. Alfven waves ([3], p.233-239; [1], p.202-237; Chen, Sec.4.18, p.136-144) We have considered two types of waves in plasma: 1. electrostatic Langmuir

More information

31545 Medical Imaging systems

31545 Medical Imaging systems 31545 Medical Imaging systems Lecture 5: Blood flow in the human body Jørgen Arendt Jensen Department of Electrical Engineering (DTU Elektro) Biomedical Engineering Group Technical University of Denmark

More information

OPAC102. The Acoustic Wave Equation

OPAC102. The Acoustic Wave Equation OPAC102 The Acoustic Wave Equation Acoustic waves in fluid Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluid. The restoring forces responsible for propagating

More information

THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS

THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS Journal of Hyperbolic Differential Equations Vol., No. 4 005 909 917 c World Scientific Publishing Company THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS VOLKER ELLING, and TAI-PING LIU, Dept.

More information

Lecture 3: 1. Lecture 3.

Lecture 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 information

Introduction to Aerospace Engineering

Introduction 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 information

Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche

Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second

More information

AA210A Fundamentals of Compressible Flow. Chapter 13 - Unsteady Waves in Compressible Flow

AA210A Fundamentals of Compressible Flow. Chapter 13 - Unsteady Waves in Compressible Flow AA210A Fundamentals of Compressible Flow Chapter 13 - Unsteady Waves in Compressible Flow The Shock Tube - Wave Diagram 13.1 Equations for irrotational, homentropic, unsteady flow ρ t + x k ρ U i t (

More information

IX. COMPRESSIBLE FLOW. ρ = P

IX. COMPRESSIBLE FLOW. ρ = P IX. COMPRESSIBLE FLOW Compressible flow is the study of fluids flowing at speeds comparable to the local speed of sound. This occurs when fluid speeds are about 30% or more of the local acoustic velocity.

More information

On Fluid Maxwell Equations

On 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 information

Turbomachinery. Hasan Ozcan Assistant Professor. Mechanical Engineering Department Faculty of Engineering Karabuk University

Turbomachinery. Hasan Ozcan Assistant Professor. Mechanical Engineering Department Faculty of Engineering Karabuk University Turbomachinery Hasan Ozcan Assistant Professor Mechanical Engineering Department Faculty of Engineering Karabuk University Introduction Hasan Ozcan, Ph.D, (Assistant Professor) B.Sc :Erciyes University,

More information

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow

AA210A 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 information

Non-relativistic flows of perfect fluids

Non-relativistic flows of perfect fluids CHAPTER IV Non-relativistic flows of perfect fluids In the previous Chapter, we have introduced the coupled dynamical equations that govern the flows of perfect fluids in the non-relativistic regime, namely

More information

1. For an ideal gas, internal energy is considered to be a function of only. YOUR ANSWER: Temperature

1. For an ideal gas, internal energy is considered to be a function of only. YOUR ANSWER: Temperature CHAPTER 11 1. For an ideal gas, internal energy is considered to be a function of only. YOUR ANSWER: Temperature 2.In Equation 11.7 the subscript p on the partial derivative refers to differentiation at

More information

Multi-D MHD and B = 0

Multi-D MHD and B = 0 CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B

More information

Shock and Expansion Waves

Shock and Expansion Waves Chapter For the solution of the Euler equations to represent adequately a given large-reynolds-number flow, we need to consider in general the existence of discontinuity surfaces, across which the fluid

More information

BERNOULLI EQUATION. The motion of a fluid is usually extremely complex.

BERNOULLI 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 information

Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015

Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous

More information

Flow rate and mass flow rate

Flow rate and mass flow rate EEN-E1040 Measurement and control of energy systems Flow measurements / 14 Sep 2017 WELCOME! v. 01 / T. Paloposki Flow rate and mass flow rate Consider the system shown here 1 Volume flow rate through

More information

Fundamentals of Acoustics

Fundamentals 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 information

Thin airfoil theory. Chapter Compressible potential flow The full potential equation

Thin airfoil theory. Chapter Compressible potential flow The full potential equation hapter 4 Thin airfoil theory 4. ompressible potential flow 4.. The full potential equation In compressible flow, both the lift and drag of a thin airfoil can be determined to a reasonable level of accuracy

More information

Partial Differential Equations II

Partial Differential Equations II Partial Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Partial Differential Equations II 1 / 28 Almost Done! Homework

More information

2 Navier-Stokes Equations

2 Navier-Stokes Equations 1 Integral analysis 1. Water enters a pipe bend horizontally with a uniform velocity, u 1 = 5 m/s. The pipe is bended at 90 so that the water leaves it vertically downwards. The input diameter d 1 = 0.1

More information

Waves in plasma. Denis Gialis

Waves in plasma. Denis Gialis Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.

More information

EGN 3353C Fluid Mechanics

EGN 3353C Fluid Mechanics Lecture 8 Bernoulli s Equation: Limitations and Applications Last time, we derived the steady form of Bernoulli s Equation along a streamline p + ρv + ρgz = P t static hydrostatic total pressure q = dynamic

More information

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

Integrodifferential 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 information

Gas Dynamics: Basic Equations, Waves and Shocks

Gas Dynamics: Basic Equations, Waves and Shocks Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks Susanne Höfner Susanne.Hoefner@fysast.uu.se Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks

More information

Compressible Flow. Professor Ugur GUVEN Aerospace Engineer Spacecraft Propulsion Specialist

Compressible Flow. Professor Ugur GUVEN Aerospace Engineer Spacecraft Propulsion Specialist Compressible Flow Professor Ugur GUVEN Aerospace Engineer Spacecraft Propulsion Specialist What is Compressible Flow? Compressible Flow is a type of flow in which the density can not be treated as constant.

More information

4 Compressible Fluid Dynamics

4 Compressible Fluid Dynamics 4 Compressible Fluid Dynamics 4. Compressible flow definitions Compressible flow describes the behaviour of fluids that experience significant variations in density under the application of external pressures.

More information

Isentropic Flow. Gas Dynamics

Isentropic Flow. Gas Dynamics Isentropic Flow Agenda Introduction Derivation Stagnation properties IF in a converging and converging-diverging nozzle Application Introduction Consider a gas in horizontal sealed cylinder with a piston

More information

ρ x + fv f 'w + F x ρ y fu + F y Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ Dv Dt + u2 tanφ + vw a a = 1 p Dw Dt u2 + v 2

ρ x + fv f 'w + F x ρ y fu + F y Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ Dv Dt + u2 tanφ + vw a a = 1 p Dw Dt u2 + v 2 Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ + uw Dt a a = 1 p ρ x + fv f 'w + F x Dv Dt + u2 tanφ + vw a a = 1 p ρ y fu + F y Dw Dt u2 + v 2 = 1 p a ρ z g + f 'u + F z Dρ Dt + ρ

More information

SOLAR MHD Lecture 2 Plan

SOLAR 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 information

Quick Recapitulation of Fluid Mechanics

Quick 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 information

Notes #6 MAE 533, Fluid Mechanics

Notes #6 MAE 533, Fluid Mechanics Notes #6 MAE 533, Fluid Mechanics S. H. Lam lam@princeton.edu http://www.princeton.edu/ lam October 1, 1998 1 Different Ways of Representing T The speed of sound, a, is formally defined as ( p/ ρ) s. It

More information

Physics 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 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 information

In this section, mathematical description of the motion of fluid elements moving in a flow field is

In 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 information

Partitioned Methods for Multifield Problems

Partitioned Methods for Multifield Problems C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible

More information

P 1 P * 1 T P * 1 T 1 T * 1. s 1 P 1

P 1 P * 1 T P * 1 T 1 T * 1. s 1 P 1 ME 131B Fluid Mechanics Solutions to Week Three Problem Session: Isentropic Flow II (1/26/98) 1. From an energy view point, (a) a nozzle is a device that converts static enthalpy into kinetic energy. (b)

More information

Section 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 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 information

PDE Solvers for Fluid Flow

PDE 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 information

SPC 407 Sheet 5 - Solution Compressible Flow Rayleigh Flow

SPC 407 Sheet 5 - Solution Compressible Flow Rayleigh Flow SPC 407 Sheet 5 - Solution Compressible Flow Rayleigh Flow 1. Consider subsonic Rayleigh flow of air with a Mach number of 0.92. Heat is now transferred to the fluid and the Mach number increases to 0.95.

More information

Fluid equations, magnetohydrodynamics

Fluid equations, magnetohydrodynamics Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics

More information

Various lecture notes for

Various 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 information

Introduction to Fluid Machines, and Compressible Flow Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur

Introduction to Fluid Machines, and Compressible Flow Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Introduction to Fluid Machines, and Compressible Flow Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 09 Introduction to Reaction Type of Hydraulic

More information

Chapter 1. Introduction to Nonlinear Space Plasma Physics

Chapter 1. Introduction to Nonlinear Space Plasma Physics Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large

More information

Fluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14

Fluid 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 information

Recapitulation: Questions on Chaps. 1 and 2 #A

Recapitulation: Questions on Chaps. 1 and 2 #A Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?

More information

III. PLASMA DYNAMICS

III. PLASMA DYNAMICS III. PLASMA DYNAMICS Prof. S. C. Brown Prof. H. H. Woodson D. T. Kyrazis Prof. O. K. Mawardi Dr. I. C. T. Nisbet A. T. Lewis Prof. D. O. Akhurst D. A. East J. K. Oddson Prof. D. C. Pridmore-Brown Z. J.

More information

Objectives. Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation

Objectives. Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation Objectives Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation Conservation of Mass Conservation of Mass Mass, like energy, is a conserved

More information

Several forms of the equations of motion

Several 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 information

Sound Waves Sound Waves:

Sound Waves Sound Waves: 3//18 Sound Waves Sound Waves: 1 3//18 Sound Waves Linear Waves compression rarefaction Inference of Sound Wave equation: Sound Waves We look at small disturbances in a compressible medium (note: compressible

More information

0.2. CONSERVATION LAW FOR FLUID 9

0.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 information

Characteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution:

Characteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution: AMath 574 January 3, 20 Today: Boundary conditions Multi-dimensional Wednesday and Friday: More multi-dimensional Reading: Chapters 8, 9, 20 R.J. LeVeque, University of Washington AMath 574, January 3,

More information

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft 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 information

Advection, Conservation, Conserved Physical Quantities, Wave Equations

Advection, Conservation, Conserved Physical Quantities, Wave Equations EP711 Supplementary Material Thursday, September 4, 2014 Advection, Conservation, Conserved Physical Quantities, Wave Equations Jonathan B. Snively!Embry-Riddle Aeronautical University Contents EP711 Supplementary

More information

Introduction to Gas Dynamics All Lecture Slides

Introduction 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 information

1/18/2011. Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis ESS227 Prof. Jin-Yi Yu

1/18/2011. Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis ESS227 Prof. Jin-Yi Yu Lecture 2: Basic Conservation Laws Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed in an inertial system = rate of change of Ua following the motion in an inertial

More information

Magnetohydrodynamic Waves

Magnetohydrodynamic Waves Magnetohydrodynamic Waves Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 17, 2016 These slides are largely based off of 4.5 and 4.8 of The Physics of

More information

Flow-Induced Vibration of Pipeline on Elastic Support

Flow-Induced Vibration of Pipeline on Elastic Support Available online at www.sciencedirect.com Procedia Engineering () 986 99 The Twelfth East Asia-Pacific Conference on Structural Engineering and Construction Flow-Induced Vibration of Pipeline on Elastic

More information

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

AA214B: 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 information

Introduction to Physical Acoustics

Introduction 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 information

THEORETICAL AND EXPERIMENTAL INVESTIGATIONS ON CHOKING PHENOMENA OF AXISYMMETRIC CONVERGENT NOZZLE FLOW

THEORETICAL AND EXPERIMENTAL INVESTIGATIONS ON CHOKING PHENOMENA OF AXISYMMETRIC CONVERGENT NOZZLE FLOW 8 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES THEORETICAL AND EXPERIMENTAL INVESTIGATIONS ON CHOKING PHENOMENA OF AXISYMMETRIC CONVERGENT NOZZLE FLOW Ryuta ISOZUMI*, Kazunori KUBO*, Daisuke

More information

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.

You 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 information

GAS DYNAMICS. M. Halük Aksel. O. Cahit Eralp. and. Middle East Technical University Ankara, Turkey

GAS DYNAMICS. M. Halük Aksel. O. Cahit Eralp. and. Middle East Technical University Ankara, Turkey GAS DYNAMICS M. Halük Aksel and O. Cahit Eralp Middle East Technical University Ankara, Turkey PRENTICE HALL f r \ New York London Toronto Sydney Tokyo Singapore; \ Contents Preface xi Nomenclature xiii

More information

Space Physics. ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen. Aalto University School of Electrical Engineering

Space Physics. ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen. Aalto University School of Electrical Engineering Space Physics ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen Aalto University School of Electrical Engineering The 6 th week: topics Last week: Examples of waves MHD: Examples

More information

Basic Fluid Mechanics

Basic Fluid Mechanics Basic Fluid Mechanics Chapter 3B: Conservation of Mass C3B: Conservation of Mass 1 3.2 Governing Equations There are two basic types of governing equations that we will encounter in this course Differential

More information

Errata for Conquering the Physics GRE, edition 2

Errata for Conquering the Physics GRE, edition 2 Errata for Conquering the Physics GRE, edition 2 March 12, 217 This document contains corrections to Conquering the Physics GRE, edition 2. It includes errata that are present in both all printings of

More information

Introduction and Basic Concepts

Introduction and Basic Concepts Topic 1 Introduction and Basic Concepts 1 Flow Past a Circular Cylinder Re = 10,000 and Mach approximately zero Mach = 0.45 Mach = 0.64 Pictures are from An Album of Fluid Motion by Van Dyke Flow Past

More information

Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion

Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Department of applied physics, Aalto University March 8, 2016 Hamiltonian versus Newtonian mechanics Newtonian mechanics:

More information

2 Equations of Motion

2 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 information

Gasdynamics 1-D compressible, inviscid, stationary, adiabatic flows

Gasdynamics 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 information

LEAST-SQUARES FINITE ELEMENT MODELS

LEAST-SQUARES FINITE ELEMENT MODELS LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows

More information

Solar Physics & Space Plasma Research Center (SP 2 RC) MHD Waves

Solar Physics & Space Plasma Research Center (SP 2 RC) MHD Waves MHD Waves Robertus vfs Robertus@sheffield.ac.uk SP RC, School of Mathematics & Statistics, The (UK) What are MHD waves? How do we communicate in MHD? MHD is kind! MHD waves are propagating perturbations

More information

MHD Modes of Solar Plasma Structures

MHD Modes of Solar Plasma Structures PX420 Solar MHD 2013-2014 MHD Modes of Solar Plasma Structures Centre for Fusion, Space & Astrophysics Wave and oscillatory processes in the solar corona: Possible relevance to coronal heating and solar

More information

Fluid Mechanics Theory I

Fluid Mechanics Theory I Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to

More information

ENGINEERING FLUID MECHANICS. CHAPTER 1 Properties of Fluids

ENGINEERING FLUID MECHANICS. CHAPTER 1 Properties of Fluids CHAPTER 1 Properties of Fluids ENGINEERING FLUID MECHANICS 1.1 Introduction 1.2 Development of Fluid Mechanics 1.3 Units of Measurement (SI units) 1.4 Mass, Density, Specific Weight, Specific Volume, Specific

More information

Finite difference method for solving Advection-Diffusion Problem in 1D

Finite difference method for solving Advection-Diffusion Problem in 1D Finite difference method for solving Advection-Diffusion Problem in 1D Author : Osei K. Tweneboah MATH 5370: Final Project Outline 1 Advection-Diffusion Problem Stationary Advection-Diffusion Problem in

More information

Physically Based Simulations (on the GPU)

Physically Based Simulations (on the GPU) Physically Based Simulations (on the GPU) (some material from slides of Mark Harris) CS535 Fall 2014 Daniel G. Aliaga Department of Computer Science Purdue University Simulating the world Floating point

More information

Why bouncing droplets are a pretty good model for quantum mechanics

Why bouncing droplets are a pretty good model for quantum mechanics Why bouncing droplets are a pretty good model for quantum mechanics Robert Brady and Ross Anderson University of Cambridge robert.brady@cl.cam.ac.uk ross.anderson@cl.cam.ac.uk Cambridge, June 2014 Robert

More information