High order cell-centered Lagrangian-type finite volume schemes with time-accurate local time stepping on unstructured triangular meshes
Introduction
In the last few years there has been a renewed interest in the development of novel accurate and robust cell-centered Lagrangian finite volume schemes for hydrodynamics. Since in the Lagrangian framework the computational mesh moves with the local fluid velocity, such schemes are regarded as the first choice in all problems presenting moving material interfaces appearing in compressible multi-phase and multi-material flows, such as in the numerical simulation of inertial confinement fusion (ICF). The vast majority of modern Lagrangian schemes adopt a cell-centered finite-volume approach, see for example [21], [29], [99], [86], [90], [91], [92], [89], where all flow variables are defined as cell-averaged quantities inside a control volume. However, also staggered Lagrangian schemes are possible, see e.g. [84], where the velocity is defined at the cell interfaces, while the other flow variables are still defined at the cell centers.
In [95], [21] Godunov-type finite volume schemes have been presented for Lagrangian hydrodynamics, while in [33], [34] the governing equations have been coupled with the equations for the evolution of the geometry and the resulting weakly hyperbolic system has been solved using a node-based finite volume solver. Unstructured multidimensional meshes have been considered by Maire in [87], [89], [88], who developed up to second order accurate cell-centered Lagrangian schemes, where the time derivatives of the fluxes have been computed with a node-centered solver. This approach may be regarded as a multi-dimensional Lagrangian extension of the Generalized Riemann problem methodology used for example in the ADER approach of Titarev and Toro [105], [107] in the Eulerian context. Arbitrary-Lagrangian–Eulerian (ALE) methods based on remeshing and remapping have also been investigated very recently for single and multi-material flows in [63], [112], [19], [100].
In [27], [81] Cheng and Shu presented the first better than second order accurate Lagrangian schemes for hydrodynamics on structured meshes, where the use of a high order Essentially Non-Oscillatory (ENO) reconstruction operator yielded high order of accuracy in space, while high order of accuracy in time was guaranteed using either a Runge–Kutta or a Lax–Wendroff-type time stepping. Arbitrary high order accurate cell-centered Lagrangian-type finite volume schemes for conservative and non-conservative hyperbolic PDEs on moving unstructured triangular and tetrahedral meshes have been considered for the first time by Boscheri et al. in a very recent series of papers [15], [48], [17], [14], [16]. A new class of meshless Lagrangian particle methods based on a high order accurate moving least-squares WENO reconstruction has been forwarded in [3].
High order accurate Lagrangian algorithms using the classical continuous finite element method (FEM) can be found, for example, in the work of Scovazzi et al. [97], [101] and Dobrev et al. [37], [38], [39], while Lagrangian discontinuous Galerkin finite elements have been recently proposed by Vilar et al. and Yu et al. in [66], [64], [65], [80]. Arbitrary-Lagrangian–Eulerian DG schemes have been developed and applied, for example, in [26], [60].
Almost all of the above mentioned algorithms use an explicit global time stepping scheme in which the timestep is computed under a classical global CFL stability condition, so that the timestep is essentially determined by the smallest control volume appearing in the mesh. In Lagrangian hydrodynamics, where the mesh follows as closely as possible the local fluid motion, very severe deformations and distortions may occur in the computational cells, especially at shocks and shear waves. As a consequence, the computational efficiency of the algorithm drastically decreases, because the smallest timestep imposed by the most deformed control volumes dictates the timestep for the entire computational grid, including those elements which are much bigger or which lie in a zone where the fluid is moving uniformly. In the Eulerian framework such a problem can be partially avoided controlling the mesh quality a priori and designing a high quality mesh at the pre-processing level, since the grid will not change anymore during the simulation. Of course, the CFL condition can be circumvented by using implicit or semi-implicit schemes, see for example [23], [25], [24], [40], [41], [42], but this approach does not yet seem to be very popular in the context of cell-centered Lagrangian-type finite volume methods. An alternative to overcome the global CFL condition consists in the development of numerical schemes that allow for time-accurate local time stepping (LTS), where each element has to obey only a less restrictive local CFL stability condition, hence using its own optimal local timestep. Therefore, many efforts have been devoted to the construction of high order accurate Eulerian schemes with time-accurate LTS, developing either discontinuous Galerkin finite element methods [62], [53], [104], [82], [67], [79], [61] or high order accurate finite volume schemes with LTS [12], [11], [111], [22], [5], [4], [20], [61], [58], [49]. The finite volume schemes with LTS adopt mainly classical adaptive mesh refinement (AMR) techniques in space and time or block-clustered local time stepping algorithms. In [70], [69] also high order accurate Runge–Kutta time integrators with local time stepping (so-called multi-rate integrators) can be found. To the knowledge of the authors, the first high order accurate Lagrangian algorithm with time accurate local time stepping on moving grids has been proposed only very recently in [46], where the equations of hydrodynamics and of classical magnetohydrodynamics (MHD) have been solved in one spatial dimension. In the present paper we extend the algorithm presented in [46] to moving unstructured triangular meshes.
The rest of the paper is structured as follows: in Section 2 the numerical scheme is described, including the details of the local time stepping algorithm on moving unstructured meshes, while numerical convergence studies as well as some classical numerical test problems for hydrodynamics are presented in Section 3. We conclude the paper giving an outlook to future research and developments in Section 4.
Section snippets
Formulation of the equations and basic set-up
In this article we consider nonlinear hyperbolic conservation laws, cast in the form where is the vector of conserved variables defined in the space of the admissible states , denotes the nonlinear flux tensor and represents a nonlinear algebraic source term, which is not allowed to be stiff. The system of equations (1) is defined in two space dimensions, hence addressing the space coordinate vector and
Test problems
In the following we solve some numerical test problems in order to validate the high order Lagrangian ADER-WENO algorithm with time accurate local time stepping (LTS) presented so far. We consider the two-dimensional Euler equations of compressible gas dynamics, which can be cast into form (1) with where the vector of conserved variables is denoted by Q and the flux tensor is addressed with . Furthermore let ρ and
Conclusions
In this article we have presented a high order Lagrangian-type finite volume scheme with time-accurate local time stepping (LTS) on moving unstructured triangular meshes. The numerical scheme is derived from [15], [48], where a classical global time stepping approach is adopted, and from the recently developed one-dimensional ALE-ADER-WENO scheme [46]. In our approach the WENO reconstruction technique is used to achieve high order of accuracy in space, while high order of accuracy in time is
Acknowledgements
The presented research has been financed by the European Research Council (ERC) under the European Union's Seventh Framework Programme (FP7/2007-2013) with the research project STiMulUs, ERC Grant agreement No. 278267. The authors acknowledge PRACE for awarding us access to the SuperMUC supercomputer of the Leibniz Rechenzentrum (LRZ) in Munich, Germany.
The authors would like to thank the two anonymous referees for their constructive comments and suggestions, which definitely helped to improve
References (116)
- et al.
A new class of moving-least-squares WENO SPH schemes
J. Comput. Phys.
(2014) - et al.
Adaptation based on interpolation errors for high order mesh refinement methods applied to conservation laws
Appl. Numer. Math.
(2012) - et al.
Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy
J. Comput. Phys.
(2000) Multidimensional HLLE Riemann solver: application to Euler and magnetohydrodynamic flows
J. Comput. Phys.
(2010)A two-dimensional HLLC Riemann solver for conservation laws: application to Euler and magnetohydrodynamic flows
J. Comput. Phys.
(2012)- et al.
Local adaptive mesh refinement for shock hydrodynamics
J. Comput. Phys.
(May 1989) - et al.
Adaptive mesh refinement for hyperbolic partial differential equations
J. Comput. Phys.
(March 1984) - et al.
Lagrangian ADER-WENO finite volume schemes on unstructured triangular meshes based on genuinely multidimensional HLL Riemann solvers
J. Comput. Phys.
(2014) - et al.
Lagrangian ADER-WENO finite volume schemes on unstructured tetrahedral meshes for conservative and nonconservative hyperbolic systems in 3D
J. Comput. Phys.
(2014) - et al.
A multi-material ReALE method with MOF interface reconstruction
Comput. Fluids
(2013)
A cell-centered Lagrangian hydrodynamics scheme on general unstructured meshes in arbitrary dimension
J. Comput. Phys.
Simulation of compressible viscous flow in time-dependent domains
Appl. Math. Comput.
A high order ENO conservative Lagrangian type scheme for the compressible Euler equations
J. Comput. Phys.
A high-order finite volume method for systems of conservation laws—multi-dimensional optimal order detection (MOOD)
J. Comput. Phys.
A new exceptional points method with application to cell-centered Lagrangian schemes and curved meshes
J. Comput. Phys.
Efficient solution algorithms for the Riemann problem for real gases
J. Comput. Phys.
Hyperbolic divergence cleaning for the MHD equations
J. Comput. Phys.
Symmetrization of Lagrangian gas dynamic in dimension two and multidimensional solvers
C. R., Méc.
A message passing benchmark for unbalanced applications
Simul. Model. Pract. Theory
Improved detection criteria for the multi-dimensional optimal order detection (MOOD) on unstructured meshes with very high-order polynomials
Comput. Fluids
High order curvilinear finite elements for axisymmetric Lagrangian hydrodynamics
Comput. Fluids
A semi-implicit discontinuous Galerkin finite element method for the numerical solution of inviscid compressible flow
J. Comput. Phys.
Analysis of semi-implicit DGFEM for nonlinear convection–diffusion problems on nonconforming meshes
Comput. Methods Appl. Mech. Eng.
Multidimensional HLLC Riemann solver for unstructured meshes – with application to Euler and MHD flows
J. Comput. Phys.
Arbitrary-Lagrangian–Eulerian ADER-WENO finite volume schemes with time-accurate local time stepping for hyperbolic conservation laws
Comput. Methods Appl. Mech. Eng.
A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes
J. Comput. Phys.
High-order unstructured Lagrangian one-step WENO finite volume schemes for non-conservative hyperbolic systems: applications to compressible multi-phase flows
Comput. Fluids
High order space–time adaptive ADER-WENO finite volume schemes for non-conservative hyperbolic systems
Comput. Methods Appl. Mech. Eng.
Efficient implementation of high order unstructured WENO schemes for cavitating flows
Comput. Fluids
Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems
J. Comput. Phys.
Quadrature-free non-oscillatory finite volume schemes on unstructured meshes for nonlinear hyperbolic systems
J. Comput. Phys.
Very high order PNPM schemes on unstructured meshes for the resistive relativistic MHD equations
J. Comput. Phys.
ADER-WENO finite volume schemes with space–time adaptive mesh refinement
J. Comput. Phys.
High accuracy numerical methods for thermally perfect gas flows with chemistry
J. Comput. Phys.
Numerical analysis of flow-induced nonlinear vibrations of an airfoil with three degrees of freedom
Comput. Fluids
Adaptive local refinement with octree load-balancing for the parallel solution of three-dimensional conservation laws
J. Parallel Distrib. Comput.
A comparative study of multimaterial Lagrangian and Eulerian methods with pressure relaxation
Comput. Fluids
Cell-centered discontinuous Galerkin discretization for two-dimensional Lagrangian hydrodynamics
Comput. Fluids
A discontinuous Galerkin discretization for solving the two-dimensional gas dynamics equations written under total Lagrangian formulation on general unstructured grids
J. Comput. Phys.
Explicit local time-stepping methods for Maxwell's equations
J. Comput. Appl. Math.
High-order explicit local time-stepping methods for damped wave equations
J. Comput. Appl. Math.
Efficient implementation of weighted ENO schemes
J. Comput. Phys.
ADER schemes on adaptive triangular meshes for scalar conservation laws
J. Comput. Phys.
An efficient local time-stepping scheme for solution of nonlinear conservation laws
J. Comput. Phys.
The cell-centered discontinuous Galerkin method for Lagrangian compressible Euler equations in two dimensions
Comput. Fluids
High order conservative Lagrangian schemes with Lax–Wendroff type time discretization for the compressible Euler equations
J. Comput. Phys.
A second-order compatible staggered Lagrangian hydrodynamics scheme using a cell-centered multidimensional approximate Riemann solver
Proc. Comput. Sci.
A high-order cell-centered Lagrangian scheme for compressible fluid flows in two-dimensional cylindrical geometry
J. Comput. Phys.
A high-order cell-centered Lagrangian scheme for two-dimensional compressible fluid flows on unstructured meshes
J. Comput. Phys.
A high-order one-step sub-cell force-based discretization for cell-centered Lagrangian hydrodynamics on polygonal grids
Comput. Fluids
Cited by (31)
A mimetic numerical scheme for multi-fluid flows with thermodynamic and geometric compatibility on an arbitrarily moving grid
2020, International Journal of Multiphase FlowAn arbitrary Lagrangian-Eulerian RKDG method for multi-material flows on adaptive unstructured meshes
2020, Computers and FluidsHigh order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes
2020, Journal of Computational PhysicsLocal time-stepping for adaptive multiresolution using natural extension of Runge–Kutta methods
2019, Journal of Computational PhysicsCitation Excerpt :This approach yields up to fourth order for compressible Euler equations in 2D. Similar and related works in this context were carried out in [15] and [2]. A detailed discussion on the stability of those MRLT schemes is presented in [14].
Reprint of: Direct Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming unstructured meshes
2018, Computers and FluidsCitation Excerpt :In the straight slip-line case no distinction between master and slave nodes is required and the mesh movement is only based on the corner-extrapolated values of the fluid velocity. At this point we also would like to refer to some recent works on high order Eulerian and ALE schemes on moving meshes with time-accurate local time stepping (LTS) presented in [55–57], where each element is allowed to run at its own optimal local time step according to a local CFL stability condition. The resulting algorithms use a conforming grid in space, but naturally produce a nonconforming mesh in time.