Vortex Dynamics & Multi-Field Applications

• Vortex is a rotating or phase-winding structure defined by local dynamics and geometry in fluids, plasmas, and quantum condensates.
• Researchers use varied criteria—including velocity gradients and eigenanalysis—to isolate vortex structures and quantify their strengths.
• Studies show vortices play crucial roles in energy transport, magnetic interactions, and wave propagation across physical and astrophysical systems.

A vortex is a rotating or phase-winding structure whose precise definition depends on the medium and observable under study. In fluid dynamics, one rigorous proposal defines a vortex as a connected region in which the magnitude of the vortex vector is larger than zero; in solar MHD, vortices are diagnosed through kinetic and magnetic swirling strengths derived from gradient-tensor eigenanalysis; in exciton-polariton condensates, a vortex carries a topological charge $M=\pm1$; and in structured light, a vortex beam is generated by a helical phase factor $\exp(i\ell\varphi)$ (Tian et al., 2017, Kannan et al., 2024, Ma et al., 2017, Rodríguez-Fajardo et al., 2024). The literature also emphasizes that there is not a consensus on the vortex definition in fluid dynamics, and that raw vorticity alone does not always distinguish rotation from shear (Tian et al., 2017).

1. Definitions, kinematics, and detection criteria

In the fluid-dynamical formulation of Tian et al., the local rotation axis is the unique unit vector $r$ for which, in the rotated frame aligned with $r$, the two velocity-gradient components normal to the axis vanish, $\partial U/\partial Z=0$ and $\partial V/\partial Z=0$. In the plane perpendicular to $r$, a point is fluid-rotational if $\alpha-\beta<0$, and the rotational strength is $\xi=2\min_\theta \Omega(\theta)$. The vortex vector is then $\omega_v\equiv R\equiv \xi r$, the total vorticity is decomposed as $\exp(i\ell\varphi)$0, and a vortex is a connected region in which $\exp(i\ell\varphi)$1 (Tian et al., 2017). This construction is explicitly local, Galilean invariant, and unique.

A different operational definition is used in solar MHD. With $\exp(i\ell\varphi)$2 the plasma velocity and $\exp(i\ell\varphi)$3 the magnetic field, one introduces the gradient tensors $\exp(i\ell\varphi)$4 and $\exp(i\ell\varphi)$5. If the eigenvalues of $\exp(i\ell\varphi)$6 are $\exp(i\ell\varphi)$7 with $\exp(i\ell\varphi)$8, then the kinetic swirling strength is $\exp(i\ell\varphi)$9 and the local rotation period is $r$0; the magnetic analogue is $r$1 (Kannan et al., 2024). In solar vortex-tube extraction, another criterion is the Instantaneous Vorticity Deviation,

$r$2

with the vortex boundary taken as the outermost closed contour around an IVD maximum subject to a convexity-deficiency threshold $r$3 (Silva et al., 2020).

These parallel definitions show that vortex identification is criterion-dependent. A plausible implication is that no single scalar diagnostic captures all aspects of vortex structure across fluids, MHD, condensates, and wave systems: some criteria isolate intrinsic rotation, some separate magnetic from kinetic twist, and some are designed to extract geometrically coherent boundaries.

2. Fluid and magnetohydrodynamic vortex dynamics

Realistic three-dimensional MURaM simulations of the Quiet Sun, Weak Plage, and Strong Plage show that similar photospheric turbulent driving does not imply similar chromospheric vortex populations. At $r$4 Mm, Quiet Sun vortices have typical diameters $r$5 Mm with area coverage $r$6, Weak Plage has $r$7 Mm and area $r$8, and Strong Plage has $r$9 Mm and area $r$0. At $r$1, the velocity-power spectra have identical inertial-range slopes $r$2 for all three cases, while magnetic-power spectra scale as $r$3 at all $r$4. Above $r$5 Mm, the observed and linear-theory magnetic swirling strengths approach each other, but at $r$6 Mm the ratio $r$7 is $r$8, $r$9, and $\partial U/\partial Z=0$0. Strong Plage also exhibits a steeper chromospheric magnetic cascade, $\partial U/\partial Z=0$1 with $\partial U/\partial Z=0$2 for $\partial U/\partial Z=0$3, compared with $\partial U/\partial Z=0$4 in Quiet Sun and Weak Plage (Kannan et al., 2024). The interpretation proposed there is a dual regime: in weaker magnetic regions, vortex-driven torsional Alfvén waves can propagate coherently upward with little damping, whereas in strongly magnetized plage, vortex interactions in narrow flux tubes promote dissipation and heating.

