Nonlinear system of mixed type and its application to steady Euler-Poisson system
|
|
- Patience Stokes
- 5 years ago
- Views:
Transcription
1 The 1st Meeting of Young Researchers in PDEs Nonlinear system of mixed type and its application to steady Euler-Poisson system Myoungjean Bae (POSTECH) -based on collaborations with- B. Duan, J. Xiao, C. Xie, H. Park, Y. Park, and H. Ryu August 18, /23
2 Introduction ρ : (electron) density, u : velocity, p : pressure, Φ : Coulomb potential Steady Euler-Poisson system div(ρu) = 0 div(ρu u+pi) = ρ Φ div(ρub) = ρu Φ Φ = ρ b Bernoulli function ( b > 0 : background charge density) B = 1 2 u 2 + γp (γ 1)ρ, γ > 1 Entropy S = ln p ρ γ Transport equations If (ρ,u,p, Φ) is a C 1 solution, then div(ρus) = div(ρu(b Φ)) = 0 div(ρu)=0 ρu S = ρu (B Φ) = 0 ** We call B Φ pseudo-bernoulli function. 2/23
3 Rankine-Hugoniot conditions Ω n Γ Ω + divq = 0 in Ω. (1) Q L 2 (Ω) C 1 (Ω ± ) C 0 (Ω ± Γ) satisfies Q Dφ dx = 0 for any φ C0 (Ω) Ω Q satisfies (1) in Ω ± and the relation [Q n] Γ := Q n x Ω Γ Q n x Ω + Γ = 0 (2) because Q Dφ = Ω + Ω Γ [Q n] φ φ divq Ω + Ω 3/23
4 Discontinuous transition [ρu n] Γ = [ρ(u n)u+pn] Γ = [ρ(u n)b] Γ = 0 Let τ be a tangential vector field on Γ. { τ [ρ(u n s )u+pn] Γ = 0 ρ(u n)[u τ] Γ = 0 n [ρ(u n)u+pn] Γ = 0 [ρ(u n) 2 +p] Γ = 0 Case 1. If u n 0 on Γ, then [ρu n] Γ = [u τ] Γ = [ρ(u n) 2 +p] Γ = [B] Γ = 0 (shock) Case 2. If u n = 0 on Γ, then [p] = 0 (contact discontinuity) 4/23
5 Example of shocks 5/23
6 Motivation: De Laval Nozzle problem for Euler system Given (ρ,u,p ) with M < 1, A A ex 0 p A t x p max p 0 x Figure: Supersonic Flow And Shock Waves by Courant-Friedrich 6/23
7 1-d solutions for E-P system References:Ascher-Markowich-Pietra-Schmeiser(1991),Luo-Xin(2012) etc. (For b 0 < ρ sonic ) E ρ = ρs b ρ Figure 1 Phase portrait for 0 < b < ρs One dimensional smooth transonic-transonic shock solution is here!! 7/23
8 Main Problem To construct multidimensional smooth transonic-transonic shock configuration for prescribed exit pressure in flat nozzle where proper electric field applies through nozzle. Γ w : u n w = 0 Γ sonic : M = 1 Γ shock M 0 < 1 M < 1 M > 1 M < 1 p = p ex x = s 0 x = s 1 Figure: Smooth transonic-transonic shock configuration Due to Φ = ρ b, a transonic shock problem is a two-phase problem with a free boundary to be determined. 8/23
9 Why is this problem difficult? Analytic property of E-P system Suppose u = 0 i.e., u = ϕ. If B Φ = 0 (pseudo-bernoulli s law), ( div (Φ 1 ) 2 ϕ 2 ) 1 γ 1 ϕ = 0, (3) Φ = (Φ 1 2 ϕ 2 ) 1 γ 1 b0 (4) Flow type M < 1 M = 1 M > 1 (3) Elliptic Degenerate Hyperbolic (3) (4) Elliptic system D.M.-E coupled system H-E coupled system W.-P. Unknown Unknown at all Unknown Note M < 1 ϕ 2 < 2(γ 1) γ+1 Φ 9/23
10 Recent results M-D subsonic flow of E-P system with u = 0 for fixed exit pressure (B.-Duan-Xie, ARMA 2016) 2-D subsonic flow of full E-P system for fixed exit pressure (B.-Duan-Xie, SIMA 2014) 2-D subsonic flow of E-P system with self-gravitation (B.-Duan-Xie, M3AS 2015) Φ = ρ 3-D subsonic flow of E-P system with nonzero swirl(=angular momentum) (B.-Weng, submitted) 2-D supersonic flow of full E-P system (B.-Duan-Xie, preprint) and so on... 10/23
11 Subsonic flow (M < 1 ϕ 2 < 2(γ 1) γ+1 Φ) Problem 1. Subsonic flow for potential flow model Given p ex with p ex p + (L) ( α,γl) 1,α,Γ L small, find a solution (ϕ,φ) to E-P system for potential flow with B.C. given as follows: x R n 1 ϕ = 0 Φ = Φ + Λ R n 1 n w p = p ex Φ = Φ + 0 nw ϕ = 0 = nw Φ L x 1 11/23
12 Theorem (B.-Duan-Xie, ARMA 2016) α (0,1), σ 0 > 0 so that if σ σ 0, then the NLBVP has a unique solution (ϕ,φ) C 1,α (N) C 2 (N) s.t. ϕ ϕ + ( 1 α,corners) 2,α,N + Φ Φ + ( 1 α,corners) 2,α,N C p ex p + (L) ( α, Λ) 1,α,Λ. Essential observation Same coupled terms ( div (Φ 1 ) 2 ϕ 2 ) 1 γ 1 ϕ = 0, Φ = (Φ 1 2 ϕ 2 ) 1 γ 1 b0 Nonlinear boundary condition p = p ex on x 1 = L Φ 1 2 ϕ 2 = p γ 1 γ ex on x 1 = L 12/23
13 Problem 2. 3-D subsonic flow with nonzero vorticity Find axi-symmetric subsonic flow of E-P system with B.C. given as follows: (r,θ) Γ en : u e r = 0 Φ = Φ 0 Λ = B 1 (0) R 2 n w Γ ex : p = p ex 0 Γ w : u n w = nw Φ = 0 L x Helmholtz decomposition for axi-symmetric flow with nonzero swirl u = = ϕ(x,r) }{{} curl free(compressibility) ( 1 r r(rψ)+ x ϕ ) + (h(x,r)e r +ψ(x,r)e θ ) }{{} div free(vorticity) e x +( x ψ + r ϕ)e r + x h(x,r)e θ }{{} nonzero swirl 13/23
14 Corresponding equations derived from E-P system ( div H(S,K+Φ 1 ) 2 u(r,ψ,dψ,dϕ,λ) 2 )u(r,ψ,dψ,dϕ,λ) = 0, }{{} (=:M) Φ = H(S,K+Φ 1 2 u(r,ψ,dψ,dϕ,λ) 2 ) b, (ψe θ ) = T(ρ,S) rs r K+ Λ r 2 r Λ 1 r e θ =: fe θ, r(rψ)+ x ϕ M (S,K,Λ) = 0. Theorem (B.-Weng, Submitted/ArXiv) Given a small perturbation of (b,p ex,s en,k en,λ en ), there exists a 3-D axisymmetric subsonic solution. The uniqueness (in small perturbation regime) is obtained under an additional assumption of boundary data. 14/23
15 Elliptic equation for ψ with a singular coefficient, W = ψ(x,r)e θ Either we find a solution ψ to ( xx + 1 r r(r r ) 1 r 2 ) ψ = f(x,r) in N = (0,L) B1 (0), s.t. ψ(x,r)e θ is C 2,α (N), or, we solve x ψ(0,r) = 0 on Γ en, ψ = 0 on Γ w, x ψ = 0 on Γ ex, ψ = 0 on N {r = 0}, W = f(x,r)e θ in N = (0,L) B 1 (0), x W = 0 on Γ en Γ ex, W = 0 on Γ w, then show that W is in the form of W = ψ(x,r)e θ. Compatibility condition for axisymmetric flow f(x, 0) = 0 15/23
16 Proposition For a fixed α (0,1), suppose that a vector field f(x,r)e θ : N R 3 is C α in N, and that f satisfies the compatibility condition f(x,0) = 0. Then the linear boundary value problem has a unique solution W : N R 3. Furthermore, W is represented as where ψ satisfies the estimate W(x) = ψ(x,r)e θ in N, (5) ψ 2,α,[0,L] [0,1] C f α,n. Note that W is well defined by (5) due to the condition ψ(x,0) = 0 for all x [0,L]. 16/23
17 Idea of Proof The unique existence of C 2,α solution W is easily obtained. 1. Let W = W 1 (x)e 1 +W 2 (x)e 2 +W 3 (x)e 3 be a C 2 solution. Then, V 1 0 in N. for W = W 2 e 2 +W 3 e 3 =: U r e r +U θ e θ in N { U r (x,r,θ) = W 2 (x,rcosθ,rsinθ)cosθ +W 3 (x,rcosθ,rsinθ)sinθ U θ (x,r,θ) = W 2 (x,rcosθ,rsinθ)sinθ +W 3 (x,rcosθ,rsinθ)cosθ. W C 2,α (U r,u θ ) C 2,α ([0,L] [0,1] T) 17/23
18 2. U r and U θ satisfy { L cyl 1 (U r,u θ ) := ( x 2 + r r r + 1 r 2 2 θ ) 1 r 2 Ur 2 r 2 θ U θ = 0, L cyl 2 (U r,u θ ) := ( x r + 1 r r + 1 r 2 2 θ ) 1 r 2 Uθ + 2 r 2 θ U r = f(x,r), 3. For each n N, set U n (x,r,θ) := 1 Then, so, U n θ 2 n 1 2 n k=0 U r ( x,r,θ + 2πk 2 n ) 1 (x,r,θ) := 2 n 1 2 n k=0 U ( ) θ x,r,θ + 2πk 2 n (U n r,u n θ )(x,r,θ) = (U n r,u n θ )(x,r,θ + 2πj 2 n ) 0 j 2n 1, n 1, (L cyl 1,Lcyl 2 )(Un r,un θ ) = (0, f). And, {(U n r,u n θ )} is uniformly bounded in C2,α ([0,L] [0,1] T). 18/23
19 4. By Arzelá-Ascoli Theorem, a subsequence {(U n k r,un k θ )} that converges to U = (Ur,U θ ), where U satisfies (L cyl 1,Lcyl 2 )U = (0, f), and U (x,r,θ) = U (x,r,θ + 2πj 2 nk ) 0 j 2 n k 1, k 1 U (x,r,θ) = U (x,r,θ +2κπ) 0 κ 1. 19/23
20 2-D supersonic flow (M > 1 ϕ 2 > 2(γ 1) γ+1 Φ) Problem 3. Given u en with u en u 0 C4 (Γ 0) small, find a supersonic potential replacements flow (ϕ,φ) of E-P system with B. C. given as follows: x 2 ϕ = 0,ϕ x1 = u en Φ x1 = E en n w Φ = Φ0 0 nw ϕ = 0 = nw Φ L x 1 Figure: D k x 2 u en = 0 on Λ for k = 1,2,3 20/23
21 Fix γ > 1, and set m 0 = J 0 S 1 γ 1 0 for J 0 = ρ 0 u 0. In N = (0,L) ( 1,1), set σ 1 := u en u 0 C 4 ([ 1,1]), σ 2 := b b 0 C 2 (N) + E en E 0 C 4 ([ 1,1]) Theorem (B.-Duan-Xiao-Xie, preprint) Given γ > 1,L > 0,m 0 > 0,u 0 > u sonic (= (γm γ 1 0 ) 1 γ+1 ),E0 > 0,b 0 > 0, S 0 > 1 depending only on γ,l,m 0,u 0,E 0 s.t. the following holds: If S 0 S 0 (J 0 small), then, σ 1, σ 2 > 0 such that if σ 1 σ 1, and σ 2 σ 2, then NLBVP II has a unique solution (ϕ,φ) [H 4 (N) C 2 (N)] 2 satisfying ϕ ϕ 0 H 4 (N) +S 1 2(γ 1) 0 Φ Φ 0 H 4 (N) C(σ 1 +S 1 2(γ 1) 0 σ 2 ). Question Is the condition S 0 S 0 ( J 0 << 1) simply technical? Or, is there a physical reason? 21/23
22 Ongoing works Tranosnic flows of E-P system M-d smooth transonic flow of E-P system (with B. Duan & C. Xie) M-d transonic shock solution of E-P system in flat nozzle (with B. Duan & C. Xie) M-d transonic shock solution of E-P system in convergent nozzle (with Y. Park) Application of Helmholtz decomposition to transonic shock solutions (with Y. Park) M-d supersonic flow of E-P system (with H. Ryu) Application of Helmholtz decomposition Transonic shocks (with Y. Park) Contact discontinuity (with H. Park) 22/23
23 Ongoing works Tranosnic flows of E-P system M-d smooth transonic flow of E-P system (with B. Duan & C. Xie) M-d transonic shock solution of E-P system in flat nozzle (with B. Duan & C. Xie) M-d transonic shock solution of E-P system in convergent nozzle (with Y. Park) Application of Helmholtz decomposition to transonic shock solutions (with Y. Park) M-d supersonic flow of E-P system (with H. Ryu) Application of Helmholtz decomposition Transonic shocks (with Y. Park) Contact discontinuity (with H. Park) THANK YOU!!! 23/23
Shock Reflection-Diffraction, Nonlinear Conservation Laws of Mixed Type, and von Neumann s Conjectures 1
Contents Preface xi I Shock Reflection-Diffraction, Nonlinear Conservation Laws of Mixed Type, and von Neumann s Conjectures 1 1 Shock Reflection-Diffraction, Nonlinear Partial Differential Equations of
More informationShock Reflection-Diffraction, Nonlinear Partial Differential Equations of Mixed Type, and Free Boundary Problems
Chapter One Shock Reflection-Diffraction, Nonlinear Partial Differential Equations of Mixed Type, and Free Boundary Problems Shock waves are steep fronts that propagate in compressible fluids when convection
More informationTHE TRANSONIC SHOCK IN A NOZZLE FOR 2-D AND 3-D COMPLETE EULER SYSTEMS
THE TRANSONIC SHOCK IN A NOZZLE FOR -D AND 3-D COMPLETE EULER SYSTEMS Zhouping Xin Department of Mathematics and IMS, CUHK, Shatin, N.T., Hong Kong Huicheng Yin Department of Mathematics and IMS, Nanjing
More informationEfficient solution of stationary Euler flows with critical points and shocks
Efficient solution of stationary Euler flows with critical points and shocks Hans De Sterck Department of Applied Mathematics University of Waterloo 1. Introduction consider stationary solutions of hyperbolic
More informationSimple waves and a characteristic decomposition of the two dimensional compressible Euler equations
Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Jiequan Li 1 Department of Mathematics, Capital Normal University, Beijing, 100037 Tong Zhang Institute
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationDirect Numerical Solution of the Steady 1D Compressible Euler Equations for Transonic Flow Profiles with Shocks
Direct Numerical Solution of the Steady 1D Compressible Euler Equations for Transonic Flow Profiles with Shocks Hans De Sterck, Scott Rostrup Department of Applied Mathematics, University of Waterloo,
More informationA Study of Transonic Flow and Airfoils. Presented by: Huiliang Lui 30 th April 2007
A Study of Transonic Flow and Airfoils Presented by: Huiliang Lui 3 th April 7 Contents Background Aims Theory Conservation Laws Irrotational Flow Self-Similarity Characteristics Numerical Modeling Conclusion
More informationON TRANSONIC SHOCKS IN A NOZZLE WITH VARIABLE END PRESSURES. Jun Li (Department of Mathematics & IMS, Nanjing University, Nanjing , P.R.
ON TRANSONIC SHOCKS IN A NOZZLE WITH VARIABLE END PRESSURES Jun Li Department of Mathematics & IMS, Nanjing University, Nanjing 210093, P.R.China Zhouping Xin Department of Mathematics IMS, CUHK, Shatin,
More informationTransonic Shocks for the Full Compressible Euler System in a General Two-Dimensional De Laval Nozzle
Arch. Rational Mech. Anal. 27 23 533 58 Digital Object Identifier DOI.7/s25-2-58-x Transonic Shocks for the Full Compressible Euler System in a General Two-Dimensional De Laval Nozzle Jun Li, Zhouping
More informationTHREE-DIMENSIONAL INTERACTION OF SHOCKS IN IRROTATIONAL FLOWS. À la mémoire de Michelle, de laquelle j ai beaucoup appris.
THREE-DIMENSIONAL INTERACTION OF SHOCKS IN IRROTATIONAL FLOWS DENIS SERRE UMPA, UMR CNRS ENS Lyon # 5669. École Normale Supérieure de Lyon 46, allée d Italie 69364 Lyon, cedex 07, France. À la mémoire
More informationFree Boundary Problems in Shock Reflection/ Diffraction and Related Transonic Flow Problems
Report no. OxPDE-15/06 Free Boundary Problems in Shock Reflection/ Diffraction and Related Transonic Flow Problems by Gui-Qiang Chen University of Oxford and Mikhail Feldman University of Wisconsin Oxford
More informationA 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 informationNonlinear stability of compressible vortex sheets in two space dimensions
of compressible vortex sheets in two space dimensions J.-F. Coulombel (Lille) P. Secchi (Brescia) CNRS, and Team SIMPAF of INRIA Futurs Evolution Equations 2006, Mons, August 29th Plan 1 2 3 Related problems
More informationShock reflection in gas dynamics
Shock reflection in gas dynamics Denis Serre École Normale Supérieure de Lyon À ma Mère Abstract This paper is about multi-dimensional shocks and their interactions. The latter take place either between
More informationSELF-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 informationConical Shock Waves for Isentropic Euler System
Conical Shock Waves for Isentropic Euler System Shuxing Chen Institute of Mathematical Research, Fudan University, Shanghai, China E-mail: sxchen@public8.sta.net.cn Dening Li Department of Mathematics,
More information2D compressible vortex sheets. Paolo Secchi
2D compressible vortex sheets Paolo Secchi Department of Mathematics Brescia University Joint work with J.F. Coulombel EVEQ 2008, International Summer School on Evolution Equations, Prague, Czech Republic,
More informationarxiv: v1 [math.ap] 21 Nov 2013
1 arxiv:1311.5596v1 [math.ap] 21 Nov 2013 SHOCK REFLECTION-DIFFRACTION, VON NEUMANN S CONJECTURES, AND NONLINEAR EQUATIONS OF MIXED TYPE GUI-QIANG CHEN AND MIKHAIL FELDMAN Abstract. Shock waves are fundamental
More informationNotes #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 informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationThin 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 informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationTHE 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 informationarxiv:astro-ph/ v1 14 Jan 2002
Accretion onto a Black Hole in the Presence of Bremsstrahlung Cooling - Parameter Space Study Sivakumar G. Manickam and Sandip K. Chakrabarti arxiv:astro-ph/020203v 4 Jan 2002 S. N. Bose National Centre
More informationPropagation of discontinuities in solutions of First Order Partial Differential Equations
Propagation of discontinuities in solutions of First Order Partial Differential Equations Phoolan Prasad Department of Mathematics Indian Institute of Science, Bangalore 560 012 E-mail: prasad@math.iisc.ernet.in
More informationMATH 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 informationA RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS
A RIEMANN PROBLEM FOR THE ISENTROPIC GAS DYNAMICS EQUATIONS KATARINA JEGDIĆ, BARBARA LEE KEYFITZ, AND SUN CICA ČANIĆ We study a Riemann problem for the two-dimensional isentropic gas dynamics equations
More informationDrag 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 information0.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É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 informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationFinal: 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 informationA Free Boundary Problem for a Quasi-linear Degenerate Elliptic Equation: Regular Reflection of Weak Shocks
A Free Boundary Problem for a Quasi-linear Degenerate Elliptic Equation: Regular Reflection of Weak Shocks SUNČICA ČANIĆ BARBARA LEE KEYFITZ AND EUN HEUI KIM University of Houston Abstract We prove the
More informationShock Waves. 1 Steepening of sound waves. We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: kˆ.
Shock Waves Steepening of sound waves We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: v u kˆ c s kˆ where u is the velocity of the fluid and k is the wave
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationCOMPARISON PRINCIPLES FOR SELF-SIMILAR POTENTIAL FLOW
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 COMPARISON PRINCIPLES FOR SELF-SIMILAR POTENTIAL FLOW GUI-QIANG G. CHEN AND MIKHAIL FELDMAN Abstract.
More informationRarefaction wave interaction for the unsteady transonic small disturbance equations
Rarefaction wave interaction for the unsteady transonic small disturbance equations Jun Chen University of Houston Department of Mathematics 4800 Calhoun Road Houston, TX 77204, USA chenjun@math.uh.edu
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 27 May, 2004 1.30 to 3.30 PAPER 64 ASTROPHYSICAL FLUID DYNAMICS Attempt THREE questions. There are four questions in total. The questions carry equal weight. Candidates
More informationThe forced response of choked nozzles and supersonic diffusers
J. Fluid Mech. (27), vol. 585, pp. 28 34. c 27 Cambridge University Press doi:.7/s22276647 Printed in the United Kingdom 28 The forced response of choked nozzles and supersonic diffusers WILLIAM H. MOASE,
More informationShock 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 informationWeak-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 informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationBlowup phenomena of solutions to the Euler equations for compressible fluid flow
J. Differential Equations 1 006 91 101 www.elsevier.com/locate/jde Blowup phenomena of solutions to the Euler equations for compressible fluid flow Tianhong Li a,, Dehua Wang b a Department of Mathematics,
More informationAST242 LECTURE NOTES PART 2
AST242 LECTURE NOTES PART 2 Contents 1. Inviscid Barotropic flow 2 1.1. Enthalpy and the Bernoulli equation 2 1.2. Dimensional analysis for Bondi accretion 3 1.3. Bondi accretion 5 1.4. The transonic point
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More information2013/5/22. ( + ) ( ) = = momentum outflow rate. ( x) FPressure. 9.3 Nozzles. δ q= heat added into the fluid per unit mass
9.3 Nozzles (b) omentum conservation : (i) Governing Equations Consider: nonadiabatic ternal (body) force ists variable flow area continuously varying flows δq f ternal force per unit volume +d δffdx dx
More informationFree Boundary Problems for Nonlinear Wave Equations: Interacting Shocks
Multiphase Fluid Flows and Multi-Dimensional Hyperbolic Problems, Newton Institute, March 31-April 4, 23 1 Free Boundary Problems for Nonlinear Wave Equations: Interacting Shocks Barbara Lee Keyfitz a
More information0.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 informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
More informationOn the Dependence of Euler Equations on Physical Parameters
On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,
More informationLEAST-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 informationQuantum Hydrodynamic Systems and applications to superfluidity at finite temperatures
Quantum Hydrodynamic Systems and applications to superfluidity at finite temperatures Paolo Antonelli 1 Pierangelo Marcati 2 1 Gran Sasso Science Institute, L Aquila 2 Gran Sasso Science Institute and
More information58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014
Professor Fred Stern Fall 04 Chapter 7 Bluff Body Fluid flows are broadly categorized:. Internal flows such as ducts/pipes, turbomachinery, open channel/river, which are bounded by walls or fluid interfaces:
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationMixed hyperbolic-elliptic systems in self-similar flows
Bol. Soc. Bras. Mat., Vol.32, No. 3, 1-23 2001, Sociedade Brasileira de Matemática Mixed hyperbolic-elliptic systems in self-similar flows Sunčica Čanić, Barbara Lee Keyfitz and Eun Heui Kim Dedicated
More informationNumerical Simulation of Supersonic Plume in a Cowl-Plug Nozzle Configuration}
Freund Publishing House Ltd. London International Journal of Turbo and Jet Engines, 16, 241-253(1999) Numerical Simulation of Supersonic Plume in a Cowl-Plug Nozzle Configuration} Tony W. H. Sheu Department
More informationQUANTUM EULER-POISSON SYSTEMS: EXISTENCE OF STATIONARY STATES. 1. Introduction
ARCHIVUM MATHEMATICUM BRNO) Tomus 4 24), 435 456 QUANTUM EULER-POISSON SYSTEMS: EXISTENCE OF STATIONARY STATES ANSGAR JÜNGEL, HAILIANG LI Abstract. A one-dimensional quantum Euler-Poisson system for semiconductors
More informationOffshore Hydromechanics Module 1
Offshore Hydromechanics Module 1 Dr. ir. Pepijn de Jong 4. Potential Flows part 2 Introduction Topics of Module 1 Problems of interest Chapter 1 Hydrostatics Chapter 2 Floating stability Chapter 2 Constant
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationChapter 9 Flow over Immersed Bodies
57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 009 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
More informationOn the Stokes semigroup in some non-helmholtz domains
On the Stokes semigroup in some non-helmholtz domains Yoshikazu Giga University of Tokyo June 2014 Joint work with Ken Abe (Nagoya U.), Katharina Schade (TU Darmstadt) and Takuya Suzuki (U. Tokyo) Contents
More informationFundamentals of Magnetic Island Theory in Tokamaks
Fundamentals of Magnetic Island Theory in Tokamaks Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX, USA Talk available at http://farside.ph.utexas.edu/talks/talks.html
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationCompressible Potential Flow: The Full Potential Equation. Copyright 2009 Narayanan Komerath
Compressible Potential Flow: The Full Potential Equation 1 Introduction Recall that for incompressible flow conditions, velocity is not large enough to cause density changes, so density is known. Thus
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationMultidimensional Conservation Laws Part II: Multidimensional Models and Problems. BCAM and UPV/EHU Courses 2010/2011
Multidimensional Conservation Laws Part II: Multidimensional Models and Problems Gui-Qiang G. Chen Oxford Centre for Nonlinear PDE Mathematical Institute, University of Oxford Website: http://people.maths.ox.ac.uk/chengq/
More informationarxiv: v1 [math.ap] 25 Sep 2017
arxiv:1709.09958v1 [math.ap] 25 Sep 2017 A variational inequality formulation for transonic compressible steady potential flows: Radially symmetric transonic shock Yung-Sze Choi Eun Heui Kim Abstract We
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationAlberto Bressan Convergence Rates for Viscous Approximations in the Presence of Linearly Degenerate Fields Gui-Qiang Chen
Day 1: June 12 Bus will leave Faculty Club for Minhang campus at 7:50 08:50-09:00 Opening ceremony Session 1 09:00--10:20 Chair: Yachun Li Alberto Bressan Convergence Rates for Viscous Approximations in
More informationOn weak solution approach to problems in fluid dynamics
On weak solution approach to problems in fluid dynamics Eduard Feireisl based on joint work with J.Březina (Tokio), C.Klingenberg, and S.Markfelder (Wuerzburg), O.Kreml (Praha), M. Lukáčová (Mainz), H.Mizerová
More informationDSMC simulations of thermal escape
DSMC simulations of thermal escape Alexey N. Volkov, R.E. Johnson, O.J. Tucker, J.T. Erwin Materials Science & Engineering University of Virginia, USA Financial support is provided by NASA through Planetary
More informationME185 Introduction to Continuum Mechanics
Fall, 0 ME85 Introduction to Continuum Mechanics The attached pages contain four previous midterm exams for this course. Each midterm consists of two pages. As you may notice, many of the problems are
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationIncompressible Navier-Stokes Equations in R 3
Incompressible Navier-Stokes Equations in R 3 Zhen Lei ( ) School of Mathematical Sciences Fudan University Incompressible Navier-Stokes Equations in R 3 p. 1/5 Contents Fundamental Work by Jean Leray
More informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationIntroduction to Gas Dynamics All Lecture Slides
Introduction to Gas Dynamics All Lecture Slides Teknillinen Korkeakoulu / Helsinki University of Technology Autumn 009 1 Compressible flow Zeroth law of thermodynamics 3 First law of thermodynamics 4 Equation
More informationAerodynamics. Lecture 1: Introduction - Equations of Motion G. Dimitriadis
Aerodynamics Lecture 1: Introduction - Equations of Motion G. Dimitriadis Definition Aerodynamics is the science that analyses the flow of air around solid bodies The basis of aerodynamics is fluid dynamics
More informationF1.9AB2 1. r 2 θ2 + sin 2 α. and. p θ = mr 2 θ. p2 θ. (d) In light of the information in part (c) above, we can express the Hamiltonian in the form
F1.9AB2 1 Question 1 (20 Marks) A cone of semi-angle α has its axis vertical and vertex downwards, as in Figure 1 (overleaf). A point mass m slides without friction on the inside of the cone under the
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationMultiscale Hydrodynamic Phenomena
M2, Fluid mechanics 2014/2015 Friday, December 5th, 2014 Multiscale Hydrodynamic Phenomena Part I. : 90 minutes, NO documents 1. Quick Questions In few words : 1.1 What is dominant balance? 1.2 What is
More informationK. Ambika and R. Radha
Indian J. Pure Appl. Math., 473: 501-521, September 2016 c Indian National Science Academy DOI: 10.1007/s13226-016-0200-9 RIEMANN PROBLEM IN NON-IDEAL GAS DYNAMICS K. Ambika and R. Radha School of Mathematics
More informationVorticity and Dynamics
Vorticity and Dynamics In Navier-Stokes equation Nonlinear term ω u the Lamb vector is related to the nonlinear term u 2 (u ) u = + ω u 2 Sort of Coriolis force in a rotation frame Viscous term ν u = ν
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationHyperbolic 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 informationTurbulent spectra generated by singularities
Turbulent spectra generated by singularities p. Turbulent spectra generated by singularities E.A. Kuznetsov Landau Institute for Theoretical Physics, Moscow Russia In collaboration with: V. Naulin A.H.
More informationEntropy 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 informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationThe Riemann problem. The Riemann problem Rarefaction waves and shock waves
The Riemann problem Rarefaction waves and shock waves 1. An illuminating example A Heaviside function as initial datum Solving the Riemann problem for the Hopf equation consists in describing the solutions
More information