r-adaptation for Shallow Water flows: conservation, well balancedness, efficiency

Luca Arpaia (1) , Mario Ricchiuto (1)

1 CARDAMOM - Certified Adaptive discRete moDels for robust simulAtions of CoMplex flOws with Moving fronts

We investigate the potential of the so-called "relocation" mesh adaptation in terms of resolution and efficiency for the simulation of free surface flows in the near shore region. Our work is developed in three main steps. First, we consider several Arbitrary Lagrangian Eulerian (ALE) formulations of the shallow water equations on moving grids, and provide discrete analogs in the Finite Volume and Residual Distribution framework. The compliance to all the physical constraints, often in competition, is taken into account. We consider different formulations allowing to combine volume conservation (DGCL) and equilibrium (Well-Balancedness), and we clarify the relations between the so-called pre-balanced form of the equations (Rogers et al., J.Comput.Phys. 192, 2003), and the classical upwiniding of the bathymetry term gradients (Bermudez and Vazquez, Computers and Fluids 235, 1994). Moreover, we propose a simple remap of the bathymetry based on high accurate quadrature on the moving mesh which, while preserving an accurate representation of the initial data, also allows to retain mass conservation within an arbitrary accuracy. Second, the coupling of the resulting schemes with a mesh partial differential equation is studied. Since the flow solver is based on genuinely explicit time stepping, we investigate the efficiency of three coupling strategies in terms of cost overhead w.r.t. the flow solver. We analyze the role of the solution remap necessary to evaluate the error monitor controlling the adaptation, and propose simplified formulations allowing a reduction in computational cost. The resulting ALE algorithm is compared with the \textit{rezoning} Eulerian approach with interpolation proposed e.g. in (Tang and Tang, SINUM 41, 2003). An alternative cost effective Eulerian approach, still allowing a full decoupling between adaptation and flow evolution steps is also proposed. Finally, a thorough numerical evaluation of the methods discussed is performed. Numerical results on propagation, and inundation problems shows that the best compromise between accuracy and CPU time is provided by a full ALE formulation. If a loose coupling with the mesh adaptation is sought, then the cheaper Eulerian approach proposed is shown to provide results quite close to the ALE.

Mots clés: Mesh adaptation; shallow water equations; run-up; mesh deformation; ALE formulation; well balanced; mass conservation; efficiency

Domaines: Modélisation et simulation; Analyse numérique [math.NA]; Mécanique des fluides [physics.class-ph]; Océan, Atmosphère; Mécanique des fluides [physics.class-ph]

Luca Arpaia, Mario Ricchiuto. r-adaptation for Shallow Water flows: conservation, well balancedness, efficiency. [Research Report] RR-8956, Inria Bordeaux Sud-Ouest. 2016. ⟨hal-01372496v2⟩