Universität zu Köln

Flad, David ORCID: 0000-0003-1155-9355 and Gassner, Gregor (2017). On the use of kinetic energy preserving DG-schemes for large eddy simulation. J. Comput. Phys., 350. S. 782 - 796. SAN DIEGO: ACADEMIC PRESS INC ELSEVIER SCIENCE. ISSN 1090-2716

Full text not available from this repository.

Abstract

Recently, element based high order methods such as Discontinuous Galerkin (DG) methods and the closely related flux reconstruction (FR) schemes have become popular for compressible large eddy simulation (LES). Element based high order methods with Riemann solver based interface numerical flux functions offer an interesting dispersion dissipation behavior for multi-scale problems: dispersion errors are very low for a broad range of scales, while dissipation errors are very low for well resolved scales and are very high for scales close to the Nyquist cutoff. In some sense, the inherent numerical dissipation caused by the interface Riemann solver acts as a filter of high frequency solution components. This observation motivates the trend that element based high order methods with Riemann solvers are used without an explicit LES model added. Only the high frequency type inherent dissipation caused by the Riemann solver at the element interfaces is used to account for the missing sub-grid scale dissipation. Due to under resolution of vortical dominated structures typical for LES type setups, element based high order methods suffer from stability issues caused by aliasing errors of the non-linear flux terms. A very common strategy to fight these aliasing issues (and instabilities) is socalled polynomial de-aliasing, where interpolation is exchanged with projection based on an increased number of quadrature points. In this paper, we start with this common no model or implicit LES (iLES) DG approach with polynomial de-aliasing and Riemann solver dissipation and review its capabilities and limitations. We find that the strategy gives excellent results, but only when the resolution is such, that about 40% of the dissipation is resolved. For more realistic, coarser resolutions used in classical LES e.g. of industrial applications, the iLES DG strategy becomes quite inaccurate. We show that there is no obvious fix to this strategy, as adding for instance a sub-grid-scale models on top doesn't change much or in worst case decreases the fidelity even more. Finally, the core of this work is a novel LES strategy based on split form DG methods that are kinetic energy preserving. The scheme offers excellent stability with full control over the amount and shape of the added artificial dissipation. This premise is the main idea of the work and we will assess the LES capabilities of the novel split form DG approach when applied to shock-free, moderate Mach number turbulence. We will demonstrate that the novel DG LES strategy offers similar accuracy as the iLES methodology for well resolved cases, but strongly increases fidelity in case of more realistic coarse resolutions. (C) 2017 Elsevier Inc. All rights reserved.

Item Type:	Journal Article
Creators:	Creators Email ORCID ORCID Put Code Flad, David UNSPECIFIED orcid.org/0000-0003-1155-9355 UNSPECIFIED Gassner, Gregor UNSPECIFIED UNSPECIFIED UNSPECIFIED
URN:	urn:nbn:de:hbz:38-210706
DOI:	10.1016/j.jcp.2017.09.004
Journal or Publication Title:	J. Comput. Phys.
Volume:	350
Page Range:	S. 782 - 796
Date:	2017
Publisher:	ACADEMIC PRESS INC ELSEVIER SCIENCE
Place of Publication:	SAN DIEGO
ISSN:	1090-2716
Language:	English
Faculty:	Unspecified
Divisions:	Unspecified
Subjects:	no entry
Uncontrolled Keywords:	Keywords Language NAVIER-STOKES EQUATIONS; DIRECT NUMERICAL-SIMULATION; FINITE-DIFFERENCE SCHEMES; SPECTRAL ELEMENT METHODS; TURBULENT FLOWS; HIGH-RESOLUTION; IMPLICIT LES; ERRORS; DISCRETIZATION; SUITABILITY Multiple languages Computer Science, Interdisciplinary Applications; Physics, Mathematical Multiple languages
Refereed:	Yes
URI:	http://kups.ub.uni-koeln.de/id/eprint/21070