Loughborough University
Leicestershire, UK
LE11 3TU
+44 (0)1509 263171
Loughborough University

Loughborough University Institutional Repository

Please use this identifier to cite or link to this item: https://dspace.lboro.ac.uk/2134/33431

Title: A new efficient implicit scheme for discretising the stiff friction terms in the shallow water equations
Authors: Xia, Xilin
Liang, Qiuhua
Keywords: Shallow water equations
Friction source terms
Implicit scheme
Stiff relaxation
Overland flow
Finite volume method
Issue Date: 2018
Publisher: © The Authors. Published by Elsevier
Citation: XIA, X. and LIANG, Q., 2018. A new efficient implicit scheme for discretising the stiff friction terms in the shallow water equations. Advances in Water Resources, 117, pp. 87-97.
Abstract: © 2018 The Authors Discretisation of the friction terms to ensure numerical stability and accuracy remains to be challenging for the development of robust numerical schemes to solve the shallow water equations (SWEs), particularly for applications involving very shallow flows (e.g. overland flows and wet/dry fronts) over complex domain topography. The key challenge is to ensure relaxation of the flow towards an equilibrium state characterised by the balance between friction and gravity in a computationally efficient way. To overcome this numerical challenge, this paper proposes a novel approach for discretising the friction source terms in the SWEs in the context of an explicit finite volume method. The overall numerical scheme adopts the HLLC Riemann solver and surface reconstruction method (SRM) to explicitly discretise the flux and bed slope source terms. Whilst a fully implicit scheme is used to handle the friction source terms, solution to the implicit formulation is analytically derived to explicitly update the flow variables. Compared with the existing approaches, the proposed scheme effectively resolves the issue associated with stiff relaxation without necessity to use an iteration method and it supports efficient simulation using time steps controlled only by the Courant–Friedrichs–Levy (CFL) condition. The current friction term discretisation scheme is not coupled with flux and bed slope calculation and therefore may be readily implemented in any other explicit finite volume SWE models. After being successfully validated against two benchmark tests with analytical solutions, the resulting new SWE model is applied to reproduce a rainfall-flooding event in the Upper Lee catchment in the UK.
Description: This is an Open Access Article. It is published by Elsevier under the Creative Commons Attribution 4.0 Unported Licence (CC BY). Full details of this licence are available at: http://creativecommons.org/licenses/by/4.0/
Sponsor: This work is funded by the NERC SINATRA and TENDERLY projects (Grant No. NE/K008781/1) and Flood-PREPARED project (Grant No. NE/P017134/1).
Version: Published
DOI: 10.1016/j.advwatres.2018.05.004
URI: https://dspace.lboro.ac.uk/2134/33431
Publisher Link: https://doi.org/10.1016/j.advwatres.2018.05.004
ISSN: 0309-1708
Appears in Collections:Published Articles (Architecture, Building and Civil Engineering)

Files associated with this item:

File Description SizeFormat
ADWR 2018.pdfPublished version2.1 MBAdobe PDFView/Open


SFX Query

Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.