High-order Gas-kinetic Scheme Research Articles

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

In this paper, the three-dimensional fully compact high-order gas-kinetic scheme (HGKS) is proposed for the direct numerical simulation of compressible turbulent flows. Because of the high-order gas evolution model, the numerical fluxes as well as the point-wise conservative variables can be evaluated from the time-accurate gas distribution function at the cell interface. As a result, both the cell-averaged variables and their cell-averaged gradients can be updated inside each cell. Based on the cell averaged values and their gradients, the compact Hermite weighted essentially non-oscillatory (HWENO) scheme is developed, in which the dimension-by-dimension reconstruction is used for three-dimensional turbulences. In both normal and tangential directions, the fifth-order HWENO reconstruction is adopted. Compared with the classical WENO scheme, the stencil for the HWENO scheme only contains 33 cells for each cell. To achieve the temporal accuracy, the two-stage fourth-order temporal discretization is used. For the evaluation of point-wise variables, the simplified third-order gas-kinetic solver is used. Several classical benchmark problems are simulated, which validate the accuracy, resolution, and robustness of compact HGKS. As a comparison, the numerical results of HGKS using non-compact WENO reconstruction are also provided. Due to the compact stencil, the compact HGKS has a favorable performance for turbulence simulation in resolving the multi-scale structures.

High-order compact gas-kinetic scheme for two-layer shallow water equations on unstructured mesh

For the two-layer shallow water equations, a high-order compact gas-kinetic scheme (GKS) on triangular mesh is proposed. The two-layer shallow water equations have complex source terms in comparison with the single layer equations. The main focus of this study is to construct a time-accurate evolution solution at a cell interface and to design a well-balanced scheme. The evolution model at a cell interface provides not only the numerical fluxes, but also the flow variables. In the development of the well-balanced scheme, the time-dependent flow variables at the cell interfaces surrounding the cell can be utilized to update the cell-averaged gradients. These time-evolving gradients are subsequently applied in the discretization of source terms inside each control volume. Based on the cell-averaged flow variables and their gradients, high-order initial data reconstruction can be achieved with compact stencils. The compact high-order GKS has advantages to simulate the flow evolution in complex domain covered by unstructured mesh. Many test cases are used to validate the accuracy and robustness of the scheme for the two-layer shallow water equations.

A one-stage high-order gas-kinetic scheme for multi-component flows with interface-sharpening technique

We propose a one-stage high-order gas-kinetic scheme (GKS) with interface-sharpening technique for the 4-equations model of compressible multi-component stiffened-gas flows. The scheme is constructed based on the framework of a one-stage efficient high-order GKS (EHGKS), which can achieve uniform high-order accuracy in both space and time. The main idea in the construction of the new scheme consists of three parts. Firstly, different from the original EHGKS, the non-oscillatory kinetic (NOK) flux is introduced to compute the leading term and the related flux evaluations, which enables the modified EHGKS suitable for stiffened gases, together with no pressure and velocity oscillations across the contact discontinuities. Then the modified scheme is utilized to discrete the conservative parts of the 4-equations model. Secondly, in order to avoid the oscillations across material interfaces, a high-order GKS is derived to solve the non-conservative parts of the 4-equations model following the idea of Abgrall [1]. Finally, to further reduce the over diffusion near a material interface, a simple interface-sharpening technique is developed, which combines the idea of interface construction and downwind scheme. The numerical results demonstrate that the proposed scheme can not only achieve the designed uniform accuracy in both space and time, but also sharply capture material interfaces.

High-order gas-kinetic schemes with non-compact and compact reconstruction for implicit large eddy simulation

