Search

In recent years, high-order discontinuous Galerkin (DG) methods have emerged as an attractive approach for numerical simulations of compressible flows. This paper presents an overview of the recent development of DG methods for compressible flows with particular focus on hypersononic flows. First, we survey state-of-the-art DG methods for computational fluid dynamics. Next, we discuss both matrix-based and matrix-free iterative methods for the solution of discrete systems stemming from the spatial DG discretizations of the compressible Navier-Stokes equations. We then describe various shock capturing methods to deal with strong shock waves in hypersonic flows. We discuss adaptivity techniques to refine high-order meshes, and synthetic boundary conditions to simulate free-stream disturbances in hypersonic boundary layers. We present a few examples to demonstrate the ability of high-order DG methods to provide accurate solutions of hypersonic laminar flows. Furthermore, we present direct numerical simulations of hypersonic transitional flow past a flared cone at Reynolds number $10.8 \times 10^6$, and hypersonic transitional shock wave boundary layer interaction flow over a flat plate at Reynolds number $3.97 \times 10^6$. These simulations run entirely on hundreds of graphics processing units (GPUs) and demonstrate the ability of DG methods to directly resolve hypersonic transitional flows, even at high Reynolds numbers, without relying on transition or turbulence models. We end the paper by offering our perspectives on error estimation, turbulence modeling, and real gas effects in hypersonic flows., Comment: 34 pages, 25 figures, and 1 table

We present an optimal transport approach for mesh adaptivity and shock capturing of compressible flows. Shock capturing is based on a viscosity regularization of the governing equations by introducing an artificial viscosity field as solution of the Helmholtz equation. Mesh adaptation is based on the optimal transport theory by formulating a mesh mapping as solution of Monge-Ampere equation. The marriage of optimal transport and viscosity regularization for compressible flows leads to a coupled system of the compressible Euler/Navier-Stokes equations, the Helmholtz equation, and the Monge-Ampere equation. We propose an iterative procedure to solve the coupled system in a sequential fashion using homotopy continuation to minimize the amount of artificial viscosity while enforcing positivity-preserving and smoothness constraints on the numerical solution. We explore various mesh monitor functions for computing r-adaptive meshes in order to reduce the amount of artificial dissipation and improve the accuracy of the numerical solution. The hybridizable discontinuous Galerkin method is used for the spatial discretization of the governing equations to obtain high-order accurate solutions. Extensive numerical results are presented to demonstrate the optimal transport approach on transonic, supersonic, hypersonic flows in two dimensions. The approach is found to yield accurate, sharp yet smooth solutions within a few mesh adaptation iterations., Comment: 41 pages, 22 figures. arXiv admin note: text overlap with arXiv:2305.00461

The accurate prediction of aerodynamic drag on satellites orbiting in the upper atmosphere is critical to the operational success of modern space technologies, such as satellite-based communication or navigation systems, which have become increasingly popular in the last few years due to the deployment of constellations of satellites in low-Earth orbit. As a result, physics-based models of the ionosphere and thermosphere have emerged as a necessary tool for the prediction of atmospheric outputs under highly variable space weather conditions. This paper proposes a high-fidelity approach for physics-based space weather modeling based on the solution of the Navier-Stokes equations using a high-order discontinuous Galerkin method, combined with a matrix-free strategy suitable for high-performance computing on GPU architectures. The approach consists of a thermospheric model that describes a chemically frozen neutral atmosphere in non-hydrostatic equilibrium driven by the external excitation of the Sun. A novel set of variables is considered to treat the low densities present in the upper atmosphere and to accommodate the wide range of scales present in the problem. At the same time, and unlike most existing approaches, radial and angular directions are treated in a non-segregated approach. The study presents a set of numerical examples that demonstrate the accuracy of the approximation and validate the current approach against observational data along a satellite orbit, including estimates of established empirical and physics-based models of the ionosphere-thermosphere system. Finally, a 1D radial derivation of the physics-based model is presented and utilized for conducting a parametric study of the main thermal quantities under various solar conditions., Comment: 26 pages, 13 figures, 5 tables

We introduce an adaptive viscosity regularization approach for the numerical solution of systems of nonlinear conservation laws with shock waves. The approach seeks to solve a sequence of regularized problems consisting of the system of conservation laws and an additional Helmholtz equation for the artificial viscosity. We propose a homotopy continuation of the regularization parameters to minimize the amount of artificial viscosity subject to positivity-preserving and smoothness constraints on the numerical solution. The regularization methodology is combined with a mesh adaptation strategy that identifies the shock location and generates shock-aligned meshes, which allows to further reduce the amount of artificial dissipation and capture shocks with increased accuracy. We use the hybridizable discontinuous Galerkin method to numerically solve the regularized system of conservation laws and the continuous Galerkin method to solve the Helmholtz equation for the artificial viscosity. We show that the approach can produce approximate solutions that converge to the exact solution of the Burgers' equation. Finally, we demonstrate the performance of the method on inviscid transonic, supersonic, hypersonic flows in two dimensions. The approach is found to be accurate, robust and efficient, and yields very sharp yet smooth solutions in a few homotopy iterations., Comment: 42 pages, 22 figures, 4 tables

