A high order least square-based finite difference-finite volume method with lattice Boltzmann flux solver for simulation of incompressible flows on unstructured grids

Abstract This paper presents a high order least square-based finite difference-finite volume (LSFD-FV) method together with lattice Boltzmann flux solver for accurate simulation of incompressible flows on unstructured grids. Within each control cell, a high order polynomial, which is based on Taylor series expansion, is applied to approximate the solution function. The derivatives in the Taylor series expansion are approximated by the functional values at the centers of neighboring cells using the mesh-free least square-based finite difference (LSFD) method developed by Ding et al. (2004) [36] . In the high order finite volume method, the recently developed lattice Boltzmann flux solver (LBFS) is applied to evaluate the inviscid and viscous fluxes physically and simultaneously at the cell interface by local reconstruction of lattice Boltzmann solution. Compared with traditional k-exact high-order finite volume method, the present method is more accurate and efficient. Various benchmark examples are tested to validate the high-order accuracy, high computational efficiency and flexibility of the proposed method on unstructured grids.

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

[2]  Chris Lacor,et al.  On the connection between the spectral volume and the spectral difference method , 2007, J. Comput. Phys..

[3]  C. Ollivier-Gooch.  Quasi-ENO Schemes for Unstructured Meshes Based on Unlimited Data-Dependent Least-Squares Reconstruction , 1997 .

[4]  Fangying Song,et al.  An Efficient Spectral Method for the Inextensible Immersed Interface in Incompressible Flows , 2019, Communications in Computational Physics.

[5]  Rémi Abgrall,et al.  On essentially non-oscillatory schemes on unstructured meshes: analysis and implementation , 1994 .

[6]  Chaowei Hu,et al.  No . 98-32 Weighted Essentially Non-Oscillatory Schemes on Triangular Meshes , 1998 .

[7]  Chaoqun Liu,et al.  Preconditioned Multigrid Methods for Unsteady Incompressible Flows , 1997 .

[8]  S. Rebay,et al.  High-Order Accurate Discontinuous Finite Element Solution of the 2D Euler Equations , 1997 .

[9]  Chiang Juay Teo,et al.  Thermal lattice Boltzmann flux solver and its application for simulation of incompressible thermal flows , 2014 .

[10]  W. X. Wu,et al.  Free vibration analysis of plates using least-square-based finite difference method , 2007 .

[11]  Fermín Navarrina,et al.  On the accuracy of finite volume and discontinuous Galerkin discretizations for compressible flow on unstructured grids , 2009 .

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

[13]  Zhi J. Wang,et al.  Spectral (finite) volume method for conservation laws on unstructured grids IV: extension to two-dimensional systems , 2004 .

[14]  Chris Lacor,et al.  An accuracy and stability study of the 2D spectral volume method , 2007, J. Comput. Phys..

[15]  Zhi J. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids. Basic Formulation , 2002 .

[16]  Y. Wang,et al.  Explicit formulations of gas-kinetic flux solver for simulation of incompressible and compressible viscous flows , 2015, J. Comput. Phys..

[17]  C. Shu,et al.  Development of least-square-based two-dimensional finite-difference schemes and their application to simulate natural convection in a cavity , 2004 .

[18]  Mahidhar Tatineni,et al.  Very high-order compact finite difference schemes on non-uniform grids for incompressible Navier-Stokes equations , 2007, J. Comput. Phys..

[19]  Fermín Navarrina,et al.  Finite volume solvers and movingleast-squares approximations for thecompressible Navier-Stokes equations onunstructured grids , 2007 .

[20]  Yuzhi Sun,et al.  Spectral (finite) volume method for conservation laws on unstructured grids VI: Extension to viscous flow , 2006, J. Comput. Phys..

[21]  Chang Shu,et al.  Numerical simulation of flows around two circular cylinders by mesh‐free least square‐based finite difference methods , 2007 .

[22]  C. Ollivier-Gooch,et al.  A high-order-accurate unstructured mesh finite-volume scheme for the advection-diffusion equation , 2002 .

[23]  Antony Jameson,et al.  Spectral Difference Method for Unstructured Grids II: Extension to the Euler Equations , 2007, J. Sci. Comput..

[24]  Hong Luo,et al.  Reconstructed Discontinuous Galerkin Methods for Hyperbolic Diffusion Equations on Unstructured Grids , 2019, Communications in Computational Physics.

[25]  Qian Wang,et al.  Compact high order finite volume method on unstructured grids III: Variational reconstruction , 2017, J. Comput. Phys..

