Mathematical balancing of flux gradient and source terms prior to using Roe's approximate Riemann solver

BD Rogers, AGL Borthwick*, PH Taylor

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

Abstract / Description of output

An algebraic technique is presented for balancing flux gradients and source terms when applying Roe's approximate Riemann solver in finite volume schemes. The numerical imbalance is eradicated by reformulating the governing matrix hyperbolic system of conservation laws in terms of deviations away from an unforced but separately specified equilibrium state. Thus, balancing is achieved by the incorporation of this extra physical information and bypasses conventional numerical treatments of the imbalance. The technique is first applied to the shallow water equations. Simulations of benchmark flows including wind-induced flow in a two-dimensional basin, transcritical flow in a one-dimensional channel and wave propagation over a two-dimensional hump are in close agreement with analytical solutions and predictions by alternative numerical schemes. The technique is then applied to a more complicated coupled pair of equation sets, the hyperbolic period- and depth-averaged ray-type wave conservation and modified shallow water equations that describe wave current interaction in the nearshore zone at the coast. Reasonable agreement is obtained with laboratory measurements of wave diffraction behind a submerged elliptical shoal [Coastal Engrg. 6 (1982) 255] and of wave-induced nearshore currents at a half-sinusoidal beach [Wave-induced nearshore currents, Ph.D. Thesis, Liverpool University, UK, 1981]. (C) 2003 Elsevier B.V. All rights reserved.

Original languageEnglish
Pages (from-to)422-451
Number of pages30
JournalJournal of Computational Physics
Issue number2
Publication statusPublished - 10 Dec 2003

Keywords / Materials (for Non-textual outputs)

  • Roe's approximate Riemann solver
  • flux gradients
  • source terms
  • shallow-water equations
  • wave-conservation equations


Dive into the research topics of 'Mathematical balancing of flux gradient and source terms prior to using Roe's approximate Riemann solver'. Together they form a unique fingerprint.

Cite this