High-order gas-kinetic scheme (HGKS) with 5th-order non-compact reconstruction has been well implemented for implicit large eddy simulation (ILES) in nearly incompressible turbulent channel flows. In this study, the HGKS with higher-order non-compact reconstruction and compact reconstruction will be validated in turbulence simulation. For higher-order non-compact reconstruction, 7th-order reconstruction in both normal and tangential directions are implemented. On the other hand, the compact HGKS with 5th-order compact reconstruction in normal direction is developed. Current work aims to show the benefits of high-order non-compact reconstruction and compact reconstruction for ILES. The accuracy of HGKS is verified by numerical simulation of three-dimensional advection of density perturbation. For the non-compact 7th-order scheme, 16 Gaussian points are required on the cell-interface to preserve the order of accuracy. Then, HGKS with non-compact and compact reconstruction is used in the three-dimensional Taylor–Green vortex (TGV) problem and turbulent channel flows. Accurate ILES solutions have been obtained from HGKS. In terms of the physical modeling underlying the numerical algorithms, the compact reconstruction has the consistent physical and numerical domains of dependence without employing additional information from cells which have no any direct physical connection with the targeted cell. The compact GKS shows a favorable performance for turbulence simulation in resolving the multi-scale structures.

Multiple-GPU accelerated high-order gas-kinetic scheme for direct numerical simulation of compressible turbulence

High-order gas-kinetic scheme (HGKS) has become a workable tool for the direct numerical simulation (DNS) of turbulence. In this paper, to accelerate the computation, HGKS is implemented with the graphical processing unit (GPU) using the compute unified device architecture (CUDA). Due to the limited available memory size, the computational scale is constrained by single GPU. For large-scale DNS of turbulence, we develop a multi-GPU HGKS simulation using message passing interface (MPI) and CUDA. The benchmark cases for compressible turbulence, including Taylor-Green vortex and turbulent channel flows, are presented to assess the numerical performance of HGKS with Nvidia TITAN RTX and Tesla V100 GPUs. For single-GPU computation, compared with the parallel central processing unit (CPU) code running on the Intel Core i7-9700 with open multi-processing (OpenMP) directives, 7x speedup is achieved by TITAN RTX and 16x speedup is achieved by Tesla V100. For multiple-GPU computation, multiple-GPU accelerated HGKS code scales properly with the increasing number of GPU. The computational time of parallel CPU code running on 1024 Intel Xeon E5-2692 cores with MPI is approximately 3 times longer than that of GPU code using 8 Tesla V100 GPUs with MPI and CUDA. Numerical results confirm the excellent performance of multiple-GPU accelerated HGKS for large-scale DNS of turbulence. Besides reducing memory access pressure, we also exploit single precision floating point arithmetic to accelerate HGKS on GPUs. Reasonably, compared to the computation with FP64 precision, the efficiency is improved and the memory cost is reduced with FP32 precision. Meanwhile, the differences in accuracy for statistical turbulent quantities appear. For turbulent channel flows, difference in long-time statistical turbulent quantities is acceptable between FP32 and FP64 precision solutions. While the obvious discrepancy in instantaneous turbulent quantities can be observed, which shows that FP32 precision is not safe for DNS in compressible turbulence. The choice of precision should be depended on the requirement of accuracy and the available computational resources.

Three-dimensional discontinuous Galerkin based high-order gas-kinetic scheme and GPU implementation

In this paper, the discontinuous Galerkin based high-order gas-kinetic schemes (DG-HGKS) are developed for the three-dimensional Euler and Navier–Stokes equations. Different from the traditional discontinuous Galerkin (DG) methods with Riemann solvers, the current method adopts a kinetic evolution process, which is provided by the integral solution of Bhatnagar–Gross–Krook (BGK) model. In the weak formulation of DG method, a time-dependent evolution function is provided, and both inviscid and viscous fluxes can be calculated uniformly. The temporal accuracy is achieved by the two-stage fourth-order discretization, and the second-order gas-kinetic solver is adopted for the fluxes over the cell interface and the fluxes inside a cell. Numerical examples, including accuracy tests and Taylor–Green vortex problem, are presented to validate the efficiency and accuracy of DG-HGKS. Both optimal convergence and super-convergence are achieved by the current scheme. The comparison between DG-HGKS and high-order gas-kinetic scheme with weighted essential non-oscillatory reconstruction (WENO-HGKS) is also given, and the numerical performances are comparable with the approximate number of degree of freedom. To accelerate the computation, the DG-HGKS is implemented with the graphics processing unit (GPU) using compute unified device architecture (CUDA). The obtained results are also compared with those calculated by the central processing units (CPU) code in terms of the computational efficiency. The speedup of GPU code suggests the potential of high-order gas-kinetic schemes for the large scale computation.

