Boundary element-finite element method for velocity-vorticity formulation of Navier-Stokes equations

A numerical method for the solution of the Navier-Stokes equations is developed using an integral representation of the conservation equations. The velocity-vorticity formulation is employed, where the kinematics is given with the Poisson equation for a velocity vector, while the kinetics is represented with the vorticity transport equation. The corresponding boundary-domain integral equations are presented along with discussions of the kinematics and kinetics of the fluid flow problem. Kinematics is solved using the boundary element method (BEM), while kinetics is solved using the finite element method (FEM). Quadratic continuous interpolation functions are used for both BEM and FEM. Two benchmark problems are considered to show the robustness and versatility of this formulation including lid driven flow in a square cavity and flow over a backward facing step.