Results 1  10
of
250
Nonoscillatory Central Schemes For Multidimensional Hyperbolic Conservation Laws
 SIAM J. Sci. Comput
, 1998
"... We construct, analyze, and implement a new nonoscillatory highresolution scheme for twodimensional hyperbolic conservation laws. The scheme is a predictorcorrector method which consists of two steps: starting with given cell averages, we first predict pointvalues which are based on nonoscillatory ..."
Abstract

Cited by 133 (11 self)
 Add to MetaCart
We construct, analyze, and implement a new nonoscillatory highresolution scheme for twodimensional hyperbolic conservation laws. The scheme is a predictorcorrector method which consists of two steps: starting with given cell averages, we first predict pointvalues which are based on nonoscillatory piecewiselinear reconstructions from the given cell averages; at the second corrector step, we use staggered averaging, together with the predicted midvalues, to realize the evolution of these averages. This results in a secondorder, nonoscillatory central scheme, a natural extension of the onedimensional secondorder central scheme of Nessyahu and Tadmor [J. Comput. Phys., 87 (1990), pp. 408448]. As in the onedimensional case, the main feature of our twodimensional scheme is simplicity. In particular, this central scheme does not require the intricate and timeconsuming (approximate) Riemann solvers which are essential for the highresolution upwind schemes; in fact, even the com...
Numerical study of timesplitting spectral discretizations of nonlinear Schrödinger equations in the semiclassical regimes
 SIAM J. SCI. COMPUT
, 2003
"... In this paper we study the performance of timesplitting spectral approximations for general nonlinear Schrödinger equations (NLS) in the semiclassical regimes, where the Planck constant ε is small. The timesplitting spectral approximation under study is explicit, unconditionally stable and conse ..."
Abstract

Cited by 91 (43 self)
 Add to MetaCart
(Show Context)
In this paper we study the performance of timesplitting spectral approximations for general nonlinear Schrödinger equations (NLS) in the semiclassical regimes, where the Planck constant ε is small. The timesplitting spectral approximation under study is explicit, unconditionally stable and conserves the position density in L 1. Moreover it is timetransverse invariant and timereversible when the corresponding NLS is. Extensive numerical tests are presented for weak/strong focusing/defocusing nonlinearities, for the Gross–Pitaevskii equation, and for currentrelaxed quantum hydrodynamics. The tests are geared towards the understanding of admissible meshing strategies for obtaining “correct” physical observables in the semiclassical regimes. Furthermore, comparisons between the solutions of the NLS and its hydrodynamic semiclassical limit are presented.
Convergence to Equilibrium for the Relaxation Approximations of Conservation Laws
, 1996
"... We study the Cauchy problem for 2\Theta2 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 ..."
Abstract

Cited by 89 (13 self)
 Add to MetaCart
