An Efficient and Accurate Two-Stage Fourth-order Gas-kinetic Scheme for the Navier-Stokes Equations

For computational fluid dynamics (CFD), the generalized Riemann problem (GRP) solver and the gas-kinetic kinetic scheme (GKS) provide a time-accurate flux function starting from a discontinuous piecewise linear flow distributions around each cell interface. With the use 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 [18]. In this paper, based on the same time-stepping method and the second-order GKS flux function [34], a fourth-order gas-kinetic scheme is constructed for the Euler and Navier-Stokes equations. In comparison with the formal one-stage time-stepping third-order gas-kinetic solver [21], the current fourth-order method not only reduces the complexity of the flux function, but also improves the accuracy of the scheme, even though the thirdand fourth-order schemes have similar computation cost. Most importantly, the robustness of the fourth-order GKS is as good as the second-order one. Perfect numerical solutions can be obtained from the high Reynolds number boundary layer solutions to the hypersonic viscous heat conducting flow computations. Many numerical tests, including many difficult ones for the Navier-Stokes solvers, have been used to validate the current fourth-order method. Following the two-stage time-stepping framework, the one-stage third-order GKS can be easily extended to 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 impact on the development of higher-order CFD methods.

[1]  Kun Xu,et al.  A Compact Third-order Gas-kinetic Scheme for Compressible Euler and Navier-Stokes Equations , 2014, 1412.4489.

[2]  V. Daru,et al.  Numerical simulation of the viscous shock tube problem by using a high resolution monotonicity-preserving scheme , 2009 .

[3]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[4]  Huazhong Tang,et al.  A High-Order Accurate Gas-Kinetic Scheme for One- and Two-Dimensional Flow Simulation , 2014 .

[5]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[6]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[7]  Sharath S. Girimaji,et al.  WENO-enhanced gas-kinetic scheme for direct simulations of compressible transition and turbulence , 2013, J. Comput. Phys..

[8]  Michael Dumbser,et al.  A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes , 2008, J. Comput. Phys..

[9]  Kun Xu,et al.  A third-order compact gas-kinetic scheme on unstructured meshes for compressible Navier-Stokes solutions , 2016, J. Comput. Phys..

[10]  Allan R. Wieting,et al.  Experimental Study of Shock Wave Interference Heating on a Cylindrical Leading Edge , 1987 .

[11]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[12]  L. Jun,et al.  A high-order multidimensional gas-kinetic scheme for hydrodynamic equations , 2013 .

[13]  Kun Xu,et al.  Comparison of the generalized Riemann solver and the gas-kinetic scheme for inviscid compressible flow simulations , 2011, J. Comput. Phys..

[14]  J. Falcovitz,et al.  A second-order Godunov-type scheme for compressible fluid dynamics , 1984 .

[15]  Domenic D'Ambrosio,et al.  Numerical Instablilities in Upwind Methods: Analysis and Cures for the “Carbuncle” Phenomenon , 2001 .

[16]  Kun Xu,et al.  A high-order multidimensional gas-kinetic scheme for hydrodynamic equations , 2013 .

[17]  Chongam Kim,et al.  Accurate, efficient and monotonic numerical methods for multi-dimensional compressible flows Part I Spatial discretization , 2005 .

[18]  R. H. Fowler The Mathematical Theory of Non-Uniform Gases , 1939, Nature.

[19]  P. Bhatnagar,et al.  A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems , 1954 .

[20]  Feng Xiao,et al.  High order multi-moment constrained finite volume method. Part I: Basic formulation , 2009, J. Comput. Phys..

[21]  Jiequan Li,et al.  A Two-Stage Fourth Order Time-Accurate Discretization for Lax-Wendroff Type Flow Solvers I. Hyperbolic Conservation Laws , 2015, SIAM J. Sci. Comput..

[22]  Gerald Warnecke,et al.  A direct Eulerian GRP scheme for compressible fluid flows , 2006, J. Comput. Phys..

[23]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[24]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[25]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[26]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[27]  Ee Han,et al.  Accuracy of the Adaptive GRP Scheme and the Simulation of 2-D Riemann Problems for Compressible Euler Equations , 2011 .

[28]  Xu-Dong Liu,et al.  Solution of Two-Dimensional Riemann Problems of Gas Dynamics by Positive Schemes , 1998, SIAM J. Sci. Comput..

[29]  Kun Xu,et al.  Microchannel flow in the slip regime: gas-kinetic BGK–Burnett solutions , 2004, Journal of Fluid Mechanics.

[30]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[31]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[32]  Shuli,et al.  The two-dimensional Riemann problem in gas dynamics , 1998 .

[33]  Kun Xu,et al.  A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method , 2001 .

[34]  Jiequan Li,et al.  The generalized Riemann problems for compressible fluid flows: Towards high order , 2014, J. Comput. Phys..

[35]  Kun Xu,et al.  Super-Burnett solutions for Poiseuille flow , 2003 .

[36]  Kun Xu,et al.  Direct modeling for computational fluid dynamics , 2015, Acta Mechanica Sinica.

[37]  Xiaodong Ren,et al.  A multi-dimensional high-order discontinuous Galerkin method based on gas kinetic theory for viscous flow computations , 2015, J. Comput. Phys..

[38]  Song Fu,et al.  A high-order gas-kinetic Navier-Stokes flow solver , 2010, J. Comput. Phys..