## Numerical Schemes For Hyperbolic Conservation Laws With Stiff Relaxation Terms (1996)

Venue: | J. Comput. Phys |

Citations: | 62 - 11 self |

### BibTeX

@ARTICLE{Jin96numericalschemes,

author = {Shi Jin and C. David Levermore},

title = {Numerical Schemes For Hyperbolic Conservation Laws With Stiff Relaxation Terms},

journal = {J. Comput. Phys},

year = {1996},

volume = {126},

pages = {449--467}

}

### Years of Citing Articles

### OpenURL

### Abstract

Hyperbolic systems often have relaxation terms that give them a partially conservative form and that lead to a long-time behavior governed by reduced systems that are parabolic in nature. In this article it is shown by asymptotic analysis and numerical examples that semidiscrete high resolution methods for hyperbolic conservation laws fail to capture this asymptotic behavior unless the small relaxation rate is resolved by a fine spatial grid. We introduce a modification of higher order Godunov methods that possesses the correct asymptotic behavior, allowing the use of coarse grids (large cell Peclet numbers). The idea is to build into the numerical scheme the asymptotic balances that lead to this behavior. Numerical experiments on 2 \Theta 2 systems verify our analysis. 1 E-mail address: jin@math.gatech.edu 2 E-mail address: lvrmr@math.arizona.edu Typeset by A M S-T E X 2 1. Introduction Hyperbolic systems of partial differential equations that arise in applications ofter have re...

### Citations

576 |
Towards the ultimate conservative difference scheme, V. A second order sequel to Godunov’s method
- Leer
- 1979
(Show Context)
Citation Context ... specifically addressing the issue of stiff relaxation. In the absence of lower order relaxation terms, these systems are just systems of conservation laws and one could employ a higher order Godunov =-=[44,10,31,21]-=- type scheme to evaluate the fluxes. These schemes properly capture the discontinuous features of the solutions. However, when the relaxation becomes stiff, it is computationally impractical to spatia... |

482 | Shock Waves and Reaction-Diffusion Equations - Smoller - 1994 |

233 | Non-oscillatory central differencing for hyperbolic conservation laws - Nessyahu, Tadmor - 1990 |

