Modeling bed erosion in free surface flows by the particle finite element method

We present a general formulation for modeling bed erosion in free surface flows using the particle finite element method (PFEM). The key feature of the PFEM is the use of an updated Lagrangian description to model the motion of nodes (particles) in domains containing fluid and solid subdomains. Nodes are viewed as material points (called particles) which can freely move and even separate from the fluid and solid subdomains representing, for instance, the effect of water drops or soil/rock particles. A mesh connects the nodes defining the discretized domain in the fluid and solid regions where the governing equations, expressed in an integral form, are solved as in the standard FEM. The necessary stabilization for dealing with the incompressibility of the fluid is introduced via the finite calculus (FIC) method. An incremental iterative scheme for the solution of the nonlinear transient coupled fluid-structure problem is described. The erosion mechanism is modeled by releasing the material adjacent to the bed surface according to the frictional work generated by the fluid shear stresses. The released bed material is subsequently transported by the fluid flow. Examples of application of the PFEM to solve a number of bed erosion problems involving large motions of the free surface and splashing of waves are presented.

[1]  Leo C. van Rijn,et al.  Mathematical Modeling of Suspended Sediment in NonUniform Flows , 1986 .

[2]  O. C. Zienkiewicz,et al.  The Finite Element Method for Fluid Dynamics , 2005 .

[3]  H. D. Vriend,et al.  Bed deformation in curved alluvial channels , 1985 .

[4]  E. Oñate,et al.  The particle finite element method. An overview , 2004 .

[5]  Eugenio Oñate,et al.  Combination of discrete element and finite element methods for dynamic analysis of geomechanics problems , 2004 .

[6]  D. M. Mcdowell.  Discussion of "Sediment Transport, Part III: Bed Forms and Alluvial Roughness" , 1986 .

[7]  James L. Smith,et al.  Compaction and water velocity effects on soil erosion in shallow flow , 1995 .

[8]  Eugenio Oñate,et al.  Polyhedrization of an arbitrary 3D point set , 2003 .

[9]  Gary Parker,et al.  A new vectorial bedload formulation and its application to the time evolution of straight river channels , 1994, Journal of Fluid Mechanics.

[10]  A. Huerta,et al.  Finite Element Methods for Flow Problems , 2003 .

[11]  Masato Sekine,et al.  Bed-load transport on transverse slope.I , 1992 .

[12]  J. Archard.  Contact and Rubbing of Flat Surfaces , 1953 .

[13]  Chi Fai Wan,et al.  Time for Development of Internal Erosion and Piping in Embankment Dams , 2003 .

[14]  Eugenio Oñate,et al.  A mesh-free finite point method for advective-diffusive transport and fluid flow problems , 1998 .

[15]  C. Fletcher.  On an alternating direction implicit finite element method for flow problems , 1982 .

[16]  A. Huerta,et al.  Finite Element Methods for Flow Problems , 2003 .

[17]  Eugenio Oñate,et al.  Derivation of stabilized equations for numerical solution of advective-diffusive transport and fluid flow problems , 1998 .

[18]  Eugenio Oñate,et al.  A Lagrangian meshless finite element method applied to fluid-structure interaction problems , 2003 .

[19]  Three-Dimensional Calculation of River Flow , 1997 .

[20]  F. Holly,et al.  Modeling of Riverbed Evolution for Bedload Sediment Mixtures , 1989 .

[21]  Eugenio Oñate,et al.  A finite element method for fluid-structure interaction with surface waves using a finite calculus formulation , 2001 .

[22]  Herbert Edelsbrunner,et al.  Three-dimensional alpha shapes , 1992, VVS.

[23]  E. Onate,et al.  An Unstructured Finite Element Solver for Ship Hydrodynamics Problems , 2003 .

[24]  O. C. Zienkiewicz,et al.  CBS versus GLS stabilization of the incompressible Navier–Stokes equations and the role of the time step as stabilization parameter , 2001 .

[25]  Eugenio Oñate,et al.  A finite point method for incompressible flow problems , 2000 .

[26]  Eugenio Oñate,et al.  Particle finite element method in fluid-mechanics including thermal convection-diffusion , 2005 .

[27]  Stephen E. Darby,et al.  Numerical Simulation of Widening and Bed Deformation of Straight Sand-Bed Rivers. II: Model Evaluation , 1996 .

[28]  Thomas J. R. Hughes,et al.  Encyclopedia of computational mechanics , 2004 .

[29]  O. C. Zienkiewicz,et al.  The Finite Element Method for Solid and Structural Mechanics , 2013 .

[30]  Alex J. Sutherland,et al.  Spatial lag effects in bed load sediment transport , 1989 .

[31]  Robin Fell,et al.  Investigation of Rate of Erosion of Soils in Embankment Dams , 2004 .

[32]  Eugenio Oñate,et al.  The meshless finite element method , 2003 .

[33]  Eugenio Oñate,et al.  Possibilities of finite calculus in computational mechanics , 2004 .

[34]  Eugenio Oñate,et al.  A stabilized finite element method for incompressible viscous flows using a finite increment calculus formulation , 2000 .

[35]  Eugenio Oñate,et al.  The particle finite element method: a powerful tool to solve incompressible flows with free‐surfaces and breaking waves , 2004 .