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.

Explicit two‐step Runge‐Kutta methods for computational fluid dynamics solvers

A RKDG Finite Element Method on Unstructured Grids for Simulating the Compressible Multimedia Fluids

Efficient implementation of complex equations of state in a high-order framework

Fast linear solvers for incompressible CFD simulations with compatible discrete operator schemes

A Numerical Investigation of Matrix-Free Implicit Time-Stepping Methods for Large CFD Simulations

A Parallel Runge–Kutta Discontinuous Galerkin Solver for Rarefied Gas Flows Based on 2D Boltzmann Kinetic Equations

Fourth-Order Runge–Kutta Schemes for Fluid Mechanics Applications

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

TWO BETTER LOOSELY COUPLED SOLUTION ALGORITHMS OF CFD BASED AEROELASTIC SIMULATION

High order implicit time integration schemes on multiresolution adaptive grids for stiff PDEs

MATRICS: The implicit matrix-free Eulerian hydrodynamics solver

On the accuracy of stiff-accurate diagonal implicit Runge-Kutta methods for finite volume based Navier-Stokes equations

A Runge-Kutta Discontinuous Galerkin Solver for 2D Boltzmann Model Equations: Verification and Analysis of Computational Performance

A PARALLEL SOLVING DESIGN FOR UNSTEADY NAVIER-STOKES EQUATION

A Mathematical Solution to the Computational Fluid Dynamics (CFD) Dilemma

A Scalable Fully Implicit Method with Adaptive Time Stepping for Unsteady Compressible Inviscid 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 Numerical Study for the Performance of the Runge-Kutta Discontinuous Galerkin Method Based on Different Numerical Fluxes

A Strongly Coupled Time-Marching Method for Solving the Navier–Stokes andk-ω Turbulence Model Equations with Multigrid

Recurrence CFD - a novel approach to simulate multiphase flows with strongly separated time scales

Efficient time stepping for reactive turbulent simulations with stiff chemistry