Vorticity Packing Effects on Long Time Turbulent Transport in Decaying Two-Dimensional Incompressible Navier-Stokes Fluids

This study demonstrates that the vorticity packing fraction in decaying two-dimensional Navier-Stokes turbulence governs the transition from point-vortex to finite-size vortex equilibria, which in turn dictates a corresponding shift in Lagrangian tracer transport from sub-diffusive orbital trapping to super-diffusive linear motion as packing increases.

The Gist

Imagine a giant, invisible swimming pool filled with water. Now, imagine you paint thin, alternating stripes of red and blue dye across the surface of this pool. The red stripes are swirling clockwise, and the blue stripes are swirling counter-clockwise. This is the starting point of the experiment described in this paper.

The scientists wanted to see what happens when these swirling stripes interact, break apart, and eventually settle down. But they didn't just watch the water; they also dropped thousands of tiny, invisible "tracers" (like tiny specks of glitter) into the water to see how they moved.

Here is the story of what they found, broken down into simple concepts:

1. The Setup: Packing the Dye

The key variable in their experiment was how tightly they packed the stripes.

• Loose Packing: Imagine just two wide stripes of red and blue. There is a lot of empty space between them.
• Tight Packing: Imagine cramping 20 narrow stripes into the same space. They are squished right next to each other.

The scientists call this the "Vorticity Packing Fraction" (VPF). It's essentially a measure of how crowded the swirling water is at the very beginning.

2. The Spark: The "Rippling" Instability

When the water starts moving, the boundary between the red and blue stripes becomes unstable. It's like when you rub your hands together quickly; the friction creates heat. Here, the friction between the opposing swirls creates a wavy, rolling motion called the Kelvin-Helmholtz instability.

Think of it like wind blowing over the ocean: the water doesn't just stay flat; it starts to curl up into little waves and eventually big whirlpools.

3. The Evolution: From Chaos to Order

As time passes, these little whirlpools crash into each other. In the world of 2D water (like a flat sheet), when two whirlpools of the same color meet, they merge to become one giant, stronger whirlpool. This is called an inverse energy cascade—small swirls combine to make big swirls.

Eventually, the chaos settles into a calm state dominated by a few massive structures. Usually, this ends up as a dipole: a giant pair of swirls (one red, one blue) locked together, drifting across the pool like a slow-moving boat.

4. The Big Discovery: How "Crowded" Changes the Journey

The paper's main finding is that how crowded the starting stripes were completely changed how the "glitter" (tracers) moved.

The "Loose" Case (Low Packing)

• The Scene: With wide gaps between stripes, the water moves slowly to start. The "glitter" gets pushed mostly in one direction (left or right) by the initial flow.
• The Movement: The glitter moves in a very predictable, straight line for a while, then gets trapped.
• The Trap: Eventually, the giant red/blue pair forms. The glitter gets stuck orbiting around these giant swirls, like a moon orbiting a planet. It doesn't go very far.
• The Result: The movement is slow and stuck (sub-diffusive). The glitter stays in a specific area and doesn't mix well.

The "Crowded" Case (High Packing)

• The Scene: With 20 tight stripes, the water goes crazy almost immediately. The instability happens fast, and the turbulence is intense and chaotic in all directions.
• The Movement: The "glitter" is thrown around violently in every direction. It mixes rapidly.
• The Result: The movement is fast and wild (super-diffusive). The glitter travels huge distances very quickly.
• The Twist: In the most crowded case (62.5% packing), the giant red/blue pair doesn't just spin in place. Instead, it shoots off in a straight line diagonally across the pool, carrying the glitter with it at high speed.

5. The Connection: The Map and the Traveler

The paper connects two different ways of looking at the water:

1. The Map (Eulerian View): Looking at the water from a fixed point (like a camera on the wall) to see the shape of the swirls.
2. The Traveler (Lagrangian View): Following the "glitter" to see where it goes.

The scientists found a perfect match between the two:

• If the water looks like a collection of distinct, separate points (loose packing), the glitter gets stuck in orbits.
• If the water looks like a dense, continuous patch of swirls (tight packing), the glitter flies freely and quickly.

Summary in a Nutshell

Think of the water as a dance floor.

• Loose Packing: The dancers are far apart. They spin slowly, and if you drop a coin on the floor, it just sits there or moves in a small circle around a dancer. It's a slow, trapped dance.
• Tight Packing: The dance floor is packed shoulder-to-shoulder. The energy is high, everyone is bumping into each other, and the coin gets thrown across the room, bouncing wildly. It's a fast, chaotic dance.

The paper proves that by simply changing how tightly you pack the initial swirls, you can switch the entire system from a slow, trapped state to a fast, explosive state. This helps scientists understand how energy and matter move in fluids, from weather patterns to plasma in stars.

Technical Summary

Technical Summary: Vorticity Packing Effects on Long Time Turbulent Transport in Decaying Two-Dimensional Incompressible Navier-Stokes Fluids

