Hostname: page-component-cd9895bd7-dzt6s Total loading time: 0 Render date: 2024-12-29T11:52:00.115Z Has data issue: false hasContentIssue false
  • English
  • Français

Comparison of Fifth-Order WENO Scheme and Finite Volume WENO-Gas-Kinetic Scheme for Inviscid and Viscous Flow Simulation

Published online by Cambridge University Press:  03 June 2015

Jun Luo ,
Lijun Xuan  and
Kun Xu
Get access

Abstract

The development of high-order schemes has been mostly concentrated on the limiters and high-order reconstruction techniques. In this paper, the effect of the flux functions on the performance of high-order schemes will be studied. Based on the same WENO reconstruction, two schemes with different flux functions, i.e., the fifth-order WENO method and the WENO-Gas-kinetic scheme (WENO-GKS), will be compared. The fifth-order finite difference WENO-SW scheme is a characteristic variable reconstruction based method which uses the Steger-Warming flux splitting for inviscid terms, the sixth-order central difference for viscous terms, and three stages Runge-Kutta time stepping for the time integration. On the other hand, the finite volume WENO-GKS is a conservative variable reconstruction based method with the same WENO reconstruction. But, it evaluates a time dependent gas distribution function along a cell interface, and updates the flow variables inside each control volume by integrating the flux function along the boundary of the control volume in both space and time. In order to validate the robustness and accuracy of the schemes, both methods are tested under a wide range of flow conditions: vortex propagation, Mach 3 step problem, and the cavity flow at Reynolds number 3200. Our study shows that both WENO-SW and WENO-GKS yield quantitatively similar results and agree with each other very well provided a sufficient grid resolution is used. With the reduction of mesh points, the WENO-GKS behaves to have less numerical dissipation and present more accurate solutions than those from the WENO-SW in all test cases. For the Navier-Stokes equations, since the WENO-GKS couples inviscid and viscous terms in a single flux evaluation, and the WENO-SW uses an operator splitting technique, it appears that the WENO-SW is more sensitive to the WENO reconstruction and boundary treatment. In terms of efficiency, the finite volume WENO-GKS is about 4 times slower than the finite differenceWENO-SW in two dimensional simulations. The current study clearly shows that besides high-order reconstruction, an accurate gas evolution model or flux function in a high-order scheme is also important in the capturing of physical solutions. In a physical flow, the transport, stress deformation, heat conduction, and viscous heating are all coupled in a single gas evolution process. Therefore, it is preferred to develop such a scheme with multi-dimensionality, and unified treatment of inviscid and dissipative terms. A high-order scheme does prefer a high-order gas evolution model. Even with the rapid advances of high-order reconstruction techniques, the first-order dynamics of the Riemann solution becomes the bottleneck for the further development of high-order schemes. In order to avoid the weakness of the low order flux function, the development of high-order schemes relies heavily on the weak solution of the original governing equations for the update of additional degree of freedom, such as the non-conservative gradients of flow variables, which cannot be physically valid in discontinuous regions.

Keywords

