Published December 2019 | Version v1
Journal article

High order direct Arbitrary-Lagrangian-Eulerian (ALE) P N P M schemes with WENO Adaptive-Order reconstruction on unstructured meshes

  • 1. Department of Mathematics and Computer Science, University of Ferrara, Via Machiavelli 30, 44121 Ferrara (Italy)
  • 2. Physics Department, University of Notre Dame, 225 Nieuwland Science Hall, Notre Dame, IN 46556 (United States)

Description

Highlights: • High order PNPM schemes on moving unstructured meshes. • PNPM schemes applied to discontinuous solutions with N>0. • WENO Adaptive-Order reconstruction for reconstructed DG schemes. • WENO Adaptive-Order strategy as DG limiter. • Shock detector for DG limiting. -- Abstract: In this work we present a conservative WENO Adaptive Order (AO) reconstruction operator applied to an explicit one-step Arbitrary-Lagrangian-Eulerian (ALE) discontinuous Galerkin (DG) method. The spatial order of accuracy is improved by reconstructing higher order piecewise polynomials of degree M>N, starting from the underlying polynomial solution of degree N provided by the DG scheme. High order of accuracy in time is achieved by the ADER approach, making use of an element-local space-time Galerkin finite element predictor that arises from a one-step time integration procedure. As a result, space-time polynomials of order M+1 are obtained and used to perform the time evolution of the numerical solution adopting a fully explicit DG scheme. To maintain algorithm simplicity, the mesh motion is restricted to be carried out using straight lines, hence the old mesh configuration at time tn is connected with the new one at time tn+1 via space-time segments, which result in space-time control volumes on which the governing equations have to be integrated in order to obtain the time evolution of the discrete solution. Our algorithm falls into the category of direct Arbitrary-Lagrangian-Eulerian (ALE) schemes, where the governing PDE system is directly discretized relying on a space-time conservation formulation and which already takes into account the new grid geometry directly during the computation of the numerical fluxes. A local rezoning strategy might be used in order to locally optimize the mesh quality and avoiding the generation of invalid elements with negative determinant. The proposed approach reduces to direct ALE finite volume schemes if N=0, while explicit direct ALE DG schemes are recovered in the case of N=M. In order to stabilize the DG solution, an a priori WENO based limiting technique is employed, that makes use of the numerical solution inside the element under consideration and its neighbor cells to find a less oscillatory polynomial approximation. By using a modal basis in a reference element, the evaluation of the oscillation indicators is very easily and efficiently carried out, hence allowing higher order modes to be properly limited, while leaving the zero-th order mode untouched for ensuring conservation. Numerical convergence rates for 2N,M4 are presented as well as a wide set of benchmark test problems for hydrodynamics on moving and fixed unstructured meshes.

Availability note (English)

Available from http://dx.doi.org/10.1016/j.jcp.2019.108899

Additional details

Identifiers

DOI
10.1016/j.jcp.2019.108899;
PII
S0021999119305972;

Publishing Information

Journal Title
Journal of Computational Physics (Print)
Journal Volume
398
Journal Page Range
vp.
ISSN
0021-9991
CODEN
JCTPAH

INIS

Country of Publication
Netherlands
Country of Input or Organization
International Atomic Energy Agency (IAEA)
INIS RN
56005746
Subject category
S71: CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS;
Descriptors DEI
ALGORITHMS; FINITE ELEMENT METHOD; LAGRANGIAN FUNCTION; PARTIAL DIFFERENTIAL EQUATIONS; POLYNOMIALS
Descriptors DEC
CALCULATION METHODS; DIFFERENTIAL EQUATIONS; EQUATIONS; FUNCTIONS; MATHEMATICAL LOGIC; MATHEMATICAL SOLUTIONS; NUMERICAL SOLUTION

Optional Information

Copyright
Copyright (c) 2019 Elsevier Inc. All rights reserved.