DIKUL - logo
E-viri
Recenzirano Odprti dostop
  • Direct Arbitrary-Lagrangian...
    Gaburro, Elena; Dumbser, Michael; Castro, Manuel J.

    Computers & fluids, 12/2017, Letnik: 159
    Journal Article

    •Novel high resolution finite volume schemes on moving nonconforming unstructured meshes.•Significantly improved robustness in shear flows compared to conforming ALE schemes.•Locally and globally conservative ALE schemes that respect the geometric conservation law (GCL).•The fully discrete one-step scheme can be derived from a space-time conservation formalism.•Well-balanced for non-trivial stationary solutions of the governing PDE. In this paper, we present a novel second-order accurate Arbitrary-Lagrangian-Eulerian (ALE) finite volume scheme on moving nonconforming polygonal grids, in order to avoid the typical mesh distortion caused by shear flowsin Lagrangian-type methods. In our new approach the nonconforming element interfaces are not defined by the user, but they are automatically detected by the algorithm if the tangential velocity difference across an element interface is sufficiently large. The grid nodes that are sufficiently far away from a shear wave are moved with a standard node solver, while at the interface we insert a new set of nodes that can slide along the interface in a nonconforming manner. In this way, the elements on both sides of the shear wave can move with a different velocity, without producing highly distorted elements. The core of the proposed method is the use of a space-time conservation formulation in the construction of the final finite volume scheme, which completely avoids the need of an additional remapping stage, hence the new method is a so-called direct ALE scheme. For this purpose, the governing PDE system is rewritten at the aid of the space-time divergence operator and then a fully discrete one-step discretization is obtained by integrating over a set of closed space-time control volumes. The nonconforming sliding of nodes along an edge requires the insertion or the deletion of nodes and edges, and in particular the space-time faces of an element can be shared between more than two cells. Due to the space-time conservation formulation, the geometric conservation law (GCL) is automatically satisfied by construction, even on moving nonconforming meshes. Moreover, the mesh quality remains high and, as a direct consequence, also the time step remains almost constant in time, even for highly sheared vortex flows. In this paper we focus mainly on logically straight slip-line interfaces, but we show also first results for general slide lines that are not logically straight. Second order of accuracy in space and time is obtained by using a MUSCL-Hancock strategy, together with a Barth and Jespersen slope limiter. The accuracy of the new scheme has been further improved by incorporating a special well balancing technique that is able to maintain particular stationary solutions of the governing PDE system up to machine precision. In particular, we consider steady vortex solutions of the shallow water equations, where the pressure gradient is in equilibrium with the centrifugal forces. A large set of different numerical tests has been carried out in order to check the accuracy and the robustness of the new method for both smooth and discontinuous problems. In particular we have compared the results for a steady vortex in equilibrium solved with a standard conforming ALE method (without any rezoning technique) and with our new nonconforming ALE scheme, to show that the new nonconforming scheme is able to avoid mesh distortion in vortex flows even after very long simulation times.