Multiple-Gpu Accelerated High-Order Gas-Kinetic Scheme for Direct Numerical Simulation of Compressible Turbulence

High-order gas-kinetic scheme (HGKS) has become a workable tool for the direct numerical simulation (DNS) of turbulence. In this paper, to accelerate the computation, HGKS is implemented with the graphical processing unit (GPU) using the compute unified device architecture (CUDA). To conduct the much large-scale DNS of turbulence, HGKS also be further upgraded with multiple GPUs using message passing interface (MPI) and CUDA architecture. The benchmark cases for compressible turbulence, including Taylor-Green vortex and turbulent channel flows, are presented to assess the numerical performance of HGKS with Nvidia TITAN RTX and Tesla V100 GPUs. For single-GPU computation, compared with the parallel central processing unit (CPU) code running on the Intel Core i7-9700 with open multi-processing (OpenMP) directives, 7x speedup is achieved by TITAN RTX and 16x speedup is achieved by Tesla V100. For multiple-GPU computation, the computational time of parallel CPU code running on 1024 Intel Xeon E5-2692 cores with MPI is approximately 3 times longer than that of GPU code using 8 Tesla V100 GPUs with MPI and CUDA. Numerical results confirm the excellent performance of multiple-GPU accelerated HGKS for large-scale DNS of turbulence. HGKS in GPU is also compiled with FP32 precision to evaluate the effect of number formats precision. Reasonably, compared to the computation with FP64 precision, the efficiency is improved and the memory cost is reduced with FP32 precision. For turbulent channel flows, difference in long-time statistical turbulent quantities is acceptable between FP32 and FP64 precision solutions. While the obvious discrepancy in instantaneous turbulent quantities can be observed, which shows that FP32 precision is not safe for DNS in compressible turbulence. The choice of precision should depended on the requirement of accuracy and the available computational resources.

High-order gas-kinetic scheme for large eddy simulation of turbulent channel flows

In this study, the high-order gas kinetic scheme (GKS) is employed for explicit large eddy simulation (hereafter referred to simply as “LES”) and implicit large eddy simulation (iLES) of turbulent channel flows. The main objective is to compare the performance of iLES and LES in the high-order finite volume framework, and study which is most suitable for turbulence simulation. The prediction abilities of different explicit LES models and iLES method on various statistics and flow structures are compared. Most results from both iLES and explicit LES are very close to those of direct numerical simulation. Moreover, iLES is generally superior to explicit LES in predicting several important flow properties, including the mean velocity profiles, Reynolds stress, and Q-criterion iso-surfaces. This superior performance of iLES may arise that the numerical dissipation of the high-order scheme is enough to replace the subgrid dissipation needed in large eddy simulation of turbulence. If the explicit LES model is further adopted, the dissipation will increase, so the results degrade. The overall satisfactory results show that the high-order GKS can provide appropriate numerical dissipation and is suitable for iLES of turbulence.

A compact high-order gas-kinetic scheme on unstructured mesh for acoustic and shock wave computations

