Vortex Dipole Evolution in Viscoelastic Media: Effects of Asymmetry, Coupling, and Transverse Shear Waves

This paper numerically investigates the dynamics of Lamb-Oseen vortex dipoles in viscoelastic fluids, revealing that while symmetric dipoles maintain translational motion, asymmetric configurations induce rotation, and increasing viscoelastic coupling generates transverse shear waves that significantly enhance vortex interaction, deformation, and dissipation in strongly coupled regimes.

The Gist

Imagine a fluid not as a simple liquid like water, but as a thick, stretchy substance—like a giant bowl of warm honey or a very dense gel. In this "viscoelastic" world, the fluid can act like a liquid (flowing) and like a solid (springing back) at the same time.

This paper explores what happens when two spinning whirlpools, called a vortex dipole, try to move through this stretchy fluid. Think of these whirlpools as two dancers holding hands, spinning in opposite directions. Usually, they push each other forward, gliding smoothly across the floor.

Here is the story of their journey, broken down into simple parts:

1. The Perfect Pair (Symmetric Dipoles)

Imagine two dancers who are identical twins. They have the same size and the same strength.

• In a normal fluid (like water): They glide in a perfectly straight line. The closer they stand to each other, the faster they move. The farther apart they are, the slower they go. It's a predictable, steady march.
• In the stretchy fluid: Things get interesting. As they move, they don't just glide; they also create ripples in the "honey" around them, like a boat creating waves. These ripples are called transverse shear waves.
  • If the fluid is only slightly stretchy, the dancers barely notice the ripples. They keep moving straight.
  • If the fluid is very stretchy (strongly coupled), the ripples become powerful. They start pushing back against the dancers. The ripples grab energy from the whirlpools, slowing them down and eventually causing them to distort and fall apart. The stronger the "stretchiness," the faster the dancers get tired and dissolve.

2. The Mismatched Pair (Asymmetric Dipoles)

Now, imagine the dancers are not twins. One is a giant, and the other is tiny. Or perhaps one is a heavyweight champion and the other is a lightweight.

• In a normal fluid: Because they are different sizes or strengths, they can't push each other in a straight line. The big one pushes the small one harder than the small one pushes back. Instead of walking straight, they start to spin in a circle. The small dancer orbits around the big one, like a moon orbiting a planet.
• In the stretchy fluid: This spinning motion makes things worse. The powerful ripples (waves) created by the stretchy fluid grab onto the smaller, weaker dancer.
  • The waves stretch the small dancer out, turning them from a round shape into a long, thin noodle.
  • Eventually, the waves completely swallow the small dancer, and they disappear. The big dancer is left alone, still spinning but now without a partner. The paper shows that the more different the two dancers are, the faster this happens.

3. The Energy Balance (The "Poynting" Rule)

The researchers also tracked the "energy budget" of this dance. They found that energy doesn't just vanish; it moves around in three specific ways:

1. The Flow (Convection): The energy moves along with the dancers as they travel.
2. The Ripples (Radiation): Energy is lost as waves shoot out into the surrounding fluid.
3. The Friction (Dissipation): Energy is lost as heat because the fluid is sticky and resists the motion.

The paper proves that these three things always balance each other out perfectly. If the dancers slow down, it's because the ripples and the stickiness took their energy. It's like a bank account where the money spent on travel, waves, and friction always equals the total money lost from the account.

The Main Takeaway

The study reveals that in complex, stretchy fluids (which are found in things like dusty space plasmas or thick gels), symmetry is key to survival.

• If the two whirlpools are perfectly matched, they can travel for a long time, even if the fluid is stretchy.
• If they are mismatched (different sizes or strengths), the stretchy fluid acts like a bully, using its own "ripples" to tear the weaker one apart.

The paper concludes that understanding how these "ripples" interact with spinning structures helps us understand how energy moves and how structures form or break down in complex fluids found in nature.

Technical Summary

Technical Summary: Vortex Dipole Evolution in Viscoelastic Media

Problem Statement
This study investigates the nonlinear evolution and energy transport of a Lamb–Oseen vortex dipole within a strongly coupled dusty plasma (SCDP), which exhibits viscoelastic (VE) fluid behavior. While symmetric dipole vortices have been extensively studied in hydrodynamic and dusty plasma contexts, systematic investigations into asymmetric counter-rotating dipoles remain limited. Realistic media often introduce asymmetries through turbulent fluctuations, spatial inhomogeneities, or variations in vortex core geometry and circulation strength. Furthermore, in SCDPs, the presence of transverse shear (TS) waves—arising from the medium's viscoelastic nature—introduces complex wave-vortex coupling mechanisms that are not present in purely inviscid or Newtonian fluids. The paper seeks to understand how asymmetry (in core radii or circulation), coupling strength, and TS waves collectively influence dipole stability, deformation, propagation speed, and energy redistribution.

