High order direct Arbitrary-Lagrangian-Eulerian (ALE) P N P M schemes with WENO Adaptive-Order reconstruction on unstructured meshes
Creators
- 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 schemes on moving unstructured meshes. • schemes applied to discontinuous solutions with . • 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 , 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 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 is connected with the new one at time 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 , while explicit direct ALE DG schemes are recovered in the case of . 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 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.108899Additional 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.