Nonlinear Wave Propagation in 1D Random Media

Size: px
Start display at page:

Download "Nonlinear Wave Propagation in 1D Random Media"

Transcription

1 Nonlinear Wave Propagation in 1D Random Media J. B. Thoo, Yuba College UCI Computational and Applied Math Seminar Winter 22 Typeset by FoilTEX

2 Background O Doherty and Anstey (Geophysical Prospecting, 1971) looked at the effects of multiple scattering by a pulse propagating through the earth s crust. With the advent of digitally-generated seismograms, they found that including multiple reflections increased the amplitude of the primary pulse enough to offset transmission losses. Idea: Multiple scattering constructive interference. incident wave directly transmitted wave sum of double-scattering ( constructive interference ) UCI Comput. Appl. Math. Seminar, Wtr 2 1

3 Many contributions since. Burridge, Papanicolaou, White (1988); Burridge and Chang (1989); Clouet and Fouque (1994); Lewicki, Burridge, and Papanicolaou (1994); Sølna (1998); and others. Some to make the O Doherty-Anstey theory rigorous and some to extend the theory to cases where fluctuations in the medium may not be weak. UCI Comput. Appl. Math. Seminar, Wtr 2 2

4 Overview We study the propagation of a weakly nonlinear wave pulse, governed by a hyperbolic system of conservation laws, through a stationary random medium in 1D. Examples Gas dynamics (interaction of sound waves with background entropy wave). Nonlinear elastic waves in stratified media. MHD (e.g., scattering of magnetoacoustic waves off Alfvén waves). UCI Comput. Appl. Math. Seminar, Wtr 2 3

5 Scalings Incident pulse (primary wave): amplitude O(ε 2 ), width O(ε 2 ) Fluctuations (background medium or scattering wave): amplitude O(ε), width O(ε 2 ) Propagation: distance O(1/ε 2 ), time O(1) UCI Comput. Appl. Math. Seminar, Wtr 2 4

6 Expansion We use the asymptotic expansion introduced by Benilov and Pelinovskĭi (1988) to get a nonlinear O Doherty-Anstey theory. The asymptotic solution includes a realization-dependent random phase shift. W.r.t. this reference frame, the pulse moves as a Brownian motion and the leading order solution is nonrandom. This self-averaging of the leading-order solution means that closure problems, u 2 u 2, do not arise in the nonlinear theory. UCI Comput. Appl. Math. Seminar, Wtr 2 5

7 In the simplest case, the normalized form of asymptotic equation for the wave profile u is u t + ( 1 2 u2) x = k(y)u xx (x y, t)dy, where k(y) = β( )β( + y) is the covariance function of scattering coefficients β (random fluctuations from the double scattering of waves). The wave profile satisfies an inviscid Burgers equation with a nonlocal, lower-order dissipative and dispersive term that describes the effects of double scattering of waves off the fluctuations in the medium or scattering wave. The fluctuations are assumed to be uniformly bounded in space. UCI Comput. Appl. Math. Seminar, Wtr 2 6

8 Outline Sketch of the derivation of the asymptotic equations Scattering by a nonuniform medium Scattering of one wave off another wave Examples: 1D gas dynamics, nonlinear elasticity Properties of the asymptotic equations (shock profiles, shock formation) UCI Comput. Appl. Math. Seminar, Wtr 2 7

9 Derivation of asymptotic equations Scattering by a nonuniform medium Quadratic nonlinearity Consider a strictly hyperbolic system of conservation laws of the form u t + [ f(u) + εb(x/ε 2 )u ] x = u = u(x,t) R m defined for < x < and t f : R m R m B = B(x/ε 2 ) R m m is a realization of a stationary random process with zero spatial mean (assumed to be uniformly bounded) ε is a small, positive parameter UCI Comput. Appl. Math. Seminar, Wtr 2 8

10 Seek an asymptotic solution u ε (ε ) of the form u ε = ε 2 u 2 (ξ,t) + ε 3 u 3 (ξ, y, t) + ε 4 u 4 (ξ, y, t) + O(ε 5 ) y = x ε 2, ξ = x λt ε3 φ(x/ε 2 ) ε 2 φ is a realization of a random process and is to be chosen u 2 is the primary wave UCI Comput. Appl. Math. Seminar, Wtr 2 9

11 Find that Primary wave is u 2 = a(ξ,t)r Random correction to the pulse location is where ( x ) εφ ε 2 = 1 x/ε 2 λ ε β(y) dy, β = l T B(y)r describes the strength of an i-wave generated by the scattering of a k-wave off the fluctuations of the medium. Remark: φ typically converges in distribution to a Brownian motion. E.g. [Lipster and Shiryaev], if β is strictly stationary with a suitable mixing condition. UCI Comput. Appl. Math. Seminar, Wtr 2 1

