Augmented Lagrangian for shallow viscoplastic flow with topography

In this paper we have developed a robust numerical algorithm for the visco-plastic Saint-Venant model with topography. For the time discretization an implicit (backward) Euler scheme was used. To solve the resulting nonlinear equations, a four steps iterative algorithm was proposed. To handle the non-differentiability of the plastic terms an iterative decomposition-coordination formulation coupled with the augmented Lagrangian method was adopted. The proposed algorithm is consistent, i.e. if the convergence is achieved then the iterative solution satisfies the nonlinear system at each time iteration. The equations for the velocity field are discretized using the finite element method, while a discontinuous Galerkin method, with an upwind choice of the flux, is adopted for solving the hyperbolic equations that describe the evolution of the thickness. The algorithm permits to solve alternatively, at each iteration, the equations for the velocity field and for the thickness. The iterative decomposition coordination formulation coupled with the augmented Lagrangian method works very well and no instabilities are present. The proposed algorithm has a very good convergence rate, with the exception of large Reynolds numbers ([email protected]?1000), not involved in the applications concerned by the shallow viscoplastic model. The discontinuous Galerkin technique assure the mass conservation of the shallow system. The model has the exact C-property for a plane bottom and an asymptotic C-property for a general topography. Some boundary value problems were selected to analyze the robustness of the numerical algorithm and the predictive capabilities of the mechanical model. The comparison with an exact rigid flow solution illustrates the accuracy of the numerical scheme in handling the non-differentiability of the plastic terms. The influence of the mesh and of the time step are investigated for the flow of a Bingham fluid in a talweg. The role of the material cohesion in stopping a viscoplastic avalanche on a talweg with barrier was analyzed. Finally, the capacities of the model to describe the flow of a Bingham fluid on a valley from the broken wall of a reservoir situated upstream were investigated.

[1]  Kolumban Hutter,et al.  The Savage‐Hutter theory: A system of partial differential equations for avalanche flows of snow, debris, and mud , 2004 .

[2]  S. Bankoff,et al.  Long-scale evolution of thin liquid films , 1997 .

[3]  F. Bouchut,et al.  On new erosion models of Savage–Hutter type for avalanches , 2008 .

[4]  Chiang C. Mei,et al.  Approximate equations for the slow spreading of a thin sheet of Bingham plastic fluid , 1990 .

[5]  R. Glowinski,et al.  Méthodes de Lagrangien augmenté : applications à la résolution numérique de problèmes aux limites , 1982 .

[6]  R. Craster,et al.  Viscoplastic flow over an inclined surface (Erratum) , 2007 .

[7]  Kolumban Hutter,et al.  Channelized free-surface flow of cohesionless granular avalanches in a chute with shallow lateral curvature , 1999, Journal of Fluid Mechanics.

[8]  B. Perthame,et al.  A new model of Saint Venant and Savage–Hutter type for gravity driven shallow water flows , 2003 .

[9]  Enrique D. Fernández-Nieto,et al.  Augmented Lagrangian Method and Compressible Visco-plastic Flows: Applications to Shallow Dense Avalanches , 2009 .

[10]  E. C. Bingham Fluidity And Plasticity , 1922 .

[11]  Luka Sopta,et al.  ENO and WENO Schemes with the Exact Conservation Property for One-Dimensional Shallow Water Equations , 2002 .

[12]  Olivier Pouliquen,et al.  A constitutive law for dense granular flows , 2006, Nature.

[13]  M. Gauthier,et al.  Rheological detection of caging and solid–liquid transitions in soft colloid–polymer mixtures , 2013 .

[14]  S. Savage,et al.  The dynamics of avalanches of granular materials from initiation to runout. Part I: Analysis , 1991 .

[15]  Amy Q. Shen,et al.  Visco-plastic models of isothermal lava domes , 2000, Journal of Fluid Mechanics.

[16]  Richard V. Craster,et al.  A consistent thin-layer theory for Bingham plastics , 1999 .

[17]  I. Ionescu Viscoplastic shallow flow equations with topography , 2013 .

[18]  D. C. Drucker,et al.  Soil mechanics and plastic analysis or limit design , 1952 .

[19]  S. Pudasaini,et al.  Rapid shear flows of dry granular masses down curved and twisted channels , 2003, Journal of Fluid Mechanics.

[20]  I. Ionescu Onset and dynamic shallow flow of a viscoplastic fluid on a plane slope , 2010 .

[21]  Kolumban Hutter,et al.  Avalanche dynamics: Dynamics of rapid flows of dense granular avalanches , 2016 .

[22]  C. Mei,et al.  Slow spreading of a sheet of Bingham fluid on an inclined plane , 1989, Journal of Fluid Mechanics.

[23]  Richard M. Iverson,et al.  Flow of variably fluidized granular masses across three‐dimensional terrain: 1. Coulomb mixture theory , 2001 .

[24]  Jean-Pierre Vilotte,et al.  Numerical modeling of avalanches based on Saint-Venant equations using a kinetic scheme , 2003 .

[25]  Roland Glowinski,et al.  On the numerical simulation of Bingham visco-plastic flow: Old and new results , 2007 .

[26]  Marina Pirulli,et al.  The effect of the earth pressure coefficients on the runout of granular material , 2007, Environ. Model. Softw..

[27]  S. Pudasaini,et al.  Gravity-Driven Rapid Shear Flows of Dry Granular Masses in Topographies with Orthogonal and Non-Orthogonal Metrics , 2003 .

[28]  Christophe Ancey,et al.  Plasticity and geophysical flows: A review , 2007 .

[29]  R. Craster,et al.  Shallow viscoplastic flow on an inclined plane , 2002, Journal of Fluid Mechanics.

[30]  R. Glowinski,et al.  Augmented Lagrangian and Operator-Splitting Methods in Nonlinear Mechanics , 1987 .

[31]  Enrique D. Fernández-Nieto,et al.  Shallow Water equations for Non-Newtonian fluids , 2010 .

[32]  Joseph B. Keller Shallow-water theory for arbitrary slopes of the bottom , 2003, Journal of Fluid Mechanics.

[33]  R. Craster,et al.  Viscoplastic flow over an inclined surface , 2006 .

[34]  Douglas R. Macayeal,et al.  Large‐scale ice flow over a viscous basal sediment: Theory and application to ice stream B, Antarctica , 1989 .

[35]  R. Dressler,et al.  NEW NONLINEAR SHALLOW-FLOW EQUATIONS WITH CURVATURE , 1978 .

[36]  Jean-Pierre Vilotte,et al.  On the use of Saint Venant equations to simulate the spreading of a granular mass , 2005 .

[37]  Kolumban Hutter,et al.  Gravity-driven free surface flow of granular avalanches over complex basal topography , 1999, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[38]  S. Savage,et al.  The motion of a finite mass of granular material down a rough incline , 1989, Journal of Fluid Mechanics.