Three-dimensional solar vortex tubes identified by IVD are approximately conical and widen with height, with mean radius $\partial U/\partial Z=0$5 km at the photosphere and $\partial U/\partial Z=0$6 km in the low chromosphere. Their tangential velocity $\partial U/\partial Z=0$7 decreases from the tube boundary toward the center at all sampled heights, which was interpreted as an eddy-viscosity effect braking the swirl. The magnetic field is $\partial U/\partial Z=0$8–$\partial U/\partial Z=0$9 stronger at the center than at the boundary, consistent with flux concentration in downdrafts. However, only one of the 17 tracked kinematic tubes achieved tangential speeds high enough that $\partial V/\partial Z=0$0 locally exceeded the tension scale $\partial V/\partial Z=0$1, thereby producing a cospatial magnetic vortex (Silva et al., 2020). This directly addresses a common misconception: kinematic vortices do not, in general, imply magnetic vortices.

In patient-specific direct numerical simulation of the left ventricle, the mitral vortex ring formed during early diastolic filling at $\partial V/\partial Z=0$2–$\partial V/\partial Z=0$3 and $\partial V/\partial Z=0$4–$\partial V/\partial Z=0$5 does not remain intact. As it propagates toward the apex, it becomes inclined and deformed by lateral straining; trailing secondary vortex tubes are shed from interactions with the ventricular wall; and helicity polarization and twisting instabilities develop in the interacting cores. By end-diastole, the ring impinges on the wall, the primary-ring circulation decays to $\partial V/\partial Z=0$6, and the global enstrophy increases by $\partial V/\partial Z=0$7 (Le et al., 2011). This demonstrates that vortex-wall and vortex-vortex interactions are central even in confined physiological flows.

3. Quantized vortices in condensates, superfluids, and magnetic matter

In incoherently driven exciton-polariton condensates, a ring-shaped off-resonant pump

$\partial V/\partial Z=0$8

with typical parameters $\partial V/\partial Z=0$9 and $r$0 nucleates crater-type vortices of topological charge $r$1 without imposing a phase profile. When two rings are separated by center-to-center distance $r$2, their overlap enforces either $r$3-locking or $r$4-locking: in practice $r$5 yields a $r$6-state and $r$7 a $r$8-state. This phase-locking enables controlled inversion of a vortex charge by temporarily switching off the target pump, as well as copying of a reference charge in triangular $r$9-locked geometries. The resulting vortex-formation time is of order $\alpha-\beta<0$0–$\alpha-\beta<0$1 ps and the switching time of order $\alpha-\beta<0$2–$\alpha-\beta<0$3 ps (Ma et al., 2017).

A related driven-dissipative polariton setting uses a homogeneous pump with an intensity groove. A circular groove of lower intensity produces dark-ring condensate states with a $\alpha-\beta<0$4-phase jump across the ring and remains stable, for example, up to $\alpha-\beta<0$5 for $\alpha-\beta<0$6, $\alpha-\beta<0$7, and $\alpha-\beta<0$8. Above that threshold, the dark ring undergoes snake instability and breaks into symmetrically arranged vortex-antivortex pairs; multiple pair-number states are stable for the same ring, i.e. true multistability. Wider grooves support higher-order dark states with multiple $\alpha-\beta<0$9-phase jumps, which act as vortex waveguides. In a circular guide, an imprinted $\xi=2\min_\theta \Omega(\theta)$0 vortex travels clockwise and $\xi=2\min_\theta \Omega(\theta)$1 travels counter-clockwise; in a U-shaped guide, each charge has only one allowed corner turn, with representative transit times of $\xi=2\min_\theta \Omega(\theta)$2 ps and $\xi=2\min_\theta \Omega(\theta)$3 ps for opposite charges (2206.12157).

In holographic superfluid, vortex-antivortex annihilation proceeds in two stages separated by the vortex core diameter $\xi=2\min_\theta \Omega(\theta)$4. The separation obeys

$\xi=2\min_\theta \Omega(\theta)$5

and

$\xi=2\min_\theta \Omega(\theta)$6

which implies effective attractive-force scalings $\xi=2\min_\theta \Omega(\theta)$7 and $\xi=2\min_\theta \Omega(\theta)$8, respectively (Lan et al., 2018). In multi-component Bose-Einstein condensates, the pairwise interaction picture is sufficient to explain dimers, trimers, and exotic lattices: intra-component interaction is always monotonic and repulsive, while inter-component interaction can be attractive or nonmonotonic depending on density-density and Rabi couplings; for trimers, the equilibrium separation obeys $\xi=2\min_\theta \Omega(\theta)$9 (Dantas et al., 2015). In ferromagnetic superconductors, magnetic-moment polarization adds a long-range attractive contribution $\omega_v\equiv R\equiv \xi r$0, stabilizing finite-size vortex clusters at low density; under current-driven motion, clusters can resonate with magnons, and above a threshold velocity $\omega_v\equiv R\equiv \xi r$1 domain walls are nucleated and the vortex configuration becomes modulated (Lin et al., 2012).

