Preprint 1998-004

Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods: The Quasi-Steady Wave-Propogation Algorithm

Randall J. LeVeque

Abstract: 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 wave-propagation 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 high-resolution 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.

Available as PostScript (1.1 Mbytes) or gzipped PostScript (248 Kbytes; uncompress using gunzip).
Balancing source terms and flux gradients in high-resolution Godunov methods: the quasi-steady wave-propogation algorithm
Randall J. LeVeque, <>
Publishing information:
J. Comput. Phys. 146 (1998) 346-365.
Submitted by:
<> January 26 1998.

[ 1996 Preprints | 1997 Preprints | 1998 Preprints | All Preprints | Preprint Server Homepage ]
© The copyright for the following documents lies with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use.

Conservation Laws Preprint Server <>
Last modified: Tue Jan 23 08:43:35 MET 2001