12 The asymptotic equation for the wave profile of the primary wave is a t + ( 1 2 Ga2) = ± Λ W(σ)a ξξ (cσ + ξ,t) dσ, ξ where W(σ) = β(y)β(y σ) and the average is given by Remarks: f(, y, ) = lim L 1 L L f(, y, )dy. G is a coefficient of quadratic nonlinearity that describes the strength of the j-wave generated by quadratic self-interaction. W is the covariance of scattering coefficients. UCI Comput. Appl. Math. Seminar, Wtr 2 11

13 Wave-wave scattering Conservative form Consider a strictly hyperbolic system of conservation laws u t + f(u) x = u = u(x,t) R m defined for < x < and t f : R m R m UCI Comput. Appl. Math. Seminar, Wtr 2 12

14 Seek asymptotic solution u ε (ε ) of the form u ε = εu 1 (ξ,t) + ε 2 u 2 (ξ,η, t) + ε 3 u 3 (ξ, η, t) + ε 4 u 4 (ξ, η, t) + O(ε 5 ) ξ = x λt ε3 φ(η, t) ε 2, η = x µt ε 2 φ is a realization of a random process to be chosen u 1 is the background scattering wave and u 2 is the primary wave UCI Comput. Appl. Math. Seminar, Wtr 2 13

15 A similar derivation as before leads us to find that the random correction to the pulse location is εφ(η, t) = 1 η λ µ Γ s(y, t) dy, where Γ describes the strength of an i-wave produced by the quadratically nonlinear interaction of a k-wave with a p-wave. UCI Comput. Appl. Math. Seminar, Wtr 2 14

16 The asymptotic equations for the wave profile of the background scattering wave s = s(η, t) and for the wave profile of the primary wave a = a(ξ, t) are s t + ( 1 3 Hs3) η =, (a + Z s 2 ) t + ( 1 2 Ga2) ξ = ± Λ W(σ, t)a ξξ (ξ cσ, t) dσ, where W(σ, t) = s(η, t)s(η σ, t). UCI Comput. Appl. Math. Seminar, Wtr 2 15

17 Remarks: s is assumed to be a stationary random process with zero spatial mean. s is not genuinely nonlinear at u =. G is a coefficient of quadratic nonlinearity that describes the strength of the j-wave generated by quadratic self-interaction. Λ, the coefficients of the nonlocal terms, describe the strength of the j-wave generated by double scattering: j p k, k p j. W is a covariance of the background scattering wave. We get the same asymptotic equations with slightly modified coefficients for the system of hyperbolic conservation laws in nonconservative form u t + A(u)u x =. UCI Comput. Appl. Math. Seminar, Wtr 2 16

18 Example: 1D gas dynamics Scattering of a nonlinear sound wave off entropy fluctuations 1D compressible Euler equations ( ρe ρv2) t + p = κρ γ exp(s/c v ), c 2 = γp ρ ρ t + (ρv) x =, (ρv) t + (ρv 2 + p) x =, ( ρev + 1 ) 2 ρv3 + pv = uniform background state with density ρ, velocity v =, entropy S, and sound speed c = γp /ρ x UCI Comput. Appl. Math. Seminar, Wtr 2 17

19 This leads to the asymptotic equations (in normalized form) for the wave profiles of the entropy wave (background scattering wave) s t = right-moving sound wave (primary wave) a t + ( 1 2 a2) x = W ( σ 2) a xx (x σ, t) dσ, where W(σ) = s(y)s(y σ). Remark: The entropy wave is not genuinely nonlinear at u = ; in fact, it is linearly degenerate. UCI Comput. Appl. Math. Seminar, Wtr 2 18

20 So, the asymptotic solution of the compressible Euler equations is where ρ v S ρ S + ε 2 a(ξ, t) + εs(η) ρ c p /c v c 2 + O(ε 3 ), η = x ε 2, ξ = x c t ε 3 φ(x/ε 2 ) ε 2 φ(x/ε 2 ) = c2 2γc v x/ε 2 s(y) dy. UCI Comput. Appl. Math. Seminar, Wtr 2 19

21 Example: Nonlinear elasticity Elastic wave propagating in 1D through an isotropic hyperelastic medium Deformation x 1 = x + u 1 (x, t), x 2 = x 2 + u 2 (x,t) Equations of motion ρu 1tt = (g(u 1x, u 2 2x)) x, ρu 2tt = (u 2x h(u 1x, u 2 2x)) x density ρ: R R strain-energy density σ: R R R constitutive functions g, h: R R R g(p, q) = σ (p, q), p h(p, q) = 2 σ(p, q) q UCI Comput. Appl. Math. Seminar, Wtr 2 2

