Results 1  10
of
35
A fast and stable wellbalanced scheme with hydrostatic reconstruction for shallow water flows
 SIAM J. Sci. Comput
"... Abstract. We consider the SaintVenant system for shallow water flows, with nonflat bottom. It is a hyperbolic system of conservation laws that approximately describes various geophysical flows, such as rivers, coastal areas, and oceans when completed with a Coriolis term, or granular flows when com ..."
Abstract

Cited by 44 (4 self)
 Add to MetaCart
Abstract. We consider the SaintVenant system for shallow water flows, with nonflat bottom. It is a hyperbolic system of conservation laws that approximately describes various geophysical flows, such as rivers, coastal areas, and oceans when completed with a Coriolis term, or granular flows when completed with friction. Numerical approximate solutions to this system may be generated using conservative finite volume methods, which are known to properly handle shocks and contact discontinuities. However, in general these schemes are known to be quite inaccurate for near steady states, as the structure of their numerical truncation errors is generally not compatible with exact physical steady state conditions. This difficulty can be overcome by using the socalled wellbalanced schemes. We describe a general strategy, based on a local hydrostatic reconstruction, that allows us to derive a wellbalanced scheme from any given numerical flux for the homogeneous problem. Whenever the initial solver satisfies some classical stability properties, it yields a simple and fast wellbalanced scheme that preserves the nonnegativity of the water height and satisfies a semidiscrete entropy inequality.
Hamiltonianpreserving schemes for the Liouville equation with discontinuous potentials
"... When numerically solving the Liouville equation with a discontinuous potential, one faces the problem of selecting a unique, physically relevant solution across the potential barrier, and the problem of a severe time step constraint due to the CFL condition. In this paper, We introduce two classes o ..."
Abstract

Cited by 24 (21 self)
 Add to MetaCart
When numerically solving the Liouville equation with a discontinuous potential, one faces the problem of selecting a unique, physically relevant solution across the potential barrier, and the problem of a severe time step constraint due to the CFL condition. In this paper, We introduce two classes of Hamiltonianpreserving schemes for such problems. By using the constant Hamiltonian across the potential barrier, we introduced a selection criterion for a unique, physically relavant solution to the underlying linear hyperbolic equation with singular coefficients. These scheme have a hyperbolic CFL condition, which is a significant improvement over a conventional discretization. We also establish the positivity, and stability in both l 1 and l ∞ norms, of these discretizations, and conducted numerical experiments to study the numerical accuracy. This work is motivated by the wellbalanced kinetic schemes by Perthame and Simeoni for the shallow water equations with a discontinuous bottom topography, and has applications to the level set methods for the computations of multivalued physical observables in the semiclassical limit of the linear Schrödinger equation with a discontinuous potential, among other applications.
From Kinetic Equations to Multidimensional Isentropic Gas Dynamics Before Shocks
, 2003
"... This article is devoted to the proof of the hydrodynamical limit from kinetic equations (including B.G.K. like equations) to multidimensional isentropic gas dynamics. It is based on a relative entropy method, hence the derivation is valid only before shocks appear on the limit system solution. Howev ..."
Abstract

Cited by 17 (2 self)
 Add to MetaCart
This article is devoted to the proof of the hydrodynamical limit from kinetic equations (including B.G.K. like equations) to multidimensional isentropic gas dynamics. It is based on a relative entropy method, hence the derivation is valid only before shocks appear on the limit system solution. However, no a priori knowledge on high velocities distributions for kinetic functions is needed. The case of the SaintVenant system with topography (where a source term is added) is included. Keywords: Hydrodynamic limit, Entropy method, B.G.K. equation, Isentropic gas dynamics, SaintVenant system.
Space Localization And WellBalanced Schemes For Discrete Kinetic Models In Diffusive Regimes
 SIAM J. Numer. Anal
, 2002
"... We derive and study WellBalanced schemes for quasimonotone discrete kinetic models. By means of a rigorous localization procedure, we reformulate the collision terms as nonconservative products and solve the resulting Riemann problem whose solution is selfsimilar. The construction of an Asymptotic ..."
Abstract

Cited by 12 (3 self)
 Add to MetaCart
