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.

Fourth-order Compact Scheme Finite Volume Method for Large Eddy Simulations of Turbulent Flows

High Accurate Finite Volume Method on Non-Uniform Meshes for Large Eddy Simulation of Turbulent Flows

High Accurate Finite Volume Method for Large Eddy Simulation of Complex Turbulent Flows

Large Eddy Simulation of Turbulent Flows Passage of Reversible Runner

Assessment of Artificial Fluid Properties for High-Order Accurate Large-Eddy Simulations of Shock-Free Compressible Turbulent Flows with Strong Temperature Gradients

Application of Artificial Fluid Properties for Stable and Accurate Large-Eddy Simulations of Compressible Turbulent Reactive Flows

Application Of High-Order Compact Difference Scheme In The Computation Of Incompressible Wall-Bounded Turbulent Flows

A Higher Order Compact Finite Difference Algorithm for Solving the Incompressible Navier–Stokes Equations

Essentially compact schemes for unsteady viscous incompressible flows

Turbulent Flow Simulation Using Large Eddy Simulation Combined with Characteristic-Based Split Scheme

Numerical Analysis of a High-Order Unstructured Overset Grid Method for Compressible LES of Turbomachinery.

A high-order upwind compact difference scheme for solving the streamfunction-velocity formulation of the unsteady incompressible Navier–Stokes equations

A Fourth‐order Compact Finite Difference Scheme for the Steady Stream Function–vorticity Formulation of the Navier–Stokes/Boussinesq Equations

A Fourth-Order Compact Finite Volume Method on Unstructured Grids for Simulation of Two Dimensional Incompressible Flow

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

Random large eddy simulation for 3-dimensional incompressible viscous flows

Development of an accurate central finite-difference scheme with a compact stencil for the simulation of unsteady incompressible flows on staggered orthogonal grids

A Nonlinear Sub-grid Scale Model for Compressible Turbulent Flow

N A ] 1 5 D ec 2 01 4 A Compact Third-order Gas-kinetic Scheme for Compressible Euler and Navier-Stokes Equations

Two-stage Fourth-Order Subcell Finite Volume Method on Hexahedral Meshes for Compressible Flows

An Upwind Compact Difference Scheme for Solving the Streamfunction-Velocity Formulation of the Unsteady Incompressible Navier-Stokes Equation