Micromagnetic vortices in confined ferromagnets exhibit another regime structure. In Permalloy dots with $\omega_v\equiv R\equiv \xi r$2 nm and $\omega_v\equiv R\equiv \xi r$3 nm, the vortex stable regime is $\omega_v\equiv R\equiv \xi r$4, the metastable regime is $\omega_v\equiv R\equiv \xi r$5, with $\omega_v\equiv R\equiv \xi r$6 Oe and $\omega_v\equiv R\equiv \xi r$7–$\omega_v\equiv R\equiv \xi r$8 Oe. At zero bias, the two lowest azimuthal spin-wave modes occur at $\omega_v\equiv R\equiv \xi r$9 GHz and $\exp(i\ell\varphi)$00 GHz, and symmetry breaking by an in-plane field produces further splitting and new metastable-state modes (0812.4954).

4. Optical, gaseous, and particle realizations

In structured-light optics, a vortex Pearcey-Gauss beam is obtained by multiplying a Gaussian-apodized Pearcey beam by the helical factor $\exp(i\ell\varphi)$01,

$\exp(i\ell\varphi)$02

A half-wave $\exp(i\ell\varphi)$03-plate with Jones operator

$\exp(i\ell\varphi)$04

converts such scalar beams into vortex Pearcey-Gauss or non-separable vector vortex Pearcey-Gauss beams, depending on the input polarization. Stokes polarimetry shows that the vector beam evolves from a near-field non-homogeneous “spider-web” polarization distribution to predominantly radial or azimuthal patterns at intermediate $\exp(i\ell\varphi)$05, and then toward a nearly uniform vector mode in the far field (Rodríguez-Fajardo et al., 2024).

Gas vortices can also function as optical elements such as lenses or waveguides. In the Burgers-vortex model, the azimuthal velocity is

$\exp(i\ell\varphi)$06

and the radial density profile determines the refractive-index well through the Gladstone-Dale relation. Mach-Zehnder interferometry provides the phase map $\exp(i\ell\varphi)$07, from which the gas number density and wavefront are reconstructed. Zernike decomposition then yields an effective focal length through

$\exp(i\ell\varphi)$08

In the reported experiments, reducing the backing pressure below $\exp(i\ell\varphi)$09 Torr drives the outlet holes into a choked-flow regime in which the mass flow rate, vortex core radius, central density depression, and defocus coefficient $\exp(i\ell\varphi)$10 all become effectively locked in (Kaganovich et al., 2020).

The vortex concept extends to matter waves in accelerator physics. A vortex particle is represented by a cylindrical-wave state

$\exp(i\ell\varphi)$11

with orbital angular momentum $\exp(i\ell\varphi)$12. In accelerator fields, radiative orbital-angular-momentum loss is slow: for electrons in $\exp(i\ell\varphi)$13–$\exp(i\ell\varphi)$14 T, $\exp(i\ell\varphi)$15 s, much longer than typical acceleration times $\exp(i\ell\varphi)$16 s, while for heavy ions $\exp(i\ell\varphi)$17 s. The non-radiative dynamics is governed by precession, and in circular accelerators the OAM tune is $\exp(i\ell\varphi)$18 for $\exp(i\ell\varphi)$19, producing resonances with spacing $\exp(i\ell\varphi)$20 MeV—far denser than conventional spin resonances. This is why linacs are identified as more feasible for vortex-particle acceleration, while Siberian snakes are proposed as tools for OAM manipulation (Karlovets et al., 11 Jul 2025).

5. Geophysical, astrophysical, and planetary contexts

In a rotating annulus with a conical bottom and oscillating upper lid, an axisymmetric inertial-wave attractor first forms in a quasi-linear regime after about 10 forcing periods. At longer times, Triadic Resonance Instability transfers energy into a slow two-dimensional manifold consisting of a regular polygonal cluster of axially oriented cyclonic vortices. In the reported case, coherent vortices emerge after $\exp(i\ell\varphi)$21 and self-organize into an $\exp(i\ell\varphi)$22 polygon with individual radius $\exp(i\ell\varphi)$23–$\exp(i\ell\varphi)$24 cm, vertical-vorticity amplitude $\exp(i\ell\varphi)$25–$\exp(i\ell\varphi)$26, cluster radius $\exp(i\ell\varphi)$27 cm, and slow prograde precession $\exp(i\ell\varphi)$28. The PDF of the vertical micro-Rossby number evolves from nearly Gaussian and symmetric to strongly cyclonic, with skewness rising to $\exp(i\ell\varphi)$29 in the saturated regime (Boury et al., 2020). The authors explicitly note analogues with the polygonal cyclones around Jupiter’s poles.

