zbMATH — the first resource for mathematics

Fourth order convergence of compact finite difference solver for 2D incompressible flow. (English) Zbl 1085.65515
Summary: We study a fourth order finite difference method for the unsteady incompressible Navier-Stokes equations in vorticity formulation. The scheme is essentially compact and can be implemented very efficiently. Either Briley’s formula, or a new higher order formula, which will be derived in this paper, can be chosen as the vorticity boundary condition. By formal Taylor expansion, the new formula for the vorticity on the boundary gives 4th order accuracy; while Briley’s formula provides only 3rd order accuracy. However, the use of either formula results in a stable method and achieves full 4th order accuracy. The convergence analysis of the scheme with our new formula will be given in this paper, while that with Briley’s formula has been established in earlier literature. The consistency analysis is easier than that of Briley’s formula, no Strang type analysis is needed. In the stability analysis part, we adopt the technique of controlling some local terms by the diffusion term via discrete elliptic regularity. Physical no-slip boundary conditions are used throughout.

65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
35Q30 Navier-Stokes equations
76D05 Navier-Stokes equations for incompressible viscous fluids
76M20 Finite difference methods applied to problems in fluid mechanics