High Order Accurate Direct Arbitrary-Lagrangian-Eulerian ADER-MOOD Finite Volume Schemes for Non-Conservative Hyperbolic Systems with Stiff Source Terms. (5th December 2016)

• Record Type: Journal Article
• Authors: Boscheri, Walter; Loubère, Raphaël
• Is Part Of: Communications in computational physics. Volume 21:Number 1(2017:Jan.)
• Page Start: 271
• Page End: 312
• Publication Date: 2016-12-05
• Subjects: 65M08, -- 65M60; Direct Arbitrary-Lagrangian-Eulerian, -- a posteriori MOOD stabilization, -- Baer-Nunziato model, -- stiff source terms, -- non-conservative products, -- unstructured mesh, -- ADER, -- high order of accuracy in space and time, -- high performance computing (HPC), -- hyperbolic conservation laws; Mathematical physics -- Data processing -- Periodicals; Physics -- Data processing -- Periodicals; 530.150285
• DOI: 10.4208/cicp.OA-2015-0024
• Languages: English
• ISSNs: 1815-2406

Abstract

Abstract: In this paper we present a 2D/3D high order accurate finite volume scheme in the context of direct Arbitrary-Lagrangian-Eulerian algorithms for general hyperbolic systems of partial differential equations with non-conservative products and stiff source terms. This scheme is constructed with a single stencil polynomial reconstruction operator, a one-step space-time ADER integration which is suitably designed for dealing even with stiff sources, a nodal solver with relaxation to determine the mesh motion, a path-conservative integration technique for the treatment of non-conservative products and an a posteriori stabilization procedure derived from the so-called Multidimensional Optimal Order Detection (MOOD) paradigm. In this work we consider the seven equation Baer-Nunziato model of compressible multi-phase flows as a representative model involving non-conservative products as well as relaxation source terms which are allowed to become stiff. The new scheme is validated against a set of test cases on 2D/3D unstructured moving meshes on parallel machines and the high order of accuracy achieved by the method is demonstrated by performing a numerical convergence study. Classical Riemann problems and explosion problems with exact solutions are simulated in 2D and 3D. The overall numerical code is also profiled to provide an estimate of the computational cost required by each component of the whole algorithm.