[26]  Bernardo Cockburn.  Discontinuous Galerkin methods , 2003 .

[27]  Chiang Juay Teo,et al.  Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows , 2014 .

[28]  Chang Shu,et al.  Mesh-free least-squares-based finite difference method for large-amplitude free vibration analysis of arbitrarily shaped thin plates , 2008 .

[29]  O. Friedrich,et al.  Weighted Essentially Non-Oscillatory Schemes for the Interpolation of Mean Values on Unstructured Grids , 1998 .

[30]  C. Shu,et al.  Lattice Boltzmann Method and Its Applications in Engineering , 2013 .

[31]  J. Douglas,et al.  Interior Penalty Procedures for Elliptic and Parabolic Galerkin Methods , 1976 .

[32]  Ignasi Colominas,et al.  High-order Finite Volume Methods and Multiresolution Reproducing Kernels , 2008 .

[33]  A. Jameson.  Time dependent calculations using multigrid, with applications to unsteady flows past airfoils and wings , 1991 .

[34]  M. Braza,et al.  Numerical study and physical analysis of the pressure and velocity fields in the near wake of a circular cylinder , 1986, Journal of Fluid Mechanics.

[35]  Weiwei Zhang,et al.  A high-order finite volume method on unstructured grids using RBF reconstruction , 2016, Comput. Math. Appl..

[36]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

[37]  Zhi Jian Wang,et al.  A unifying lifting collocation penalty formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids , 2009, J. Comput. Phys..

[38]  Weiwei Zhang,et al.  A Novel Multi-Dimensional Limiter for High-Order Finite Volume Methods on Unstructured Grids , 2017 .

[39]  Yong Wang,et al.  A simplified circular function–based gas kinetic scheme for simulation of incompressible flows , 2017 .

[40]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[41]  P. Frederickson,et al.  Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction , 1990 .

[42]  Qian Wang,et al.  Compact high order finite volume method on unstructured grids II: Extension to two-dimensional Euler equations , 2016, J. Comput. Phys..

[43]  Xiaoyi He,et al.  Lattice Boltzmann Method on Curvilinear Coordinates System , 1997 .

[44]  S. Dennis,et al.  Numerical solutions for steady flow past a circular cylinder at Reynolds numbers up to 100 , 1970, Journal of Fluid Mechanics.

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

[46]  Ilaria Perugia,et al.  Superconvergence of the Local Discontinuous Galerkin Method for Elliptic Problems on Cartesian Grids , 2001, SIAM J. Numer. Anal..

[47]  R. Hartmann.  Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element methods , 2008 .

[48]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[49]  Michael Dumbser,et al.  Efficient implementation of high order unstructured WENO schemes for cavitating flows , 2013 .

[50]  C. Shu,et al.  Simulation of incompressible viscous flows past a circular cylinder by hybrid FD scheme and meshless least square-based finite difference method , 2004 .

[51]  Zhi Jian Wang,et al.  A conservative correction procedure via reconstruction formulation with the Chain-Rule divergence evaluation , 2013, J. Comput. Phys..

[52]  J. Wu,et al.  Implicit velocity correction-based immersed boundary-lattice Boltzmann method and its applications , 2009 .

[53]  Ghader Ghorbaniasl,et al.  A stability analysis for the spectral volume method on tetrahedral grids , 2009, J. Comput. Phys..

[54]  Fan Zhang,et al.  A Simplified Artificial Compressibility Flux for the Discontinuous Galerkin Solution of the Incompressible Navier-Stokes Equations , 2019, Communications in Computational Physics.

[55]  Timothy J. Barth,et al.  Recent developments in high order K-exact reconstruction on unstructured meshes , 1993 .

[56]  Yu Wang,et al.  An improved multiphase lattice Boltzmann flux solver for three-dimensional flows with large density ratio and high Reynolds number , 2015, J. Comput. Phys..

[57]  Taohua Liu,et al.  A Consistent Fourth-Order Compact Finite Difference Scheme for Solving Vorticity-Stream Function Form of Incompressible Navier-Stokes Equations , 2019, Numerical Mathematics: Theory, Methods and Applications.

[58]  Marcelo Reggio,et al.  Solving incompressible fluid flows on unstructured meshes with the lattice Boltzmann flux solver , 2017 .

[59]  Chiang Juay Teo,et al.  Development of LBGK and incompressible LBGK‐based lattice Boltzmann flux solvers for simulation of incompressible flows , 2014 .