A Unified Framework for the Solution of Hyperbolic PDE Systems Using High Order Direct Arbitrary-Lagrangian–Eulerian Schemes on Moving Unstructured Meshes with Topology Change (original) (raw)
Related papers
Computers & Fluids
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 flows in 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.
FORCE schemes on moving unstructured meshes for hyperbolic systems
Computers & Mathematics with Applications, 2019
The aim of this paper is to propose a new simple and robust numerical flux of the centered type in the context of Arbitrary-Lagrangian-Eulerian (ALE) finite volume schemes. The work relies on the FORCE flux of Toro and Billet and is concerned with the solution of general hyperbolic systems of nonlinear equations involving both conservative and nonconservative terms as well as sources which might become stiff. The proposed approach is formulated in a general way using a path-conservative method and the Roe-type system matrix is computed numerically in order to provide a numerical flux function that can be applied to any given hyperbolic system. Furthermore, one great advantage of the FORCE flux is that no information about the eigenstructure of the system is needed, not even eigenvalues, but only information regarding the geometry of the control volumes are required, which are automatically available in the moving mesh framework. Our method is of the finite volume type, high order accurate in space, thanks to a WENO reconstruction operator, and even in time, due to a fully-discrete ADER one-step discretization. The algorithm applies to moving multidimensional unstructured meshes composed by triangles and tetrahedra. Both accuracy and robustness of the scheme are assessed on a series of test problems for the Euler equations of compressible gas dynamics, for the magnetohydrodynamics equations as well as for the Baer-Nunziato model of compressible multi-phase flows.
High Order ADER Schemes for Continuum Mechanics
Frontiers in Physics
In this paper we first review the development of high order ADER finite volume and ADER discontinuous Galerkin schemes on fixed and moving meshes, since their introduction in 1999 by Toro et al. We show the modern variant of ADER based on a space-time predictor-corrector formulation in the context of ADER discontinuous Galerkin schemes with a posteriori subcell finite volume limiter on fixed and moving grids, as well as on space-time adaptive Cartesian AMR meshes. We then present and discuss the unified symmetric hyperbolic and thermodynamically compatible (SHTC) formulation of continuum mechanics developed by Godunov, Peshkov and Romenski (GPR model), which allows to describe fluid and solid mechanics in one single and unified first order hyperbolic system. In order to deal with free surface and moving boundary problems, a simple diffuse interface approach is employed, which is compatible with Eulerian schemes on fixed grids as well as direct Arbitrary-Lagrangian-Eulerian methods on moving meshes. We show some examples of moving boundary problems in fluid and solid mechanics.
ADER-WENO finite volume schemes with space–time adaptive mesh refinement
Journal of Computational Physics, 2013
We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space-time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.
Computers & Fluids, 2009
We develop a new family of well-balanced path-conservative quadrature-free one-step ADER finite volume and discontinuous Galerkin finite element schemes on unstructured meshes for the solution of hyperbolic partial differential equations with non-conservative products and stiff source terms. The fully discrete formulation is derived using the recently developed framework of explicit one-step P N P M schemes of arbitrary high order of accuracy in space and time for conservative hyperbolic systems [Dumbser M, Balsara D, Toro EF, Munz CD. A unified framework for the construction of one-step finitevolume and discontinuous Galerkin schemes. J Comput Phys 2008;227:8209-53]. The two key ingredients of our high order approach are: first, the high order accurate P N P M reconstruction operator on unstructured meshes, using the WENO strategy presented in [Dumbser M, Käser M, Titarev VA Toro EF. Quadrature-free non-oscillatory finite volume schemes on unstructured meshes for nonlinear hyperbolic systems. J Comput Phys 2007;226:204-43] to ensure monotonicity at discontinuities, and second, a local space-time Galerkin scheme to predict the evolution of the reconstructed polynomial data inside each element during one time step to obtain a high order accurate one-step time discretization. This approach is also able to deal with stiff source terms as shown in [Dumbser M, Enaux C, Toro EF. Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws. J Comput Phys 2008;227:3971-4001]. These two key ingredients are combined with the recently developed path-conservative methods of Parés [Parés C. Numerical methods for nonconservative hyperbolic systems: a theoretical framework. SIAM J Numer Anal 2006;44:300-21] and Castro et al. [Castro MJ, Gallardo JM, Parés C. High-order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems. Math Comput 2006;75:1103-34] to treat the non-conservative products properly. We show applications of our method to the twolayer shallow water equations as well as applications to the recently published depth-averaged two-fluid flow model of Pitman and Le [Pitman EB, Le L. A two-fluid model for avalanche and debris flows. Philos Trans Roy Soc A 2005;363:1573-601].
Computer Methods in Applied Mechanics and Engineering, 2014
We present a class of high order finite volume schemes for the solution of nonconservative hyperbolic systems that combines the one-step ADER-WENO finite volume approach with space-time adaptive mesh refinement (AMR). The resulting algorithm, which is particularly well suited for the treatment of material interfaces in compressible multi-phase flows, is based on: (i) high order of accuracy in space obtained through WENO reconstruction, (ii) a high order one-step time discretization via a local space-time discontinuous Galerkin predictor method, and (iii) the use of a path conservative scheme for handling the non-conservative terms of the equations. The AMR property with time accurate local time stepping, which has been treated according to a cell-by-cell strategy, strongly relies on the high order one-step time discretization, which naturally allows a high order accurate and consistent computation of the jump terms at interfaces between elements using different time steps. The new scheme has been succesfully validated on some test problems for the Baer-Nunziato model of compressible multiphase flows.
Journal of Computational Physics, 2018
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. Highlights • A velocity filter is presented to reduce spurious mesh motion. • A new multidirectional approximate Riemann problem is presented. • More robust mesh motion with higher-order and polygonal meshes.
16th AIAA Computational Fluid Dynamics Conference, 2003
We consider the solution of inviscid as well as viscous unsteady flow problems with moving boundaries by the arbitrary Lagrangian Eulerian (ALE) method. We present two computational approaches for achieving formal second-order time-accuracy on moving grids. The first approach is based on flux time-averaging, and the second one on mesh configuration time-averaging. In both cases, we show that formally second-order time-accurate ALE schemes can be designed. We illustrate our theoretical findings and highlight their impact on practice with the solution of inviscid as well as viscous, unsteady, nonlinear flow problems associated with the AGARD Wing 445.6 and a complete F-16 configuration.