Problem Statement
This study investigates the long-time Lagrangian transport dynamics of passive tracer particles in decaying, two-dimensional (2D), incompressible Navier–Stokes turbulence. Building upon recent findings by Biswas et al. (2022), which established that the initial total circulation—quantified by the vorticity packing fraction (VPF)—determines the late-time Eulerian statistical equilibria (transitioning from point-vortex to finite-size patch-vortex regimes), this paper examines how these initial conditions influence the evolution of tracer transport. The central objective is to correlate the Eulerian statistical states with the underlying Lagrangian transport across three distinct evolutionary stages: the early linear onset, the intermediate turbulence development, and the late-time coherent dipole evolution.

Methodology
The research employs high-resolution, high-Reynolds-number direct numerical simulations (DNS) using a custom-developed, GPU-accelerated framework.

• Fluid Solver: The in-house solver GHD2D utilizes a pseudospectral method with 2/3 dealiasing and a second-order Adams–Bashforth time integration scheme to solve the dimensionless 2D incompressible Navier–Stokes equations in a doubly periodic domain ($L=2\pi$).
• Initial Conditions: Turbulence is triggered via the Kelvin–Helmholtz instability (KHI) using parallel, alternating strips of positive and negative vorticity. The study systematically varies the number of strips ($n = 2, 4, 8, 16, 20$), corresponding to initial VPFs of 6.25% to 62.5%. Two initial vorticity profiles were tested (tapered and sharp-edged); the sharp-edged profile was selected for subsequent analysis due to its slightly faster instability growth and robustness.
• Tracer Solver: A dedicated particle-tracking solver (GHD2D-TP) was developed and benchmarked against a kinematic 2D chaotic flow (Folgia et al., 2023) and analytical KHI growth rates. It integrates the Lagrangian equations of motion using a fourth-order Runge–Kutta (RK4) algorithm with linear interpolation of the velocity/stream-function fields.
• Analysis: Transport is quantified via the time evolution of the mean-square displacement (MSD), time-dependent diffusion coefficients, and the probability distribution functions (PDFs) of particle positions and velocities. The study analyzes the transition from ballistic to diffusive/super-diffusive regimes and the emergence of Gaussianity in PDFs.

Key Results
The simulations reveal a strong dependence of turbulent transport characteristics on the initial vorticity packing fraction (VPF):

1. Turbulence Onset and Development: Increasing the VPF accelerates the transition from the linear shear-driven regime to the nonlinear turbulent regime.

   • Low VPF (6.25%–12.5%): Turbulence onset is delayed ($T_{TO} \approx 43\text{--}45$). The flow evolves into a sparse, highly anisotropic state with few large-scale vortex rolls.
   • High VPF (50%–62.5%): Turbulence onset occurs earlier ($T_{TO} \approx 33\text{--}36$). The flow rapidly develops into a dense, nearly isotropic background with numerous interacting vortex rolls.
   • Transport Evolution: In the intermediate stage, increasing VPF drives a transition from sub-diffusive, anisotropic transport (dominated by shear in the x-direction) to super-diffusive, isotropic transport (driven by instability-induced fluctuations in the y-direction).

2. Late-Time Coherent Structures: The long-time dynamics are governed by the largest-scale coherent dipolar vortices.

   • Low to Moderate VPF: Dipoles form and undergo orbital rotation along large-radius trajectories. This results in sub-diffusive trapping of tracer particles.
   • Highest VPF (62.5%): The dipoles exhibit linear translational motion along diagonal paths rather than orbital rotation. This leads to super-diffusive transport, where the total MSD grows faster than linearly ($\alpha > 1$).

3. Statistical Distributions:

   • Position PDFs: Low VPF cases show a slow evolution toward Gaussianity and persistent anisotropy (unequal $n(x)$ and $n(y)$). High VPF cases rapidly achieve Gaussianity and isotropy. Notably, the 62.5% case develops heavy tails at very long times, signaling a departure from normality consistent with super-diffusion.
   • Velocity PDFs: Initially bimodal due to shear, velocity distributions evolve toward Gaussian forms as turbulence develops. At late times, they regain bimodality corresponding to the counter-rotating dipole populations. High VPF cases show faster Gaussianization and sustained overlap between $x$ and $y$ components.

4. Eulerian-Lagrangian Correlation: The study establishes a direct correspondence between the Eulerian statistical equilibrium and Lagrangian transport.

   • Low VPF cases align with point-vortex theory (sinh–Poisson) and exhibit sub-diffusive, anisotropic transport.
   • High VPF cases align with the finite-size patch-vortex (KMRS) theory and exhibit super-diffusive, isotropic transport.

Significance and Claims
The paper claims that the initial vorticity packing fraction is a critical parameter governing not only the statistical equilibrium of the fluid field but also the fundamental nature of particle transport in decaying 2D turbulence. The primary contribution is the demonstration of a strong correlation between the transition from point-vortex-dominated to finite-size-vortex-dominated Eulerian equilibria and the concurrent shift from sub-diffusive to super-diffusive Lagrangian transport.

Specifically, the authors highlight that while increasing VPF generally reduces overall transport in the intermediate stages, the extreme case of the highest packing fraction (62.5%) results in the largest total transport at late times due to anomalous super-diffusive behavior driven by linear dipole translation. This finding underscores the role of the classical exclusion principle (incompressibility) and total circulation in determining the final dynamical state of the fluid and the associated mixing efficiency. The work provides a unified framework linking statistical mechanical models of 2D turbulence with observable Lagrangian transport phenomena.