We derive and study WellBalanced schemes for quasimonotone discrete kinetic models. By means of a rigorous localization procedure, we reformulate the collision terms as nonconservative products and solve the resulting Riemann problem whose solution is selfsimilar. The construction of an Asymptotic Preserving (AP) Godunov scheme is straightforward and various compactness properties are established within different scalings. At last, some computational results are supplied to show that this approach is realizable and ecient on concrete 2 × 2 models.
Asymptotic preserving (AP) schemes for multiscale kinetic and hyperbolic equations: a review, Lecture Notes for Summer School on ”Methods and Models of Kinetic Theory
, 2010
"... 2. Hyperbolic systems with stiff relaxations 3 3. Kinetic equations: the Euler regime 8 4. Linear transport equations: the diffusion regime 15 ..."
Abstract

Cited by 10 (5 self)
 Add to MetaCart
2. Hyperbolic systems with stiff relaxations 3 3. Kinetic equations: the Euler regime 8 4. Linear transport equations: the diffusion regime 15
Godunovtype Approximation for a General Resonant Balance Law With Large Data
, 2003
"... We consider the Cauchy problem for the 2 2 nonstrictly hyperbolic system u t + f(a; u) x g(a; u)a x = 0 (a; u)(0; ) = (a o ; u o ): For possibly large, discontinuous and resonant data, the generalized solution to the Riemann problem is introduced, interaction estimates are carried out using a ne ..."
Abstract

Cited by 9 (3 self)
 Add to MetaCart
We consider the Cauchy problem for the 2 2 nonstrictly hyperbolic system u t + f(a; u) x g(a; u)a x = 0 (a; u)(0; ) = (a o ; u o ): For possibly large, discontinuous and resonant data, the generalized solution to the Riemann problem is introduced, interaction estimates are carried out using a new change of variables and the convergence of Godunov approximations is shown. Uniqueness is addressed relying on a suitable extension of Kruzkov's techniques. 1
Relaxation Schemes for the Shallow Water Equations
 INT. J. NUMER. METH. FLUIDS
, 2003
"... We present a class of first and second order in space and time relaxation schemes for the shallow water (SW) equations. A new approach of incorporating the geometrical source term in the relaxation model is also presented. The schemes are based on classical relaxation models combined with RungeKut ..."
Abstract

Cited by 7 (0 self)
 Add to MetaCart
We present a class of first and second order in space and time relaxation schemes for the shallow water (SW) equations. A new approach of incorporating the geometrical source term in the relaxation model is also presented. The schemes are based on classical relaxation models combined with RungeKutta time stepping mechanisms. Numerical results are presented for several benchmark test problems with or without the source term present.
High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallowwater systems
 Math. Comp
"... Abstract. This paper is concerned with the development of high order methods for the numerical approximation of onedimensional nonconservative hyperbolic systems. In particular, we are interested in high order extensions of the generalized Roe methods introduced by I. Toumi in 1992, based on WENO r ..."
Abstract

Cited by 6 (2 self)
 Add to MetaCart
Abstract. This paper is concerned with the development of high order methods for the numerical approximation of onedimensional nonconservative hyperbolic systems. In particular, we are interested in high order extensions of the generalized Roe methods introduced by I. Toumi in 1992, based on WENO reconstruction of states. We also investigate the wellbalanced properties of the resulting schemes. Finally, we will focus on applications to shallowwater systems. 1.
The L 1 stability of a Hamiltonianpreserving scheme for the Liouville equation with discontinuous potentials, preprint
"... equation with discontinuous potentials ∗ ..."
Computation of transmissions and reflections in geometrical optics via the reduced Liouville equation
 WAVE MOTION
, 2006
"... We develop a numerical scheme for the wave front computation of complete transmissions and reflections in geometrical optics. Such a problem can be formulated by a reduced Liouville equation with a discontinuous local wave speed or index of refraction, arising in the high frequency limit of linear w ..."
Abstract

Cited by 5 (5 self)
 Add to MetaCart
We develop a numerical scheme for the wave front computation of complete transmissions and reflections in geometrical optics. Such a problem can be formulated by a reduced Liouville equation with a discontinuous local wave speed or index of refraction, arising in the high frequency limit of linear waves through inhomogeneous media. The key idea is to incorporate Snell’s Law of Refraction into the numerical flux for the reduced Liouville equation. This scheme allows a hyperbolic CFL condition, under which positivity, and stabilities in both l∞ and l¹ norms, are established. Numerical experiments are carried out to demonstrate the validity and accuracy of this new scheme.