Following the development of a third-order compact gas-kinetic scheme (GKS) for the Euler and Navier-Stokes equations (Ji et al. (2020) [40]), in this paper a higher-order compact GKS up to sixth order of accuracy will be constructed for the shock and acoustic wave computation on unstructured mesh. The compactness is defined by the physical domain of dependence for an unstructured triangular cell, which may involve the closest neighbors of neighboring cells. The compactness and high-order accuracy of the scheme are coming from the consistency between the high-order initial reconstruction and the high-order gas evolution model under GKS framework. The high-order evolution solution at a cell interface provides not only a time-accurate flux function, but also the time-evolving flow variables. Therefore, the cell-averaged flow variables and their gradients can be explicitly updated at the next time level from the moments of the same time-dependent gas distribution function. Based on the cell averaged flow variables and their cell-averaged derivatives, the nonlinear high-order compact reconstruction can be obtained for these variables in the evaluation of local equilibrium and non-equilibrium states. The current nonlinear reconstruction is a combination of WENO and ENO methodology, which is specifically suitable for compact GKS on unstructured mesh with a high-order (≥4) accuracy. The initial piecewise discontinuous reconstruction is used for the determination of non-equilibrium state and an evolved smooth reconstruction for the equilibrium state. The evolution model in GKS is based on a relaxation process from non-equilibrium to equilibrium state. The time-accurate gas distribution function provides the Navier-Stokes (NS) flux function directly without separating the inviscid and viscous terms, which simplifies the numerical discretization of the NS equations directly on unstructured mesh. Based on the time-accurate flux solver, the two-stage fourth-order time discretization can be applied to get a fourth-order time-accurate solution with only two stages, which reduces two reconstructions in comparison with the same order Runge-Kutta time stepping method. The current high-order GKS can uniformly capture acoustic and shock waves without identifying troubled cells and implementing additional limiting procedure. In addition, the fourth- up to sixth-order compact GKS can use almost the same time step as a second-order shock capturing scheme in most flow problems. The fourth-order GKS on unstructured mesh will be used in the computations from low speed incompressible viscous flow to the high Mach number shock interaction. The accuracy, efficiency, and robustness of the scheme have been validated. The main conclusion of the paper is that beyond the first-order Riemann solver, the use of high-order gas evolution model seems have advantages in the development of high-order schemes.

High-order gas-kinetic scheme with parallel computation for direct numerical simulation of turbulent flows

The performance of high-order gas-kinetic scheme (HGKS) has been investigated for the direct numerical simulation (DNS) of isotropic compressible turbulence up to the supersonic regime [9]. Due to the multi-scale nature and coupled temporal-spatial evolution process, HGKS provides a valid tool for the numerical simulation of compressible turbulent flow. Based on the domain decomposition and message passing interface (MPI), a parallel HGKS code is developed for large-scale computation in this paper. The standard tests from the nearly incompressible flow to the supersonic one, including Taylor-Green vortex problem, turbulent channel flow and isotropic compressible turbulence, are presented to validate the parallel scalability, efficiency, accuracy and robustness of parallel implementation. The performance of HGKS for the nearly incompressible turbulence is comparable with the high-order finite difference scheme, including the resolution of flow structure and efficiency of computation. Based on the accuracy of the numerical solution, the numerical dissipation of the scheme in the turbulence simulation is quantitatively evaluated. Meanwhile, based on the kinetic formulation HGKS shows advantage for supersonic turbulent flow simulation with its accuracy and robustness. The current work demonstrates the capability of HGKS as a powerful DNS tool from the low speed to supersonic turbulence study, which is less reported under the framework of finite volume scheme.

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

In this paper, a high-order gas-kinetic scheme is developed on three-dimensional unstructured meshes for compressible Euler and Navier–Stokes equations. To achieve the high-order spatial accuracy, the three-dimensional weighted essentially non-oscillatory (WENO) reconstruction is extended to the unstructured tetrahedral and hexahedral meshes. A simple strategy is adopted for the selection of candidate stencils, and the topologically independent linear weights are used for the spatial reconstruction. The efficiency and robustness of the classical WENO reconstruction are improved. In addition to the two-stage fourth-order temporal discretization and lower–upper symmetric Gauss–Seidel method, the explicit and implicit high-order gas-kinetic schemes are developed for unsteady and steady problems. Accuracy tests on hexahedral and tetrahedral grids validate the third-order of accuracy, and various three-dimensional incompressible and compressible numerical experiments are also presented. The results validate the accuracy and robustness of the proposed scheme for both inviscid and viscous flows. In the future, the current scheme will be extended to the hybrid unstructured meshes and Reynolds-averaged Navier–Stokes simulation with high Reynolds numbers.

