Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas
|
|
- Joella Dixon
- 5 years ago
- Views:
Transcription
1 Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas Pedro González-Casanova Henríquez Unidad de Investigación en Cómputo Aplicado DGSCA, UNAM Elvira Torondel y Ricard Garrido Departamento de Matemática Aplicada, Universidad de Valencia, España Seminarios de Modelación Matemática y Computacional Instituto de Geofísica, junio 2008
2 Contents: Nonlinear Hyperbolic Systems of conservation laws. Restricted relativity: Conservative representation. Hybrid Compact-WENO method. Numerical results. Conclusions.
3 Hyperbolic Systems of Conservation Laws: Basics 1D Let u = (u 1,..., u p ) and x IR u t + f x = 0 Integral representation: udx + t D D The Jacobian A matrix of f is A(u) = f u (= s(u)) fds = 0
4 Quasi-linear representation -u smooth enough- is given by u t + A(u) u x = 0 The system is hyperbolic if the Jacobian A has d real eigenvalues λ 1 (u)... λ p (u) and p linearly independent eigenvectors If the eigenvalues are all different, the system is called strictly hyperbolic. Initial value problem: u(x, t = 0) = u 0 (x)
5 Hyperbolic Systems of Conservation Laws: Characteristic Form Consider a complete set of left eigenvectors of A, l j, with i components l i j, then l T j A = λ jl T j T, denotes transpose. Hyperbolic systems can be written in characteristic form as l T j u t + λ jl T j f x = 0 Can we uncouple the system of equations?: nonlinear case the eigenvalues and eigenvectors depend on the solution u.
6 Hyperbolic Systems of Conservation Laws: Linear Case Suppose A is a constant matrix, i.e.. eigenvectors and eigenvalues constant. Then l T j u t + λ j l T j f x = 0 this is an uncoupled system of equations: 1D hyperbolic scalar conservation laws. Initial conditions: lj T u(x, t = 0) = lt j u 0(x) for j = 1,..., d. We can now solve d scalar problems in the characteristic space, and then, using the right eigenvectors go back to the physical space.
7 In the non linear case the Jacobian and thus the eigenvalues and eigenvectors depends on u. Solution: frozen the Jacobian at some u. In the numerical algorithm, this implies that we can locally solve the problem in the characteristic space.
8 Hyperbolic Systems of Conservation Laws: Back to reality, IR d Let u = (u 1,..., u p ) and x = (x 1,..., x d ) IR d and f j (u) = (f 1j, f 2j,,..., f pj ) is the vector of fluxes u t + d j=1 f j x j = 0 (= s(u)) Integral representation: t D d udx + j=1 D f j (u) nds = 0 For all j = 1,..., d the Jacobian A j matrix of f j is A j (u) = f j u
9 Summary: We have one Jacobian for each spatial direcction. For each Jacobian apply the characteristic algorithm: Use left eigenvalues to uncouple the system. Solve the 1D problem in each characteristic direction. Use the right eigenvalues to return to the physical space. Joint the solutions corresponding to each spatial dimension.
10 Special Relativistic Fluid Hydrodynamics Relativistic simulations are more difficult than Newtonian simulations (Norman and Winkler 1986) The characteristic wave speeds in the relativistic case not only depends on the fluid velocity components in the wave propagation direction, but also depends on the normal velocity components. This coupling complicates the solution of the Riemann problem severely, (Donat et al. 1998).
11 Relativistic Fluid Dynamics: Local conservation laws Restricted relativity: Relativistic flux local conservation laws: Local conservation of rest mass, µ (ρu µ ) = 0 Local conservation of stress-enery of a fluid, µ T µν = 0 ρ is the rest mass density measured in the fluid frame, U µ = W (c, u) the fluid 4-velocity vector, W is the Lorentz factor, c the speed of light, u is the classical three-velocity.
12 Also, T µν is the energy-momentum tensor, which for a perfect fluid can be written as T µν = ρhu µ u ν + pg µν where p is the pressure and h = 1 + ɛ + p ρ is the specific enthalpy, with ɛ being the specific internal energy. The tensor g µν defines the metric of the space-time where the fluid evolves.
13 Relativistic Fluid Dynamics: System of conservation Laws The first step is to rewrite the local conservation laws in conservative form, (see Martí and Müller, 1999): Consider the following variables: D = ρw S j = ρhw 2 v j τ = ρhw 2 p ρw which are, respectively, the rest mass, momentum and total energy densities, measured in the laboratory frame. Where v i are the components of the threevelocity of the fluid and W is the Lorentz factor.
14 Definining now the conserved quantities u = (D, S j, τ ) T and the fluxes f i (u) = (Dv i, S j v i + pδ ij, S i Dv i) T where, δ ij is the Kroneker symbol and v i is the velocity. We obtain the desired conservative representation u t + f i = 0 x i i This system, which is closed by an equation of state: p = p(ρ, ɛ), is hyperbolic.
15 Thus using the spectral decomposition of the three 5 5 Jacobian matices. (Donat, et al.1998) B i = f i (u) u numerical methods can be applied.
16 Compact-WENO schemes Compact schemes very accurate in smooth regions with spectral-like resolution (see Lele 1992), however produce non-physical oscillations (Gibbs phenomena) when applied to flow with discontinuities. WENO schemes are less accurate in smooth regions but are good to approximate discontinuous flows. Hybrid schemes: combines the advantages of compact schemes in smooth regions with the sharp WENO technique near the discontinuities.
17 The finite difference equation: 1D Problem Scalar hyperbolic conservation law: u t + f x = 0 where f = f(u) Semi-discrete conservative finite difference scheme u j t + 1 ( ) ˆf h j+1/2 ˆf j 1/2 = 0 ˆf j+1/2 is a numerical flux function 1 ( ˆfj+1/2 ˆf ) ( ) f h j 1/2 = x j + O(h k )
18 The time integration: three-stage, TVD Runge- Kutta scheme. Let L j (u) = 1 h ( ˆfj+1/2 ˆf j 1/2 ) the TVD Runge-Kutta scheme is given by u (1) j = u n j + tl j(u n ) (1) u (2) j = 3 4 un j u(1) j tl j(u (1) ) (2) u n+1 j = 1 3 un j u(2) j tl j(u (2) ) (3)
19 Compact-WENO scheme in 1D: Numerical fluxes Step 1.- Compute the numerical flux function using fifth-order conservative upwind compact scheme proposed by Pirozzoli, φ j+1/2 ˆfj 1/2 + ˆf j+1/2 + ψ j+1/2 ˆfj+3/2 = ˆb j+1/2 Step 2.- Compute the numerical flux function of fifth-order WENO scheme ˆf W ENO j+1/2 = 2 γ=0 ω γ f γ j+1/2 Step 3.- Combine both fluxes: a wight σ j+1/2, depending on a threshold value r c, measure the smoothness of the numerical solution. If
20 the solution is smooth the compact scheme is used, otherwise the WENO method is utilized. Step 4.- The scalar Compact-WENO scheme, is given by σ j+1/2 φ j+1/2 ˆf j 1/2 + ˆf j+1/2 +σ j+1/2 ψ j+1/2 ˆf j+3/2 = ĉ j+1/2 The weight smoother indicator is σ j+1/2 = min(1, r j+1/2 r c ) the threshold value r c is chosen in an appropriated form.
21 Compact-WENO scheme for RFD 1. Upwind method. At each fixed x j+1/2, the average state u j+1/2 is computed by the simple mean: u j+1/2 = 1 2 (u j + u j+1 ) 2. Freeze the Jacobian. Compute right λ (i) j+1/2 (i = 1, 2, 3) and left eigenvectors l (i) at u j+1/2. j+1/2 (i = 1, 2, 3) 3. Characteristic projection. Compute local characteristic projected fluxes. w (i) m = l (i) j+1/2 f m, i = 1, 2, 3; m = j 1,..., j+2.
22 4. 1D compact-weno method For i = 1, 2, and 3, hybrid compact-weno scalar scheme applied to the local characteristic variables w (i) (i = 1, 2, 3) σ j+1/2 φ (i) j+1/2ŵ(i) j 1/2 + ŵ(i) j+1/2 + where σ j+1/2 ψ (i) j+1/2ŵ(i) j+3/2 = ĉ(i) j+1/2 5. For i = 1, 2, and 3, the characteristic compact- WENO is Φ j+1/2 ˆfj 1/2 + L j+1/2 ˆfj+1/2 + Ψ j+1/2 ˆfj+3/2 = ĉ j+1/2
23 where Φ j+1/2 = Ψ j+1/2 = L j+1/2 = σ (1) j+1/2 φ(1) j+1/2 l(1) j+1/2 σ (2) j+1/2 φ(2) j+1/2 l(2) j+1/2 σ (3) j+1/2 φ(3) j+1/2 l(3) j+1/2 l (1) j+1/2 l (2) j+1/2 l (3) j+1/2 σ (1) j+1/2 ψ(1) j+1/2 l(1) j+1/2 σ (2) j+1/2 ψ(2) j+1/2 l(2) j+1/2 σ (3) j+1/2 ψ(3) j+1/2 l(3) j+1/2
24 and ĉ j+1/2 = ĉ (1) j+1/2 ĉ (2) j+1/2 ĉ (3) j+1/2
25 Numerical test: Relativistic shock tube problem, 1D 1. Hot dense fluid on one side is separated by a membrane from a cool rarefied gas on the other. 2. The membrane is removed and the high density fluid pushes into the low density fluid.
26 Numerical test: Relativistic shock tube problem, 1D Sketch of the shock tube solution: 1. (1) the undisturbed high density fluid; 2. (2) the rarefraction wave; 3. (3) a region of constant velocity and pressure which features a contact discontinuity separating regions of different density; 4. (4) the shock itself; 5. (5) the undisturbed low density fluid.
27
28 Numerical test: 1D, problem, No. 1 Riemann Problem: The initial condition for the first problem is: {p L = 13.3, ρ L = 10, v L = 0} x [0, 0.5[ {p R = 0, ρ R = 1, v R = 0} x [0.5, 1] The results are showed in the following tables.
29 Numerical test: 1D, Problem 1 Schemes: Hybrid and WENO for problem 1 in 1D with N = 625 points, T = 0.4, t x = 0.8 HYB r c CPU 2 Loc. err. % WENO WEN
30
31 Numerical test: 1D, Problem 2 Riemann Problem: The second problem initial condition is: {D(x) = 10, S(x) = 0, τ (x) = 19.95} L x [0, 0.5[ {D(x) = sin(20πx), S(x) = 0, τ (x) = } R x [0.5, 1] Turbulent right side conditions.
32 Numerical test: 1D, Problem 2 Schemes: Hybrid and WENO for problem 2 in 1D with N = 625 points, T = 0.4, t x = 0.8 HYB r c CPU 2 Loc. err. % WENO WEN
33
34 Numerical test: 2D Problem. A square [0, 1] 2, initially there is a uniformly distributed gas whose density is everywhere. At t = 0, an inflow flux is injected into the square from the region 0.4 < x < 0.6, y = 0, while other boundary conditions are set to be absorbent.
35 Schemes for problem in 2D: Hybrid ( x t = 0.5), WENO ( x t = 0.5), N N = points, T = 0.8, v y 0 = 0.9c HYBRID r c CPU ERROR % WENO WENO
36 Comparative density ρ between hybrid compact- WENO (λ = 0.5, r c = 0.3) and WENO (λ = 0.5) for the 2D problem with points
37 Schemes for problem 1 in 2D: Hybrid ( x t = 0.5), WENO ( x t = 0.5), N N = points, T = 0.8, v y 0 = 0.9c HYBRID r c CPU-Time (s) ERROR % WENO WENO
38 Comparative density ρ between hybrid compact- WENO (λ = 0.5, r c = 0.3) and WENO (λ = 0.5) for the 2D problem with points.
39 Conclusions Hybrid compact-weno method has been formulated which accurately solve relativistic flow problems characterized by having strong shocks and turbulence.
40 Conclusions Hybrid compact-weno method has been formulated which accurately solve relativistic flow problems characterized by having strong shocks and turbulence. Numerical experiments in 1D and 2D for two methods applied to relativistic problems were performed.
41 Conclusions Hybrid compact-weno method has been formulated which accurately solve relativistic flow problems characterized by having strong shocks and turbulence. Numerical experiments in 1D and 2D for two methods applied to relativistic problems were performed. We show that the hybrid method gives, a better resolution in 1D and 2D, owing to its spectral like resolution property.
42 Conclusions Hybrid compact-weno method has been formulated which accurately solve relativistic flow problems characterized by having strong shocks and turbulence. Numerical experiments in 1D and 2D for two methods applied to relativistic problems were performed. We show that the hybrid method gives, a better resolution in 1D and 2D, owing to its spectral like resolution property. The 2D experiments were performed for a v y 0 = 0.99c.
43 Conclusions Hybrid compact-weno method has been formulated which accurately solve relativistic flow problems characterized by having strong shocks and turbulence. Numerical experiments in 1D and 2D for two methods applied to relativistic problems were performed. We show that the hybrid method gives, a better resolution in 1D and 2D, owing to its spectral like resolution property. The 2D experiments were performed for a v y 0 = 0.99c.
44 Further work is in progress to calculate ultrarelativistic problems for v y 0 = 0.999c.
Improvement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationA New Relativistic Hydrodynamic Code
A New Relativistic Hydrodynamic Code Dongsu Ryu, Indranil Chattopadhyay Department of Astronomy and Space Science, Chungnam National University, Daejeon 305-764, Korea: ryu@canopus.cnu.ac.kr, indra@canopus.cnu.ac.kr
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 informationNumerical Hydrodynamics in Special Relativity
Numerical Hydrodynamics in Special Relativity J. M. Martí Departamento de Astronomía y Astrofísica Universidad de Valencia 46100 Burjassot (Valencia), Spain martij@godunov.daa.uv.es E. Müller Max-Planck-Institut
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationNumerical 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 informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
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 informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
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. 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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationarxiv: v4 [physics.flu-dyn] 23 Mar 2019
High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows Man Long Wong a, Sanjiva K. Lele a,b arxiv:7.895v4 [physics.flu-dyn] 3 Mar 9
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 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 informationSMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES
Proceedings of ALGORITMY 2016 pp. 155 164 SMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES TAMER H. M. A. KASEM AND FRANÇOIS G. SCHMITT Abstract. The weighted essentially non-oscillatory
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationFundamental equations of relativistic fluid dynamics
CHAPTER VI Fundamental equations of relativistic fluid dynamics When the energy density becomes large as may happen for instance in compact astrophysical objects, in the early Universe, or in high-energy
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
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 informationMathematical foundations of relativistic hydrodynamics. José Antonio Font
Mathematical foundations of relativistic hydrodynamics José Antonio Font University of Valencia NewCompStar School 2017, University of Sofia Outline Lecture 1: Relativistic Hydrodynamics Lecture 2: Numerical
More informationConstrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics
Project work at the Department of Mathematics, TUHH Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics Katja Baumbach April 4, 005 Supervisor:
More informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
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 informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationRunge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter
Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter Jun Zhu, inghui Zhong, Chi-Wang Shu 3 and Jianxian Qiu 4 Abstract In this paper, we propose a new type of weighted
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
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 informationMulti-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws
Multi-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws Bruno Costa a, Wai Sun Don b, a Departamento de Matemática Aplicada, IM-UFRJ, Caia Postal 6853, Rio de Janeiro, RJ, C.E.P. 945-97,
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 informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationTwo-waves PVM-WAF method for non-conservative systems
Two-waves PVM-WAF method for non-conservative systems Manuel J. Castro Díaz 1, E.D Fernández Nieto, Gladys Narbona Reina and Marc de la Asunción 1 1 Departamento de Análisis Matemático University of Málaga
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
More informationEntropy stable high order discontinuous Galerkin methods. for hyperbolic conservation laws
Entropy stable high order discontinuous Galerkin methods for hyperbolic conservation laws Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Tianheng Chen, and with Yong Liu
More informationState of the Art MHD Methods for Astrophysical Applications p.1/32
State of the Art MHD Methods for Astrophysical Applications Scott C. Noble February 25, 2004 CTA, Physics Dept., UIUC State of the Art MHD Methods for Astrophysical Applications p.1/32 Plan of Attack Is
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More informationMulti-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws
Multi-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws Bruno Costa Wai Sun Don February 5, 6 Abstract In this article we introduce the multi-domain hybrid Spectral-WENO method aimed
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationarxiv:gr-qc/ v1 18 Jun 2004
Code Development of Three-Dimensional General Relativistic Hydrodynamics with AMR(Adaptive-Mesh Refinement) and Results From Special and General Relativistic Hydrodynamic Orhan Dönmez Nigde University
More informationEulerian Vortex Motion in Two and Three Dimensions
Eulerian Vortex Motion in Two and Three Dimensions Igor Yanovsky June 2006 1 Objective An Eulerian approach is used to solve the motion of an incompressible fluid, in two and three dimensions, in which
More informationAnomalous wave structure in magnetized materials described by non-convex equations of state
This is a preprint of: Anomalous wave structure in magnetized materials described by non-convex equations of state, Susana Serna, Antonio Marquina, Phys. Fluids, vol. 6, 6, 4. DOI: [.63/.48545] Anomalous
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More informationA Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1
A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationArt Checklist. Journal Code: Article No: 6959
Art Checklist Journal Code: JCPH Article No: 6959 Disk Recd Disk Disk Usable Art # Y/N Format Y/N Remarks Fig. 1 Y PS Y Fig. 2 Y PS Y Fig. 3 Y PS Y Fig. 4 Y PS Y Fig. 5 Y PS Y Fig. 6 Y PS Y Fig. 7 Y PS
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationA Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws
J Sci Comput (011) 46: 359 378 DOI 10.1007/s10915-010-9407-9 A Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws Fei Teng Li Yuan Tao Tang Received: 3 February 010 / Revised:
More informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationImproved Seventh-Order WENO Scheme
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Eposition 4-7 January 2, Orlando, Florida AIAA 2-45 Improved Seventh-Order WENO Scheme Yiqing Shen Gecheng Zha Dept.
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationA General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations
A General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations Ronald P. Fedkiw Xu-Dong Liu Stanley Osher September, 2 Abstract Standard
More informationA method for avoiding the acoustic time step restriction in compressible flow
A method for avoiding the acoustic time step restriction in compressible flow Nipun Kwatra Jonathan Su Jón T. Grétarsson Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationA re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws
A re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws Chieh-Sen Huang a,, Todd Arbogast b,, Chen-Hui Hung c,3 a Department of Applied Mathematics and National
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationFinite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime
Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math 2017 41: 1027 1041 c TÜBİTAK doi:10.906/mat-1602-8 Finite volume approximation of the relativistic
More informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationImproved constrained scheme for the Einstein equations: An approach to the uniqueness issue
Improved constrained scheme for the Einstein equations: An approach to the uniqueness issue Jérôme Novak (Jerome.Novak@obspm.fr) Laboratoire Univers et Théories (LUTH) CNRS / Observatoire de Paris / Université
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
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 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 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 informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationPROBLEM SET. Heliophysics Summer School. July, 2013
PROBLEM SET Heliophysics Summer School July, 2013 Problem Set for Shocks and Particle Acceleration There is probably only time to attempt one or two of these questions. In the tutorial session discussion
More informationHybrid DG/FV schemes for magnetohydrodynamics and relativistic hydrodynamics
Hybrid DG/FV schemes for magnetohydrodynamics and relativistic hydrodynamics Jonatan Núñez-De La Rosa, Claus-Dieter Munz To cite this version: Jonatan Núñez-De La Rosa, Claus-Dieter Munz. Hybrid DG/FV
More informationNumber-Flux Vector and Stress-Energy Tensor
Massachusetts Institute of Technology Department of Physics Physics 8.962 Spring 2002 Number-Flux Vector and Stress-Energy Tensor c 2000, 2002 Edmund Bertschinger. All rights reserved. 1 Introduction These
More informationSolving Einstein s Equations: PDE Issues
Solving Einstein s Equations: PDE Issues Lee Lindblom Theoretical Astrophysics, Caltech Mathematical and Numerical General Relativity Seminar University of California at San Diego 22 September 2011 Lee
More informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
More informationGravitation: Tensor Calculus
An Introduction to General Relativity Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Notes based on textbook: Spacetime and Geometry by S.M. Carroll Spring 2013
More informationGoing with the flow: A study of Lagrangian derivatives
1 Going with the flow: A study of Lagrangian derivatives Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc/ 12 February
More informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
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 informationExtremum-Preserving Limiters for MUSCL and PPM
arxiv:0903.400v [physics.comp-ph] 7 Mar 009 Extremum-Preserving Limiters for MUSCL and PPM Michael Sekora Program in Applied and Computational Mathematics, Princeton University Princeton, NJ 08540, USA
More informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
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 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 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 informationStable Semi-Discrete Schemes for the 2D Incompressible Euler Equations
Stable Semi-Discrete Schemes for the D Incompressible Euler Equations Doron Levy Department of Mathematics Stanford University dlevy@math.stanford.edu http://math.stanford.edu/ dlevy Maryland, May 004
More informationNumerical Methods for Modern Traffic Flow Models. Alexander Kurganov
Numerical Methods for Modern Traffic Flow Models Alexander Kurganov Tulane University Mathematics Department www.math.tulane.edu/ kurganov joint work with Pierre Degond, Université Paul Sabatier, Toulouse
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationS hw 2 v and hw 2 p D are the relativistic densities
THE ASTROPHYSICAL JOURNAL, 448 : L105 L108, 1995 August 1 1995. The American Astronomical Society. All rights reserved. Printed in U.S.A. MORPHOLOGY AND DYNAMICS OF HIGHLY SUPERSONIC RELATIVISTIC JETS
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 informationHIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF STATES FOR SOLVING HYPERBOLIC SYSTEMS WITH NONCONSERVATIVE PRODUCTS
MATHEMATICS OF COMPUTATION Volume 75, Number 255, July 26, Pages 113 1134 S 25-5718(6)1851-5 Article electronically published on March 21, 26 HIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF
More information