22 Scattering of longitudinal waves in a nonuniform medium Assume the Lamé moduli and density have the expansions (ε ) λ(x) = λ + ελ 1 (x/ε 2 ) + O(ε 2 ), µ(x) = µ + εµ 1 (x/ε 2 ) + O(ε 2 ), ρ(x) = ρ + ερ 1 (x/ε 2 ) + O(ε 2 ). Then we can write the equations of motion in the form u t + [ f(u) + εb(x/ε 2 )u ] x =, where B = (2µ 1 + λ 1 ) µ 1 ρ 1 /ρ 2 ρ 1 /ρ 2. UCI Comput. Appl. Math. Seminar, Wtr 2 21

23 This leads to the asymptotic equation (in normalized form) for the wave profile of the right-moving longitudinal wave a t + ( 1 2 a2 ) x = W ( σ 2) a xx (x σ, t) dσ, where W(σ) = β + (y)β (y σ), β ± (y) = 1 2ρ c l [ 2µ1 (y) + λ 1 (y) ± c 2 l ρ 1 (y) ] and c l is the longitudinal wave speed. Note: The nonlocal term on the RHS originates from the scattering of a right-moving longitudinal wave into a left-moving longitudinal wave. The scattering of a longitudinal wave does not produce any transverse waves, and vice versa, in an isotropic medium. UCI Comput. Appl. Math. Seminar, Wtr 2 22

24 So, the asymptotic solution of the equations of motion is where m 1 m 2 w 1 w 2 ε2 a(ξ, t) φ(x/ε 2 ) = 1 2ρ c l c t 1 1/(ρ c l ) m k = ρu kt, w k = u kx, ξ = x c tt ε 3 φ(x/ε 2 ) ε 2, x/ε 2 and c t is the transverse wave speed. + O(ε3 ), [ 2µ1 (y) + λ 1 (y) c 2 l ρ 1 (y) ] dy UCI Comput. Appl. Math. Seminar, Wtr 2 23

25 Scattering of transverse waves in a nonuniform medium The asymptotic equation (in normalized form) for the wave profile of the right-moving transverse wave is a t ( 1 3 a3) x = V ( σ 2) a xx (x σ, t) dσ, where V (σ) = β + (y)β (y σ), β ± (y) = 1 2ρ c t [ µ1 (y) ± c 2 tρ 1 (y) ]. UCI Comput. Appl. Math. Seminar, Wtr 2 24

26 The asymptotic solution of the equations of motion is where m 1 m 2 w 1 w 2 εa(ξ, t) φ(x/ε 2 ) = 1 2ρ c l c t 1 1/(ρ c t ) m k = ρu kt, w k = u kx, ξ = x c lt ε 3 φ(x/ε 2 ) ε 2, x/ε 2 + O(ε2 ), [ µ1 (y) c 2 tρ 1 (y) ] dy. UCI Comput. Appl. Math. Seminar, Wtr 2 25

27 Scattering of longitudinal wave off transverse wave Density, ρ = ρ, Lamé moduli λ = λ and µ = µ The asymptotic equations (in normalized form) for the wave profiles of the right-moving transverse wave (background scattering wave) s t + ( 1 3 Hs3) η =. right-moving longitudinal wave (primary wave) a t + ( 1 2 a2) ξ = W ( σ c, t ) a ξξ (ξ σ, t)dσ W(σ, t) = s(η, t)s(η σ, t). Remark: The transverse waves are not genuinely nonlinear at u =. UCI Comput. Appl. Math. Seminar, Wtr 2 26

28 The asymptotic solution of the equations of motion is m 1 m 2 w 1 w 2 εs(η, t) 1 1/(ρ c t ) + ε 2 a(ξ,t) 1 1/(ρ c l ) + O(ε3 ), where m k = ρu kt, w k = u kx, η = x c tt ε 2, ξ = x c lt ε 2. UCI Comput. Appl. Math. Seminar, Wtr 2 27

29 Properties of the asymptotic equations Shock profiles In the simplest case, the normalized form of asymptotic equation for the wave profile u is u t + ( 1 2 u2) x = k(y)u xx (x y, t) dy, where k is the covariance function of the random fluctuations. The right-hand side describes the effects of double scattering of the incident pulse wave. is a lower-order dissipative and dispersive term, and its dissipative effects may prevent the formation of shocks. UCI Comput. Appl. Math. Seminar, Wtr 2 28