A compact and efficient high-order gas-kinetic scheme

In this paper, a compact fifth-order gas-kinetic scheme (GKS) is proposed to solve the Euler equations of compressible flows. It is based on the framework of a one-stage efficient high-order GKS (EHGKS), which can achieve arbitrary high-order accuracy in both space and time, and the newly-developed compact reconstruction technique, combining both the ideas of the simple weighted essentially non-oscillatory scheme (WENO) and Hermite WENO (HWENO). The evolution model of EHGKS provides the time-dependent variables at cell interfaces, which can directly compute the cell-average slope. Thus the necessary information for a fifth-order compact HWENO reconstruction technique is obtained, and also making it possible to improve simple WENO as a compact reconstruction by determining the quartic polynomial based on the compact stencil. The one-stage time discretization technique of EHGKS makes it possible to utilize the reconstruction only once in each computing time step, which is very different from the present compact multi-stage GKS schemes. The newly-developed compact reconstruction provides the polynomial, not only the pointwise variables but also their spatial derivatives for the interface variables evaluations directly. Several numerical examples are given to validate the accuracy and effectiveness of the proposed compact one-stage scheme. And the numerical results also demonstrate that the new scheme is more efficient than the original EHGKS method with WENO reconstruction technique.

The study of shallow water flow with bottom topography by high-order compact gas-kinetic scheme on unstructured mesh

A well-balanced compact high-order gas-kinetic scheme (GKS) on unstructured mesh is first developed for solving the shallow water equations with source terms. The distinguishable feature of the finite volume GKS is that based on the gas-kinetic formulation, a time-accurate gas distribution function can be constructed, from which both the fluxes and the flow variables can be explicitly evaluated at the cell interface. As a result, besides the update of cell-averaged conservative variables, the cell-averaged slopes of the flow variables can be updated as well. Equipped with both flow variables and their slopes, a fourth-order compact spatial reconstruction on unstructured mesh can be obtained as the initial condition at the beginning of each time step. For the shallow water flow, in order to preserve the well-balanced property, the advection and the source terms in the flux function have to be balanced properly. The current compact GKS achieves high-order accuracy, keeps the well-balanced property, and has super-robustness in the simulation of bore waves. The scheme is used in the shallow water flow studies, such as dam breaking and bore wave propagation. In addition, the pollution transport, morphodynamics, and bottom friction in the shallow water flow have been included in the scheme. In the end, the water discharge in the Pearl River estuary and the dam-break experiment with movable bed topography have been simulated.

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

In the previous study [J. Comput. Phys. 417 (2020) 109558], under arbitrary-Lagrangian-Eulerian (ALE) formulation a high-order gas-kinetic scheme has been developed for the computation of two-dimensional flows. For the three-dimensional flows, due to the distorted mesh, it becomes more difficult to develop robust high-order ALE methods with the precise preservation of geometric conservation law. In this paper, the high-order gas-kinetic ALE scheme will be constructed for three-dimensional flows. The key ingredients of the scheme are the use of weighted essentially non-oscillatory (WENO) scheme for spatial reconstruction and the two-stage fourth-order discretization for temporal evolution. In the ALE formulation, in order to release the problems associated with mesh distortion and non-coplanar vertexes of a control volume, in the spatial reconstruction the selection of candidate stencils and the topologically independent linear weights have to be carefully designed. In the surface integrals for the flux transport, a bilinear interpolation is used to parameterize both grid coordinates and grid moving velocity with the preservation of the geometric conservation law. In the computation, the grid velocity is determined by the variational formulation based Lagrangian nodal solver. Numerical examples are presented to evaluate the accuracy, robustness, and the preservation of geometric conservation law of the current scheme.

Three dimensional high-order gas-kinetic scheme for supersonic isotropic turbulence II: Coarse-graining analysis of compressible Ksgs budget

