Convergence of relaxation schemes for hyperbolic conservation laws with stiff source terms

We focus in this study on the convergence of a class of relaxation numerical schemes for hyperbolic scalar conservation laws including stiff source terms. Following Jin and Xin, we use as approximation of the scalar conservation law, a semi-linear hyperbolic system with a second stiff source term. This allows us to avoid the use of a Riemann solver in the construction of the numerical schemes. The convergence of the approximate solution toward a weak solution is established in the cases of first and second order accurate MUSCL relaxed methods.

     

On the Cauchy problem for macroscopic model of pedestrian flows

In this paper, we discuss the limit behavior of hyperbolic systems of conservation laws with stiff relaxation terms to the local systems as the relaxation time tends to zero. The prototype is crowd models derived from crowd dynamics according to macroscopic scaling when the flow of crowds is supposed to satisfy the paradigms of continuum mechanics. Under an appropriate structural stability condition, the asymptotic expansion is obtained when one assumes the existence of a smooth solution to the equilibrium system. In this case, the local existence of a classical solution is also shown.     

Hyperbolic conservation laws with stiff relaxation terms and entropy

We study the limiting behavior of systems of hyperbolic conservation laws with stiff relaxation terms. Reduced systems, inviscid and viscous local conservation laws, and weakly nonlinear limits are derived through asymptotic expansions. An entropy condition is introduced for N × N systems that ensures the hyperbolicity of the reduced inviscid system. The resulting characteristic speeds are shown to be interlaced with those of the original system. Moreover, the first correction to the reduced system is shown to be dissipative. A partial converse is proved for 2 × 2 systems. This structure is then applied to study the convergence to the reduced dynamics for the 2 × 2 case. © 1994 John Wiley & Sons, Inc.     

The relaxation schemes for systems of conservation laws in arbitrary space dimensions

We present a class of numerical schemes (called the relaxation schemes) for systems of conservation laws in several space dimensions. The idea is to use a local relaxation approximation. We construct a linear hyperbolic system with a stiff lower order term that approximates the original system with a small dissipative correction. The new system can be solved by underresolved stable numerical discretizations without using either Riemann solvers spatially or a nonlinear system of algebraic equations solvers temporally. Numerical results for 1-D and 2-D problems are presented. The second-order schemes are shown to be total variation diminishing (TVD) in the zero relaxation limit for scalar equations. ©1995 John Wiley & Sons, Inc.     

SOME RESULTS ON HYPERBOLIC SYSTEMS WITH RELAXATION

In this article, first, the authors prove that there exists a unique global smooth solution for the Cauthy problem to the hyperbolic conservation laws systems with relaxation; second, in the large time station, they prove that the global smooth solutions of the hyperbolic conservation laws systems with relaxation converge to rarefaction waves solution at a determined L^P（p ≥ 2） decay rate.     

ON THE CENTRAL RELAXING SCHEME Ⅱ: SYSTEMS OF HYPERBOLIC CONSERVATION LAWS

1. IntroductionWe are interested in construction of the central reltalng sChemes for system of noIilinearhyperbolic conservation lawswith initial data U(0, x) = Uo(x), x = (x1 ? ...! xd), based on the local relaJxation approkimationof Eq.(1.1) [2, 3, 6, 8, 9, 12].To i11ustrate the basic idea of the relaalng schemes, for the sake of simplicity in the presentation, we restrict our attention to onedimensional scalar conservaioll lawsFirst, introduce a linear hyperbollc system with a stiff sourc…     

Relaxation limit for a symmetrically hyperbolic system

In this paper, we apply the invariant region theory to get an a prioriL^∞ estimate of the relaxation approximated solutions to the Cauchy problem of a symmetrically hyperbolic system with stiff relaxation and dominant diffusion, and then obtain that the relaxation approximated solutions converge almost everywhere to the equilibrium state of the symmetrical system with the aid of the compactness framework about the scalar equation.     

A roe-type Riemann solver for hyperbolic systems with relaxation based on time-dependent wave decomposition

Summary. The paper is devoted to the construction of a higher order Roe-type numerical scheme for the solution of hyperbolic systems with relaxation source terms. It is important for applications that the numerical scheme handles both stiff and non stiff source terms with the same accuracy and computational cost and that the relaxation variables are computed accurately in the stiff case. The method is based on the solution of a Riemann problem for a linear system with constant coefficients: a study of the behavior of the solutions of both the nonlinear and linearized problems as the relaxation time tends to zero enables to choose a convenient linearization such that the numerical scheme is consistent with both the hyperbolic system when the source terms are absent and the correct relaxation system when the relaxation time tends to zero. The method is applied to the study of the propagation of sound waves in a two-phase medium. The comparison between our numerical scheme, usual fractional step methods, and numerical simulation of the relaxation system shows the necessity of using the solutions of a fully coupled hyperbolic system with relaxation terms as the basis of a numerical scheme to obtain accurate solutions regardless of the stiffness. Received October 7, 1994 / Revised version received September 27, 1995     

