EP711 Supplementary Material Thursday, September 4, 2014 Advection, Conservation, Conserved Physical Quantities, Wave Equations Jonathan B. Snively!Embry-Riddle Aeronautical University
Contents EP711 Supplementary Material Thursday, September 4, 2014 Lagrangian vs. Eulerian Coordinates Advection and Conservation Conserved Physical Quantities Wave Equations
Lagrangian Description of Fluid Motion Describes temporal evolution of a quantity Q in a parcel of gas or fluid. Assume Q(x(t),t), where x(t) may be a vector position and t is time, tracking the evolution of a single parcel s properties. Effectively, Q is a function of time only. Conservation of Q within the parcel is expressed by: dq dt =0 Which, simply, states that some quantity Q in a parcel does not vary in time the total derivative of Q is zero.
Total Derivative Since Q is a function of x(t) and t, the Lagrangian total derivative, may be expanded with help from the chain rule, yielding the following expression: dq dt = @Q @t + @~x @t rq Considering that x/ t=v (velocity), dq dt = @Q + ~v rq @t Thus, the material derivative can be expanded as: d dt = @ @t + ~v r
Total Derivative In case you are having doubts about the total derivative and the chain rule, let s write out the full expression, where x=(x(t), y(t), z(t)): dq(x, y, z, t) dt = @Q @t + @Q @x @x @t + @Q @y @y @t + @Q @z @z @t Or, Where dq dt = @Q @t + @~x @t rq rq = @Q @x î + @Q @y ĵ + @Q @z ˆk
Eulerian Description of Fluid Motion Describes temporal evolution of quantity Q in a gas or fluid at fixed points in space. Conservation of Q in space and time is expressed:!! @Q @t + ~v rq =0 This states that a quantity Q varies in space and time, while it is advected (transported) by a fluid velocity field v. This is also known as an advection equation, which here is not expressed in its most-general form.
Contents EP711 Supplementary Material Thursday, September 4, 2014 Lagrangian vs. Eulerian Coordinates Advection and Conservation Conserved Physical Quantities Wave Equations
Advection Equation for Non-Divergent Flow A more general conservative advection equation, describing conservation of a quantity in a divergent flow, is given by: where @Q @t + r (Q~v) =0 r (Q~v) =~v rq + Qr ~v In contrast, the previous expressions, given by: dq dt = @Q + ~v rq =0 @t have assumed non-divergent velocity: r ~v =0
Conservation Laws and Continuity Equations This more general advection equation can be considered a continuity equation or conservation law for a quantity Q that is conserved under all conditions for example, mass density, momentum, or energy. @Q @t + r (Q~v) =0 Here, the quantity in parentheses is known as the flux f, such that the expression can be written as: @Q @t = r ~f Although Q is here a scalar, it does not need to be!
Conservation Laws and Continuity Equations fx Q dz fx + ( fx/ x) dx dx dy @Q @t = r ~f
Flux Form (for Conservative Advection) For a scalar quantity (such as density), the flux f depends on Q(x,t), such that the expression can be written as: @Q @t = r ~f Where Q(x,t) is a vector, this expression effectively contains one equation for each spatial dimension; f is then a tensor. An example we will discuss is the momentum conservation equation, here shown using somewhat mixed notation: @ @t ( ~v) = r ( v iv j + p ij )
Just a note... The advection equation is a form of a one-way wave equation, with solutions traveling at velocity v. Combine two of them to obtain the familiar second-order wave equation!
Contents EP711 Supplementary Material Thursday, September 4, 2014 Lagrangian vs. Eulerian Coordinates Advection and Conservation Conserved Physical Quantities Wave Equations
Conservation of Mass For a scalar quantity (such as density), the flux f depends on Q(x,t), such that the expression can be written as: @Q @t = r ~f The best-known conservation law is the continuity equation for neutral mass density: @ @t = r ( ~v) In 1-Dimension: @ @t = @ @x ( v) = v @ @x @v @x
Conservation of Charge (Same!) A vector current density is indicative of a flux of scalar charge density at a certain drift velocity: The resulting conservation law is the continuity equation for electric charge density: @ @t = r ( ~v) Integral form:
Conservation of Momentum Momentum equations for each spatial direction (x, y, z for Cartesian coordinates) here contained in one expression as a vector conservation law: @ @t ( ~v) = r ( v iv j + p ij ) This expression effectively describes Newton s law. In one dimension, momentum conservation can be expressed: @ @t ( v) = @ @x ( v2 + p)
Conservation of Energy Energy is a scalar quantity (like density): @E @t = r {(E + p)~v} Here, it is important to define how E relates to other quantities via an equation of state We will assume an ideal single-constituent gas. E = + 1 2 (~v ~v) Internal + Kinetic The internal energy is given by: = c v T = p ( 1)
Euler Equations of Motion: Example equations for a simple atmosphere supporting acoustic and gravity waves, but (lazily) excluding viscosity and thermal conduction (e.g., Snively s 2003 MS thesis): We model the atmosphere as a non-rotating, fully-nonlinear, compressible gas: Mass Density: Momentum: Energy: + ( v) t = 0 (1) ( v)+ ( vv) t = p g (2) E t + { (E + p)v} = gv z (3) where the energy equation and the equation of state for an ideal gas are defined as: p State: E = ( 1) + 1 (v v) (4) 2 where is density, is pressure, is the fluid velocity, along with energy density
Contents EP711 Supplementary Material Thursday, September 4, 2014 Lagrangian vs. Eulerian Coordinates Advection and Conservation Conserved Physical Quantities Wave Equations
Wave Equations The hyperbolic PDEs that we will investigate in this course yield wave equation solutions when manipulated in the appropriate context: @ 2 u @t 2 c 2 r 2 u =0 @ 2 u @t 2 c 2 @2 u @x 2 =0 The wave equation can be described instead in terms of two advection equations with opposite speeds: @u @t c @u @x =0 @u @t + c@u @x =0 A system involving two advection equations and two variables can be constructed...
Idealized System of Two Advection Equations A system of two coupled PDEs describing advection of quantities u and v can be specified: @u @t + c @v @x =0 @v @t + c@u @x =0 Combined, they form a wave equations with characteristic speed c. This form is consistent with that describing onedimensional EM wave propagation, acoustic wave propagation, or wave propagation on a string.
Recall from Undergrad Maxwell s Equations: Note: Potter demonstrates calculation of Lax method stability for a linear system of PDEs based on Maxwell s equations! The same can be applied for acoustics (hint-hint)...
Linear Acoustics Equations The equations of linear acoustics describe conservation of mass density (continuity) and conservation of linear momentum, which in the following system have been coupled via an ideal gas equation of state: @p @t + p @u 0 @x =0 @u @t + 1 @p 0 @x =0 Here, p is pressure, u is velocity, rho0 is density, and gamma is the ratio of specific heats.
Wave Equation Solution Like general wave equation solutions, the wave equation can be obtained as follows: @ 2 p @t 2 + p 0 @ @t @u @x =0 Substitute from momentum equation: @ 2 p @t 2 p 0 0 @ 2 p @x 2 =0 @u @t = The speed of sound can be expressed: c 2 s = p 0 @ 2 p @t 2 = c2 s @ 2 p @x 2 0 1 0 @p @x
Dispersion Relations Wave dispersion relations relate the effective wavenumber and frequency of waves within a system under the assumption of Fourier time-harmonic solutions of the form: u(x, t) =U o e j(!t kx) @ 2 u @t 2 c 2 @2 u @x 2 =0! 2 + c 2 k 2 =0 The resulting expression can be rewritten as:! = ck Which is an ideal dispersion (non-dispersive!) relation for waves with phase velocity and group velocity equal to c. NOTE: c is not the fluid velocity, but a characteristic speed of the physical system