Linear dispersion-diffusion analysis and its application to under-resolved turbulence simulations using discontinuous Galerkin spectral/hp methods
File(s)accepted manuscript.pdf (2.28 MB)
Accepted version
Author(s)
Moura, RC
Sherwin, SJ
Peiro, J
Type
Journal Article
Abstract
We investigate the potential of linear dispersion–diffusion analysis in providing direct guidelines for turbulence simulations through the under-resolved DNS (sometimes called implicit LES) approach via spectral/hp methods. The discontinuous Galerkin (DG) formulation is assessed in particular as a representative of these methods. We revisit the eigensolutions technique as applied to linear advection and suggest a new perspective to the role of multiple numerical modes, peculiar to spectral/hp methods. From this new perspective, “secondary” eigenmodes are seen to replicate the propagation behaviour of a “primary” mode, so that DG's propagation characteristics can be obtained directly from the dispersion–diffusion curves of the primary mode. Numerical dissipation is then appraised from these primary eigencurves and its effect over poorly-resolved scales is quantified. Within this scenario, a simple criterion is proposed to estimate DG's effective resolution in terms of the largest wavenumber it can accurately resolve in a given hp approximation space, also allowing us to present points per wavelength estimates typically used in spectral and finite difference methods. Although strictly valid for linear advection, the devised criterion is tested against (1D) Burgers turbulence and found to predict with good accuracy the beginning of the dissipation range on the energy spectra of under-resolved simulations. The analysis of these test cases through the proposed methodology clarifies why and how the DG formulation can be used for under-resolved turbulence simulations without explicit subgrid-scale modelling. In particular, when dealing with communication limited hardware which forces one to consider the performance for a fixed number of degrees of freedom, the use of higher polynomial orders along with moderately coarser meshes is shown to be the best way to translate available degrees of freedom into resolution power.
Date Issued
2015-10-01
Date Acceptance
2015-06-26
Citation
Journal of Computational Physics, 2015, 298 (1), pp.695-710
ISSN
0021-9991
Publisher
Elsevier
Start Page
695
End Page
710
Journal / Book Title
Journal of Computational Physics
Volume
298
Issue
1
Copyright Statement
© 2015, Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International http://creativecommons.org/licenses/by-nc-nd/4.0/
Sponsor
Royal Academy Of Engineering
McLaren Racing Limited
Engineering & Physical Science Research Council (EPSRC)
Identifier
https://www.sciencedirect.com/science/article/pii/S0021999115004179
Grant Number
AEDZ_P40009
AEDZ_P42726
EP/L000407/1
Subjects
Science & Technology
Technology
Physical Sciences
Computer Science, Interdisciplinary Applications
Physics, Mathematical
Computer Science
Physics
Under-resolved DNS
Implicit LES
Dispersion-diffusion analysis
Discontinuous Galerkin formulation
Spectral/hp methods
LARGE-EDDY SIMULATION
DRIVEN BURGERS-EQUATION
ELEMENT METHODS
VANISHING VISCOSITY
NONUNIFORM GRIDS
FLOWS
SCHEMES
01 Mathematical Sciences
02 Physical Sciences
09 Engineering
Applied Mathematics
Publication Status
Published
Date Publish Online
2015-06-30