The direct numerical simulation (DNS) of compressible isotropic turbulence up to the supersonic regime with Mat=1.2 has been investigated by the high-order gas-kinetic scheme (HGKS) (Cao et al. (2019) [8]). In this study, the DNS on a much higher initial turbulent Mach number up to Mat=2.0 is performed by HGKS, which confirms the super robustness of HGKS. The coarse-graining analysis of subgrid-scale (SGS) turbulent kinetic energy Ksgs budget is fully analyzed for constructing one-equation SGS model in the compressible large eddy simulation (LES). The exact compressible SGS turbulent kinetic energy Ksgs transport equation is derived with density weighted filtering process. Based on the compressible Ksgs transport equation, the coarse-graining processes are implemented on three sets of unresolved grids with the Box filter. The coarse-graining analysis of compressible Ksgs budgets shows that all unresolved source terms are dominant in the current system. Especially, the magnitude of SGS pressure-dilation term is on the order of SGS solenoidal dissipation term within the initial acoustic time scale. Therefore, the SGS pressure-dilation term cannot be neglected as that in the previous work. The delicate coarse-graining analysis of SGS diffusion terms in compressible Ksgs equation also confirms that both the fluctuation velocity triple correlation term and the pressure-velocity correlation term are dominant terms. The current analysis provides an indication on the order of magnitude of all SGS terms in compressible Ksgs budget, which provides a solid basis for compressible LES modeling of high Mach number turbulent flow.

Fourth-order gas-kinetic scheme for turbulence simulation with multi-dimensional WENO reconstruction

In this paper, a fourth-order multi-dimensional weighted essentially non-oscillatory (WENO) reconstruction is developed, which is aiming at the complicated flows. In corporate with two-stage fourth-order temporal discretization, a high-order gas-kinetic scheme is proposed for the direct numerical simulation (DNS) and large eddy simulation (LES) of turbulence. In such WENO scheme, a simple strategy of selecting stencils is adopted and the topology independent linear weights are used. The fourth-order reconstruction improves the order of accuracy of the previous third-order WENO reconstruction [Computers and Fluids 198 (2020) 104401], which is too dissipative for the simulation of turbulent flows. Compared with the classical dimension-by-dimension WENO reconstruction, the order of accuracy can be also achieved for any non-uniform and curvilinear meshes in the finite volume framework. Numerical examples, including the classical turbulence cases, are provided to illustrate the good performance of such WENO scheme for turbulence simulation. More importantly, the fourth-order WENO scheme is robust and works well from the subsonic to hypersonic flows. In the future, the current scheme will be extended to the genuinely unstructured meshes for more complicated turbulent flows.

An Acoustic and Shock Wave Capturing Compact High-Order Gas-Kinetic Scheme with Spectral-Like Resolution

A compact high-order gas-kinetic scheme (GKS) with spectral resolution will be presented and used in the simulation of acoustic and shock waves. Based on high-order gas evolution model, the GKS provides a time accurate solution at a cell interface, from which both time accurate flux function and flow variables can be obtained for updating cell averaged conservative flow variables and gradients. Based on cell-averaged flow variables and gradients, a reconstruction with compact stencil up to eighth-order spacial accuracy can be achieved, which include the nonlinear reconstruction for initial non-equilibrium state and linear reconstruction for the equilibrium one. The GKS unifies the nonlinear and linear reconstruction through a time evolution process at a cell interface from the non-equilibrium state to an equilibrium one. In the smooth acoustic wave region, the linear reconstruction will contribute mainly to the evolution. In the non-equilibrium shock region, the nonlinear reconstruction will play a dominant role.

Performance Enhancement for High-Order Gas-Kinetic Scheme Based on WENO-Adaptive-Order Reconstruction

