Fast relaxation of a viscous vortex in an external flow

This paper analyzes the evolution of concentrated vortices in two-dimensional external flows, deriving an accurate approximation for high Reynolds numbers and proving that sharply peaked Gaussian vortices rapidly relax to this solution on a timescale much shorter than the diffusive time due to enhanced dissipation.

The Gist

Imagine you are watching a drop of ink swirl in a glass of water. If the water is perfectly still, the ink drop slowly spreads out into a fuzzy, round cloud, getting wider and wider as time goes on. This is what happens to a "vortex" (a spinning whirlpool) in a calm fluid.

But what happens if you start stirring the glass? What if the water itself is flowing, stretching, and twisting around that ink drop?

This is the question Martin Donati and Thierry Gallay tackle in their paper. They study how a concentrated, spinning vortex behaves when it gets swept along by a smooth, flowing current. They look at two specific scenarios: one where the vortex starts as a perfect, mathematical point, and another where it starts as a slightly "messy" blob.

Here is the story of their findings, explained with some everyday analogies.

1. The Perfectly Prepared Vortex (The "Ideal" Start)

Imagine you drop a single, perfect grain of sand into a river. In the world of fluid dynamics, this is a "Dirac mass"—a point with no size at all.

The authors show that even though this grain is infinitely small, it doesn't just get washed away randomly. Instead, it behaves like a surfer riding a wave.

• The Surf: The vortex center moves exactly where the river current tells it to go.
• The Shape: As it moves, the river's current stretches and squeezes the vortex. It stops being a perfect circle and turns into an ellipse (like a squashed circle).
• The Result: The authors created a mathematical "map" that predicts exactly where this vortex will be and how squashed it will get. They proved that if the water is very thick (viscous) but not too thick, this map is incredibly accurate. The vortex stays tight and follows the path predicted by their formula for a very long time.

2. The Ill-Prepared Vortex (The "Messy" Start)

Now, imagine you don't drop a perfect grain of sand. Instead, you drop a small, round puff of smoke that is perfectly circular, but the river is flowing in a way that wants to stretch it into an oval immediately.

This is the "ill-prepared" scenario. The vortex starts as a round blob, but the river is trying to turn it into an oval.

• The Struggle: At first, the vortex fights back. It wobbles. It oscillates. It's like a rubber band that you stretch and then let go; it snaps back and forth a few times before settling into the shape the river wants.
• The Surprise: The authors discovered that this wobbly phase is extremely short. Much shorter than you would expect.
• The "Magic" Dissipation: Usually, we think of viscosity (thickness) as a slow, lazy force that just makes things spread out. But here, the authors found that the internal friction inside the vortex core acts like a super-charged shock absorber. It quickly damps out those wobbles.
• The Outcome: The vortex rapidly "relaxes." It stops fighting the current, snaps into the perfect oval shape predicted by the "ideal" map, and then happily surfs along.

The "Enhanced Dissipation" Analogy

Why does the messy vortex fix itself so fast? The paper uses a concept called enhanced dissipation.

Think of a spinning top. If you spin it on a rough table, it slows down. But if you spin it on a table that is also vibrating or shearing (moving sideways), the top loses its energy much faster. The motion of the table helps the friction work harder.

In the paper, the external flow (the river) shears the vortex. This shearing action mixes the fluid inside the vortex so efficiently that the viscosity can "eat up" the irregularities (the roundness) much faster than it would in still water. It's like the river is helping the vortex clean itself up in record time.

The Big Picture

The paper is essentially a rigorous proof of two things:

1. Prediction: We can accurately predict how a tiny, concentrated whirlpool will move and deform in a flowing river, even if the river is twisting and turning.
2. Resilience: Even if the whirlpool starts with the wrong shape (too round for the current), it doesn't stay that way. It quickly "learns" the shape the current demands and settles down, thanks to a hidden speed-up in how fast it loses energy.