30 Example 1 k(x) = exp( x ) random telegraph process Ornstein-Uhlenbeck process Applying the operator x + 1 to the normalized equation yields the relaxing gas equation, [ ( x + 1) u t + u x + ( 1 2 u2 ) x ] u x =. Traveling wave solution u = u(x ct, t) satisfies (c 1 u)u u2 cu =. lim x = u L, lim x u = c = u L /2 UCI Comput. Appl. Math. Seminar, Wtr 2 29

31 If u L > 2 (c > 1), then the solution is smooth with an embedded subshock with shock speed c = 1 + u. u = 2c L 2(c 1) c = 1 + <u> u = R If < u L < 2 ( < c < 1), then the solution has a smooth wave profile. UCI Comput. Appl. Math. Seminar, Wtr 2 3

32 Example 2 k(x) = { 1 x if x 1, if x > 1. This leads to the partial differential-difference equation u t + ( u u2) x = u(x + 1, t) u(x,t). The linearized equation has solution u t + u x = u(x + 1, t) u(x, t) u(x, t) = n= t n n! e t f(x t + n), u(x,) = f(x). So, an initial wave profile splits into multiple pulses as the wave propagates. UCI Comput. Appl. Math. Seminar, Wtr 2 31

33 1.5 q(1) at time t = 1.5 q(1) at time t = q(1) at time t = q(1) at time t = u t + ( 1 2 u2) x = u(x + 1, t) u(x, t), u(x,) = exp( 25x 2 ) UCI Comput. Appl. Math. Seminar, Wtr 2 32

34 Shock formation For appropriate covariance functions k, if a smooth pulse satisfying u t + ( 1 2 u2 ) x = k(y)u xx (x y, t) dy is nowhere too steep when t =, then it will not break; and if the pulse is anywhere too steep at t =, then it will break in finite time. Of importance is the condition that k ( ) >, for it provides a damping term, evident after integration by parts twice: u t + (k( ) + u)u x = k ( )u + k (y)u(x y, t) dy. The integral term alone on the RHS may not be enough to prevent wave breaking. UCI Comput. Appl. Math. Seminar, Wtr 2 33

35 Theorem. Assume that 1. k H 3 (R) is a covariance function and k ( ) >. 2. u = u(x, t) is a smooth solution of u t + ( 1 2 u2) x = k(y)u xx (x y, t)dy that decays rapidly as x ±. Let α = k ( ) and C = k ( ) u + ( 1 α k ( ) k 2 + k 2 ) u If C < α 2 /4 and u x > (α + α 2 4C)/2 for all x R, then the solution u does not break. 2. If u x < (α + α 2 + 4C)/2 for some x R, then the solution breaks in finite time. In fact, u x as t [ (α + α 2 + 4C)/2 u x ] 1. UCI Comput. Appl. Math. Seminar, Wtr 2 34

36 Sketch of proof energy estimate u 2 u 2 sup-norm estimate u u e αt + 1 α k 2 u 2 α, C R + and C dv dt + v2 + αv C v is globally smooth if α 2 4C > and v() > (α + α 2 4C)/2 v blows up in finite time if v() < (α + α 2 + 4C)/2 UCI Comput. Appl. Math. Seminar, Wtr 2 35

From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction

From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy

More information

Alexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011

Alexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011 Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,

More information

Travelling waves. Chapter 8. 1 Introduction

Travelling waves. Chapter 8. 1 Introduction Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part

More information

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

Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.

Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,

More information

Math Partial Differential Equations 1

Math Partial Differential Equations 1 Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)

More information

Introduction LECTURE 1

Introduction LECTURE 1 LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in

More information

Fission of a longitudinal strain solitary wave in a delaminated bar

Fission of a longitudinal strain solitary wave in a delaminated bar Fission of a longitudinal strain solitary wave in a delaminated bar Karima Khusnutdinova Department of Mathematical Sciences, Loughborough University, UK K.Khusnutdinova@lboro.ac.uk and G.V. Dreiden, A.M.

More information

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso

More information

Waves and characteristics: Overview 5-1

Waves and characteristics: Overview 5-1 Waves and characteristics: Overview 5-1 Chapter 5: Waves and characteristics Overview Physics and accounting: use example of sound waves to illustrate method of linearization and counting of variables

More information

Ray Theory for a Locally Layered Random Medium

Ray Theory for a Locally Layered Random Medium Ray Theory for a Locally Layered Random Medium K. Sølna G. Papanicolaou August 4, 24 Abstract We consider acoustic pulse propagation in inhomogeneous media over relatively long propagation distances. Our

More information

Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine

Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine Lecture 2 The wave equation Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.0 28/09/2018 1 Learning objectives of this lecture Understand the fundamental properties of the wave equation

More information