In protoplanetary discs, large anticyclonic vortices formed by the Rossby Wave Instability at a stationary Gaussian pressure bump can shepherd low-mass planets. Hydrodynamical simulations show that for the standard bump amplitude $\exp(i\ell\varphi)$30, a planet locks to the vortex in a 1:1 configuration at $\exp(i\ell\varphi)$31 with $\exp(i\ell\varphi)$32 behind the vortex center. For a stronger vortex, $\exp(i\ell\varphi)$33, migration is halted farther out at $\exp(i\ell\varphi)$34. Trapping persists down to $\exp(i\ell\varphi)$35, and embryos formed inside the vortex are expelled within $\exp(i\ell\varphi)$36 orbits and then re-lock at the same orbital distance (Ataiee et al., 2014). This suggests that a single vortex may act as a repetitive site of planet formation, provided it survives planetary growth.

Solar and stellar atmospheres provide another astrophysical setting in which vortex geometry, magnetic topology, and wave transport are intertwined. The solar studies summarized above indicate that flux-tube expansion, vortex packing, and vortex-induced torsional Alfvén waves jointly determine whether energy is transported upward or deposited locally as chromospheric heating (Kannan et al., 2024), while the conical solar vortex tubes detected with IVD quantify the local force balance between pressure gradients, Lorentz forces, and eddy-viscosity braking (Silva et al., 2020). A plausible implication is that vortices in magnetized astrophysical plasmas must be analyzed simultaneously as flow structures and as wave-launching or wave-dissipating elements.

6. Mathematical formulations and reduced models

In two-dimensional incompressible Euler flow on $\exp(i\ell\varphi)$37, the vorticity $\exp(i\ell\varphi)$38 is transported by the Biot-Savart velocity

$\exp(i\ell\varphi)$39

Approximating $\exp(i\ell\varphi)$40 by a sum of Dirac masses yields the classical point-vortex ODE system. In exterior domains, the impermeable boundary is replaced by a boundary vortex sheet or by a collection of boundary point vortices chosen so that the flow remains tangent at midpoints between adjacent vortices and the total circulation is conserved. For smooth exterior domains, the corresponding discretization converges on compact subsets away from the boundary at rate $\exp(i\ell\varphi)$41. The same work also introduces a fluid-charge method based on inverting $\exp(i\ell\varphi)$42 rather than $\exp(i\ell\varphi)$43; for strictly convex boundaries, the resulting matrix is strictly diagonally dominant and its condition number remains $\exp(i\ell\varphi)$44, which yields significant numerical improvements (Arsénio et al., 2017).

For the exterior of the unit disk, the construction becomes explicit. The harmonic correction is represented by fixed boundary vortices at uniformly distributed angular points, and the discrete no-penetration condition leads to a discrete Hilbert-transform system. For each $\exp(i\ell\varphi)$45, the system has a unique solution, and for smooth data the approximation converges with rate

$\exp(i\ell\varphi)$46

while the numerical tests reported there show the optimal $\exp(i\ell\varphi)$47 rate in the smooth case (Arsénio et al., 2017).

Reduced point-vortex Hamiltonians provide a complementary analytical framework. For a symmetric pair of counter-rotating point vortices around a circular cylinder in a uniform stream, the dynamics admits Föppl equilibria that are centers in the symmetric subspace. The phase portrait also contains a saddle on the normal line and nilpotent saddles at infinity whose homoclinic loops delimit the region of nonlinear stability. Under antisymmetric perturbations, the same Föppl equilibrium becomes a saddle, and the unstable eigendirections correspond to the alternating pattern associated with vortex shedding (Vasconcelos et al., 2012). This reduced model is important because it shows that bounded vortex motion, escape to upstream infinity, and symmetry-breaking instability can all emerge from a minimal Hamiltonian description.

Taken together, these results suggest that “vortex” is less a single object than a family of rotational or phase-winding structures whose operational definition, interaction law, and dynamical role are fixed by constitutive physics, geometry, and boundary conditions. In some settings vortices transport energy or information coherently upward or along a guide; in others they cluster, annihilate, shed, or cascade toward dissipative scales. The unifying theme across these literatures is not one universal diagnostic, but the recurrent appearance of organized rotation as a dynamically consequential degree of freedom.