Results 1  10
of
278
Balancing Source Terms and Flux Gradients in HighResolution Godunov Methods: The QuasiSteady WavePropogation Algorithm
 J. Comput. Phys
, 1998
"... . Conservation laws with source terms often have steady states in which the flux gradients are nonzero but exactly balanced by source terms. Many numerical methods (e.g., fractional step methods) have difficulty preserving such steady states and cannot accurately calculate small perturbations of suc ..."
Abstract

Cited by 78 (5 self)
 Add to MetaCart
. Conservation laws with source terms often have steady states in which the flux gradients are nonzero but exactly balanced by source terms. Many numerical methods (e.g., fractional step methods) have difficulty preserving such steady states and cannot accurately calculate small perturbations of such states. Here a variant of the wavepropagation algorithm is developed which addresses this problem by introducing a Riemann problem in the center of each grid cell whose flux difference exactly cancels the source term. This leads to modified Riemann problems at the cell edges in which the jump now corresponds to perturbations from the steady state. Computing waves and limiters based on the solution to these Riemann problems gives highresolution results. The 1D and 2D shallow water equations for flow over arbitrary bottom topography are use as an example, though the ideas apply to many other systems. The method is easily implemented in the software package clawpack. Keywords: Godunov meth...
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
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.
Some approximate Godunov schemes to compute shallowwater equatons with topography
, 2003
"... WestEC here te comput tmp of shallowwath equat ons wi ttN:LE&BN y byFinit Volume metmeN , in a onedimensional framework(tNL&3 allmetB: sint oduced may benatEfiLLN extEfiLL in t o dimensions) . AllmetA:3 are based on a discretcrNB on of tN tBCEAEN(B by a piecewisefunctE n constEC on each ce ..."
Abstract

Cited by 54 (5 self)
 Add to MetaCart
WestEC here te comput tmp of shallowwath equat ons wi ttN:LE&BN y byFinit Volume metmeN , in a onedimensional framework(tNL&3 allmetB: sint oduced may benatEfiLLN extEfiLL in t o dimensions) . AllmetA:3 are based on a discretcrNB on of tN tBCEAEN(B by a piecewisefunctE n constEC on each cell of tN mesh, from an original idea of Le Rouxet al. Whereaste WellBalanced scheme of Le Roux is based on tN exact resol utol of each Riemann problem, we consider here approximat Riemann solvers. Several singlestg metleN are derived from tom formalism, and numerical result are comparedt a fract ionalsta metN4 . Some tme cases arepresent& : convergencetn ardsstsN: stsN: in subcritECC and supercriterc configuratfiE s, occurrence of dry area by a drain over a bump and occurrence of vacuum by a double rarefactNL wave over astAA Numerical schemes, combined wi t an appropriat highorder extNfi ion, provideaccurat e and convergent approxim atximN # 2003 Elsevier ScienceLtc Allright s reserved. 1. I5964V139 We stE4 intNL paper some approximat Godunov schemest comput shallowwatN equatlow wit a sourcetur oftBAL4fiN(LA in a onedimensional framework. Allmetfi3L presentL may beextB3L4 natB3L4N t tt tB3L4N(&CL4 al model.
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 53 (6 self)
 Add to MetaCart
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 51 (14 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 46 (13 self)
 Add to MetaCart
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.
D.: A thirdorder semidiscrete central scheme for conservation laws and convectiondiffusion equations
 SIAM J. Sci. Comput
, 2000
"... We present a new thirdorder, semidiscrete, central method for approximating solutions to multidimensional systems of hyperbolic conservation laws, convectiondiffusion equations, and related problems. Our method is a highorder extension of the recently proposed secondorder, semidiscrete method ..."
Abstract

Cited by 45 (4 self)
 Add to MetaCart
We present a new thirdorder, semidiscrete, central method for approximating solutions to multidimensional systems of hyperbolic conservation laws, convectiondiffusion equations, and related problems. Our method is a highorder extension of the recently proposed secondorder, semidiscrete method in [16]. The method is derived independently of the specific piecewise polynomial reconstruction which is based on the previously computed cellaverages. We demonstrate our results, by focusing on the new thirdorder CWENO reconstruction presented in [21]. The numerical results we present, show the desired accuracy, high resolution and robustness of our method. Key words. Hyperbolic systems, convectiondiffusion equations, central difference schemes, highorder accuracy, nonoscillatory schemes, WENO reconstruction. AMS(MOS) subject classification. Primary 65M10; secondary 65M05.
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 45 (11 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.
Fully adaptive multiresolution finite volume schemes for conservation laws
 Math. Comp
, 2003
"... Abstract. The use of multiresolution decompositions in the context of finite volume schemes for conservation laws was first proposed by A. Harten for the purpose of accelerating the evaluation of numerical fluxes through an adaptive computation. In this approach the solution is still represented at ..."
Abstract

Cited by 41 (13 self)
 Add to MetaCart
Abstract. The use of multiresolution decompositions in the context of finite volume schemes for conservation laws was first proposed by A. Harten for the purpose of accelerating the evaluation of numerical fluxes through an adaptive computation. In this approach the solution is still represented at each time step on the finest grid, resulting in an inherent limitation of the potential gain in memory space and computational time. The present paper is concerned with the development and the numerical analysis of fully adaptive multiresolution schemes, in which the solution is represented and computed in a dynamically evolved adaptive grid. A crucial problem is then the accurate computation of the flux without the full knowledge of fine grid cell averages. Several solutions to this problem are proposed, analyzed, and compared in terms of accuracy and complexity. 1.
Numerical Solution Of Reservoir Flow Models Based On Large Time Step Operator Splitting Algorithms
 FILTRATION IN POROUS MEDIA AND INDUSTRIAL APPLICATIONS, LECTURE NOTES IN MATHEMATICS
, 1999
"... During recent years the authors and collaborators have been involved in an activity related to the construction and analysis of large time step operator splitting algorithms for the numerical simulation of multiphase flow in heterogeneous porous media. The purpose of these lecture notes is to revie ..."
Abstract

Cited by 35 (20 self)
 Add to MetaCart
During recent years the authors and collaborators have been involved in an activity related to the construction and analysis of large time step operator splitting algorithms for the numerical simulation of multiphase flow in heterogeneous porous media. The purpose of these lecture notes is to review some of this activity. We illustrate the main ideas behind these novel operator splitting algorithms for a basic twophase flow model. Special focus is posed on the numerical solution algorithms for the saturation equation, which is a convection dominated, degenerate convectiondiffusion equation. Both theory and applications are discussed. The general background for the reservoir flow model is reviewed, and the main features of the numerical algorithms are presented. The basic mathematical results supporting the numerical algorithms are also given. In addition, we present some results from the BV solution theory for quasilinear degenerate parabolic equations, which provides the correct ...