Methodology
The authors employ the incompressible generalized hydrodynamic (i-GHD) model, which suppresses longitudinal (compressional) modes to isolate the influence of transverse shear waves. The model consists of coupled continuity, momentum, and Poisson equations, reformulated into a set of convective equations for vorticity ($\xi_z$) and strain ($\vec{\psi}$).

• Numerical Approach: Simulations are conducted in a two-dimensional domain ($24\pi \times 24\pi$) using the LCPFCT (Flux-Corrected Transport) method to solve the coupled evolution equations.
• Configuration: The study examines three distinct dipole configurations:
  • Case A: Symmetric dipoles ($\Omega_1 = \Omega_2, a_1 = a_2$) with varying initial separation distances ($b_0$) and circulation strengths ($\Omega_0$).
  • Case B: Asymmetric dipoles with unequal core radii ($a_1 \neq a_2$) but equal circulation.
  • Case C: Asymmetric dipoles with unequal circulation strengths ($\Omega_1 \neq \Omega_2$) but equal core radii.
• Regimes: Simulations cover inviscid fluids and viscoelastic fluids across weak, moderate, and strong coupling regimes, defined by the ratio of viscosity to relaxation time ($\eta/\tau_m$).
• Conservation Analysis: A Poynting-like conservation theorem is derived and applied to quantify the relative contributions of radiative flux (TS wave emission), convective transport, and dissipative processes during dipole evolution.

Key Contributions and Results

1. Symmetric Dipole Dynamics:

   • In inviscid fluids, symmetric dipoles exhibit sustained translational motion with a speed inversely proportional to the initial separation distance and directly proportional to circulation strength, consistent with point-vortex theory.
   • In viscoelastic media, the emission of TS waves modifies this behavior. In weakly coupled regimes, the dipole remains largely intact with minor structural distortion. However, as coupling strength increases, TS waves accelerate energy and momentum extraction from the dipole. In strongly coupled regimes, this leads to rapid structural deformation and the eventual dissipation of the vortex pair, with higher circulation strengths providing greater stability against this decay.

2. Asymmetric Dipole Dynamics:

   • Inviscid Limit: Asymmetry in core radii or circulation breaks the balance of induced velocities. Instead of linear translation, the dipole undergoes rotational motion where the weaker (or smaller) vortex orbits the stronger (or larger) one, resulting in curved trajectories.
   • Viscoelastic Effects: In VE media, the interplay between asymmetry and TS waves is critical. In weakly coupled regimes, the dynamics resemble the inviscid case with slight modifications. In moderately and strongly coupled regimes, TS waves significantly enhance the strain interaction between the vortices. This leads to the accelerated deformation of the weaker vortex into elliptical or filamentary structures, followed by its rapid engulfment and dissipation by the wave field. The stronger vortex often persists longer, behaving approximately as a rotating monopole.

3. Energy Transport and Conservation:

   • The study validates a Poynting-like conservation theorem, demonstrating that the time evolution of the domain-integrated energy density is governed by the balance of three terms: radiative flux ($S$) due to TS wave emission, convective transport ($T$) due to dipole motion, and dissipative loss ($P$) due to viscoelastic relaxation.
   • Numerical results confirm that for symmetric dipoles, convective transport dominates the energy balance during propagation. For asymmetric dipoles, the curved trajectory leads to continuous interaction with the boundary, resulting in persistent radiative flux contributions. In all cases, the sum of these terms dynamically compensates to maintain global energy balance.

Significance and Claims
The paper claims to provide insight into wave-vortex coupling in complex fluids, specifically within the context of strongly coupled dusty plasmas. The findings suggest that transverse shear waves act as a primary mechanism for enhancing vortex-vortex coupling in strongly coupled regimes, driving the transition from coherent dipole motion to asymmetric, dissipative dynamics.

The authors note that these results have implications for understanding transport processes and structure formation in strongly coupled plasmas. They also suggest that similar mechanisms involving elastic stresses and wave propagation may be relevant to other viscoelastic media, such as polymeric fluids, geophysical flows, and biological systems, where vortex stability and mixing are influenced by wave modes. The study explicitly states that future work will extend these findings to fully compressible regimes, include higher-order nonlinear effects, and explore three-dimensional configurations to better represent realistic conditions.