Singular Limits of Stiff Relaxation and Dominant Diffusion for Nonlinear Systems

We are concerned with singular limits of stiff relaxation and dominant diffusion for general 2×2 nonlinear systems of conservation laws, that is, the relaxation time τ tends to zero faster than the diffusion parameter ε, τ=o(ε), ε→0. We establish the following general framework: If there exists an a priori L^∞ bound that is uniformly with respect to ε for the solutions of a system, then the solution sequence converges to the corresponding equilibrium solution of this system. Our results indicate that the convergent behavior of such a limit is independent of either the stability criterion or the hyperbolicity of the corresponding inviscid quasilinear systems, which is not the case for other type of limits. This framework applies to some important nonlinear systems with relaxation terms, such as the system of elasticity, the system of isentropic fluid dynamics in Eulerian coordinates, and the extended models of traffic flows. The singular limits are also considered for some physical models, without L^∞ bounded estimates, including the system of isentropic fluid dynamics in Lagrangian coordinates and the models of traffic flows with stiff relaxation terms. The convergence of solutions in L^p to the equilibrium solutions of these systems is established, provided that the relaxation time τ tends to zero faster than ε.     

Relaxation WENO schemes for multidimensional hyperbolic systems of conservation laws

We present a class of high‐order weighted essentially nonoscillatory (WENO) reconstructions based on relaxation approximation of hyperbolic systems of conservation laws. The main advantage of combining the WENO schemes with relaxation approximation is the fact that the presented schemes avoid solution of the Riemann problems due to the relaxation approach and high‐resolution is obtained by applying the WENO approach. The emphasis is on a fifth‐order scheme and its performance for solving a wide class of systems of conservation laws. To show the effectiveness of these methods, we present numerical results for different test problems on multidimensional hyperbolic systems of conservation laws. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Zero relaxation and dissipation limits for hyperbolic conservation laws

We are interested in hyperbolic systems of conservation laws with relaxation and dissipation, particularly the zero relaxation limit. Such a limit is of interest in several physical situations, including gas flow near thermo-equilibrium, kinetic theory with small mean free path, and viscoelasticity with vanishing memory. In this article we study hyperbolic systems of two conservation laws with relaxation. For the stable case where the equilibrium speed is subcharacteristic with respect to the frozen speeds, we illustrate for a model in viscoelasticity that no oscillation develops for the nonlinear system in the zero relaxation limit. For the marginally stable case where the equilibrium speed may equal one of the frozen speeds, we show for a model in phase transitions that no oscillation arises when the dissipation is present and goes to zero more slowly than the relaxation. Our analysis includes the construction of suitable entropy pairs to derive energy estimates. We need such energy estimates not only for the compactness properties but also for the deviation from the equilibrium of the solutions for the relaxation systems. The theory of compensated compactness is then applied to study the oscillation in the zero relaxation limit. © 1993 John Wiley & Sons, Inc.     

Generalization of a Relaxation Scheme for Systems of Forced Nonlinear Hyperbolic Conservation Laws with Spatially Dependent Flux Functions

A generalization of a finite difference method for calculating numerical solutions to systems of nonlinear hyperbolic conservation laws in one spatial variable is investigated. A previously developed numerical technique called the relaxation method is modified from its initial application to solve initial value problems for systems of nonlinear hyperbolic conservation laws. The relaxation method is generalized in three ways herein to include problems involving any combination of the following factors: systems of nonlinear hyperbolic conservation laws with spatially dependent flux functions, nonzero forcing terms, and correctly posed boundary values. An initial value problem for the forced inviscid Burgers' equation is used as an example to show excellent agreement between theoretical solutions and numerical calculations. An initial boundary value problem consisting of a system of four partial differential equations based on the two-layer shallow-water equations is solved numerically to display a more general applicability of the method than was previously known.     

Critical thresholds in hyperbolic relaxation systems

Critical threshold phenomena in one-dimensional 2×2 quasi-linear hyperbolic relaxation systems are investigated. Assuming both the subcharacteristic condition and genuine nonlinearity of the flux, we prove global in time regularity and finite-time singularity formation of solutions simultaneously by showing the critical threshold phenomena associated with the underlying relaxation systems. Our results apply to the well-known isentropic Euler system with damping. Within the same framework it is also shown that the solution of the semi-linear relaxation system remains smooth for all time, provided the subcharacteristic condition is satisfied.     

