A comparison of interpolation techniques for non-conformal high-order discontinuous Galerkin methods

Abstract

The capability to incorporate moving geometric features within models for complex simulations is a common requirement in many fields. Fluid mechanics within aeronautical applications, for example, routinely feature rotating (e.g. turbines, wheels and fan blades) or sliding components (e.g. in compressor or turbine cascade simulations). With an increasing trend towards the high-fidelity modelling of these cases, in particular combined with the use of high-order discontinuous Galerkin methods, there is therefore a requirement to understand how different numerical treatments of the interfaces between the static mesh and the sliding/rotating part impact on overall solution quality. In this article, we compare two different approaches to handle this non-conformal interface. The first is the so-called mortar approach, where flux integrals along edges are split according to the positioning of the non-conformal grid. The second is a less-documented point-to-point interpolation method, where the interior and exterior quantities for flux evaluations are interpolated from elements lying on the opposing side of the interface. Although the mortar approach has significant advantages in terms of its numerical properties, in that it preserves the local conservation properties of DG methods, in the context of complex 3D meshes it poses notable implementation difficulties which the point-to-point method handles more readily. In this paper we examine the numerical properties of each method, focusing not only on observing convergence orders for smooth solutions, but also how each method performs in under-resolved simulations of linear and nonlinear hyperbolic problems, to inform the use of these methods in implicit large-eddy simulations.

arXiv comment37 pages, 15 figures, 5 tables, submitted to Computer Methods in Applied Mechanics and Engineering, revision 1

References in corpus (8)

• Nektar++: enhancing the capability and application of high-fidelity spectral/$hp$ element methods
• On the use of kinetic energy preserving DG-schemes for large eddy simulation
• A Simple, Efficient, High-order Accurate Sliding-Mesh Interface Approach to the Spectral Difference Method on Coupled Rotating and Stationary Domains
• Geophysical-astrophysical spectral-element adaptive refinement (GASpAR): Object-oriented h-adaptive code for geophysical fluid dynamics simulation
• A matrix-free high-order discontinuous Galerkin compressible Navier-Stokes solver: A performance comparison of compressible and incompressible formulations for turbulent incompressible flows
• Nektar++: Design and implementation of an implicit, spectral/$hp$ element, compressible flow solver using a Jacobian-free Newton Krylov approach
• An Improved Roe Scheme for All Mach-Number Flows Simultaneously Curing Known Problems
• Fast Barycentric-Based Evaluation Over Spectral/hp Elements

Cited by in corpus (5)

• A unifying algebraic framework for discontinuous Galerkin and flux reconstruction methods based on the summation-by-parts property
• A Time-Accurate Inflow Coupling for Zonal LES
• Stable, entropy-pressure compatible subsonic Riemann boundary condition for embedded DG compressible flow simulations
• Fast Barycentric-Based Evaluation Over Spectral/hp Elements
• An Efficient Sliding Mesh Interface Method for High-Order Discontinuous Galerkin Schemes