Abstract:For computational fluid dynamics (CFD), the generalized Riemann problem (GRP) solver and the second-order gas-kinetic scheme (GKS) provide a time-accurate flux function starting from a discontinuous piecewise linear flow distributions around a cell interface. With the adoption of time derivative of the flux function, a two-stage Lax–Wendroff-type (L-W for short) time stepping method has been recently proposed in the design of a fourth-order time accurate method for inviscid flow [21]. In this paper, based on the same time-stepping method and the second-order GKS flux function [42], a fourth-order gas-kinetic scheme is constructed for the Euler and Navier–Stokes (NS) equations. In comparison with the formal one-stage time-stepping third-order gas-kinetic solver [24], the current fourth-order method not only reduces the complexity of the flux function, but also improves the accuracy of the scheme. In terms of the computational cost, a two-dimensional third-order GKS flux function takes about six times of the computational time of a second-order GKS flux function. However, a fifth-order WENO reconstruction may take more than ten times of the computational cost of a second-order GKS flux function. Therefore, it is fully legitimate to develop a two-stage fourth order time accurate method (two reconstruction) instead of standard four stage fourth-order Runge–Kutta method (four reconstruction). Most importantly, the robustness of the fourth-order GKS is as good as the second-order one. In the current computational fluid dynamics (CFD) research, it is still a difficult problem to extend the higher-order Euler solver to the NS one due to the change of governing equations from hyperbolic to parabolic type and the initial interface discontinuity. This problem remains distinctively for the hypersonic viscous and heat conducting flow. The GKS is based on the kinetic equation with the hyperbolic transport and the relaxation source term. The time-dependent GKS flux function provides a dynamic process of evolution from the kinetic scale particle free transport to the hydrodynamic scale wave propagation, which provides the physics for the non-equilibrium numerical shock structure construction to the near equilibrium NS solution. As a result, with the implementation of the fifth-order WENO initial reconstruction, in the smooth region the current two-stage GKS provides an accuracy of O((Δx)5,(Δt)4) for the Euler equations, and O((Δx)5,τ2Δt) for the NS equations, where τ is the time between particle collisions. Many numerical tests, including difficult ones for the Navier–Stokes solvers, have been used to validate the current method. Perfect numerical solutions can be obtained from the high Reynolds number boundary layer to the hypersonic viscous heat conducting flow. Following the two-stage time-stepping framework, the third-order GKS flux function can be used as well to construct a fifth-order method with the usage of both first-order and second-order time derivatives of the flux function. The use of time-accurate flux function may have great advantages on the development of higher-order CFD methods.

High-order Gas-Kinetic Scheme with Three-Dimensional WENO Reconstruction for the Euler and Navier-Stokes Solutions

Two-stage fourth-order gas-kinetic scheme for three-dimensional Euler and Navier-Stokes solutions

Gas-Kinetic Scheme Coupled with Turbulent Kinetic Energy Equation for Computing Hypersonic Turbulent and Transitional Flows

A Third-Order Gas-Kinetic Scheme for Three-Dimensional Inviscid and Viscous Flow Computations

Compact and Three-dimensional Extension of Two-Stage Fourth-Order Gas-Kinetic Scheme

High-order gas-kinetic scheme on three-dimensional unstructured meshes for compressible flows

A High-Order Multidimensional Gas-Kinetic Scheme for Hydrodynamic Equations

Three-dimensional Third-Order Gas-Kinetic Scheme on Hybrid Unstructured Meshes for Euler and Navier-Stokes Equations

An efficient high-order gas-kinetic scheme with hybrid WENO-AO method for the Euler and Navier-Stokes solutions

A HWENO Reconstruction Based High-Order Compact Gas-Kinetic Scheme on Unstructured Mesh.

Implicit high-order gas-kinetic schemes for compressible flows on three-dimensional unstructured meshes I: Steady flows

High-Order Unified Gas-kinetic Scheme

Compact High-Order Gas-Kinetic Scheme for Three-Dimensional Flow Simulations

High-Order Gas-Kinetic Scheme in General Curvilinear Coordinate for iLES of Compressible Wall-Bounded Turbulent Flows

An efficient and accurate two-stage fourth-order gas-kinetic scheme for the Euler and Navier–Stokes equations

Performance Enhancement for High-order Gas-kinetic Scheme Based on WENO-adaptive-order Reconstruction

Two-step Multi-Resolution Reconstruction-Based Compact Gas-Kinetic Scheme on Tetrahedral Mesh

Compact high-order gas-kinetic scheme for direct numerical simulation of compressible turbulent flows

An Arbitrary-Lagrangian-Eulerian High-Order Gas-Kinetic Scheme for Three-Dimensional Computations.

A three-dimensional compact high-order gas-kinetic scheme on structured mesh