Why Does This Matter?

While this sounds like abstract math, it helps us understand real-world phenomena:

• Weather: How do hurricanes or tornadoes move and change shape when they get caught in larger wind currents?
• Oceanography: How do swirling eddies in the ocean interact with ocean currents?
• Engineering: How do fuel injectors in engines mix fuel (which often starts as a concentrated drop) with air?

The authors have provided a "user manual" for these swirling things, showing us that no matter how you start them, they quickly find their rhythm and follow the flow, provided the fluid is thick enough to keep them from flying apart.

Technical Summary

1. Problem Statement

The paper investigates the long-time evolution of a concentrated vortex advected by a smooth, divergence-free, external velocity field $f(x,t)$ in two spatial dimensions. The fluid dynamics are governed by the 2D incompressible Navier-Stokes equations in vorticity form:
$$\partial_t \omega + (u + f) \cdot \nabla \omega = \nu \Delta \omega$$
where $\omega$ is the vorticity, $u$ is the self-induced velocity (via the Biot-Savart law), and $\nu$ is the kinematic viscosity.

The study focuses on the high Reynolds number regime ($\nu \to 0$) and addresses two distinct initial data scenarios:

1. Well-prepared data: The initial vorticity is a Dirac mass ($\omega_0 = \Gamma \delta_{z_0}$). This represents an idealized point vortex.
2. Ill-prepared data: The initial vorticity is a sharply peaked Gaussian (or radially symmetric blob) that does not initially match the deformation induced by the external strain.

The Core Question: How does the vortex core deform under external shear, and how quickly does an "ill-prepared" (symmetric) vortex relax into the "well-prepared" (deformed) state? The authors aim to provide rigorous error estimates for asymptotic approximations of the solution.

2. Methodology

The authors employ a combination of perturbative asymptotic analysis and energy estimates in weighted spaces, leveraging recent advances in enhanced dissipation.

A. Self-Similar Variables and Perturbation Expansion

To handle the singularity of the Dirac mass and the small viscosity limit, the authors introduce self-similar variables centered at the vortex position $z(t)$:
$$\xi = \frac{x - z(t)}{\sqrt{\nu t}}, \quad \Omega(\xi, t) = \frac{\nu t}{\Gamma} \omega(x, t)$$
They construct an approximate solution $\Omega_{\text{app}}$ as a perturbation expansion in the aspect ratio $\varepsilon(t) = \sqrt{\nu t}/d$ (where $d = \sqrt{\Gamma T_0}$ is the effective vortex size):
$$\Omega_{\text{app}} = \Omega_0 + \varepsilon^2 \Omega_2 + \varepsilon^3 \Omega_3 + \varepsilon^4 \Omega_4$$

• $\Omega_0$ is the standard Lamb-Oseen vortex (Gaussian).
• Higher-order terms ($\Omega_2, \dots$) account for the deformation of streamlines due to the external strain rates $a(t)$ and $b(t)$.
• The vortex position $z(t)$ is determined by an ODE that includes a viscous correction term: $z'(t) = f(z(t), t) + \nu t \Delta f(z(t), t)$.

B. Linearized Stability and Enhanced Dissipation

For the "ill-prepared" case (Gaussian initial data), the initial error is large because the symmetric profile does not match the external strain. The authors decompose the error $w$ into:

1. A linear component $w_0$ solving the linearized equation $\partial_\tau w_0 = (L - \delta^{-1}\Lambda)w_0$.
2. A nonlinear correction $\tilde{w}$.

They utilize enhanced dissipation estimates (from Li, Wei, and Zhang) which show that the linear operator $L - \delta^{-1}\Lambda$ (where $\Lambda$ is the advection operator) generates a semigroup that decays faster than the standard diffusive rate. Specifically, the decay rate is proportional to $\exp(-c \tau / \delta^{1/3})$, where $\delta = \nu/\Gamma$ is the inverse Reynolds number. This implies relaxation occurs on a time scale much shorter than the diffusive time.