0.3.4 Burgers Equation and Nonlinear Wave

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 information

1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.

1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω. 1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the

More information

Simple 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 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

SPHERICAL SHOCK WAVES IN MAGNETO-GAS-DYNAMICS

SPHERICAL SHOCK WAVES IN MAGNETO-GAS-DYNAMICS CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 15 Number 1 Spring 27 SPHERICAL SHOCK WAVES IN MAGNETO-GAS-DYNAMICS RAJAN ARORA ABSTRACT. A Spherically symmetric model describing a blast wave in the solar

More information

ASYMPTOTIC THEORY FOR WEAKLY NON-LINEAR WAVE EQUATIONS IN SEMI-INFINITE DOMAINS

ASYMPTOTIC THEORY FOR WEAKLY NON-LINEAR WAVE EQUATIONS IN SEMI-INFINITE DOMAINS Electronic Journal of Differential Equations, Vol. 004(004), No. 07, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ASYMPTOTIC

More information

(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves

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

Propagation of Solitons Under Colored Noise

Propagation of Solitons Under Colored Noise Propagation of Solitons Under Colored Noise Dr. Russell Herman Departments of Mathematics & Statistics, Physics & Physical Oceanography UNC Wilmington, Wilmington, NC January 6, 2009 Outline of Talk 1

More information

Stochastic Stokes drift of a flexible dumbbell

Stochastic Stokes drift of a flexible dumbbell Stochastic Stokes drift of a flexible dumbbell Kalvis M. Jansons arxiv:math/0607797v4 [math.pr] 22 Mar 2007 Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK January

More information

Équation de Burgers avec particule ponctuelle

Équation de Burgers avec particule ponctuelle Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin

More information

Lecture Notes on Hyperbolic Conservation Laws

Lecture Notes on Hyperbolic Conservation Laws Lecture Notes on Hyperbolic Conservation Laws Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa. 16802, USA. bressan@math.psu.edu May 21, 2009 Abstract These notes provide

More information

THE STOKES SYSTEM R.E. SHOWALTER

THE STOKES SYSTEM R.E. SHOWALTER THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes

More information

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16. CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo

More information

MATH 220: Problem Set 3 Solutions

MATH 220: Problem Set 3 Solutions MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )

More information

MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20

MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20 MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20 MA 201 (2016), PDE 2 / 20 Vibrating string and the wave equation Consider a stretched string of length

More information

Entropy and Relative Entropy

Entropy and Relative Entropy Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for

More information

1 Energy dissipation in astrophysical plasmas

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

Math 124A October 11, 2011

Math 124A October 11, 2011 Math 14A October 11, 11 Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This corresponds to a string of infinite length. Although

More information

Numerical Solutions to Partial Differential Equations

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

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

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

Applications of the compensated compactness method on hyperbolic conservation systems

Applications of the compensated compactness method on hyperbolic conservation systems Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,

More information

Existence and exponential stability of the damped wave equation with a dynamic boundary condition and a delay term.

Existence and exponential stability of the damped wave equation with a dynamic boundary condition and a delay term. Existence and exponential stability of the damped wave equation with a dynamic boundary condition and a delay term. Stéphane Gerbi LAMA, Université de Savoie, Chambéry, France Fachbereich Mathematik und

More information

1 Introduction to PDE MATH 22C. 1. Introduction To Partial Differential Equations Recall: A function f is an input-output machine for numbers:

1 Introduction to PDE MATH 22C. 1. Introduction To Partial Differential Equations Recall: A function f is an input-output machine for numbers: 1 Introduction to PDE MATH 22C 1. Introduction To Partial Differential Equations Recall: A function f is an input-output machine for numbers: y = f(t) Output y 2R Input t 2R Name of function f t=independent

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

Entropy stable schemes for degenerate convection-diffusion equations

Entropy stable schemes for degenerate convection-diffusion equations Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.

More information

Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System

Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System Joshua Ballew University of Maryland College Park Applied PDE RIT March 4, 2013 Outline Description of the Model Relative Entropy Weakly

More information

Quasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths

Quasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths in the presence of plasma sheaths Department of Mathematical Sciences Ulsan National Institute of Science and Technology (UNIST) joint work with Masahiro Suzuki (Nagoya) and Chang-Yeol Jung (Ulsan) The

More information

Solitons : An Introduction

Solitons : An Introduction Solitons : An Introduction p. 1/2 Solitons : An Introduction Amit Goyal Department of Physics Panjab University Chandigarh Solitons : An Introduction p. 2/2 Contents Introduction History Applications Solitons

More information

Asymptotic Behavior of Waves in a Nonuniform Medium

Asymptotic Behavior of Waves in a Nonuniform Medium Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform

