Fourth‐order finite difference scheme for the incompressible Navier–Stokes equations in a disk

We develop an efficient fourth-order finite difference method for solving the incompressible Navier–Stokes equations in the vorticity-stream function formulation on a disk. We use the fourth-order Runge–Kutta method for the time integration and treat both the convection and diffusion terms explicitly. Using a uniform grid with shifting a half mesh away from the origin, we avoid placing the grid point directly at the origin; thus, no pole approximation is needed. Besides, on such grid, a fourth-order fast direct method is used to solve the Poisson equation of the stream function. By Fourier filtering the vorticity in the azimuthal direction at each time stage, we are able to increase the time step to a reasonable size. The numerical results of the accuracy test and the simulation of a vortex dipole colliding with circular wall are presented. Copyright © 2003 John Wiley & Sons, Ltd.

[1]  John D. Anderson,et al.  A computational fluid dynamic technique valid at the centerline for non-axisymmetric problems in cylindrical coordinates , 1979 .

[2]  T. Colonius,et al.  Numerical Treatment of Polar Coordinate Singularities , 2000 .

[3]  Jian‐Guo Liu,et al.  Vorticity Boundary Condition and Related Issues for Finite Difference Schemes , 1996 .

[4]  Ming-Chih Lai,et al.  Fast direct solvers for Poisson equation on 2D polar and spherical geometries , 2002 .

[5]  Parviz Moin,et al.  An Efficient Method for Temporal Integration of the Navier-Stokes Equations in Confined Axisymmetric Geometries , 1996 .

[6]  W. Roger Briley,et al.  A numerical study of laminar separation bubbles using the Navier-Stokes equations , 1971, Journal of Fluid Mechanics.

[7]  Bengt Fornberg,et al.  Comparison of finite difference‐ and pseudospectral methods for convective flow over a sphere , 1997 .

[8]  Evangelos A. Coutsias,et al.  Pseudospectral Solution of the Two-Dimensional Navier-Stokes Equations in a Disk , 1999, SIAM J. Sci. Comput..

[9]  Linear Stability of Finite-Difference Approximations on a Uniform Latitude-Longitude Grid with Fourier Filtering , 1976 .

[10]  A. Thom,et al.  The flow past circular cylinders at low speeds , 1933 .

[11]  M. Lai.  A note on finite difference discretizations for Poisson equation on a disk , 2001 .

[12]  Andrea Vacca,et al.  An Efficient Method for the Solution of the Incompressible Navier-Stokes Equations in Cylindrical Geometries , 1999 .

[13]  Arieh Iserles.  A First Course in the Numerical Analysis of Differential Equations: Multigrid techniques , 1996 .

[14]  C. Canuto.  Spectral methods in fluid dynamics , 1991 .

[15]  L. Quartapelle,et al.  Numerical solution of the incompressible Navier-Stokes equations , 1993, International series of numerical mathematics.