Purpose: This study aims to assess the robustness and accuracy of the face-centred finite volume (FCFV) method for the simulation of compressible laminar flows in different regimes, using numerical benchmarks. Design/methodology/approach: The work presents a detailed comparison with reference solutions published in the literature -- when available -- and numerical results computed using a commercial cell-centred finite volume software. Findings: The FCFV scheme provides first-order accurate approximations of the viscous stress tensor and the heat flux, insensitively to cell distortion or stretching. The strategy demonstrates its efficiency in inviscid and viscous flows, for a wide range of Mach numbers, also in the incompressible limit. In purely inviscid flows, non-oscillatory approximations are obtained in the presence of shock waves. In the incompressible limit, accurate solutions are computed without pressure correction algorithms. The method shows its superior performance for viscous high Mach number flows, achieving physically admissible solutions without carbuncle effect and predictions of quantities of interest with errors below 5%. Originality/value: The FCFV method accurately evaluates, for a wide range of compressible laminar flows, quantities of engineering interest, such as drag, lift and heat transfer coefficients, on unstructured meshes featuring distorted and highly stretched cells, with an aspect ratio up to ten thousand. The method is suitable to simulate industrial flows on complex geometries, relaxing the requirements on mesh quality introduced by existing finite volume solvers and alleviating the need for time-consuming manual procedures for mesh generation to be performed by specialised technicians., Comment: 39 pages, 18 figures, 12 tables

This paper presents an overview of the functionalities and applications of Exasim, an open-source code for generating high-order discontinuous Galerkin codes to numerically solve parametrized partial differential equations (PDEs). The software combines high-level and low-level languages to construct parametrized PDE models via Julia, Python or Matlab scripts and produce high-performance C++ codes for solving the PDE models on CPU and Nvidia GPU processors with distributed memory. Exasim provides matrix-free discontinuous Galerkin discretization schemes together with scalable reduced basis preconditioners and Newton-GMRES solvers, making it suitable for accurate and efficient approximation of wide-ranging classes of PDEs., Comment: 19 pages, 4 figures, 3 tables

This work presents the face-centred finite volume (FCFV) paradigm for the simulation of compressible flows. The FCFV method defines the unknowns at the face barycentre and uses a hybridisation procedure to eliminate all the degrees of freedom inside the cells. In addition, Riemann solvers are defined implicitly within the expressions of the numerical fluxes. The resulting methodology provides first-order accurate approximations of the conservative quantities, i.e. density, momentum and energy, as well as of the viscous stress tensor and of the heat flux, without the need of any gradient reconstruction procedure. Hence, the FCFV solver preserves the accuracy of the approximation in presence of distorted and highly stretched cells, providing a solver insensitive to mesh quality. In addition, FCFV is capable of constructing non-oscillatory approximations of sharp discontinuities without resorting to shock capturing or limiting techniques. For flows at low Mach number, the method is robust and is capable of computing accurate solutions in the incompressible limit without the need of introducing specific pressure correction strategies. A set of 2D and 3D benchmarks of external flows is presented to validate the methodology in different flow regimes, from inviscid to viscous laminar flows, from transonic to subsonic incompressible flows, demonstrating its potential to handle compressible flows in realistic scenarios., Comment: 38 pages, 22 figures

This work presents a review of high-order hybridisable discontinuous Galerkin (HDG) methods in the context of compressible flows. Moreover, an original unified framework for the derivation of Riemann solvers in hybridised formulations is proposed. This framework includes, for the first time in an HDG context, the HLL and HLLEM Riemann solvers as well as the traditional Lax-Friedrichs and Roe solvers. HLL-type Riemann solvers demonstrate their superiority with respect to Roe in supersonic cases due to their positivity preserving properties. In addition, HLLEM specifically outstands in the approximation of boundary layers because of its shear preservation, which confers it an increased accuracy with respect to HLL and Lax-Friedrichs. A comprehensive set of relevant numerical benchmarks of viscous and inviscid compressible flows is presented. The test cases are used to evaluate the competitiveness of the resulting high-order HDG scheme with the aforementioned Riemann solvers and equipped with a shock treatment technique based on artificial viscosity., Comment: 60 pages, 31 figures. arXiv admin note: substantial text overlap with arXiv:1912.00044

This work proposes a high-order hybridised discontinuous Galerkin (HDG) formulation of the Harten-Lax-Van Leer (HLL) Riemann solver for compressible flows. A unified framework is introduced to present Lax-Friedrichs, Roe and HLL Riemann solvers via appropriate definitions of the HDG numerical fluxes. The resulting high-order HDG method with HLL Riemann solver is evaluated through a set of numerical simulations of inviscid compressible flows in different regimes, from subsonic isentropic flows to transonic and supersonic problems with shocks. The accuracy of the proposed method is comparable with the one of Lax-Friedrichs and Roe numerical fluxes in subsonic and transonic flows. The superior performance of HLL is highlighted in supersonic cases, where the method provides extra robustness, being able to produce positivity preserving approximations without the need of any user-defined entropy fix., Comment: An extended and updated version of the paper covering both inviscid and viscous compressible flows is available at arXiv:2009.06396