In this article we present a new class of high order accurate Arbitrary-Eulerian-Lagrangian (ALE) one-step WENO finite volume schemes for solving nonlinear hyperbolic systems of conservation laws on moving two dimensional unstructured triangular meshes. A WENO reconstruction algorithm is used to achieve high order accuracy in space and a high order one-step time discretization is achieved by using the local space-time Galerkin predictor proposed in [25]. For that purpose, a new element-local weak formulation of the governing PDE is adopted on moving space-time elements. The space-time basis and test functions are obtained considering Lagrange interpolation polynomials passing through a predefined set of nodes. Moreover, a polynomial mapping defined by the same local space-time basis functions as the weak solution of the PDE is used to map the moving physical space-time element onto a space-time reference element. To maintain algorithmic simplicity, the final ALE one-step finite volume scheme uses moving triangular meshes with straight edges. This is possible in the ALE framework, which allows a local mesh velocity that is different from the local fluid velocity. We present numerical convergence rates for the schemes presented in this paper up to sixth order of accuracy in space and time and show some classical numerical test problems for the two-dimensional Euler equations of compressible gas dynamics.

In this paper we present a new family of high order accurate Arbitrary-LagrangianEulerian (ALE) one-step ADER-WENO finite volume schemes for the solution of nonlinear systems of conservative and non-conservative hyperbolic partial differential equations with stiff source terms on moving tetrahedral meshes in three space dimensions. A WENO reconstruction technique is used to achieve high order of accuracy in space, while an element-local space-time Discontinuous Galerkin finite element predictor on moving meshes is used to obtain a high order accurate one-step time discretization. Within the space-time predictor the physical element is mapped onto a reference element using an isoparametric approach, where the space-time basis and test functions are given by the Lagrange interpolation polynomials passing through a predefined set of space-time nodes. Since our algorithm is cell-centered, the final mesh motion is computed by using a suitable node solver algorithm. A rezoning step as well as a flattener strategy are used in some of the test problems to avoid mesh tangling or excessive element deformations that may occur when the computation involves strong shocks or shear waves. The ALE algorithm presented in this article belongs to the so-called direct ALE methods because the final Lagrangian finite volume scheme is based directly on a space-time conservation formulation of the governing PDE system, with the rezoned geometry taken already into account during the computation of the fluxes.We apply our new high order unstructured ALE schemes to the 3D Euler equations of compressible gas dynamics, for which a set of classical numerical test problems has been solved and for which convergence rates up to sixth order of accuracy in space and time have been obtained. We furthermore consider the equations of classical ideal magnetohydrodynamics (MHD) as well as the non-conservative seven-equation Baer-Nunziato model of compressible multi-phase flows with stiff relaxation source terms.

In this article we present the first better than second order accurate unstructured Lagrangian-type one-step WENO finite volume scheme for the solution of hyperbolic partial differential equations with non-conservative products. The method achieves high order of accuracy in space together with essentially non-oscillatory behaviour using a nonlinear WENO reconstruction operator on unstructured triangular meshes. High order accuracy in time is obtained via a local Lagrangian space-time Galerkin predictor method that evolves the spatial reconstruction polynomials in time within each element. The final one-step finite volume scheme is derived by integration over a moving space-time control volume, where the nonconservative products are treated by a path-conservative approach that defines the jump terms on the element boundaries. The entire method is formulated as an Arbitrary-Lagrangian-Eulerian (ALE) method, where the mesh velocity can be chosen independently of the fluid velocity.The new scheme is applied to the full seven-equation Baer-Nunziato model of compressible multi-phase flows in two space dimensions. The use of a Lagrangian approach allows an excellent resolution of the solid contact and the resolution of jumps in the volume fraction. The high order of accuracy of the scheme in space and time is confirmed via a numerical convergence study. Finally, the proposed method is also applied to a reduced version of the compressible Baer-Nunziato model for the simulation of free surface water waves in moving domains. In particular, the phenomenon of sloshing is studied in a moving water tank and comparisons with experimental data are provided.

SUMMARYIn this paper we present a class of high order accurate cell-centered Arbitrary-Eulerian-Lagrangian (ALE) one-step ADER-WENO finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two-dimensional unstructured triangular meshes along the principles laid out in [1]. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space-time Galerkin predictor allows the schemes to be high order accurate also in time by using an element-local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained i) either as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu [2], or, ii) as a solution of multiple one-dimensional half-Riemann problems around a vertex, as suggested by Maire [3], or, iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional HLL Riemann solver recently proposed by Balsara et al. [4]. Once the vertex velocity and thus the new node location has been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level t n with the new ones at the next time level t n+1 . If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space-time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unneccesary, since the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic MHD equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration.

No abstract