A FREE BOUNDARY PROBLEM FOR TWO-DIMENSIONAL GAS DYNAMICS EQUATIONS
|
|
- Alexander Rodgers
- 5 years ago
- Views:
Transcription
1 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 Lee Keyfitz, Ohio State University Texas PDE Seminar University of Texas Pan American March 26-27, 2011
2 Reflection of a shock by a wedge INTRODUCTION two-dimensional Riemann problems for systems of conservation laws steady / unsteady transonic small disturbance equations nonlinear wave system gas dynamics equations (isentropic, adiabatic) shock reflection regular reflection (strong, weak) Mach reflection Guderley reflection self-similar coordinates 8 < : mixed type (hyperbolic, elliptic) system free boundary value problem nonlinear problem theory of second order elliptic equations with mixed boundary conditions (Dirichlet and oblique) and fixed point theory Result Riemann problem for the isentropic gas dynamics equations strong regular reflection existence of a local solution in weighted Holder spaces A free boundary problem for two-dimensional gas dynamics equations 1
3 SHOCK REFLECTION PSfrag replacements flow reflected shock incident shock t < 0 t = 0 t > 0 ments motivating problem: airflow past object types of reflection strong regular reflection weak regular reflection Mach reflection inc. shock refl. shock refl. shock inc. shock refl. shock subsonic inc. shock subsonic subsonic region region region Mach stem supersonic region A free boundary problem for two-dimensional gas dynamics equations 2
4 ements TWO-DIMENSIONAL RIEMANN PROBLEMS U t + F (U) x + G(U) y = 0, U : R 2 [0, ) R m y similarity variables ξ = x/t and η = y/t U 1 reduced system: U 0 U1 y = kx ξ (F ξu) + η (G ηu) = 2U = kx initial data becomes boundary data at infinity x the reduced system changes type: hyperbolic in the far field, mixed type near the origin resolve one-dimensional discontinuities in the far field formulate a free boundary problem for the reflected shock and the subsonic state theory of second order elliptic equations with mixed boundary conditions (Gilbarg, Trudinger, Lieberman) fixed point theorems A free boundary problem for two-dimensional gas dynamics equations 3
5 adiabatic gas dynamics equations E = 1 γ 1 p + u2 +v 2 ρ 2 ρ t + (ρu) x + (ρv) y = 0 (ρu) t + (ρu 2 + p) x + (ρuv) y = 0 (ρv) t + (ρuv) x + (ρv 2 + p) y = 0 (ρe) t + (u(ρe + p)) x + (v(ρe + p)) y = 0 isentropic gas dynamics equations nonlinear wave system p = p(ρ) ρ t + (ρu) x + (ρv) y = 0 (ρu) t + (ρu 2 + p) x + (ρuv) y = 0 (ρv) t + (ρuv) x + (ρv 2 + p) y = 0 p = p(ρ) ρ t + (ρu) x + (ρv) y = 0 (ρu) t + p x = 0 (ρv) t + p y = 0 unsteady transonic small disturbance equation u t + uu x + v y = 0, v x + u y = 0 Euler equations for potential flow, pressure gradient system A free boundary problem for two-dimensional gas dynamics equations 4
6 ISENTROPIC GAS DYNAMICS EQUATIONS lacements ρ t + (ρu) x + (ρv) y = 0 y (ρu) t + (ρu 2 + p) x + (ρuv) y = 0 x = ky (ρv) t + (ρuv) x + (ρv 2 + p) y = 0 U 0 U 1 x p = p(ρ), c 2 (ρ) := p (ρ) > 0 and increasing x = ky U = (ρ, u, v) Initial data: U 0 = (ρ 0, u 0, 0), U 1 = (ρ 1, 0, 0) - fix ρ 0 > ρ 1 > 0 - fix k (k C (ρ 0, ρ 1 ), ) - set u 0 = 1 + k 2 q (ρ 0 ρ 1 )(p(ρ 0 ) p(ρ 1 )) ρ 0 ρ 1 A free boundary problem for two-dimensional gas dynamics equations 5
7 Formulate the problem in self-similar coordinates: ξ = x/t, η = y/t (ρu) ξ + (ρv ) η + 2ρ = 0 (U, V ) U + U + p ξ /ρ = 0 (U, V ) V + V + p η /ρ = 0 - here, PSfragU replacements := u ξ and V := v η - when linearized about a constant state U = (ρ, u, v ), the system is hyperbolic outside the circle C U : (u ξ) 2 + (v η) 2 = c 2 (ρ ) η U 0 l 1 C Ub U a S 1 U 1 ξ U 0 = (ρ 0, 0, n 0 ), U 1 = (ρ 1, 0, 0) S a : ξ = kη + χ a, l a : ξ = kη S b : ξ = kη + χ a, l b : ξ = kη ρ 0, m a, n a ), U b = (ρ 0, m a, n a ) m a, n a, χ a - functions of ρ 0, ρ 1, k C U 1 C U a U b l 2 S 2 A free boundary problem for two-dimensional gas dynamics equations 6
8 Sfrag replacements Solution in the hyperbolic part of the domain η U 0 l 1 PSfrag replacements 1.4 k C Ub U a U 1 1.2, U 1 = (ρ 1, 0, 0) χ a, l a : ξ = kη a, l b : ξ = kη = (ρ 0, m a, n a ) nctions of ρ 0, ρ 1, k C U1 C Ua S 1 U 1 U a 0.6 Ξ s U b C U1 C 0.4 S U a 2 C U Ub b U 0 = (ρ 0, 0, n 0), U 1 = (ρ 1, 0, 0) 0.2 S a : ξ = kη + χ a, l a : ξ = kη S l b : ξ = kη + χ a, l b : ξ = kη 2 U a = (ρ 0, m a, n a), U b = (ρ 0, m a, n a) m a, n a, χ a - functions of ρ 0, ρ 1, k ξ U 0ξ η S 1 1 S 2 l1 0.8 l 2 Ξ s Ρ 0 Ρ 1 If k (k C (ρ 0, ρ 1 ), ), then: 1) Ξ s is outside the sonic circles C U 1, C U a, C U b 2) quasi-1-d Riemann problem at Ξ s has a solution (shock shock) - two solutions: U F = (ρ F, u F, 0) and U R = (ρ R, u R, 0) Ξ s is subsonic w.r.t. U F when k (k C, ) Ξ s is subsonic w.r.t. U R when k (k C, k ) Ξ s is supersonic w.r.t. U R when k (k, ) A free boundary problem for two-dimensional gas dynamics equations 7
9 Two types of regular reflection η ρ 0 η ρ 0 ements ρ ρ 1 ρ ρ 1 Ξ s ξ Ξ s ξ C U F C U R Case 1: Ξ s is subsonic w.r.t. the solution at Ξ s strong regular reflection Case 2: Ξ s is supersonic w.r.t. the solution at Ξ s weak regular reflection A free boundary problem for two-dimensional gas dynamics equations 8
10 Numerical simulations for the adiabatic gas dynamics equations using CLAWPACK software transonic (strong) and supersonic (weak) regular reflection 1 q(1) at time q(1) at time y 0.5 y x x A free boundary problem for two-dimensional gas dynamics equations 9
11 S b l a Free boundary lb problem describing transonic regular reflection η reduced system in Ω Ω σ Ω Σ 0 Σ U a Ξ s U 1 ξ jump conditions along Σ symmetry conditions along Σ 0 Dirichlet condition for ρ on σ Dirichlet conditions at Ξ s Theorem: Main idea. There exists a local solution ρ H ( γ) 1+α, U, V H( γ) 1+ɛ and Σ H 1+α of the above free boundary problem. (0 < ɛ < α < γ < 1) j 2nd order equation for ρ in Ω reduced system in Ω two 1st order transport equations for U and V in Ω 8 < oblique derivative boundary condition for ρ along Σ jump conditions on Σ Dirichlet conditions for U and V along Σ : shock evolution equation dη/dξ = Ψ(ρ, ξ, η) A free boundary problem for two-dimensional gas dynamics equations 10
12 Holder spaces Let S R 2, u : S R. Define u 0;S := sup u(x) [u] α;s := sup x y u(x) u(y) x y α u α;s := u 0;S + [u] α;s supermum norm α-holder seminorm α-holder norm u k+α;s := P k j=0 Dj u 0;S + [D k u] α;s (k + α)-holder norm Let T S and δ > 0. Let S δ;t := {x S : dist(x, T ) > δ}, and for a > 0 and b such that a b 0, define u ( b) a;s\t := sup δ a b u a;sδ;t δ>0 weighted interior Holder norm A free boundary problem for two-dimensional gas dynamics equations 11
13 Theorem: There exists a local solution ρ H ( γ) of the free boundary problem 1+α, U, V H( γ) 1+ɛ and Σ H 1+α a ij (ρ, U, V )D ij ρ + b i (ρ, U, V )D i ρ + c ij (ρ, U, V )D i ρd j ρ = 0 (U, V ) U + U + p ξ /ρ = 0 (U, V ) V + V + p η /ρ = 0 β(ρ, U, V ) ρ = F (ρ, U, V ) U = G(ρ, ξ, η) V = H(ρ, ξ, η) dη dξ = Ψ(ρ, ξ, η) 9 >= >; ρ η = U η = V = 0 on Σ 0 ρ = f on σ U(Ξ s ) = U F, η(ξ s ) = 0 Remark. We introduce cut-off functions to ensure that: on Σ : η = η(ξ) the 2nd order equation for ρ is uniformly elliptic in Ω the condition for ρ on Σ is uniformly oblique the shock evolution equation is well-defined 9 = ; in Ω A free boundary problem for two-dimensional gas dynamics equations 12
14 Step 1. given Σ m K H 1+α, find ρ m, U m, V m solving the fixed BP in Ω m ff a ij (ρ, U, V )D ij ρ + b i (ρ, U, V )D i ρ + c ij (ρ, U, V )D i ρd j ρ = 0 (U, V ) U + U + p ξ /ρ = 0, (U, V ) V + V + p η /ρ = 0 in Ω β(ρ, U, V ) ρ = F (ρ, U, V ) on Σ, ρ η = 0 on Σ 0, ρ σ = f, ρ Ξs = ρ s U Σ = G(ρ), U Ξs = u F ξ s, U η = 0 on Σ 0 V Σ = H(ρ), V Ξs = 0, V = 0 on Σ 0 fix ω, W, Z H ( γ) 1+ɛ linearize the second order problem for density using ω, W, Z in the coefficients show that there exists a solution ρ H ( γ) 1+α show that the map ω ρ has a fixed point ρ[w, Z] of this linear problem linearize the first order problem for pseudo-velocities using W, Z, ρ[w, Z] show that there exists a solution U, V H ( γ) 1+ɛ of this linear problem show that the map (W, Z) (U, V ) is a contraction Step 2. find Σ m+1 using the shock evolution equation dη/dξ = Ψ(ρ m, ξ, η) show that the map Σ m Σ m+1 has a fixed point Σ A free boundary problem for two-dimensional gas dynamics equations 13
15 RELATED WORK steady transonic small disturbance equation shock perturbation: Čanić-Keyfitz-Lieberman unsteady transonic small disturbance equation strong regular reflection: Čanić-Keyfitz-Kim weak regular reflection: Čanić-Keyfitz-Kim Guderley reflection: Tesdall-Hunter nonlinear wave system Mach reflection: Čanić-Keyfitz-Kim, Sever strong regular reflection: Jegdić-Keyfitz-Čanić weak regular reflection: Jegdić Guderley reflection: Tesdall-Sanders-Keyfitz pressure-gradient system weak regular reflection: Y.Zheng-D.Wang Euler equations for potential flow weak regular reflection: G.Q.Chen-M.Feldman adiabatic/isentropic gas dynamics equations T.Chang-G.Q.Chen, S.X.Chen et al., T.Zhang-Y.Zheng, T.P.Liu-V.Elling, A.Tesdall-R.Sanders A free boundary problem for two-dimensional gas dynamics equations 14
16 REFERENCES S. Čanić, B. L. Keyfitz, G. M. Lieberman, A proof of existence of perturbed steady transonic shocks via a free boundary problem, Communications on Pure and Applied Mathematics, LIII (2000), S. Čanić, B. L. Keyfitz, E. H. Kim, Free boundary problems for the unsteady transonic small disturbance equation: transonic regular reflection, Methods and Applications of Analysis, 7 (2000), S. Čanić, B. L. Keyfitz, E. H. Kim, A free boundary problem for a quasi-linear degenerate elliptic equation: regular reflection of weak shocks, Communications on Pure and Applied Mathematics, LV (2002), S. Čanić, B. L. Keyfitz, E. H. Kim, Free boundary problems for nonlinear wave systems: Mach stems for interacting shocks, SIAM Journal on Mathematical Analysis, 37 (2006), K. Jegdić, B. L. Keyfitz, S. Čanić, Transonic regular reflection for the nonlinear wave system, Journal of Hyperbolic Differential Equations, 3 (2006), K. Jegdić, Weak regular reflection for the nonlinear wave system, Journal of Hyperbolic Differential Equations, 5 (2008), K. Jegdić, B. L. Keyfitz, S. Čanić, A free boundary problem for the isentropic gas dynamics equations transonic regular reflection, submitted. A free boundary problem for two-dimensional gas dynamics equations 15
A 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 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 informationHyperbolic 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 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 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 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 informationSELF-SIMILAR SOLUTIONS FOR THE TRIPLE POINT PARADOX IN GASDYNAMICS
SIAM J. APPL. MATH. Vol. 0, No. 0, pp. 000 000 c XXXX Society for Industrial and Applied Mathematics SELF-SIMILAR SOLUTIONS FOR THE TRIPLE POINT PARADOX IN GASDYNAMICS ALLEN M. TESDALL, RICHARD SANDERS,
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 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 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 informationShock 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 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 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 informationSimple waves and characteristic decompositions of quasilinear hyperbolic systems in two independent variables
s and characteristic decompositions of quasilinear hyperbolic systems in two independent variables Wancheng Sheng Department of Mathematics, Shanghai University (Joint with Yanbo Hu) Joint Workshop on
More information03ER25575, and the National Science Foundation, grant a Research supported by the Department of Energy, grant DE-FG02-
The XXIst Annual Courant Lectures, Courant Institute, NYU, November, 3 1 What Studying Quasi-Steady Problems Can Tell Us About Steady Transonic Flow November, 3 Barbara Lee Keyfitz a Department of Mathematics
More informationc 2002 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 63, No., pp. 4 6 c 00 Society for Industrial and Applied Mathematics SELF-SIMILAR SOLUTIONS FOR WEAK SHOCK REFLECTION ALLEN M. TESDALL AND JOHN K. HUNTER Abstract. We present numerical
More informationChapter 6 WEAK SHOCK REFLECTION. 1. Introduction. John K. Hunter. Allen M. Tesdall
Chapter 6 WEAK SHOCK REFLECTION John K. Hunter University of California at Davis Davis, CA 95616, USA jkhunter@ucdavis.edu Allen M. Tesdall University of California at Davis Davis, CA 95616, USA amtesdal@ucdavis.edu
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 informationCHAPTER 1 The Compressible Euler System in Two Space Dimensions. Introduction
1 CHAPTER 1 The Compressible Euler System in Two Space Dimensions Yuxi Zheng Department of Mathematics The Pennsylvania State University E-mail: yzheng@math.psu.edu Abstract We analyze the self-similar
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 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 informationTwo-Dimensional Riemann Solver for Euler Equations of Gas Dynamics
Journal of Computational Physics 167, 177 195 (2001) doi:10.1006/jcph.2000.6666, available online at http://www.idealibrary.com on Two-Dimensional Riemann Solver for Euler Equations of Gas Dynamics M.
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 informationTwo-Dimensional Regular Shock Reflection for the Pressure Gradient System of Conservation Laws
Acta Mathematicae Applicatae Sinica, English Series Vol. 22, No. 2 (2006) 1 34 Two-Dimensional Regular Shock Reflection for the Pressure Gradient System of Conservation Laws Yuxi Zheng Department of Mathematics,
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 informationWorkshop on Multi-Dimensional Euler Equations and Conservation Laws. Department of Mathematics, University of Pittsburgh. November 6-9, 2003.
Workshop on Multi-Dimensional Euler Equations and Conservation Laws Department of Mathematics, University of Pittsburgh November 6-9, 2003 Schedule All talks are in Room 704, Thackeray Hall, 139 University
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 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 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 informationNonlinear system of mixed type and its application to steady Euler-Poisson system
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,
More informationApplications 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 informationTWO-DIMENSIONAL RIEMANN PROBLEMS: FROM SCALAR CONSERVATION LAWS TO COMPRESSIBLE EULER EQUATIONS
Acta Mathematica Scientia 2009,29B(4):777 802 http://actams.wipm.ac.cn TWO-DIMENSIONAL RIEMANN PROBLEMS: FROM SCALAR CONSERVATION LAWS TO COMPRESSIBLE EULER EQUATIONS Dedicated to Professor Wu Wenjun on
More informationIntroduction to Fluid Mechanics. Chapter 13 Compressible Flow. Fox, Pritchard, & McDonald
Introduction to Fluid Mechanics Chapter 13 Compressible Flow Main Topics Basic Equations for One-Dimensional Compressible Flow Isentropic Flow of an Ideal Gas Area Variation Flow in a Constant Area Duct
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationAA210A Fundamentals of Compressible Flow. Chapter 13 - Unsteady Waves in Compressible Flow
AA210A Fundamentals of Compressible Flow Chapter 13 - Unsteady Waves in Compressible Flow The Shock Tube - Wave Diagram 13.1 Equations for irrotational, homentropic, unsteady flow ρ t + x k ρ U i t (
More 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 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 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 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 informationOptimal Transportation. Nonlinear Partial Differential Equations
Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007
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 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 informationMATANIA BEN-ARTZI, JIEQUAN LI AND GERALD WARNECKE
A DIRECT EULERIAN GRP SCHEME FOR COMPRESSIBLE FLUID FLOWS MATANIA BEN-ARTZI JIEQUAN LI AND GERALD WARNECKE Abstract. A direct Eulerian generalized Riemann problem (GRP) scheme is derived for compressible
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More information1 PART1: Bratu problem
D9: Advanced Numerical Analysis: 3 Computer Assignment 3: Recursive Projection Method(RPM) Joakim Möller PART: Bratu problem. Background The Bratu problem in D is given by xx u + yy u + λe u =, u Γ = ()
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 informationREGULARITY RESULTS FOR THE EQUATION u 11 u 22 = Introduction
REGULARITY RESULTS FOR THE EQUATION u 11 u 22 = 1 CONNOR MOONEY AND OVIDIU SAVIN Abstract. We study the equation u 11 u 22 = 1 in R 2. Our results include an interior C 2 estimate, classical solvability
More informationNonlinear Wave Propagation in 1D Random Media
Nonlinear Wave Propagation in 1D Random Media J. B. Thoo, Yuba College UCI Computational and Applied Math Seminar Winter 22 Typeset by FoilTEX Background O Doherty and Anstey (Geophysical Prospecting,
More informationAlexei 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 informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationActa Mathematicae Applicatae Sinica, 19(2003), A global solution to a two-dimensional Riemann problem involving shocks as free boundaries
Acta Mathematicae Alicatae Sinica, 19(003), 559 57 A global solution to a two-dimensional Riemann roblem involving shocks as free boundaries Yuxi Zheng 1 Deartment of Mathematics The Pennsylvania State
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 informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationAeroelasticity. Lecture 9: Supersonic Aeroelasticity. G. Dimitriadis. AERO0032-1, Aeroelasticity and Experimental Aerodynamics, Lecture 9
Aeroelasticity Lecture 9: Supersonic Aeroelasticity G. Dimitriadis AERO0032-1, Aeroelasticity and Experimental Aerodynamics, Lecture 9 1 Introduction All the material presented up to now concerned incompressible
More informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationNotes #6 MAE 533, Fluid Mechanics
Notes #6 MAE 533, Fluid Mechanics S. H. Lam lam@princeton.edu http://www.princeton.edu/ lam October 1, 1998 1 Different Ways of Representing T The speed of sound, a, is formally defined as ( p/ ρ) s. It
More information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More 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 informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Martin Dindos Sukjung Hwang University of Edinburgh Satellite Conference in Harmonic Analysis Chosun University, Gwangju,
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationVariational 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 informationSTABILITY OF STEADY MULTI-WAVE CONFIGURATIONS FOR THE FULL EULER EQUATIONS OF COMPRESSIBLE FLUID FLOW. GUI-QIANG G. CHEN University of Oxford.
Report no. PDE-CDE-18/03 STABILITY OF STEADY MULTI-WAVE CONFIGURATIONS FOR THE FULL EULER EQUATIONS OF COMPRESSIBLE FLUID FLOW by GUI-QIANG G. CHEN University of Oxford and MATTHEW RIGBY University of
More informationCoupling conditions for transport problems on networks governed by conservation laws
Coupling conditions for transport problems on networks governed by conservation laws Michael Herty IPAM, LA, April 2009 (RWTH 2009) Transport Eq s on Networks 1 / 41 Outline of the Talk Scope: Boundary
More informationComputation of entropic measure-valued solutions for Euler equations
Computation of entropic measure-valued solutions for Euler equations Eitan Tadmor Center for Scientific Computation and Mathematical Modeling (CSCAMM) Department of Mathematics, Institute for Physical
More informationConservation 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 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 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 informationAnalysis on Linear Stability of Oblique Shock Waves in Steady Supersonic Flow
Analysis on Linear Stability of Oblique Shock Waves in Steady Supersonic Flow Dening Li Department of Mathematics, West Virginia University, USA Abstract An attached oblique shock wave is generated when
More informationBoundary regularity of solutions of degenerate elliptic equations without boundary conditions
Boundary regularity of solutions of elliptic without boundary Iowa State University November 15, 2011 1 2 3 4 If the linear second order elliptic operator L is non with smooth coefficients, then, for any
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 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 informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationContinuity Equation for Compressible Flow
Continuity Equation for Compressible Flow Velocity potential irrotational steady compressible Momentum (Euler) Equation for Compressible Flow Euler's equation isentropic velocity potential equation for
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationOn a simple model of isothermal phase transition
On a simple model of isothermal phase transition Nicolas Seguin Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie Paris 6 France Micro-Macro Modelling and Simulation of Liquid-Vapour Flows
More informationRegularity 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 informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 225 Estimates of the second-order derivatives for solutions to the two-phase parabolic problem
More informationGLOBAL 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 informationOblique derivative problems for elliptic and parabolic equations, Lecture II
of the for elliptic and parabolic equations, Lecture II Iowa State University July 22, 2011 of the 1 2 of the of the As a preliminary step in our further we now look at a special situation for elliptic.
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 informationGoverning Equations of Fluid Dynamics
Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present
More informationDeforming 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 informationFLUID MECHANICS 3 - LECTURE 4 ONE-DIMENSIONAL UNSTEADY GAS
FLUID MECHANICS 3 - LECTURE 4 ONE-DIMENSIONAL UNSTEADY GAS Consider an unsteady 1-dimensional ideal gas flow. We assume that this flow is spatially continuous and thermally isolated, hence, it is isentropic.
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 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 informationA 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 informationElliptic Problems for Pseudo Differential Equations in a Polyhedral Cone
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 9, Number 2, pp. 227 237 (2014) http://campus.mst.edu/adsa Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone
More informationElliptic PDEs of 2nd Order, Gilbarg and Trudinger
Elliptic PDEs of 2nd Order, Gilbarg and Trudinger Chapter 2 Laplace Equation Yung-Hsiang Huang 207.07.07. Mimic the proof for Theorem 3.. 2. Proof. I think we should assume u C 2 (Ω Γ). Let W be an open
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationSOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD. 1. Dirichlet integral
SOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD CRISTIAN E. GUTIÉRREZ FEBRUARY 3, 29. Dirichlet integral Let f C( ) with open and bounded. Let H = {u C ( ) : u = f on } and D(u) = Du(x) 2 dx.
More informationASYMPTOTIC 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 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 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 informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More information