Results 1  10
of
110
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 56 (4 self)
 Add to MetaCart
(Show Context)
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.
Numerical Solution of an Inverse Problem in SizeStructured Population Dynamics, in "Inverse Problems
 n o 4, april 2009 BR . Activity Report INRIA 2009
"... We consider a sizestructured model for cell division and address the question of determining the division (birth) rate from the measured stable size distribution of the population. We propose a new regularization technique based on a filtering approach. We prove convergence of the algorithm and val ..."
Abstract

Cited by 15 (4 self)
 Add to MetaCart
(Show Context)
We consider a sizestructured model for cell division and address the question of determining the division (birth) rate from the measured stable size distribution of the population. We propose a new regularization technique based on a filtering approach. We prove convergence of the algorithm and validate the theoretical results by implementing numerical simulations, based on classical techniques. We compare the results for direct and inverse problems, for the filtering method and for the quasireversibility method proposed in [1]. 1
2005), On the use of Saint Venant equations to simulate the spreading of a granular mass
 J. Geophys. Res
"... [1] Cliff collapse is an active geomorphological process acting at the surface of the Earth and telluric planets. Recent laboratory studies have investigated the collapse of an initially cylindrical granular mass along a rough horizontal plane for different initial aspect ratios a = Hi/Ri, where Hi ..."
Abstract

Cited by 14 (8 self)
 Add to MetaCart
(Show Context)
[1] Cliff collapse is an active geomorphological process acting at the surface of the Earth and telluric planets. Recent laboratory studies have investigated the collapse of an initially cylindrical granular mass along a rough horizontal plane for different initial aspect ratios a = Hi/Ri, where Hi and Ri are the initial height and radius, respectively. A numerical simulation of these experiments is performed using a minimal depthintegrated model based on a longwave approximation. A dimensional analysis of the equations shows that such a model exhibits the scaling laws observed experimentally. Generic solutions are independent of gravity and depend only on the initial aspect ratio a and an effective friction angle. In terms of dynamics, the numerical simulations are consistent with the experiments for a 1. The experimentally observed saturation of the final height of the deposit, when normalized with respect to the initial radius of the cylinder, is accurately reproduced numerically. Analysis of the results sheds light on the correlation between the area overrun by the granular mass and its initial potential energy. The extent of the deposit, the final height, and the arrest time of the front can be directly estimated from the ‘‘generic solution’ ’ of the model for terrestrial and extraterrestrial avalanches. The effective friction, a parameter classically used to describe the mobility of gravitational flows, is shown to depend on the initial aspect ratio a. This dependence should be taken into account when interpreting the high mobility of large volume events.
APPROXIMATION OF HYPERBOLIC MODELS FOR CHEMOSENSITIVE MOVEMENT
, 2005
"... Numerical methods with different orders of accuracy are proposed to approximate hyperbolic models for chemosensitive movements. On the one hand, first and secondorder wellbalanced finite volume schemes are presented. This approach provides exact conservation of the steady state solutions. On th ..."
Abstract

Cited by 11 (2 self)
 Add to MetaCart
(Show Context)
Numerical methods with different orders of accuracy are proposed to approximate hyperbolic models for chemosensitive movements. On the one hand, first and secondorder wellbalanced finite volume schemes are presented. This approach provides exact conservation of the steady state solutions. On the other hand, a highorder finite difference weighted essentially nonoscillatory (WENO) scheme is constructed and the wellbalanced reconstruction is adapted to this scheme to exactly preserve steady states and to retain highorder accuracy. Numerical simulations are performed to verify accuracy and the wellbalanced property of the proposed schemes and to observe the formation of networks in the hyperbolic models similar to those observed in the experiments.
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 11 (3 self)
 Add to MetaCart
(Show Context)
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.
Asymptotic preserving (AP) schemes for multiscale kinetic and hyperbolic equations: a review
, 2012
"... Kinetic and hyperbolic equations contain small scales (mean free path/time, Debye length, relaxation or reaction time, etc.) that lead to various different asymptotic regimes, in which the classical numerical approximations become prohibitively expensive. Asymptoticpreserving (AP) schemes are schem ..."
Abstract

Cited by 10 (5 self)
 Add to MetaCart