We study the Cauchy problem for 2\Theta2 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 zero. This research was strongly motivated by the recent numerical investigations of S. Jin and Z. Xin on the relaxation schemes for conservation laws. 1. Introduction In this paper we are interested to the relaxation behaviour of the following system of hyperbolic conservation laws with a singular perturbation source (1.1) ae @ t u + @ x v = 0 ; @ t v + @ x oe(u) = \Gamma 1 " (v \Gamma f(u)) (" ? 0); for (x; t) 2 IR \Theta (0; 1). Here oe, f are some given smooth functions such that oe 0 (u) ( ? 0), f(0) = 0. The system (1.1) is equivalent to the onedimensional perturbed wave equation (1.2) @ tt w \Gamma @ x oe(@ x...
Uniformly accurate schemes for hyperbolic systems with relaxations
 SIAM J. Numer. Anal
, 1997
"... Abstract. We develop highresolution shockcapturing numerical schemes for hyperbolic systems with relaxation. In such systems the relaxation time may vary from order1 to much less than unity. When the relaxation time is small, the relaxation term becomes very strong and highly stiff, and underreso ..."
Abstract

Cited by 75 (23 self)
 Add to MetaCart
(Show Context)
Abstract. We develop highresolution shockcapturing numerical schemes for hyperbolic systems with relaxation. In such systems the relaxation time may vary from order1 to much less than unity. When the relaxation time is small, the relaxation term becomes very strong and highly stiff, and underresolved numerical schemes may produce spurious results. Usually one cannot decouple the problem into separate regimes and handle different regimes with different methods. Thus it is important to have a scheme that works uniformly with respect to the relaxation time. Using the Broadwell model of the nonlinear Boltzmann equation we develop a secondorder scheme that works effectively, with a fixed spatial and temporal discretization, for all ranges of the mean free path. Formal uniform consistency proof for a firstorder scheme and numerical convergence proof for the secondorder scheme are also presented. We also make numerical comparisons of the new scheme with some other schemes. This study is motivated by the reentry problem in hypersonic computations.
Third Order Nonoscillatory Central Scheme For Hyperbolic Conservation Laws
"... . A thirdorder accurate Godunovtype scheme for the approximate solution of hyperbolic systems of conservation laws is presented. Its two main ingredients include: #1. A nonoscillatory piecewisequadratic reconstruction of pointvalues from their given cell averages; and #2. A central differencing ..."
Abstract

Cited by 69 (12 self)
 Add to MetaCart
. A thirdorder accurate Godunovtype scheme for the approximate solution of hyperbolic systems of conservation laws is presented. Its two main ingredients include: #1. A nonoscillatory piecewisequadratic reconstruction of pointvalues from their given cell averages; and #2. A central differencing based on staggered evolution of the reconstructed cell averages. This results in a thirdorder central scheme, an extension along the lines of the secondorder central scheme of Nessyahu and Tadmor [NT]. The scalar scheme is nonoscillatory (and hence  convergent), in the sense that it does not increase the number of initial extrema (as does the exact entropy solution operator). Extension to systems is carried out by componentwise application of the scalar framework. In particular, we have the advantage that, unlike upwind schemes, no (approximate) Riemann solvers, fieldbyfield characteristic decompositions, etc., are required. Numerical experiments confirm the highresolution content of...
Relaxation of energy and approximate Riemann solvers for general pressure laws in dynamics
 SIAM J. Num. Anal
, 1998
"... Abstract. We consider the Euler equations for a compressible inviscid fluid with a general pressure law p(ρ, ε), where ρ represents the density of the fluid and ε its specific internal energy. We show that it is possible to introduce a relaxation of the nonlinear pressure law introducing an energy d ..."
Abstract

Cited by 62 (7 self)
 Add to MetaCart
(Show Context)
Abstract. We consider the Euler equations for a compressible inviscid fluid with a general pressure law p(ρ, ε), where ρ represents the density of the fluid and ε its specific internal energy. We show that it is possible to introduce a relaxation of the nonlinear pressure law introducing an energy decomposition under the form ε = ε1 + ε2. The internal energy ε1 is associated with a (simpler) pressure law p1(ρ, ε1); the energy ε2 is advected by the flow. These two energies are also subject to a relaxation process and in the limit of an infinite relaxation rate, we recover the initial pressure law p. We show that, under some conditions of subcharacteristic type, for any convex entropy associated with the pressure p, we can find a global convex and uniform entropy for the relaxation system. From our construction, we also deduce the extension to general pressure laws of classical approximate Riemann solvers for polytropic gases, which only use a single call to the pressure law (per mesh point and time step). For the Godunov scheme, we show that this extension satisfies stability, entropy, and accuracy conditions.
Central WENO Schemes for Hyperbolic Systems of Conservation Laws
 MATH. MODEL. NUMER. ANAL
, 2001
"... We present a family of highorder, essentially nonoscillatory, central schemes for approximating solutions of hyperbolic systems of conservation laws. These schemes are based on a new centered version of the Weighed Essentially NonOscillatory (WENO) reconstruction of pointvalues from cellaverages ..."
Abstract

Cited by 59 (13 self)
 Add to MetaCart
We present a family of highorder, essentially nonoscillatory, central schemes for approximating solutions of hyperbolic systems of conservation laws. These schemes are based on a new centered version of the Weighed Essentially NonOscillatory (WENO) reconstruction of pointvalues from cellaverages, which is then followed by an accurate approximation of the fluxes via a natural continuous extension of RungeKutta solvers. We explicitly construct the third and fourthorder scheme and demonstrate their highresolution properties in several numerical tests.
Discrete Kinetic Schemes For Multidimensional Systems Of Conservation Laws
 SIAM J. Numer. Anal
, 2000
"... We present here some numerical schemes for general multidimensional systems of conservation laws based on a class of discrete kinetic approximations, which includes the relaxation schemes by S. Jin and Z. Xin. These schemes have a simple formulation even in the multidimensional case and do not need ..."
Abstract

Cited by 56 (13 self)
 Add to MetaCart
(Show Context)
We present here some numerical schemes for general multidimensional systems of conservation laws based on a class of discrete kinetic approximations, which includes the relaxation schemes by S. Jin and Z. Xin. These schemes have a simple formulation even in the multidimensional case and do not need the solution of the local Riemann problems. For these approximations we give a suitable multidimensional generalization of the Whitham's stability subcharacteristic condition. In the scalar multidimensional case we establish the rigorous convergence of the approximated solutions to the unique entropy solution of the equilibrium Cauchy problem.
An Evans function approach to spectral stability of smallamplitude viscous shock profiles
, 2002
"... In recent work, the second author and various collaborators have shown using Evans function/refined semigroup techniques that, under very general circumstances, the problems of determining one or multidimensional nonlinear stability of a smooth shock profile may be reduced to that of determining ..."
Abstract

Cited by 55 (41 self)
 Add to MetaCart
In recent work, the second author and various collaborators have shown using Evans function/refined semigroup techniques that, under very general circumstances, the problems of determining one or multidimensional nonlinear stability of a smooth shock profile may be reduced to that of determining spectral stability of the corresponding linearized operator about the wave. It is expected that this condition should in general be analytically verifiable in the case of small amplitude profiles, but this has so far been shown only on a casebycase basis using clever (and difficult to generalize) energy estimates. Here, we describe how the same set of Evans function tools that were used to accomplish the original reduction can be used to show also smallamplitude spectral stability by a direct and readily generalizable procedure. This approach both recovers the results obtained by energy methods, and yields new results not previously obtainable. In particular, we establish onedimensional stability of small amplitude relaxation profiles, completing the Evans function program set out in Mascia&Zumbrun [MZ.1]. Multidimensional stability of small amplitude viscous profiles will be addressed in a companion paper [PZ], completing the program of Zumbrun [Z.3].
Approximate Solutions of Nonlinear Conservation Laws and Related Equations
, 1997
"... During the recent decades there was an enormous amount of activity related to the construction and analysis of modern algorithms for the approximate solution of nonlinear hyperbolic conservation laws and related problems. To present some aspects of this successful activity, we discuss the analytical ..."
Abstract

Cited by 52 (10 self)
 Add to MetaCart
During the recent decades there was an enormous amount of activity related to the construction and analysis of modern algorithms for the approximate solution of nonlinear hyperbolic conservation laws and related problems. To present some aspects of this successful activity, we discuss the analytical tools which are used in the development of convergence theories for these algorithms. These include classical compactness arguments (based on BV a priori estimates), the use of compensated compactness arguments (based on H^1compact entropy production), measure valued solutions (measured by their negative entropy production), and finally, we highlight the most recent addition to this bag of analytical tools  the use of averaging lemmas which yield new compactness and regularity results for nonlinear conservation laws and related equations. We demonstrate how these analytical tools are used in the convergence analysis of approximate solutions for hyperbolic conservation laws and related equations. Our discussion includes examples of Total Variation Diminishing (TVD) finitedifference schemes; error estimates derived from the onesided stability of Godunovtype methods for convex conservation laws (and their multidimensional analogue  viscosity solutions of demiconcave HamiltonJacobi equations); we outline, in the onedimensional case, the convergence proof of finiteelement streamlinediffusion and spectral viscosity schemes based on the divcurl lemma; we also address the questions of convergence and error estimates for multidimensional finitevolume schemes on nonrectangular grids; and finally, we indicate the convergence of approximate solutions with underlying kinetic formulation, e.g., finitevolume and relaxation schemes, once their regularizing effect is quantified in terms of the averaging lemma.