A higher order compact finite difference algorithm for solving the incompressible Navier–Stokes equations

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 (C) 2011 John Wiley & Sons, Ltd.

[1]  M. Ben-Artzi,et al.  A pure-compact scheme for the streamfunction formulation of Navier-Stokes equations , 2005 .

[2]  Jie Shen,et al.  An Efficient Spectral-Projection Method for the Navier-Stokes Equations in Cylindrical Geometries , 2002 .

[3]  Nicola Parolini,et al.  Numerical investigation on the stability of singular driven cavity flow , 2002 .

[4]  Alexandre J. Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations , 1969 .

[5]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[6]  B. Fornberg,et al.  A compact fourth‐order finite difference scheme for the steady incompressible Navier‐Stokes equations , 1995 .

[7]  Shiyi Chen,et al.  Simulation of Cavity Flow by the Lattice Boltzmann Method , 1994, comp-gas/9401003.

[8]  M. Minion,et al.  Accurate projection methods for the incompressible Navier—Stokes equations , 2001 .

[9]  Jian‐Guo Liu,et al.  Projection method I: convergence and numerical boundary layers , 1995 .

[10]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[11]  F. DUCROSb,et al.  COMPACT FINITE DIFFERENCE SCHEMES ON NON-UNIFORM MESHES . APPLICATION TO DIRECT NUMERICAL SIMULATIONS OF COMPRESSIBLE FLOWS , 1998 .

[12]  R. Knikker,et al.  Study of a staggered fourth‐order compact scheme for unsteady incompressible viscous flows , 2009 .

[13]  Nicola Parolini,et al.  A mixed-basis spectral projection method , 2002 .

[14]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[15]  Jean-Pierre Croisille,et al.  A Compact Difference Scheme for the Biharmonic Equation in Planar Irregular Domains , 2009, SIAM J. Numer. Anal..

[16]  Stéphane Abide,et al.  A 2D compact fourth-order projection decomposition method , 2005 .

[17]  Tao Tang,et al.  A Compact Fourth-Order Finite Difference Scheme for Unsteady Viscous Incompressible Flows , 2001, J. Sci. Comput..

[18]  Fourth‐ and tenth‐order compact finite difference solutions of perturbed circular vortex flows , 2005 .

[19]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[20]  Michel Visonneau,et al.  Incompressible flow calculations with a consistent physical interpolation finite volume approach , 1994 .

[21]  Olivier Botella,et al.  On the solution of the Navier-Stokes equations using Chebyshev projection schemes with third-order accuracy in time , 1997 .

[22]  E Weinan,et al.  Projection Method II: Godunov--Ryabenki Analysis , 1996 .

[23]  Alexandre Joel Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations* , 1989 .

[24]  Yongbin Ge,et al.  A fourth‐order compact finite difference scheme for the steady stream function–vorticity formulation of the Navier–Stokes/Boussinesq equations , 2003 .

[25]  Toshio Kobayashi,et al.  Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme , 1999 .

[26]  A. Chorin.  Numerical solution of the Navier-Stokes equations , 1968 .

[27]  G. Carey,et al.  High‐order compact scheme for the steady stream‐function vorticity equations , 1995 .

[28]  Charles-Henri Bruneau,et al.  An efficient scheme for solving steady incompressible Navier-Stokes equations , 1990 .

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

[30]  M. Sahin,et al.  A novel fully implicit finite volume method applied to the lid‐driven cavity problem—Part I: High Reynolds number flow calculations , 2003 .

[31]  MA YANWENa,et al.  NUMERICAL SOLUTION OF THE INCOMPRESSIBLE NAVIER–STOKES EQUATIONS WITH AN UPWIND COMPACT DIFFERENCE SCHEME , 1999 .

[32]  Murli M. Gupta,et al.  A new paradigm for solving Navier-Stokes equations: streamfunction-velocity formulation , 2005 .

[33]  Y. V. S. S. Sanyasiraju,et al.  Higher order semi compact scheme to solve transient incompressible Navier-Stokes equations , 2005 .

[34]  H. B. Keller,et al.  Driven cavity flows by efficient numerical techniques , 1983 .