More information

Fractional Laplacian

Fractional Laplacian Fractional Laplacian Grzegorz Karch 5ème Ecole de printemps EDP Non-linéaire Mathématiques et Interactions: modèles non locaux et applications Ecole Supérieure de Technologie d Essaouira, du 27 au 30 Avril

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

Computing High Frequency Waves By the Level Set Method

Computing High Frequency Waves By the Level Set Method Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard

More information

MATH 205C: STATIONARY PHASE LEMMA

MATH 205C: STATIONARY PHASE LEMMA MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)

More information

M.Sc. in Meteorology. Numerical Weather Prediction

M.Sc. in Meteorology. Numerical Weather Prediction M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section

More information

Nonlinear stabilization via a linear observability

Nonlinear stabilization via a linear observability via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result

More information

A stochastic particle system for the Burgers equation.

A stochastic particle system for the Burgers equation. A stochastic particle system for the Burgers equation. Alexei Novikov Department of Mathematics Penn State University with Gautam Iyer (Carnegie Mellon) supported by NSF Burgers equation t u t + u x u

More information

Derivation of Generalized Camassa-Holm Equations from Boussinesq-type Equations

Derivation of Generalized Camassa-Holm Equations from Boussinesq-type Equations Derivation of Generalized Camassa-Holm Equations from Boussinesq-type Equations H. A. Erbay Department of Natural and Mathematical Sciences, Faculty of Engineering, Ozyegin University, Cekmekoy 34794,

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

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

Drag of a thin wing and optimal shape to minimize it

Drag of a thin wing and optimal shape to minimize it Drag of a thin wing and optimal shape to minimize it Alejandro Pozo December 21 st, 211 Outline 1 Statement of the problem 2 Inviscid compressible flows 3 Drag for supersonic case 4 Example of optimal

More information

Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation

Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Peter Howard, Texas A&M University University of Louisville, Oct. 19, 2007 References d = 1: Commun. Math. Phys. 269 (2007) 765

More information

A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations

A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with

More information

Chapter 16 Waves in One Dimension

Chapter 16 Waves in One Dimension Chapter 16 Waves in One Dimension Slide 16-1 Reading Quiz 16.05 f = c Slide 16-2 Reading Quiz 16.06 Slide 16-3 Reading Quiz 16.07 Heavier portion looks like a fixed end, pulse is inverted on reflection.

More information

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant

More information

MATH 220 solution to homework 5