143 | Hyperbolic conservation laws with stiff relaxation terms and entropy
- Chen, Levermore, et al.
- 1994
(Show Context)
Citation Context ...atisfied then a simple linear analysis shows that some perturbations can grow at rates on the order of 1=�� . The relation between solutions of system (1.3) and those of equation (1.4) was studied=-= in [6,7,36]. A crit-=-ical nondimensional parameter ffl is the ratio of the distance a typical sound wave travels (at speed c = p p 0 (h)) over the relaxation time �� to a typical gradient length L: ffl = c�� L : (... |

136 |
Hyperbolic conservation laws with relaxation
- Liu
- 1987
(Show Context)
Citation Context ...atisfied then a simple linear analysis shows that some perturbations can grow at rates on the order of 1=�� . The relation between solutions of system (1.3) and those of equation (1.4) was studied=-= in [6,7,36]. A crit-=-ical nondimensional parameter ffl is the ratio of the distance a typical sound wave travels (at speed c = p p 0 (h)) over the relaxation time �� to a typical gradient length L: ffl = c�� L : (... |

126 | A finite difference method for the numerical computation and discontinuous solutions of the equations of fluid namics - Godunov - 1959 |

126 |
Moment closure hierarchies for kinetic theory
- Levermore
- 1996
(Show Context)
Citation Context ...46], multiphase flow [15] and turbulence modeling [13,30], where the relaxation terms describe thermal coupling, phase coupling or turbulent kinematics respectively. They also arise in kinetic theory =-=[2,5,33,34] whe-=-re the relaxation terms model the interaction of particles. One of the simplest such systems is the p-system with damping: @ t h + @ x w = 0 ; @ t w + @ x p(h) = \Gamma 1 �� w : (1.1) Here the par... |

122 |
Introduction to physical gas dynamics
- Vincenti, Kruger
- 1965
(Show Context)
Citation Context ...ation terms that give the system a partially conservative form and lead to a long-time behavior governed by a reduced system that is parabolic in nature. Such systems are found in relaxing gas theory =-=[8,46]-=-, multiphase flow [15] and turbulence modeling [13,30], where the relaxation terms describe thermal coupling, phase coupling or turbulent kinematics respectively. They also arise in kinetic theory [2,... |

96 |
Numerical Methods for Conservation Laws, Birkhäuser
- LeVeque
- 1990
(Show Context)
Citation Context ... specifically addressing the issue of stiff relaxation. In the absence of lower order relaxation terms, these systems are just systems of conservation laws and one could employ a higher order Godunov =-=[44,10,31,21]-=- type scheme to evaluate the fluxes. These schemes properly capture the discontinuous features of the solutions. However, when the relaxation becomes stiff, it is computationally impractical to spatia... |

77 |
Fluid Dynamics of Viscoelastic Liquids
- Joseph
- 1990
(Show Context)
Citation Context ...strative nondimensional linear model. When p(h) = h and �� = ffl in system (1.1), it becomes @ t h + @ x w = 0 ; @ t w + @ x h = \Gamma 1 ffl w ; (2.1) which are also known as the telegraph equati=-=ons [27]-=-. This system has characteristic speeds \Sigma1. Accordingly, its solutions exhibit an approximate parabolic behavior (1.2) governed by @ t h \Gamma ffl @ xx h = 0 ; w = \Gammaffl @ x h : (2.2) Here h... |

68 |
Uniformly high-order accurate essentially non-oscillatory schemes
- Harten, Enquist, et al.
- 1987
(Show Context)
Citation Context ... specifically addressing the issue of stiff relaxation. In the absence of lower order relaxation terms, these systems are just systems of conservation laws and one could employ a higher order Godunov =-=[44,10,31,21]-=- type scheme to evaluate the fluxes. These schemes properly capture the discontinuous features of the solutions. However, when the relaxation becomes stiff, it is computationally impractical to spatia... |

61 | Numerical Methods for Hyperbolic Conservation laws with stiff relaxation - Pember |

57 |
A study of numerical methods for hyperbolic conservation laws with sti source terms
- LeVeque, Yee
- 1990
(Show Context)
Citation Context ...ptotic regimes. Other early literature studied stiff source problems in reactive flows where an incorrect numerical shock speed was found when the reaction time scale ffl was not resolved numerically =-=[1,9,11,12,18,32]-=-. These works focused on source terms that take on a more complicated form than those considered in this article (for example, exhibiting multiple equilibria, or the characteristic speeds of the origi... |

50 | Runge-Kutta methods for hyperbolic conservation laws with stiff relaxation terms
- Jin
- 1995
(Show Context)
Citation Context ...cretization. We will only present the numerical results for H. The behavior of W is similar under the resolved time discretization. Fully-discrete numerical schemes are studied in a forthcoming paper =-=[24]. He-=-re we test the above five schemes for the linear telegraph equations (2.1), with c = 1, over the periodic interval 0sxs2 with initial condition h(0; x) = 2 + sin(��x) and w(0; x) = \Gamma0:1. We t... |

47 |
Theoretical and numerical structure for reacting shock waves
- Colella, Majda, et al.
- 1986
(Show Context)
Citation Context ...ptotic regimes. Other early literature studied stiff source problems in reactive flows where an incorrect numerical shock speed was found when the reaction time scale ffl was not resolved numerically =-=[1,9,11,12,18,32]-=-. These works focused on source terms that take on a more complicated form than those considered in this article (for example, exhibiting multiple equilibria, or the characteristic speeds of the origi... |

41 |
Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes
- Larsen, Morel
- 1989
(Show Context)
Citation Context ...rical schemes arise in linear transport theory where the density of particles undergoing collisions with a background medium will solve the diffusion equation in the limit of vanishing mean free path =-=[28,29]-=-. In this case much progress has been made in developing schemes that incorporate boundary layers and material interfaces [17,25,26]. However little attention has been given to this problem for the ca... |

36 |
Convergence to Nonlinear Diffusion Waves for Solutions of a System of Hyperbolic Conservation Laws with
- Hasio, Liu
- 1992
(Show Context)
Citation Context ...tem, which is hyperbolic provided p 0 (h) ? 0. The system can be viewed as the isentropic Euler equations in Lagrangian form with a drag term in the momentum equation. It was studied by Hsiao and Liu =-=[22] who sho-=-wed that its solutions exhibit a long-time behavior governed by @ t h \Gamma �� @ xx p(h) = 0 ; w = \Gamma�� @ x p(h) : (1.2) The first equation above is sometimes referred to as the porous me... |

36 |
On the Relation Between the Upwind-Differencing Schemes of Godunov, Enguist-Osher and Roe
- Leer
- 1985
(Show Context)
Citation Context ...ly build the proper asymptotic balances into the spatial difference scheme. We began such an investigation in [23]. In particular, in spirit close to the idea of modeling the flow as piecewise steady =-=[14,40,45]-=-, a semidiscrete numerical scheme will be said to have the correct parabolic behavior if the asymptotics that leads from the hyperbolic system to the parabolic approximation in continuous space is mim... |

31 |
Zero relaxation and dissipation limits for hyperbolic conservation laws
- Chen, Liu
- 1993
(Show Context)
Citation Context ...atisfied then a simple linear analysis shows that some perturbations can grow at rates on the order of 1=�� . The relation between solutions of system (1.3) and those of equation (1.4) was studied=-= in [6,7,36]. A crit-=-ical nondimensional parameter ffl is the ratio of the distance a typical sound wave travels (at speed c = p p 0 (h)) over the relaxation time �� to a typical gradient length L: ffl = c�� L : (... |

30 |
Quasilinear hyperbolic systems
- Liu
- 1979
(Show Context)
Citation Context ...nserved quantity v is asymptotically governed by a system of convection-diffusion equations, analogous to (1.4a). There have been many theoretical investigations of such systems in this regime (e.g., =-=[6,7,34,35]-=-). In particular, the general framework of the parabolic approximation for N \ThetaN systems in any spatial dimension was laid down in [7]. In this article we are concerned with the numerical methods ... |

29 |
The Fluid-Dynamical Limit of a Nonlinear Model Boltzmann Equation
- Caflish, Papanicolaou
- 1979
(Show Context)
Citation Context ...of the introduction. In particular, these ideas may be helpful in numerically understanding kinetic equations close to their macroscopic limit described by either the Euler or Navier-Stokes equations =-=[2,3,4,33]-=-. We have only examined semidiscrete schemes in this article; the time steps in our numerical experiments being small enough so no new errors were introduced to affect our asymptotic analysis. However... |

27 |
Upwind differencing schemes for hyperbolic conservation laws with source terms
- Roe
- 1986
(Show Context)
Citation Context ...ly build the proper asymptotic balances into the spatial difference scheme. We began such an investigation in [23]. In particular, in spirit close to the idea of modeling the flow as piecewise steady =-=[14,40,45]-=-, a semidiscrete numerical scheme will be said to have the correct parabolic behavior if the asymptotics that leads from the hyperbolic system to the parabolic approximation in continuous space is mim... |

26 |
Grounds effects on pressure fluctuations in the atmospheric boundary layer
- Gibson, Launder
- 1978
(Show Context)
Citation Context ...tive form and lead to a long-time behavior governed by a reduced system that is parabolic in nature. Such systems are found in relaxing gas theory [8,46], multiphase flow [15] and turbulence modeling =-=[13,30]-=-, where the relaxation terms describe thermal coupling, phase coupling or turbulent kinematics respectively. They also arise in kinetic theory [2,5,33,34] where the relaxation terms model the interact... |

26 | The Convergence of Numerical Transfer Schemes in Diffusive Regimes I: The Discrete-Ordinate Method
- Golse, Jin, et al.
- 1999
(Show Context)
Citation Context ... solve the diffusion equation in the limit of vanishing mean free path [28,29]. In this case much progress has been made in developing schemes that incorporate boundary layers and material interfaces =-=[17,25,26]-=-. However little attention has been given to this problem for the case of genuinely nonlinear systems; in such cases the difficulties of properly capturing the development and evolution of discontinui... |

23 |
The discrete-ordinate method in diffusive regime, Transport Theory and Statistical
- Jin, Levermore
- 1991
(Show Context)
Citation Context ... solve the diffusion equation in the limit of vanishing mean free path [28,29]. In this case much progress has been made in developing schemes that incorporate boundary layers and material interfaces =-=[17,25,26]-=-. However little attention has been given to this problem for the case of genuinely nonlinear systems; in such cases the difficulties of properly capturing the development and evolution of discontinui... |

19 |
The piecewise-parabolic method (PPM) for gas-dynamical simulations
- Colella, Woodward
- 1984
(Show Context)
Citation Context |

18 |
The generalized Riemann problem for reactive flows
- Ben-Artzi
- 1989
(Show Context)
Citation Context ...ptotic regimes. Other early literature studied stiff source problems in reactive flows where an incorrect numerical shock speed was found when the reaction time scale ffl was not resolved numerically =-=[1,9,11,12,18,32]-=-. These works focused on source terms that take on a more complicated form than those considered in this article (for example, exhibiting multiple equilibria, or the characteristic speeds of the origi... |

17 |
The asymptotic analysis of wave interactions and numerical calculations of transonic nozzle flow”, Adv
- Glaz, Liu
- 1984
(Show Context)
Citation Context ...ly build the proper asymptotic balances into the spatial difference scheme. We began such an investigation in [23]. In particular, in spirit close to the idea of modeling the flow as piecewise steady =-=[14,40,45]-=-, a semidiscrete numerical scheme will be said to have the correct parabolic behavior if the asymptotics that leads from the hyperbolic system to the parabolic approximation in continuous space is mim... |

15 |
Fully-discrete numerical transfer in diffusive regimes, Transport Theory Statist. Phys
- JIN, LEVERMORE
- 1993
(Show Context)
Citation Context ...realized by properly modifying their numerical fluxes. The main idea we utilize is to simply build the proper asymptotic balances into the spatial difference scheme. We began such an investigation in =-=[23]-=-. In particular, in spirit close to the idea of modeling the flow as piecewise steady [14,40,45], a semidiscrete numerical scheme will be said to have the correct parabolic behavior if the asymptotics... |

12 |
E.: Navier-Stokes and Boltzmann shock profiles for a model of gas dynamics
- Caflisch
- 1979
(Show Context)
Citation Context ...46], multiphase flow [15] and turbulence modeling [13,30], where the relaxation terms describe thermal coupling, phase coupling or turbulent kinematics respectively. They also arise in kinetic theory =-=[2,5,33,34] whe-=-re the relaxation terms model the interaction of particles. One of the simplest such systems is the p-system with damping: @ t h + @ x w = 0 ; @ t w + @ x p(h) = \Gamma 1 �� w : (1.1) Here the par... |

12 |
Numerical Wave Propagation in an Advection Equation with a Nonlinear Source Term
- Griffiths, Stuart, et al.
- 1992
(Show Context)
Citation Context |

10 |
The asymptotic diffusion limit of discretized transport problems
- Larsen
- 1992
(Show Context)
Citation Context ...rical schemes arise in linear transport theory where the density of particles undergoing collisions with a background medium will solve the diffusion equation in the limit of vanishing mean free path =-=[28,29]-=-. In this case much progress has been made in developing schemes that incorporate boundary layers and material interfaces [17,25,26]. However little attention has been given to this problem for the ca... |

8 |
Numerical approximation of hyperbolic conservation laws with stiff terms, in
- Engquist, Sjogreen
- 1991
(Show Context)
Citation Context |

7 |
The continuous structure of discontinuities
- Glimm
- 1989
(Show Context)
Citation Context ...e system a partially conservative form and lead to a long-time behavior governed by a reduced system that is parabolic in nature. Such systems are found in relaxing gas theory [8,46], multiphase flow =-=[15]-=- and turbulence modeling [13,30], where the relaxation terms describe thermal coupling, phase coupling or turbulent kinematics respectively. They also arise in kinetic theory [2,5,33,34] where the rel... |

7 | Characteristic-based schemes for dispersive waves. I. The method of characteristics for smooth solutions - Roe, Arora - 1993 |

5 |
The Boltzmann Equation and its Applications
- Cerecignani
- 1988
(Show Context)
Citation Context ...46], multiphase flow [15] and turbulence modeling [13,30], where the relaxation terms describe thermal coupling, phase coupling or turbulent kinematics respectively. They also arise in kinetic theory =-=[2,5,33,34] whe-=-re the relaxation terms model the interaction of particles. One of the simplest such systems is the p-system with damping: @ t h + @ x w = 0 ; @ t w + @ x p(h) = \Gamma 1 �� w : (1.1) Here the par... |

4 |
Second Moment Closure for the Near-Wall Sublayer: Development and Application
- Launder, Shima
- 1989
(Show Context)
Citation Context ...tive form and lead to a long-time behavior governed by a reduced system that is parabolic in nature. Such systems are found in relaxing gas theory [8,46], multiphase flow [15] and turbulence modeling =-=[13,30]-=-, where the relaxation terms describe thermal coupling, phase coupling or turbulent kinematics respectively. They also arise in kinetic theory [2,5,33,34] where the relaxation terms model the interact... |

3 |
A subcell resolution method for viscous systems and computations of combustion waves
- Harabetian
- 1992
(Show Context)
Citation Context ...ristic speeds of the original system and the local equilibrium system do not interlace), and the resulting long-time behavior can be quite different. Also related to our work is those of Harabetian's =-=[19,20]-=-, in which the dissipative effects of numerical schemes for viscous hyperbolic systems have been studied without specifically addressing the issue of stiff relaxation. In the absence of lower order re... |

3 |
Fluid Dynamical Limits of Discrete Kinetic Theories, in "Macroscopic Simulations of Complex Hydrodynamic Phenomena
- Levermore
- 1992
(Show Context)
Citation Context |

2 |
Uniformly accurate numerical schemes for hyperbolic systems with relaxation
- Jin, Russo
- 1997
(Show Context)
Citation Context ...of the introduction. In particular, these ideas may be helpful in numerically understanding kinetic equations close to their macroscopic limit described by either the Euler or Navier-Stokes equations =-=[2,3,4,33]-=-. We have only examined semidiscrete schemes in this article; the time steps in our numerical experiments being small enough so no new errors were introduced to affect our asymptotic analysis. However... |

1 | Gas Dynamics with Relaxation Effects, Rep. Prog. Phys - Clark - 1978 |

1 |
Robust Difference
- Engquist, Sjogreen
(Show Context)
Citation Context |

1 |
A Numerical Method for Viscous Perturbations of Hyperbolic Conservation
- Harabetian
- 1990
(Show Context)
Citation Context ...ristic speeds of the original system and the local equilibrium system do not interlace), and the resulting long-time behavior can be quite different. Also related to our work is those of Harabetian's =-=[19,20]-=-, in which the dissipative effects of numerical schemes for viscous hyperbolic systems have been studied without specifically addressing the issue of stiff relaxation. In the absence of lower order re... |

1 |
The Modified Equation Approach to the Sability and Accuracy of Finite Difference Methods
- Warming, Hyett
- 1973
(Show Context)
Citation Context ...Gamma 2H j +H j \Gamma1 2ffi + W j+1 \Gamma W j \Gamma1 2ffi = 0 ; @ t W j \Gamma W j+1 \Gamma 2W j +W j \Gamma1 2ffi + H j+1 \Gamma H j \Gamma1 2ffi = \Gamma 1 ffl W j : (2.8) The modified equations =-=[47]-=- corresponding to this scheme are @ t h + @ x w = ffi 2 @ xx h \Gamma ffi 2 6 @ xxx w ; @ t w + @ x h = \Gamma 1 ffl w + ffi 2 @ xx w \Gamma ffi 2 6 @ xxx h : These show the first order nature of the ... |