This paper is aimed at eliciting consistency conditions for the existence of unsteady incompressible axisymmetric swirling viscous Beltrami vortices and explicitly constructing solutions that obey the conditions as well as the Navier–Stokes equations. By Beltrami flow, it is meant that vorticity, i.e., the curl of velocity, is proportional to velocity at any local point in space and time. The consistency conditions are derived for the proportionality coefficient, the velocity field and external force. The coefficient, whose dimension is of [length−1], is either constant or nonconstant. In the former case, the well-known exact nondivergent three-dimensional unsteady vortex solutions are obtained by solving the evolution equations for the stream function directly. In the latter case, the consistency conditions are given by nonlinear equations of the stream function, one of which corresponds to the Bragg–Hawthorne equation for steady inviscid flow. Solutions of a novel type are found by numerically solving the nonlinear constraint equation at a fixed time. Time dependence is recovered by taking advantage of the linearity of the evolution equation of the stream function. The proportionality coefficient is found to depend on space and time. A phenomenon of partial restoration of the broken scaling invariance is observed at short distances.