High-order gas-kinetic scheme (HGKS) has been well-developed in the past years. Abundant numerical tests including hypersonic flow, turbulence, and aeroacoustic problems, have been used to validate its accuracy, efficiency, and robustness. However, there are still rooms for its further improvement. Firstly, the reconstruction in the previous scheme mainly achieves a third-order accuracy for the initial non-equilibrium states. At the same time, the equilibrium state in space and time in HGKS has to be reconstructed separately. Secondly, it is complicated to get reconstructed data at Gaussian points from the WENO-type method in high dimensions. For HGKS, besides the point-wise values at the Gaussian points it also requires the slopes in both normal and tangential directions of a cell interface. Thirdly, there exists visible spurious overshoot/undershoot at weak discontinuities from the previous HGKS with the standard WENO reconstruction. In order to overcome these difficulties, in this paper we use an improved reconstruction for HGKS. The WENO with adaptive order (WENO-AO) method is implemented for reconstruction.A whole polynomial inside each cell is provided in WENO-AO reconstruction. The HGKS becomes simpler than the previous one with the direct implementation of cell interface values and their slopes from WENO-AO. The additional reconstruction of equilibrium state at the beginning of each time step can be avoided as well by dynamically merging the reconstructed non-equilibrium slopes. The new HGKS essentially releases or totally removes the above existing problems in previous HGKS. The accuracy of the scheme from 1D to 3D from the new HGKS can recover the theoretical order of accuracy of the WENO reconstruction.In the two- and three-dimensional simulations, the new HGKS shows better robustness and efficiency than the previous scheme in all test cases.

High-Order Gas-Kinetic Scheme in Curvilinear Coordinates for the Euler and Navier-Stokes Solutions

The high-order gas-kinetic scheme (HGKS) has achieved success in simulating compressible flow in Cartesian mesh. To study the flow problem in general geometry, such as the flow over a wing-body configuration, the development of a three-dimensional HGKS in general curvilinear coordinates becomes necessary. In this paper, a two-stage fourth-order gas-kinetic scheme is developed for the Euler and Navier-Stokes solutions in the curvilinear coordinates. Based on the coordinate transformation, the kinetic equation is transformed first to the computational space, and the flux function in the gas-kinetic scheme is obtained there and is transformed back to the physical domain for the update of conservative flow variables inside each control volume. To achieve the expected order of accuracy, the dimension-by-dimension reconstruction based on the WENO scheme is adopted in the computational domain, where the reconstructed variables are the cell averaged Jacobian and the Jacobian-weighted conservative variables, and the conservative variables are obtained by ratio of the above reconstructed data at Gaussian quadrature points of each cell interface. In the two-stage fourth-order gas kinetic scheme (GKS), similar to the generalized Riemann solver (GRP), the initial spatial derivatives of conservative variables have to be used in the evaluation of the time dependent flux function in GKS, which are reconstructed as well through orthogonalization in physical space and chain rule. A variety of numerical examples from the order tests to the solutions with strong discontinuities are presented to validate the accuracy and robustness of the current scheme. The precise satisfaction of the geometrical conservation law in non-orthogonal mesh is also demonstrated through the numerical example.

An efficient high-order gas-kinetic scheme (I): Euler equations

In this paper, an efficient high-order gas-kinetic scheme (EHGKS) is proposed to solve the Euler equations for compressible flows. We re-investigate the underlying mechanism of the high-order gas-kinetic scheme (HGKS) and find a new strategy to improve its efficiency. The main idea of the new scheme contains two parts. Firstly, inspired by the state-of-art simplifications on the third-order HGKS, we extend the HGKS to the case of arbitrary high-order accuracy and eliminate its unnecessary high-order dissipation terms. Secondly, instead of computing the derivatives of particle distribution function and their complex moments, we introduce a Lax-Wendroff procedure to compute the high-order derivatives of macroscopic quantities directly. The new scheme takes advantage of both HGKS and the Lax-Wendroff procedure, so that it can be easily extended to the case of arbitrary high-order accuracy with practical significance. Typical numerical tests are carried out by EHGKS, with the third, fifth and seventh order of accuracy. The presence of good resolution of the discontinuities and flow details, together with the optimal CFL numbers, validates the high accuracy and strong robustness of EHGKS. To compare the efficiency, we present the results computed by the EHGKS, the original HGKS and Runge-Kutta-WENO-GKS. This further demonstrates the advantages of EHGKS.