In this work, a high-order weighted essentially nonoscillatory (WENO) finite-difference lattice Boltzmann method (WENOLBM) is developed and assessed for an accurate simulation of incompressible flows. To handle curved geometries with nonuniform grids, the incompressible form of the discrete Boltzmann equation with the Bhatnagar-Gross-Krook (BGK) approximation is transformed into the generalized curvilinear coordinates and the spatial derivatives of the resulting lattice Boltzmann equation in the computational plane are solved using the fifth-order WENO scheme. The first-order implicit-explicit Runge-Kutta scheme and also the fourth-order Runge-Kutta explicit time integrating scheme are adopted for the discretization of the temporal term. To examine the accuracy and performance of the present solution procedure based on the WENOLBM developed, different benchmark test cases are simulated as follows: unsteady Taylor-Green vortex, unsteady doubly periodic shear layer flow, steady flow in a two-dimensional (2D) cavity, steady cylindrical Couette flow, steady flow over a 2D circular cylinder, and steady and unsteady flows over a NACA0012 hydrofoil at different flow conditions. Results of the present solution are compared with the existing numerical and experimental results which show good agreement. To show the efficiency and accuracy of the solution methodology, the results are also compared with the developed second-order central-difference finite-volume lattice Boltzmann method and the compact finite-difference lattice Boltzmann method. It is shown that the present numerical scheme is robust, efficient, and accurate for solving steady and unsteady incompressible flows even at high Reynolds number flows.
SUMMARYThis paper uses a fourth-order compact finite-difference scheme for solving steady incompressible flows. The high-order compact method applied is an alternating direction implicit operator scheme, which has been used by Ekaterinaris for computing two-dimensional compressible flows. Herein, this numerical scheme is efficiently implemented to solve the incompressible Navier-Stokes equations in the primitive variables formulation using the artificial compressibility method. For space discretizing the convective fluxes, fourth-order centered spatial accuracy of the implicit operators is efficiently obtained by performing compact space differentiation in which the method uses block-tridiagonal matrix inversions. To stabilize the numerical solution, numerical dissipation terms and/or filters are used. In this study, the highorder compact implicit operator scheme is also extended for computing three-dimensional incompressible flows. The accuracy and efficiency of this high-order compact method are demonstrated for different incompressible flow problems. A sensitivity study is also conducted to evaluate the effects of grid resolution and pseudocompressibility parameter on accuracy and convergence rate of the solution. The effects of filtering and numerical dissipation on the solution are also investigated. Test cases considered herein for validating the results are incompressible flows in a 2-D backward facing step, a 2-D cavity and a 3-D cavity at different flow conditions. Results obtained for these cases are in good agreement with the available numerical and experimental results. The study shows that the scheme is robust, efficient and accurate for solving incompressible flow problems.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.