Some flows such as the wakes of rotating devices often display helical symmetry. We present an original DNS code for the dynamics of such helically symmetric systems. We show that, by enforcing helical symmetry, the three-dimensional Navier–Stokes equations can be reduced to a two-dimensional unsteady problem. The numerical method is a generalisation of the vorticity/streamfunction formulation in a circular domain, with finite differences in the radial direction and spectral decomposition along the azimuth. When compared to a standard three-dimensional code, this allows us to reach larger Reynolds numbers and to compute quasi-steady patterns. We illustrate the importance of helical pitch by some physical cases: the dynamics of several helical vortices and a quasi-steady vortex flow. We also study the linear dynamics and nonlinear saturation in the Batchelor vortex basic flow, a paradigmatic example of trailing vortex instability. We retrieve the behaviour of inviscid modes and present new results concerning the saturation of viscous centre modes.