A steady-state capturing method for hyperbolic systems with geometrical source terms
Department of Mathematics, University of Wisconsin, Madison,
WI 53706, USA. (email@example.com)
Revised: 7 February 2001
We propose a simple numerical method for capturing the steady state solution of hyperbolic systems with geometrical source terms. We use the interface value, rather than the cell-averages, for the source terms that balance the nonlinear convection at the cell interface, allowing the numerical capturing of the steady state with a formal high order accuracy. This method applies to Godunov or Roe type upwind methods but requires no modification of the Riemann solver. Numerical experiments on scalar conservation laws and the one dimensional shallow water equations show much better resolution of the steady state than the conventional method, with almost no new numerical complexity.
Mathematics Subject Classification: 35L65 / 65M06 / 76B15
Key words: Hyperbolic systems / source terms / steady state solution / shallow water equations / shock capturing methods.
© EDP Sciences, SMAI, 2001