C. Weighted Energy Estimates

To control the remainder terms globally in time, the authors construct sophisticated weighted $L^2$ energy functionals.

• They define a weight function $p_\varepsilon(\xi, t)$ that adapts to the elliptical deformation of the vortex core.
• The domain is split into three regions: an inner region (near the core), an intermediate region, and a far-field region.
• This allows them to derive differential inequalities that control the growth of the error, ensuring the approximation remains valid for long times (up to $T \sim T_0 \ln(1/\delta)$).

3. Key Contributions and Results

Result 1: Well-Prepared Initial Data (Dirac Mass)

Theorem 1.6 establishes that if the initial vorticity is a Dirac mass, the solution remains close to the constructed approximation $\omega_{\text{app}}$ for all $t \in (0, T)$.

• Approximation: The solution is a Lamb-Oseen vortex with a deformed core (elliptical streamlines) and a position $z(t)$ corrected by viscosity.
• Error Estimate: The $L^1$ error is bounded by $O(\varepsilon^2(\varepsilon + \delta))$, which is significantly sharper than previous $O(\varepsilon)$ estimates.
• Significance: This proves that the "well-prepared" state (deformed vortex) is the correct leading-order description for point vortices in external flows, including viscous corrections to the trajectory.

Result 2: Ill-Prepared Initial Data (Gaussian Vortex)

Theorem 1.9 addresses the case where the initial data is a radially symmetric Gaussian.

• Fast Relaxation: The solution rapidly relaxes to the "well-prepared" approximation derived in Theorem 1.6.
• Relaxation Rate: The convergence rate is governed by the enhanced dissipation effect, scaling as $\delta^{1/6} (\log(1/\delta))^{1/2}$.
• Time Scale: The relaxation happens on a time scale $t \sim t_0 (1 + \tau_\delta)$ where $\tau_\delta \sim \delta^{1/3} \log(1/\delta)$. This is much faster than the diffusive time scale ($t \sim 1/\delta$).
• Significance: This rigorously justifies the physical observation that symmetric vortices quickly adapt their shape to external shear, regardless of their initial symmetry.

Result 3: Vortex Trajectory

The paper clarifies the motion of the vortex center.

• The center of vorticity $\bar{z}(t)$ follows the ODE $z' = f(z) + \nu t \Delta f(z) + O(\varepsilon^3)$.
• The naive trajectory $z' = f(z)$ (ignoring viscosity) is only accurate to $O(\varepsilon^2)$, whereas the corrected trajectory is accurate to $O(\varepsilon^3)$.

4. Significance and Impact

1. Rigorous Justification of Asymptotics: The paper provides the first rigorous proof that the "metastable" deformed vortex state (often used in physics literature) is indeed the attractor for concentrated vortices in external flows.
2. Enhanced Dissipation Application: It successfully applies enhanced dissipation estimates (previously used for axisymmetrization) to the problem of shape relaxation in an external strain field, demonstrating the universality of this mechanism in vortex dynamics.
3. Improved Trajectory Estimates: By identifying the viscous correction term $\nu t \Delta f$ in the vortex trajectory, the authors refine the classical point-vortex model, showing that viscosity plays a non-negligible role in the drift of concentrated vortices even at high Reynolds numbers.
4. Mathematical Technique: The construction of the specific non-radial weight function $p_\varepsilon$ allows for precise control of the solution in the presence of time-dependent external flows, overcoming limitations of previous methods that were restricted to short times or specific flow configurations.

In summary, Donati and Gallay demonstrate that a concentrated vortex in a smooth external flow quickly forgets its initial shape (if ill-prepared) and settles into a specific deformed profile determined by the external strain, while its center follows a trajectory slightly modified by viscosity. This relaxation occurs on a "fast" time scale due to the interplay between advection and diffusion (enhanced dissipation).