An accurate and efficient semi‐implicit method for section‐averaged free‐surface flow modelling

Abstract

SUMMARYAn accurate, efficient and robust numerical method for the solution of the section-averaged De St. Venant equations of open channel flow is presented and discussed. The method consists in a semi-implicit, finitevolume discretization of the continuity equation capable to deal with arbitrary cross-section geometry and in a semi-implicit, finite-difference discretization of the momentum equation. By using a proper semiLagrangian discretization of the momentum equation, a highly efficient scheme that is particularly suitable for subcritical regimes is derived. Accurate solutions are obtained in all regimes, except in presence of strong unsteady shocks as in dam-break cases. By using a suitable upwind, Eulerian discretization of the same equation, instead, a scheme capable of describing accurately also unsteady shocks can be obtained, although this scheme requires to comply with a more restrictive stability condition. The formulation of the two approaches allows a unified implementation and an easy switch between the two. The code is verified in a wide range of idealized test cases, highlighting its accuracy and efficiency characteristics, especially for long time range simulations of subcritical river flow. Finally, a model validation on field data is presented, concerning simulations of a flooding event of the Adige river.

A semi‐implicit numerical model for urban drainage systems

“…At each pressure point, a semi‐implicit finite volume approximation of the free surface Equation is taken to be where and σ i , j is a sign function associated with the orientation of the j th axial velocity. Specifically, …”

A semi-implicit, semi-Lagrangian, p-adaptive discontinuous Galerkin method for the shallow water equations

“…As discussed in [53], methods employing the non conservative formulation of the momentum equation satisfy automatically the so called C-property and are naturally well balanced, so that no problems arise due to the variable bathymetry in the still water case.…”

Modelling the transition between fixed and mobile bed conditions in two-phase free-surface flows: The Composite Riemann Problem and its numerical solution

“…In the literature there is a broad set of different depth-averaged mathematical models which concern the mobile-bed case (see, among many others, [5,7,16,23,28,30]). This is essentially due to the lack of a unified approach which is able to describe any type of sediment-laden flows, independently from the concentration of the granular phase.…”