A finite element method for localized erosion in porous media with applications to backward piping in levees

This paper presents a computational method able to effectively model both the simultaneous processes typically observed in backward erosion piping, ie, the pipe tip propagation and the conduit cross‐section enlargement. The numerical method is based on the novel formulation of a problem of localized erosion along a line propagating in a multidimensional porous medium. In this line, a conduit with evolving transverse size is embedded, which conveys a multiphase flow. The two systems, porous medium and pipe, are bridged by exchange terms of multiphase fluid mass and by a shared fluid pressure field. On the contrary, different fields are considered to describe flows, which are assumed as Darcian in the porous medium and turbulent in the conduit. These two flows drive pipe propagation and enlargement, respectively, as modeled by means of proper erosion kinetic laws. The corresponding numerical formulation is based on the combination between one‐ and multidimensional finite elements, to model the erosion conduit and the porous medium, respectively. Several simulations are proposed to demonstrate the ability of the proposed approach in reproducing available experimental data of real‐scale tests on levees. Our results point out the crucial role played by the combined influence of pipe propagation and enlargement, as well as of three‐dimensional (3D) effects. We also assess the mesh independence of the proposed numerical solution, particularly as concerns the calculated pipe propagation history.