Inclusion of no-slip boundary conditions in the MEEVC scheme

This work presents three methods for enforcing tangential velocity boundary conditions for the MEEVC scheme, which was shown to be mass, enstrophy, energy and vorticity conserving scheme in the case of inviscid flow [1]. While the normal velocity component can be strongly imposed in a div-conforming formulation for the velocity field, inclusion of the tangential velocity needs to be set through an appropriate choice of vorticity boundary conditions. Three methods to impose the tangential velocity boundary condition will be discussed: The kinematic Dirichlet formulation, the kinematic Neumann formulation and the dynamic Neumann formulation. The conservation properties of each of the resulting schemes are analyzed and numerical results are shown for the Taylor–Green vortex and for the dipole collision test cases. These confirm that kinematic Neumann vorticity boundary conditions perform best.