Two-dimensional fluids via matrix hydrodynamics

Two-dimensional (2-D) incompressible, inviscid fluids produce fascinating patterns of swirling motion. How and why the patterns emerge are long-standing questions, first addressed in the 19th century by Helmholtz, Kirchhoff, and Kelvin. Countless researchers have since contributed to innovative techniques and results. Yet, the overarching problem of swirling 2-D motion and its long-time behavior remains largely open. Here we shed light on this problem via a link to isospectral matrix flows. The link is established through V. Zeitlin’s beautiful model for the numerical discretization of Euler’s equations in 2-D. When considered on the sphere, Zeitlin’s model offers deep connections between 2-D hydrodynamics and unitary representations of the rotation group. Consequently, it provides a dictionary that maps hydrodynamical concepts to matrix Lie theory, which in turn gives connections to matrix factorizations, random matrices, and integrability theory, for example. Results about finite-dimensional matrices can then be transferred to infinite-dimensional fluids via quantization theory, which is here used as an analysis tool (albeit traditionally describing the limit between quantum and classical physics). We demonstrate how the dictionary is constructed and how it unveils techniques for 2-D hydrodynamics. We also give accompanying convergence results for Zeitlin’s model on the sphere.

💡 Research Summary

The paper investigates the long‑standing problem of vortex formation and long‑time dynamics in two‑dimensional (2‑D) incompressible, inviscid fluids by establishing a precise correspondence between the infinite‑dimensional Euler equations on the sphere and a finite‑dimensional isospectral matrix flow introduced by V. Zeitlin. The authors begin by recalling the geometric formulation of the 2‑D Euler equations: the vorticity ω evolves according to ∂ₜω + {ω,ψ}=0 with ψ = (−Δ)⁻¹ω, and the system possesses a Lie–Poisson structure on the dual of the Poisson algebra C∞(S²)/ℝ. This structure yields infinitely many Casimir invariants C_f(ω)=∫_S² f(ω) dμ, which encode the continuous distribution of vorticity and are essential for the dynamics but are invisible in the classical point‑vortex model.

To retain these invariants while obtaining a tractable numerical scheme, the authors adopt Zeitlin’s discretisation, which is a Lie‑Poisson preserving quantisation of the fluid on the sphere. Using the representation theory of the rotation group SO(3), they map spherical harmonics Y_{ℓm} to the matrix basis T_{ℓm} of the Lie algebra 𝔲(N) with N=2ℓ+1. The Poisson bracket of functions is replaced by the commutator of matrices scaled by iℏ, where ℏ≈1/N. In this way the continuous vorticity field ω is replaced by a Hermitian matrix ω_N(t)∈𝔲(N) that evolves according to an isospectral Lax equation
  ˙ω_N =