02.60Cb47.11.Df47.45.AbWENO schemegas-kinetic schemeEuler equationsNavier-Stokes equationshighorder methods
Type
Research Article
Information
Communications in Computational Physics , Volume 14 , Issue 3 , September 2013 , pp. 599 - 620
Copyright
Copyright © Global Science Press Limited 2013

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Ben-Artzi, M., Li, J. and Warnecke, G., A direct Eulerian GRP scheme for compressible fluid flows, J. Comput. Phys., 218 (2006), pp. 19–43.CrossRefGoogle Scholar
[2]Bhatnagar, P.L., Gross, E.P., and Krook, M., A Model for Collision Processes in Gases I: Small Amplitude Processes in Charged and Neutral One-Component Systems, Phys. Rev., 94 (1954), pp. 511–525.Google Scholar
[3]Casper, J., Finite-volume implementation of high-order essentially non-oscillatory schemes in two dimensions, AIAA Journal, 30 (1992), pp. 2829–2835.Google Scholar
[4]Deng, X.G., Mao, M.L., Tu, G.H., Zhang, H.X. and Zhang, Y.F., High-order and high accurate CFD methods and their applications for complex grid problems, Commun. Comput. Phys., 11 (2012), pp. 1081–1102.Google Scholar
[5]Ghia, U, Ghia, K.N, Shin, C.T, High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method, J. Comput. Phys., 48 (1982), pp. 387–411.Google Scholar
[6]Guo, Z.L., Liu, H.W., Luo, L.S., and Xu, K., A comparative study of the LBM and GKS methods for 2D near incompressible flows, J. Comput. Phys. 227 (2008), pp. 4955–4976.CrossRefGoogle Scholar
[7]Harten, A., Engquist, B., Osher, S., Chakravarthy, S., Uniformly high order Essentially Non-Oscillatory schemes III, J. Comput. Phys. 71 (1987) 231C303.CrossRefGoogle Scholar
[8]Huang, J.C., Linb, H., Yang, J.Y., Implicit preconditioned WENO scheme for steady viscous flow computation, J. Comput. Phys., 228 (2009), pp. 420–438.CrossRefGoogle Scholar
[9]Huang, J.C., Xu, K. and Yu, P.B., A unified gas-kinetic scheme for continuum and rarefied flows II: multi-dimensional cases., Commun. Comput. Phys., 12 (2012), pp. 662–690.CrossRefGoogle Scholar
[10]Jiang, G.S., Shu, C.W., Efficient implementation of Weighted ENO schemes, J. Comput. Phys. 126 (1996) 202C228.Google Scholar
[11]Li, J.Q., Li, Q.B., Xu, K., Comparison of the Generalized Riemann Solver and the Gas-Kinetic Scheme for Inviscid Compressible Flow Simulations, J. Comput. Phys., 230 (2011), pp. 5080–5099.Google Scholar
[12]Li, Q.B., Xu, K., and Fu, S., A high-order gas-kinetic Navier-Stokes solver, J. Comput. Phys., 229 (2010), pp. 6715–6731.Google Scholar
[13]Li, W., Ren, Y.X., High-order k-exact WENO finite volume schemes for solving gas dynamic Euler equations on unstructured grids, Int. J.Numer. Meth. Fluids, (2011).Google Scholar
[14]Liu, X.D., Osher, S., Weighted essentially non-oscillatory schemes, J. Comput. Phys., 115 (1994), pp. 200–212.Google Scholar
[15]Luo, J., Xu, K., A high-order WENO-gas-kinetic scheme for hydrodynamic equations, preprint (2012).Google Scholar
[16]Ohwada, T. and Fukata, S., Simple derivation of high-resolution schemes for compressible flows by kinetic approach, J. Comput. Phys. 211 (2006), pp. 424.CrossRefGoogle Scholar
[17]Qiu, J.M. and Shu, C.W., Conservative semi-Lagrangian finite difference WENO formulations with applications to the Vlasov equation, Commun. Comput. Phys., 10 (2011), pp. 979–1000.Google Scholar
[18]Ren, Y.X., Liu, M., Zhang, H.X., A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws, J. Comput. Phys., 192 (2003), pp. 365–386.Google Scholar
[19]Shen, Y.Q., Zha, G.C.,Low diffusion E-CUSP scheme with implicit high order WENO scheme for preconditioned Navier-Stokes equations, Comput Fluids, 50(2012), pp. 13–23.Google Scholar
[20]Shu, C.W., Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, Lecture Notes in Mathematics, Springer, 1998.Google Scholar
[21]Shu, C.W., High-order finite difference and finite volume WENO schemes and discontinuous Galerkin methods for CFD, International Journal of Computational Fluid Dynamics, 17 (2003), pp. 107–118.CrossRefGoogle Scholar
[22]Shu, C.W. and Osher, S., Efficient implementation of essentially nonoscillatory shock-capturing schemes II, J. Comput. Phys., 83 (1989), pp. 32–78.CrossRefGoogle Scholar
[23]Steger, J.L. and Warming, R.,Flux vector splitting of the inviscid gas dynamic equation with application to finite difference methods, J Compt. Phys, 40(1981), pp. 263–293.Google Scholar
[24]Su, M.D., Xu, K., Ghidaoui, M., Low Speed Flow Simulation by the Gas-kinetic Scheme, J. Comput. Phys., 150 (1999), pp. 17–39.CrossRefGoogle Scholar
[25]Woodward, P. and Colella, P., Numerical simulations of two-dimensional fluid flow with strong shocks, J. Comput. Phys., 54 (1984), pp. 115–173.Google Scholar
[26]Xu, K., A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method, J. Comput. Phys., 171 (2001), pp. 289–335.Google Scholar
[27]Xu, K. and He, X.Y., Lattice Boltzmann method and Gas-kinetic BGK scheme in the low Mach number viscous flow simulations, J. Comput. Phys. 190 (2003), pp. 100–117.Google Scholar
[28]Xu, K. and Huang, J.C., A unified gas-kinetic scheme for continuum and rarefied flows, J. Comput. Phys. 229 (2010), pp. 7747–7764.Google Scholar
[29]Xu, K. and Josyula, E., Continuum formulation for non-equilibrium shock structure calculation, Commun.Comput. Phys. 1 (2006), pp. 425–448.Google Scholar