The need to numerically represent a free vortex system arises frequently in fundamental and applied research. Many possible techniques for realizing this vortex system exist but most tend to prioritize accuracy either inside or outside of the vortex core, which therefore makes them unsuitable for a stability analysis considering the entire flow field. In this article, a simple method is presented that is shown to yield an accurate representation of the flow inside and outside of the vortex core. The method is readily implemented in any incompressible Navier–Stokes solver using primitive variables and Cartesian coordinates. It can potentially be used to model a wide range of vortices but is here applied to the case of two helices, which is of renewed interest due to its relevance for wind turbines and helicopters. Three-dimensional stability analysis is performed in both a rotating and a translating frame of reference, which yield eigenvalue spectra that feature both mutual inductance and elliptic instabilities. Comparison of these spectra with available theoretical predictions is used to validate the proposed baseflow model, and new insights into the elliptic instability of curved Batchelor vortices are presented. Furthermore, it is shown that the instabilities in the rotating and the translating reference frames have the same structure and growth rate, but different frequency. A relation between these frequencies is provided.