A matrix-free high-order discontinuous Galerkin compressible Navier-Stokes solver: A performance comparison of compressible and incompressible formulations for turbulent incompressible flows

Both compressible and incompressible Navier-Stokes solvers can be used and are used to solve incompressible turbulent flow problems. In the compressible case, the Mach number is then considered as a solver parameter that is set to a small value, $\mathrm{M}\approx 0.1$, in order to mimic incompressible flows. This strategy is widely used for high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations. The present work raises the question regarding the computational efficiency of compressible DG solvers as compared to a genuinely incompressible formulation. Our contributions to the state-of-the-art are twofold: Firstly, we present a high-performance discontinuous Galerkin solver for the compressible Navier-Stokes equations based on a highly efficient matrix-free implementation that targets modern cache-based multicore architectures. The performance results presented in this work focus on the node-level performance and our results suggest that there is great potential for further performance improvements for current state-of-the-art discontinuous Galerkin implementations of the compressible Navier-Stokes equations. Secondly, this compressible Navier-Stokes solver is put into perspective by comparing it to an incompressible DG solver that uses the same matrix-free implementation. We discuss algorithmic differences between both solution strategies and present an in-depth numerical investigation of the performance. The considered benchmark test cases are the three-dimensional Taylor-Green vortex problem as a representative of transitional flows and the turbulent channel flow problem as a representative of wall-bounded turbulent flows

Linear stability, transient energy growth and the role of viscosity stratification in compressible plane Couette flow

Linear stability and the non-modal transient energy growth in compressible plane Couette flow are investigated for two prototype mean flows: (a) the {\it uniform shear} flow with constant viscosity, and (b) the {\it non-uniform shear} flow with {\it stratified} viscosity. Both mean flows are linearly unstable for a range of supersonic Mach numbers ($M$). For a given $M$, the critical Reynolds number ($Re$) is significantly smaller for the uniform shear flow than its non-uniform shear counterpart. An analysis of perturbation energy reveals that the instability is primarily caused by an excess transfer of energy from mean-flow to perturbations. It is shown that the energy-transfer from mean-flow occurs close to the moving top-wall for ``mode I'' instability, whereas it occurs in the bulk of the flow domain for ``mode II''. For the non-modal analysis, it is shown that the maximum amplification of perturbation energy, $G_{\max}$, is significantly larger for the uniform shear case compared to its non-uniform counterpart. For $\alpha=0$, the linear stability operator can be partitioned into ${\cal L}\sim \bar{\cal L} + Re^2{\cal L}_p$, and the $Re$-dependent operator ${\cal L}_p$ is shown to have a negligibly small contribution to perturbation energy which is responsible for the validity of the well-known quadratic-scaling law in uniform shear flow: $G(t/{\it Re}) \sim {\it Re}^2$. A reduced inviscid model has been shown to capture all salient features of transient energy growth of full viscous problem. For both modal and non-modal instability, it is shown that the {\it viscosity-stratification} of the underlying mean flow would lead to a delayed transition in compressible Couette flow

Implicit High-Order Flux Reconstruction Solver for High-Speed Compressible Flows

The present paper addresses the development and implementation of the first high-order Flux Reconstruction (FR) solver for high-speed flows within the open-source COOLFluiD (Computational Object-Oriented Libraries for Fluid Dynamics) platform. The resulting solver is fully implicit and able to simulate compressible flow problems governed by either the Euler or the Navier-Stokes equations in two and three dimensions. Furthermore, it can run in parallel on multiple CPU-cores and is designed to handle unstructured grids consisting of both straight and curved edged quadrilateral or hexahedral elements. While most of the implementation relies on state-of-the-art FR algorithms, an improved and more case-independent shock capturing scheme has been developed in order to tackle the first viscous hypersonic simulations using the FR method. Extensive verification of the FR solver has been performed through the use of reproducible benchmark test cases with flow speeds ranging from subsonic to hypersonic, up to Mach 17.6. The obtained results have been favorably compared to those available in literature. Furthermore, so-called super-accuracy is retrieved for certain cases when solving the Euler equations. The strengths of the FR solver in terms of computational accuracy per degree of freedom are also illustrated. Finally, the influence of the characterizing parameters of the FR method as well as the the influence of the novel shock capturing scheme on the accuracy of the developed solver is discussed

Institute for Computational Mechanics in Propulsion (ICOMP)

The Institute for Computational Mechanics in Propulsion (ICOMP) is a combined activity of Case Western Reserve University, Ohio Aerospace Institute (OAI) and NASA Lewis. The purpose of ICOMP is to develop techniques to improve problem solving capabilities in all aspects of computational mechanics related to propulsion. The activities at ICOMP during 1991 are described

Computational methods for internal flows with emphasis on turbomachinery

Current computational methods for analyzing flows in turbomachinery and other related internal propulsion components are presented. The methods are divided into two classes. The inviscid methods deal specifically with turbomachinery applications. Viscous methods, deal with generalized duct flows as well as flows in turbomachinery passages. Inviscid methods are categorized into the potential, stream function, and Euler aproaches. Viscous methods are treated in terms of parabolic, partially parabolic, and elliptic procedures. Various grids used in association with these procedures are also discussed

Spatial three-dimensional secondary instability compressible boundary-layer flows

Three-dimensional linear secondary instability theory is extended for compressible and high Mach number boundary layer flows. The small but finite amplitude compressible Tollmien-Schlichting wave effect on the growth of 3-D perturbations is investigated. The focus is on principal parametric resonance responsible for the strong growth of subharmonic in low disturbance environment. The effect of increasing Mach number on the onset, growth, the shape of eigenfunctions of the subharmonic is assessed, and the resulting vortical structure is examined