(Show Context)
Kinetic and hyperbolic equations contain small scales (mean free path/time, Debye length, relaxation or reaction time, etc.) that lead to various different asymptotic regimes, in which the classical numerical approximations become prohibitively expensive. Asymptoticpreserving (AP) schemes are schemes that are efficient in these asymptotic regimes. The designing principle of AP schemes is to preserve, at the discrete level, the asymptotic limit that drives one (usually the microscopic) equation to its asymptotic (macroscopic) equation. An AP scheme is based on solving the microscopic equation, instead of using a multiphysics approach that couples different physical laws at different scales. When the small scale is not numerically resolved, an AP scheme automatically becomes a macroscopic solver for the limiting equation. The AP methodology offers simple, robust and efficient computational methods for a large class of multiscale kinetic, hyperbolic and other physical problems. This
Why many theories of shock waves are necessary. Convergence error in formally pathconsistent schemes
, 2008
"... ..."
An asymptotic high order masspreserving scheme for a hyperbolic model of chemotaxis
 SIAM J. Num. Anal
"... Abstract. We introduce a new class of finite difference schemes for approximating the solutions to an initialboundary value problem on a bounded interval for a one dimensional dissipative hyperbolic system with an external source term, which arises as a simple model of chemotaxis. Since the solutio ..."
Abstract

Cited by 8 (2 self)
 Add to MetaCart
(Show Context)
Abstract. We introduce a new class of finite difference schemes for approximating the solutions to an initialboundary value problem on a bounded interval for a one dimensional dissipative hyperbolic system with an external source term, which arises as a simple model of chemotaxis. Since the solutions to this problem may converge to non constant asymptotic states for large times, standard schemes usually fail to yield a good approximation. Therefore, we propose a new class of schemes, which use an asymptotic higher order correction, second and third order in our examples, to balance the effects of the source term and the influence of the asymptotic solutions. A special care is needed to deal with boundary conditions, to avoid harmful loss of mass. Convergence results are proven for these new schemes, and several numerical tests are presented and discussed to verify the effectiveness of their behavior.
A MULTIWAVE APPROXIMATE RIEMANN SOLVER FOR IDEAL MHD BASED ON RELAXATION II NUMERICAL IMPLEMENTATION WITH 3 AND 5 WAVES
"... Abstract. In the first part of this work ([5]), we introduced an approximate Riemann solver for onedimensional ideal MHD derived from a relaxation system. We gave sufficient conditions for the solver to satisfy discrete entropy inequalities, and to preserve positivity of density and internal energy ..."
Abstract

Cited by 8 (4 self)
 Add to MetaCart
(Show Context)
Abstract. In the first part of this work ([5]), we introduced an approximate Riemann solver for onedimensional ideal MHD derived from a relaxation system. We gave sufficient conditions for the solver to satisfy discrete entropy inequalities, and to preserve positivity of density and internal energy. In this paper we consider the practical implementation, and derive explicit wave speed estimates satisfying the stability conditions of [5]. We present a 3wave solver that well resolves fast waves and material contacts, and a 5wave solver that accurately resolves the cases when two eigenvalues coincide. A full 7wave solver, which is highly accurate on all types of waves, will be described in a followup paper. We test the solvers on onedimensional shock tube data and smooth shear waves. (1.1) (1.2)
A SecondOrder WellBalanced Positivity Preserving CentralUpwind Scheme for the SaintVenant System
 Communications in Mathematical Sciences
"... Abstract. A family of Godunovtype centralupwind schemes for the SaintVenant system of shallow water equations has been first introduced in [A. Kurganov and D. Levy, M2AN Math. Model. Numer. Anal., 36 (2002), pp. 397–425]. Depending on the reconstruction step, the secondorder versions of the sche ..."
Abstract

Cited by 8 (0 self)
 Add to MetaCart
(Show Context)
Abstract. A family of Godunovtype centralupwind schemes for the SaintVenant system of shallow water equations has been first introduced in [A. Kurganov and D. Levy, M2AN Math. Model. Numer. Anal., 36 (2002), pp. 397–425]. Depending on the reconstruction step, the secondorder versions of the schemes there could be made either wellbalanced or positivity preserving, but fail to satisfy both properties simultaneously. Here, we introduce an improved secondorder centralupwind scheme which, unlike its forerunners, is capable to both preserve stationary steady states (lake at rest) and to guarantee the positivity of the computed fluid depth. Another novel property of the proposed scheme is its applicability to models with discontinuous bottom topography. We demonstrate these features of the new scheme, as well as its high resolution and robustness in a number of one and twodimensional examples. Key words. Hyperbolic systems of conservation and balance laws, semidiscrete centralupwind schemes, SaintVenant system of shallow water equations. AMS subject classifications. 65M99, 35L65 1.