Sullivan vortex

In fluid dynamics, the Sullivan vortex is an exact solution of the Navier–Stokes equations describing a two-celled vortex in an axially strained flow, that was discovered by Roger D. Sullivan in 1959.

[1][2] At large radial distances, the Sullivan vortex resembles a Burgers vortex, however, it exhibits a two-cell structure near the center, creating a downdraft at the axis and an updraft at a finite radial location.

[3] Specifically, in the outer cell, the fluid spirals inward and upward and in the inner cell, the fluid spirals down at the axis and spirals upwards at the boundary with the outer cell.

[4] Due to its multi-celled structure, the vortex is used to model tornadoes[5] and large-scale complex vortex structures in turbulent flows.

[6] Consider the velocity components

of an incompressible fluid in cylindrical coordinates in the form[7] where

η = α

is the strain rate of the axisymmetric stagnation-point flow.

The Burgers vortex solution is simply given by

Sullivan showed that there exists a non-trivial solution for

from the Navier-Stokes equations accompanied by a function

is the exponential integral.

being is the Euler–Mascheroni constant, whereas for large values of

, which is obtained by solving the equation

Within the inner cell, the transition between the downdraft and the updraft occurs at

, which is obtained by solving the equation

The vorticity components of the Sullivan vortex are given by The pressure field

with respect to its central value

is the fluid density.

The first term on the right-hand side corresponds to the potential flow motion, i.e.,

, whereas the remaining two terms originates from the motion associated with the Sullivan vortex.

Explicit solution of the Navier–Stokes equations for the Sullivan vortex in stretched cylindrical stagnation surfaces was solved by P. Rajamanickam and A. D. Weiss and is given by[8] where

η = α

, Note that the location of the stagnation cylindrical surface is not longer given by

is the principal branch of the Lambert W function.

here should be interpreted as the measure of the volumetric source strength

= 2 π α

and not the location of the stagnation surface.

Here, the vorticity components of the Sullivan vortex are given by