Abstract
The aim of this work is to develop an original semi-coupled approach for solving the Saint-Venant and Exner equation system. The shallow-water equations are solved first to calculate the full evolution of flow fields, and then the Exner equation is solved to model the morphodynamic response. The numerical method is based on a projection method, which consists in combining the momentum and continuity equations in order to establish a Poisson-type equation for water surface levels. The present formulation is identified as a semi-coupled approach because the projection method incorporates the change in the bed evolution in its approximation. A second-order numerical scheme has been proposed using an unstructured finite-volume technique and an implicit time integration method. Several benchmarks are used to demonstrate the capabilities, accuracy and performance of the model.
Notation
= | constant coefficient (s m) | |
= | Chézy coefficient (–) | |
= | sediment grain size (m) | |
= | vector of source terms (–) | |
f | = | Coriolis parameter (–) |
F | = | bed friction force (N) |
g | = | gravity acceleration (m s) |
h | = | total water depth (m) |
H | = | depth of reference (m) |
= | Manning's coefficient (–) | |
= | vector of unit discharge (m s) | |
= | components of in x and y directions (–) | |
= | vector of sediment transport discharge (m s) | |
= | components of in x and y directions (–) | |
t | = | time (s) |
= | vector of depth-averaged velocity (–) | |
u,v | = | components of in x and y directions (–) |
= | friction velocity (m ) | |
= | bed thickness of non-erodible bottom (m) | |
= | water surface (m) | |
δ | = | time difference of water surface elevation (m) |
ν | = | horizontal dispersion coefficient (m s) |
ρ | = | water density (kg m) |
= | sediment density (kg m) | |
σ | = | porosity of the bed material (–) |
= | bed friction force (N) | |
= | Shields parameter (–) | |
= | critical shear stress (–) | |
ξ | = | model porosity (–) |