We investigate the potential of the so-called "relocation" mesh adaptation in terms of resolution and e ciency 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 di↵erent 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 di↵erential equation is studied. Since the flow solver is based on genuinely explicit time stepping, we investigate the e ciency 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 rezoning Eulerian approach with interpolation proposed e.g. in (Tang and Tang, SINUM 41, 2003). An alternative cost e↵ective 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.

No abstract

In this paper we consider the solution of hyperbolic conservation laws on moving meshes by means of an Arbitrary Lagrangian Eulerian (ALE) formulation of the Runge-Kutta Residual Distribution (RD) schemes of Ricchiuto and Abgrall (J Comput Phys 229(16):5653-5691, 2010). Up to the authors knowledge, the problem of recasting RD schemes into ALE framework has been solved with first order explicit schemes and with second order implicit schemes. Our resulting scheme is explicit and second order accurate when computing discontinuous solutions.

We consider the numerical approximation of the Shallow Water Equations (SWEs) in spherical geometry for oceanographic applications. To provide enhanced resolution of moving fronts present in the flow we consider adaptive discrete approximations on moving triangulations of the sphere. To this end, we restate all Arbitrary Lagrangian Eulerian (ALE) transport formulas, as well as the volume transformation laws, for a 2D manifold. Using these results, we write the set of ALE-SWEs on the sphere. We then propose a Residual Distribution discrete approximation of the governing equations. Classical properties as the DGCL and the C-property (well balancedness) are reformulated in this more general context. An adaptive mesh movement strategy is proposed. The discrete framework obtained is thoroughly tested on standard benchmarks in large scale oceanography to prove their potential as well as the advantage brought by the adaptive mesh movement.