Three dimensional modeling of free surface flow and sediment transport with bed deformation using automatic mesh motion

This study presents the development of a 3D numerical code for flow-sediment interactions with associated bed changes in free surface flows. To capture the water-air interface, a novel volume of fluid (VOF) formulation is implemented using the ghost fluid method with one-side extrapolation for dynamic pressure. Equations for fluid flow and free surface motion are decoupled due to separation of time scales. Accordingly, time step size can be increased by a factor of 3-orders of magnitude and still preserve numerical stability and computational efficiency. A novel finite area method (FAM) is utilized to discretize Exner equation on irregular bed surface providing a 3D finite volume-like discretization on curved surfaces. The evolution of the water-sediment interface is captured by a novel vertex-based unstructured mesh dynamic motion solve using Laplace operator with variable diffusivity. The code is implemented in foam-extend and tested against two classical experiments. Good results are obtained with correct trends and lower absolute error compared to previous mobile bed models. This shows that the developed code has a good potential of being applied to mobile-bed hydraulic real problems.

[1]  N. Olsen,et al.  Three-dimensional numerical modeling of water and sediment flow in a sand trap , 1994 .

[2]  R. T. Cheng,et al.  SEMI-IMPLICIT FINITE DIFFERENCE METHODS FOR THREE-DIMENSIONAL SHALLOW WATER FLOW , 1992 .

[3]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[4]  Qiuhua Liang,et al.  A two-dimensional hydro-morphological model for river hydraulics and morphology with vegetation , 2017, Environ. Model. Softw..

[5]  W. Rodi,et al.  3D NUMERICAL MODELING OF FLOW AND SEDIMENT TRANSPORT IN OPEN CHANNELS , 2000 .

[6]  A. Probst,et al.  Modelling trace metal transfer in large rivers under dynamic hydrology: A coupled hydrodynamic and chemical equilibrium model , 2017, Environ. Model. Softw..

[7]  Marian Muste,et al.  Two-phase formulation of suspended sediment transport , 1999 .

[8]  Paolo Burlando,et al.  Coupling spatially distributed river and groundwater transport models to investigate contaminant dynamics at river corridor scales , 2016, Environ. Model. Softw..

[9]  Jan S. Ribberink,et al.  The validity of a depth-integrated model for suspended sediment transport , 1986 .

[10]  Alessandro Dal Palù,et al.  A non-uniform efficient grid type for GPU-parallel Shallow Water Equations models , 2017, Environ. Model. Softw..

[11]  Roger Alexander Falconer,et al.  Numerical modelling of three-dimensional suspended sediment for estuarine and coastal waters , 1996 .

[12]  H. Jasak,et al.  Automatic mesh motion for the unstructured finite volume method , 2007 .

[13]  Jørgen Fredsøe,et al.  A Sediment Transport Model for Straight Alluvial Channels , 1976 .

[14]  O. Ubbink.  Numerical prediction of two fluid systems with sharp interfaces , 1997 .

[15]  Weiming Wu.  Computational River Dynamics , 2007 .

[16]  Hrvoje Jasak,et al.  Error analysis and estimation for the finite volume method with applications to fluid flows , 1996 .

[17]  Michele Guidolin,et al.  A weighted cellular automata 2D inundation model for rapid flood analysis , 2016, Environ. Model. Softw..

[18]  F. Menter.  Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[19]  H. Jasak,et al.  Implementation of the Ghost Fluid Method for Free Surface Flows in Polyhedral Finite Volume Framework , 2017 .

[20]  A. D. Gosman,et al.  AUTOMATIC RESOLUTION CONTROL FOR THE FINITE-VOLUME METHOD, PART 1: A-POSTERIORI ERROR ESTIMATES , 2000 .

[21]  Hrvoje Jasak,et al.  Application of the finite volume method and unstructured meshes to linear elasticity , 2000 .

[22]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[23]  Hrvoje Jasak,et al.  Rapid Free Surface Simulation for Steady-State Hull Resistance with FVM using OpenFOAM , 2014 .

[24]  H. Jasak,et al.  A moving mesh finite volume interface tracking method for surface tension dominated interfacial fluid flow , 2012 .

[25]  Forrest M. Holly,et al.  Three-dimensional numerical simulation of mobile-bed hydrodynamics , 1994 .

[26]  L. Rijn.  Sediment Transport, Part II: Suspended Load Transport , 1984 .

[27]  Zhixian Cao,et al.  Sediment-Laden Flow in Open Channels from Two-Phase Flow Viewpoint , 1995 .

[28]  Xiaofeng Liu,et al.  Three-Dimensional Numerical Model with Free Water Surface and Mesh Deformation for Local Sediment Scour , 2008 .

[29]  B. Brørs.  Numerical Modeling of Flow and Scour at Pipelines , 1999 .

[30]  H. Rusche.  Computational fluid dynamics of dispersed two-phase flows at high phase fractions , 2003 .

[31]  H. Fang,et al.  Three-dimensional calculations of flow and suspended sediment transport in the neighborhood of the dam for the Three Gorges Project (TGP) reservoir in the Yangtze River , 2003 .