MATH 220 solution to homework 5 MATH 220 solution to homework 5 Problem. (i Define E(t = k(t + p(t = then E (t = 2 = 2 = 2 u t u tt + u x u xt dx u 2 t + u 2 x dx, u t u xx + u x u xt dx x [u tu x ] dx. Because f and g are compactly

More information

J10M.1 - Rod on a Rail (M93M.2)

J10M.1 - Rod on a Rail (M93M.2) Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight

More information

Hyperbolic Systems of Conservation Laws. I - Basic Concepts

Hyperbolic Systems of Conservation Laws. I - Basic Concepts Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The

More information

Introduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0

Introduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0 Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods

More information

Hyperbolic Conservation Laws Past and Future

Hyperbolic Conservation Laws Past and Future Hyperbolic Conservation Laws Past and Future Barbara Lee Keyfitz Fields Institute and University of Houston bkeyfitz@fields.utoronto.ca Research supported by the US Department of Energy, National Science

More information

Conditional Symmetry Reduction and Invariant Solutions of Nonlinear Wave Equations

Conditional Symmetry Reduction and Invariant Solutions of Nonlinear Wave Equations Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 1, 229 233 Conditional Symmetry Reduction and Invariant Solutions of Nonlinear Wave Equations Ivan M. TSYFRA Institute of Geophysics

More information

Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions

Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Tom Elsden 1 Andrew Wright 1 1 Dept Maths & Stats, University of St Andrews DAMTP Seminar - 8th May 2017 Outline Introduction Coordinates

More information

Shock formation in the compressible Euler equations and related systems

Shock formation in the compressible Euler equations and related systems Shock formation in the compressible Euler equations and related systems Geng Chen Robin Young Qingtian Zhang Abstract We prove shock formation results for the compressible Euler equations and related systems

More information

A FREE BOUNDARY PROBLEM FOR TWO-DIMENSIONAL GAS DYNAMICS EQUATIONS

A FREE BOUNDARY PROBLEM FOR TWO-DIMENSIONAL GAS DYNAMICS EQUATIONS A FREE BOUNDARY PROBLEM FOR TWO-DIMENSIONAL GAS DYNAMICS EQUATIONS Katarina Jegdić Department of Computer and Mathematical Sciences University of Houston Downtown Sunčica Čanić, University of Houston Barbara

More information

Regularity of Weak Solution to Parabolic Fractional p-laplacian

Regularity of Weak Solution to Parabolic Fractional p-laplacian Regularity of Weak Solution to Parabolic Fractional p-laplacian Lan Tang at BCAM Seminar July 18th, 2012 Table of contents 1 1. Introduction 1.1. Background 1.2. Some Classical Results for Local Case 2

More information

Ray equations of a weak shock in a hyperbolic system of conservation laws in multi-dimensions

Ray equations of a weak shock in a hyperbolic system of conservation laws in multi-dimensions Proc. Indian Acad. Sci. (Math. Sci.) Vol. 126, No. 2, May 2016, pp. 199 206. c Indian Academy of Sciences Ray equations of a weak in a hyperbolic system of conservation laws in multi-dimensions PHOOLAN

More information

Hyperbolic Systems of Conservation Laws

Hyperbolic Systems of Conservation Laws Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/

More information

x ct x + t , and the characteristics for the associated transport equation would be given by the solution of the ode dx dt = 1 4. ξ = x + t 4.

x ct x + t , and the characteristics for the associated transport equation would be given by the solution of the ode dx dt = 1 4. ξ = x + t 4. . The solution is ( 2 e x+ct + e x ct) + 2c x+ct x ct sin(s)dx ( e x+ct + e x ct) + ( cos(x + ct) + cos(x ct)) 2 2c 2. To solve the PDE u xx 3u xt 4u tt =, you can first fact the differential operat to

More information

Math 311, Partial Differential Equations, Winter 2015, Midterm

Math 311, Partial Differential Equations, Winter 2015, Midterm Score: Name: Math 3, Partial Differential Equations, Winter 205, Midterm Instructions. Write all solutions in the space provided, and use the back pages if you have to. 2. The test is out of 60. There

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

Are Solitary Waves Color Blind to Noise?

Are Solitary Waves Color Blind to Noise? Are Solitary Waves Color Blind to Noise? Dr. Russell Herman Department of Mathematics & Statistics, UNCW March 29, 2008 Outline of Talk 1 Solitary Waves and Solitons 2 White Noise and Colored Noise? 3

More information

Lecture17: Generalized Solitary Waves

Lecture17: Generalized Solitary Waves Lecture17: Generalized Solitary Waves Lecturer: Roger Grimshaw. Write-up: Andrew Stewart and Yiping Ma June 24, 2009 We have seen that solitary waves, either with a pulse -like profile or as the envelope

More information

Maxwell s equations derived from minimum assumptions

Maxwell s equations derived from minimum assumptions Maxwell s equations derived from minimum assumptions Valery P.Dmitriyev (Dated: December 4, 2011) Maxwell s equations, under disguise of electromagnetic fields occurred in empty space, describe dynamics

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

SELF-SIMILAR SOLUTIONS FOR THE 2-D BURGERS SYSTEM IN INFINITE SUBSONIC CHANNELS

SELF-SIMILAR SOLUTIONS FOR THE 2-D BURGERS SYSTEM IN INFINITE SUBSONIC CHANNELS Bull. Korean Math. oc. 47 010, No. 1, pp. 9 37 DOI 10.4134/BKM.010.47.1.09 ELF-IMILAR OLUTION FOR THE -D BURGER YTEM IN INFINITE UBONIC CHANNEL Kyungwoo ong Abstract. We establish the existence of weak

More information

Answers to Problem Set Number 04 for MIT (Spring 2008)

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

Hyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan

Hyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global

More information

Hysteresis rarefaction in the Riemann problem

Hysteresis rarefaction in the Riemann problem Hysteresis rarefaction in the Riemann problem Pavel Krejčí 1 Institute of Mathematics, Czech Academy of Sciences, Žitná 25, 11567 Praha 1, Czech Republic E-mail: krejci@math.cas.cz Abstract. We consider

More information

Variational formulation of entropy solutions for nonlinear conservation laws

Variational formulation of entropy solutions for nonlinear conservation laws Variational formulation of entropy solutions for nonlinear conservation laws Eitan Tadmor 1 Center for Scientific Computation and Mathematical Modeling CSCAMM) Department of Mathematics, Institute for

More information

GLOBAL EXISTENCE FOR THE ONE-DIMENSIONAL SEMILINEAR TRICOMI-TYPE EQUATIONS

GLOBAL EXISTENCE FOR THE ONE-DIMENSIONAL SEMILINEAR TRICOMI-TYPE EQUATIONS GLOBAL EXISTENCE FOR THE ONE-DIMENSIONAL SEMILINEAR TRICOMI-TYPE EQUATIONS ANAHIT GALSTYAN The Tricomi equation u tt tu xx = 0 is a linear partial differential operator of mixed type. (For t > 0, the Tricomi

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme

More information

Cut-on, cut-off transition of sound in slowly varying flow ducts

Cut-on, cut-off transition of sound in slowly varying flow ducts Cut-on, cut-off transition of sound in slowly varying flow ducts Sjoerd W. Rienstra 19-3-1 Department of Mathematics and Computing Science, Eindhoven University of Technology, The Netherlands, s.w.rienstra@tue.nl

More information

Stochastic Volatility and Correction to the Heat Equation

Stochastic Volatility and Correction to the Heat Equation Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century

More information

BLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF TYPE WITH ARBITRARY POSITIVE INITIAL ENERGY

BLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF TYPE WITH ARBITRARY POSITIVE INITIAL ENERGY Electronic Journal of Differential Equations, Vol. 6 6, No. 33, pp. 8. ISSN: 7-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu BLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF

More information

New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows

New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows Plenary talk on the conference Stochastic and Analytic Methods in Mathematical Physics, Yerevan, Armenia,

More information

Deforming Composite Grids for Fluid Structure Interactions

Deforming Composite Grids for Fluid Structure Interactions Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,

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

2 u 1-D: 3-D: x + 2 u

2 u 1-D: 3-D: x + 2 u c 2013 C.S. Casari - Politecnico di Milano - Introduction to Nanoscience 2013-14 Onde 1 1 Waves 1.1 wave propagation 1.1.1 field Field: a physical quantity (measurable, at least in principle) function

More information

Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms

Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms Future Generation Computer Systems () 65 7 Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms R. Naidoo a,b, S. Baboolal b, a Department

More information

Modeling & Control of Hybrid Systems Chapter 4 Stability

Modeling & Control of Hybrid Systems Chapter 4 Stability Modeling & Control of Hybrid Systems Chapter 4 Stability Overview 1. Switched systems 2. Lyapunov theory for smooth and linear systems 3. Stability for any switching signal 4. Stability for given switching

More information

Blowup for Hyperbolic Equations. Helge Kristian Jenssen and Carlo Sinestrari

Blowup for Hyperbolic Equations. Helge Kristian Jenssen and Carlo Sinestrari Blowup for Hyperbolic Equations Helge Kristian Jenssen and Carlo Sinestrari Abstract. We consider dierent situations of blowup in sup-norm for hyperbolic equations. For scalar conservation laws with a

More information

Reading: P1-P20 of Durran, Chapter 1 of Lapidus and Pinder (Numerical solution of Partial Differential Equations in Science and Engineering)

Reading: P1-P20 of Durran, Chapter 1 of Lapidus and Pinder (Numerical solution of Partial Differential Equations in Science and Engineering) Chapter 1. Partial Differential Equations Reading: P1-P0 of Durran, Chapter 1 of Lapidus and Pinder (Numerical solution of Partial Differential Equations in Science and Engineering) Before even looking

More information

n v molecules will pass per unit time through the area from left to

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

Resolvent Estimates and Quantification of Nonlinear Stability

Resolvent Estimates and Quantification of Nonlinear Stability Resolvent Estimates and Quantification of Nonlinear Stability Heinz Otto Kreiss Department of Mathematics, UCLA, Los Angeles, CA 995 Jens Lorenz Department of Mathematics and Statistics, UNM, Albuquerque,

More information

NONLOCAL DIFFUSION EQUATIONS

NONLOCAL DIFFUSION EQUATIONS NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,

More information

On some weighted fractional porous media equations

On some weighted fractional porous media equations On some weighted fractional porous media equations Gabriele Grillo Politecnico di Milano September 16 th, 2015 Anacapri Joint works with M. Muratori and F. Punzo Gabriele Grillo Weighted Fractional PME

More information

Singularity formation for compressible Euler equations

Singularity formation for compressible Euler equations Singularity formation for compressible Euler equations Geng Chen Ronghua Pan Shengguo Zhu Abstract In this paper, for the p-system and full compressible Euler equations in one space dimension, we provide

More information

THE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK

THE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK THE INTERACTION OF TURBULENCE WITH THE HELIOSPHERIC SHOCK G.P. Zank, I. Kryukov, N. Pogorelov, S. Borovikov, Dastgeer Shaikh, and X. Ao CSPAR, University of Alabama in Huntsville Heliospheric observations

More information

Receiver. Johana Brokešová Charles University in Prague

Receiver. Johana Brokešová Charles University in Prague Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate

More information

Waves in a Shock Tube

Waves in a Shock Tube Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid

More information