The third-order relaxation schemes for hyperbolic conservation laws

A class of semi-discrete third-order relaxation schemes are presented for relaxation systems which approximate systems of hyperbolic conservation laws. These schemes for the scalar conservation law are shown to satisfy the property of total variation diminishing (TVD) in the zero relaxation limit. A third-order TVD Runge–Kutta splitting method is developed for the temporal discretization of the semi-discrete schemes. Numerical results are given illustrating these schemes on one-dimensional nonlinear problems.     

Asymptotic stability of relaxation shock profiles for hyperbolic conservation laws

This paper studies the asymptotic stability of traveling relaxation shock profiles for hyperbolic systems of conservation laws. Under a stability condition of subcharacteristic type the large time relaxation dynamics on the level of shocks is shown to be determined by the equilibrium conservation laws. The proof is due to the energy principle, using the weighted norms, the interaction of waves from various modes is treated by imposing suitable weight matrix.     

MULTIDIMENSIONAL RELAXATION APPROXIMATIONS FOR HYPERBOLIC SYSTEMS OF CONSERVATION LAWS

We construct and implement a non-oscillatory relaxation scheme for multidimensional hyperbolic systems of conservation laws. The method transforms the nonlinear hyperbolic system to a semilinear model with a relaxation source term and linear characteristics which can be solved numerically without using either Riemann solver or linear iterations. To discretize the relaxation system we consider a high-resolution reconstruction in space and a TVD Runge-Kutta time integration. Detailed formulation of the scheme is given for problems in three space dimensions and numerical experiments are implemented in both scalar and system cases to show the effectiveness of the method.     

Convergence to equilibrium for the relaxation approximations of conservation laws

We study the Cauchy problem for 2 × 2 semilinear and quasilinear hyperbolic systems with a singular relaxation term. Special comparison and compactness properties are established by assuming the subcharacteristic condition. Therefore we can prove the convergence to equilibrium of the solutions of these problems as the singular perturbation parameter tends to 0. This research was strongly motivated by the recent numerical investigations of S. Jin and Z. Xin on the relaxation schemes for conservation laws. © 1996 John Wiley & Sons, Inc.     

Stability of contact discontinuity for Jin-Xin relaxation system

In this paper, we first construct a contact wave for 1-dimensional Jin-Xin relaxation system [S. Jin, Z.P. Xin, The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Comm. Pure Appl. Math. 48 (1995) 235-277]. This wave serves as the relaxation version of contact discontinuity of the corresponding hyperbolic system at equilibrium. Such a contact wave is shown to be nonlinearly stable under small initial perturbation. The time-decay rate is also obtained by weighted energy estimates.     

Generalized Glimm scheme to the initial boundary value problem of hyperbolic systems of balance laws

In this paper we provide a generalized version of the Glimm scheme to establish the global existence of weak solutions to the initial-boundary value problem of 2×2 hyperbolic systems of conservation laws with source terms. We extend the methods in [J.B. Goodman, Initial boundary value problem for hyperbolic systems of conservation laws, Ph.D. Dissertation, Stanford University, 1982; J.M. Hong, An extension of Glimm’s method to inhomogeneous strictly hyperbolic systems of conservation laws by “weaker than weak” solutions of the Riemann problem, J. Differential Equations 222 (2006) 515-549] to construct the approximate solutions of Riemann and boundary Riemann problems, which can be adopted as the building block of approximate solutions for our initial-boundary value problem. By extending the results in [J. Glimm, Solutions in the large for nonlinear hyperbolic systems of equations, Comm. Pure Appl. Math. 18 (1965) 697-715] and showing the weak convergence of residuals, we obtain stability and consistency of the scheme.     

A Hyperbolic Heat Transfer Problem with Phase Changes

There is a large literature on hyperbolic heat-transfer problems,and a few investigators have studied phase-change or melt problemsfor materials governed by such systems. Most formulations ofmelt problems for hyperbolic models insist on continuity ofthe temperature across the melt interface, and a number of investigatorshave observed that this insistence leads to mathematical difficulties.In this paper an alternative model is explored, where continuityof the temperature at the melt interface is not imposed. Instead,we insist that the relaxation process describing the relationbetween heat flux and temperature gradient be interpreted asa conservation equation that must hold across a melt interface.With this formulation, we are led to a solvable problem withdesirable asymptotic properties.