An efficient high-order least square-based finite difference-finite volume method for solution of compressible Navier-Stokes equations on unstructured grids. (30th May 2021)

• Record Type: Journal Article
• Authors: Liu, Y.Y.; Shu, C.; Zhang, H.W.; Yang, L.M.; Lee, Cunbiao
• Is Part Of: Computers & fluids. Volume 222(2021)
• Publication Date: 2021-05-30
• Subjects: High-order finite volume method -- Least square-based finite difference -- Gas-kinetic flux solver -- Compressible viscous flows -- Unstructured grids; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2021.104926
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: A high-order finite volume method is developed for compressible viscous flows. A meshless method is applied to approximate derivatives. Unstructured grids can be easily applied. Inviscid and viscous fluxes are evaluated physically and simultaneously. Higher accuracy than the second-order method is achieved. Abstract: This paper presents an efficient high-order finite volume method for solution of compressible Navier-Stokes equations on unstructured grids. In this method, a high-order polynomial which is based on Taylor series expansion is applied to approximate the solution function within each control cell. The derivatives in the Taylor series expansion are approximated by the functional values at the cell centers of the considered cell and its neighboring cells using the mesh-free least square-based finite difference (LSFD) scheme. Naturally, this least square-based finite difference-finite volume (LSFD-FV) method inherits appealing characteristics of the LSFD scheme, i.e., the simple algorithm and easy implementation. In addition, since the LSFD scheme is mesh-free in nature, the developed high-order LSFD-FV method is endowed with the flexibility to handle the multi-dimensional problems with complex geometries on arbitrary grids. Different from other high-order methods, this LSFD-FV method applies a novel and effective strategy, i.e., the discrete gas-kinetic flux solver (DGKFS), to compute numerical fluxes at the cell interface. In this way, the inviscid and viscous fluxes are simultaneously and effectively evaluated by local reconstruction of solutions for the Boltzmann equation. The efficient time marching strategy coupled with the high-order LSFD-FV method is adopted to solve the resultant ordinary differential equations with the matching accuracy. A series of numerical examples are presented to validate the accuracy, robustness and flexibility of the developed method on unstructured grids. Numerical results demonstrate the superior performance of the developed high-order method on simulating compressible viscous flows, in comparison with the low-order counterpart and the k -exact method of the same order of accuracy.