A novel vorticity-velocity formulation of the Navier-Stokes equations -the MassConserving, Smooth (MC-Smooth) vorticity-velocity formulation -is developed in this work. The governing equations of the MC-Smooth formulation include a new secondorder Poisson-like elliptic velocity equation, along with the vorticity transport equation, the energy conservation equation, and N spec species mass balance equations. In this study, the MC-Smooth formulation is compared to two pre-existing vorticity-velocity formulations by applying each formulation to confined and unconfined axisymmetric laminar diffusion flame problems. For both applications, very good to excellent agreement for the simulation results of the three formulations has been obtained. The MC-Smooth formulation requires the least CPU time and can overcome the limitations of the other two pre-existing vorticity-velocity formulations by ensuring mass conservation and solution smoothness over a broader range of flow conditions. In addition to these benefits, other important features of the MC-Smooth formulation include: (1) it does not require the use of a staggered grid, and (2) it does not require excessive grid refinement to ensure mass conservation. The MC-Smooth formulation is a computationally attractive approach that can effectively extend the applicability of the vorticity-velocity formulation.