David Flad

Simulation of underresolved turbulent flows by adaptive filtering using the high order discontinuous Galerkin spectral element method

Scale-resolving simulations of turbulent flows in complex domains demand accurate and efficient numerical schemes, as well as geometrical flexibility. For underresolved situations, the avoidance of aliasing errors is a strong demand for stability. For continuous and discontinuous Galerkin schemes, an effective way to prevent aliasing errors is to increase the quadrature precision of the projection operator to account for the non-linearity of the operands (polynomial dealiasing, overintegration). But this increases the computational costs extensively. In this work, we present a novel spatially and temporally adaptive dealiasing strategy by projection filtering. We show this to be more efficient for underresolved turbulence than the classical overintegration strategy. For this novel approach, we discuss the implementation strategy and the indicator details, show its accuracy and efficiency for a decaying homogeneous isotropic turbulence and the transitional Taylor-Green vortex and compare it to the original overintegration approach and a state of the art variational multi-scale eddy viscosity formulation.

A Discontinuous Galerkin Spectral Element Method for the direct numerical simulation of aeroacoustics

The discontinuous Galerkin Spectral Element Method (DG-SEM) is highly attractive for both DNS and LES of turbulent flows due to its low dispersion and dissipation errors, but also because of its good parallel scaling property. We show that especially for underresolved simulations the method has highly beneficial properties for LES, as well as for the direct treatment of the acoustic propagation. We also discuss different approaches for non-reflecting boundary conditions for DG-SEM and show the behavior of the methods for airfoil flows. Our main intend is to directly simulate trailing edge noise for airfoil flows at medium Reynolds numbers.

High-Order Discontinuous Galerkin Schemes for Large-Eddy Simulations of Moderate Reynolds Number Flows

In this article, we describe the capabilities of high order discontinuous Galerkin methods at the Institute for Aerodynamics and Gasdynamics for the Large-Eddy Simulation of wall-bounded flows at moderate Reynolds numbers. In these scenarios, the prediction of laminar regions, flow transition and developed turbulence poses a great challenge to the numerical scheme, as overprediction of numerical dissipation can significantly influence the accuracy of the integral quantities. While this increases the burden on the numerical scheme and the LES subgrid model, the moderate Reynolds numbers prevent the occurrence of thin wall boundary layers and allows the resolution of the boundary layer without the need for wall modelling strategies. We take full advantage of the low numerical errors and associated superior scale resolving capabilities of high order spectral approximations by using high order ansatz functions up to 12th order, which allows us to resolve the significant features of these flows at a very low number of degrees of freedom. Without the need for any additional filtering, explicit or implicit modelling or artificial dissipation, the high order scheme capture the turbulent flow at the considered Reynolds number range very well.

Application and Development of the High Order Discontinuous Galerkin Spectral Element Method for Compressible Multiscale Flows

This paper summarizes our progress in the application of a high-order discontinuous Galerkin (DG) method for scale resolving fluid dynamics simulations on the Cray XC40 Hazel Hen cluster at HLRS. We present the large eddy simulation (LES) of flow around a wall mounted cylinder, a LES of flow around an airfoil at realistic Reynolds number using a recently introduced kinetic energy preserving flux formulation and a simulation of transitional flow in a low pressure turbine. Furthermore, it provides an overview over the parallel efficiency reached by our code when using up to 49, 152 CPUs and the latest developments of our DG framework.

High Fidelity Scale-Resolving Computational Fluid Dynamics Using the High Order Discontinuous Galerkin Spectral Element Method

In this report we give an overview of our high-order simulations of turbulent flows carried out on the HLRS systems. The simulation framework is built around a highly scalable solver based on the discontinuous Galerkin spectral element method (DGSEM). It has been designed to support large scale simulations on massively parallel architectures and at the same time enabling the use of complex geometries with unstructured, nonconforming meshes. We are thus capable of fully exploiting the performance of HLRS Cray XE6 (Hermit) and XC40 (Hornet) systems not just for academic benchmark problems but also industrial applications. We exemplify the capabilities of our framework at three recent simulations, where we have performed direct numerical and large eddy simulations of turbulent compressible flows. The test cases include a high-speed turbulent boundary layer flow utilizing close to 94,000 physical cores, a DNS of a NACA 0012 airfoil at Re = 100, 000 and direct aeroacoustic simulations of a close-to-production car mirror at Re c = 100, 000.