This paper describes a novel vortex element method for simulating incompressible laminar flow over a two-dimensional backward-facing step. The model employs an operator-splitting technique to compute the evolution of the vorticity field downstream of abrupt changes in flow geometry. During the advective stage of the computation, a semi-Lagrangian scheme is used to update the positions of the vortex elements, whilst an analytical diffusion algorithm employing Oseen vortices is implemented during the diffusive time step. Redistributing the vorticity analytically instead of using the more traditional random-walk method enables the numerical model to simulate steady flows directly and avoids the need to filter the results to remove the oscillations created by the random-walk procedure.

Model validation has been achieved by comparing the length of the recirculating eddy behind a confined backward-facing step against data from experimental and alternative numerical investigations. In addition, results from the vortex element method are compared against predictions obtained using the commercial finite-volume computational fluid dynamics code, CFD-ACE+. The results show that the vortex element scheme marginally overpredicts the length of the downstream recirculating eddy, implying that the method may be associated with an artificial reduction in the vorticity diffusion rate. Nevertheless the results demonstrate that the proposed vortex redistribution scheme provides a practical alternative to traditional random-walk discrete vortex algorithms.