Linearized implicit time advancing and defect correction applied to sediment transport simulations

The numerical simulation of sediment transport problems is considered in this paper. The physical problem is modeled through the shallow-water equations coupled with the Exner equation to describe the time evolution of the bed profile. The spatial discretization of the governing equations is carried out by a finite-volume method and a modified Roe scheme designed for non-conservative systems. As for the time advancing, starting from an explicit method, a linearized implicit scheme is generated, in which the flux Jacobian is computed through automatic differentiation. Second-order accuracy in space is then obtained through MUSCL reconstruction and in time through a backward differentiation formula associated with a defect-correction approach. The implicit time advancing is compared in terms of accuracy and computational time with the explicit approach for one-dimensional and two-dimensional sediment transport problems, characterized by different time scales for the evolution of the bed and of the water flow. It is shown that, whenever the use of large time steps is compatible with the capture of the water flow dynamics and of the bedload evolution, the implicit scheme is far more efficient than its explicit counterpart with a CPU reduction up to more than two orders of magnitude. This makes implicit time differencing an attractive option for complex real life applications in this area.

[1]  T. Gallouët,et al.  Some approximate Godunov schemes to compute shallow-water equations with topography , 2003 .

[2]  Hervé Guillard,et al.  A Second Order Defect Correction Scheme for Unsteady Problems , 1996 .

[3]  Enrique D. Fernández-Nieto,et al.  Sediment transport models in Shallow Water equations and numerical approach by high order finite volume methods , 2008 .

[4]  M. Vázquez-Cendón.  Improved Treatment of Source Terms in Upwind Schemes for the Shallow Water Equations in Channels with Irregular Geometry , 1999 .

[5]  D. Causon,et al.  The surface gradient method for the treatment of source terms in the shallow-water equations , 2001 .

[6]  Luigi Fraccarollo,et al.  A well-balanced approach for flows over mobile-bed with high sediment-transport , 2006, J. Comput. Phys..

[7]  Dafydd Gibbon,et al.  1 User’s guide , 1998 .

[8]  Fayssal Benkhaldoun,et al.  Solution of the Sediment Transport Equations Using a Finite Volume Method Based on Sign Matrix , 2009, SIAM J. Sci. Comput..

[9]  T. C. Rebollo,et al.  A flux‐splitting solver for shallow water equations with source terms , 2003 .

[10]  Carlos Parés,et al.  Two-dimensional sediment transport models in shallow water equations. A second order finite volume approach on unstructured meshes , 2009 .

[11]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[12]  Enrique D. Fernández Nieto,et al.  A family of stable numerical solvers for the shallow water equations with source terms , 2003 .

[13]  Javier Murillo,et al.  Generalized Roe schemes for 1D two-phase, free-surface flows over a mobile bed , 2008, J. Comput. Phys..

[14]  A. I. Delis,et al.  Relaxation approximation to bed-load sediment transport , 2008 .

[15]  Manuel Jesús Castro Díaz,et al.  High Order Extensions of Roe Schemes for Two-Dimensional Nonconservative Hyperbolic Systems , 2009, J. Sci. Comput..

[16]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[17]  Justin Hudson,et al.  Formulations for Numerically Approximating Hyperbolic Systems Governing Sediment Transport , 2003, J. Sci. Comput..

[18]  Luca Bonaventura,et al.  An accurate and efficient semi‐implicit method for section‐averaged free‐surface flow modelling , 2011 .

[19]  Fayssal Benkhaldoun,et al.  A two‐dimensional finite volume morphodynamic model on unstructured triangular grids , 2009 .