A compact monotonic discretization scheme for solving second-order vorticity–velocity equations Original Research Article

This paper presents a numerical method for solving the steady-state Navier–Stokes equations for incompressible fluid flows using velocities and vorticity as working variables. The method involves solving a second-order differential equation for the velocity and a convection–diffusion equation for the vorticity in Cartesian grids. The key to the success of the numerical simulation of this class of flow equations depends largely on proper simulation of vorticity transport equation subject to proper boundary vorticity. In this paper, we present a monotonic advection–diffusion multi-dimensional scheme and a theoretically rigorous implementation of vorticity boundary conditions. While the derivation of the proposed integral vorticity boundary condition is more elaborate and is more difficult to solve than conventional local approaches, the present approach offers significant advantages. In this study, both lid-driven and backward-facing step problems have been selected for comparison and validation purposes.