High resolution numerical schemes for solving kinematic wave equation

Chunshui Yu, Guohong Duan

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

This paper compares the stability, accuracy, and computational cost of several numerical methods for solving the kinematic wave equation. The numerical methods include the second-order MacCormack finite difference scheme, the MacCormack scheme with a dissipative interface, the second-order MUSCL finite volume scheme, and the fifth-order WENO finite volume scheme. These numerical schemes are tested against several synthetic cases and an overland flow experiment, which include shock wave, rarefaction wave, wave steepening, uniform/non-uniform rainfall generated overland flows, and flow over a channel of varying bed slope. The results show that the MacCormack scheme is not a Total Variation Diminishing (TVD) scheme because oscillatory solutions occurred at the presence of shock wave, rarefaction wave, and overland flow over rapidly varying bed slopes. The MacCormack scheme with a dissipative interface is free of oscillation but with considerable diffusions. The Godunov-type schemes are accurate and stable when dealing with discontinuous waves. Furthermore the Godunov-type schemes, like MUSCL and WENO scheme, are needed for simulating surface flow from spatially non-uniformly distributed rainfalls over irregular terrains using moderate computing resources on current personal computers.

Original languageEnglish (US)
Pages (from-to)823-832
Number of pages10
JournalJournal of Hydrology
Volume519
Issue numberPA
DOIs
StatePublished - Nov 7 2014

Fingerprint

wave equation
overland flow
kinematics
shock wave
numerical method
rainfall
oscillation
resource
cost
experiment

Keywords

  • Godunov-type scheme
  • Kinematic wave equation
  • MacCormack scheme
  • Rainfall-runoff overland flow
  • Rarefaction wave
  • Shock wave

ASJC Scopus subject areas

  • Water Science and Technology

Cite this

High resolution numerical schemes for solving kinematic wave equation. / Yu, Chunshui; Duan, Guohong.

In: Journal of Hydrology, Vol. 519, No. PA, 07.11.2014, p. 823-832.

Research output: Contribution to journalArticle

@article{1c4fb3f2318040828e7a7e9d4111ee01,
title = "High resolution numerical schemes for solving kinematic wave equation",
abstract = "This paper compares the stability, accuracy, and computational cost of several numerical methods for solving the kinematic wave equation. The numerical methods include the second-order MacCormack finite difference scheme, the MacCormack scheme with a dissipative interface, the second-order MUSCL finite volume scheme, and the fifth-order WENO finite volume scheme. These numerical schemes are tested against several synthetic cases and an overland flow experiment, which include shock wave, rarefaction wave, wave steepening, uniform/non-uniform rainfall generated overland flows, and flow over a channel of varying bed slope. The results show that the MacCormack scheme is not a Total Variation Diminishing (TVD) scheme because oscillatory solutions occurred at the presence of shock wave, rarefaction wave, and overland flow over rapidly varying bed slopes. The MacCormack scheme with a dissipative interface is free of oscillation but with considerable diffusions. The Godunov-type schemes are accurate and stable when dealing with discontinuous waves. Furthermore the Godunov-type schemes, like MUSCL and WENO scheme, are needed for simulating surface flow from spatially non-uniformly distributed rainfalls over irregular terrains using moderate computing resources on current personal computers.",
keywords = "Godunov-type scheme, Kinematic wave equation, MacCormack scheme, Rainfall-runoff overland flow, Rarefaction wave, Shock wave",
author = "Chunshui Yu and Guohong Duan",
year = "2014",
month = "11",
day = "7",
doi = "10.1016/j.jhydrol.2014.08.003",
language = "English (US)",
volume = "519",
pages = "823--832",
journal = "Journal of Hydrology",
issn = "0022-1694",
publisher = "Elsevier",
number = "PA",

}

TY - JOUR

T1 - High resolution numerical schemes for solving kinematic wave equation

AU - Yu, Chunshui

AU - Duan, Guohong

PY - 2014/11/7

Y1 - 2014/11/7

N2 - This paper compares the stability, accuracy, and computational cost of several numerical methods for solving the kinematic wave equation. The numerical methods include the second-order MacCormack finite difference scheme, the MacCormack scheme with a dissipative interface, the second-order MUSCL finite volume scheme, and the fifth-order WENO finite volume scheme. These numerical schemes are tested against several synthetic cases and an overland flow experiment, which include shock wave, rarefaction wave, wave steepening, uniform/non-uniform rainfall generated overland flows, and flow over a channel of varying bed slope. The results show that the MacCormack scheme is not a Total Variation Diminishing (TVD) scheme because oscillatory solutions occurred at the presence of shock wave, rarefaction wave, and overland flow over rapidly varying bed slopes. The MacCormack scheme with a dissipative interface is free of oscillation but with considerable diffusions. The Godunov-type schemes are accurate and stable when dealing with discontinuous waves. Furthermore the Godunov-type schemes, like MUSCL and WENO scheme, are needed for simulating surface flow from spatially non-uniformly distributed rainfalls over irregular terrains using moderate computing resources on current personal computers.

AB - This paper compares the stability, accuracy, and computational cost of several numerical methods for solving the kinematic wave equation. The numerical methods include the second-order MacCormack finite difference scheme, the MacCormack scheme with a dissipative interface, the second-order MUSCL finite volume scheme, and the fifth-order WENO finite volume scheme. These numerical schemes are tested against several synthetic cases and an overland flow experiment, which include shock wave, rarefaction wave, wave steepening, uniform/non-uniform rainfall generated overland flows, and flow over a channel of varying bed slope. The results show that the MacCormack scheme is not a Total Variation Diminishing (TVD) scheme because oscillatory solutions occurred at the presence of shock wave, rarefaction wave, and overland flow over rapidly varying bed slopes. The MacCormack scheme with a dissipative interface is free of oscillation but with considerable diffusions. The Godunov-type schemes are accurate and stable when dealing with discontinuous waves. Furthermore the Godunov-type schemes, like MUSCL and WENO scheme, are needed for simulating surface flow from spatially non-uniformly distributed rainfalls over irregular terrains using moderate computing resources on current personal computers.

KW - Godunov-type scheme

KW - Kinematic wave equation

KW - MacCormack scheme

KW - Rainfall-runoff overland flow

KW - Rarefaction wave

KW - Shock wave

UR - http://www.scopus.com/inward/record.url?scp=84906734941&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84906734941&partnerID=8YFLogxK

U2 - 10.1016/j.jhydrol.2014.08.003

DO - 10.1016/j.jhydrol.2014.08.003

M3 - Article

AN - SCOPUS:84906734941

VL - 519

SP - 823

EP - 832

JO - Journal of Hydrology

JF - Journal of Hydrology

SN - 0022-1694

IS - PA

ER -