Premium
A higher order compact finite difference algorithm for solving the incompressible Navier–Stokes equations
Author(s) -
Tian Zhenfu,
Liang Xian,
Yu Peixiang
Publication year - 2011
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.3184
Subject(s) - stencil , compact finite difference , mathematics , discretization , finite difference , poisson's equation , navier–stokes equations , mathematical analysis , finite difference method , order of accuracy , pressure correction method , central differencing scheme , compressibility , projection method , temporal discretization , finite difference coefficient , numerical analysis , finite element method , mathematical optimization , numerical stability , physics , mechanics , dykstra's projection algorithm , computational science , mixed finite element method , thermodynamics
On the basis of the projection method, a higher order compact finite difference algorithm, which possesses a good spatial behavior, is developed for solving the 2D unsteady incompressible Navier–Stokes equations in primitive variable. The present method is established on a staggered grid system and is at least third‐order accurate in space. A third‐order accurate upwind compact difference approximation is used to discretize the non‐linear convective terms, a fourth‐order symmetrical compact difference approximation is used to discretize the viscous terms, and a fourth‐order compact difference approximation on a cell‐centered mesh is used to discretize the first derivatives in the continuity equation. The pressure Poisson equation is approximated using a fourth‐order compact difference scheme constructed currently on the nine‐point 2D stencil. New fourth‐order compact difference schemes for explicit computing of the pressure gradient are also developed on the nine‐point 2D stencil. For the assessment of the effectiveness and accuracy of the method, particularly its spatial behavior, a problem with analytical solution and another one with a steep gradient are numerically solved. Finally, steady and unsteady solutions for the lid‐driven cavity flow are also used to assess the efficiency of this algorithm. Copyright © 2011 John Wiley & Sons, Ltd.