On the stability of viscous Riemann ellipsoids

This study investigates the linear stability of viscous Riemann ellipsoids by deriving a generalized Poincaré equation for inviscid oscillations and applying boundary-layer analysis to quantify first-order viscous corrections, ultimately providing comprehensive stability diagrams that elucidate the roles of rotation, strain, and diffusion in geophysical and astrophysical flows.

The Gist

Imagine a giant, spinning ball of fluid floating in space. It's not a perfect sphere; it's squashed into an egg shape (an ellipsoid) because it's spinning so fast. Now, imagine that inside this spinning ball, the fluid isn't just rotating as a solid block; it's also sloshing around with its own internal currents. This is what scientists call a Riemann ellipsoid.

For over a century, physicists have tried to figure out: Is this spinning, sloshing ball stable, or will it eventually tear itself apart?

This paper by Joris Labarbe is like a new, high-tech manual for answering that question, looking at the problem in two different scenarios: when the fluid is perfectly slippery (no friction) and when it has a tiny bit of stickiness (viscosity).

Here is the breakdown of what the paper does, using simple analogies:

1. The "Perfectly Slippery" Scenario (Inviscid Limit)

First, the author looks at the ball as if the fluid were like water with zero friction. In this world, the fluid can slide past itself without any resistance.

• The Old Way vs. The New Way: Previously, scientists tried to solve this using a method called the "virial tensor method." Think of this like trying to solve a complex puzzle by moving huge, heavy blocks around. It gets incredibly difficult and slow if you want to look at tiny, detailed ripples on the surface. Another method was like using a telescope that only sees things far away (short-wavelength approximations), missing the details up close.
• The New Tool: Labarbe invents a new mathematical "lens" (a generalized Poincaré equation). Imagine this as a super-smart calculator that can instantly tell you how any size of ripple—from a tiny pebble-sized wave to a massive ocean swell—will behave on this spinning ball.
• The Discovery: Using this new tool, the author confirms that almost all of these spinning, sloshing balls are actually unstable. They are like a spinning top that is wobbling so much it's about to fall over. The paper maps out exactly when and why they become unstable, showing that the internal sloshing (strain) and the spinning (rotation) work together to make the shape wobble and eventually break apart.

2. The "Sticky" Scenario (Viscosity)

Next, the author adds a tiny bit of "honey" to the fluid. In the real world, fluids have viscosity (thickness/friction). Usually, we think of friction as a stabilizer—like how a brake slows a car down to stop it from crashing.

• The Counter-Intuitive Twist: The paper finds something surprising. In these spinning balls, adding a tiny bit of friction doesn't just slow the wobble down; it can actually make the instability worse.
• The Analogy: Imagine a child on a swing. If you push them at the wrong time, they go higher. Friction in this specific spinning system acts like a mischievous friend who pushes the swing at the exact wrong moment, causing the wobble to grow faster than it would have without the friction.
• The Boundary Layer: To figure this out, the author looks at a very thin layer of fluid right against the surface of the ball (the "boundary layer"). It's like looking at the very thin skin of an orange to understand how the whole fruit reacts to being squeezed. By analyzing this thin skin, the author calculated exactly how much the "stickiness" changes the stability.

3. The Big Picture

The paper doesn't just say "it's unstable." It draws a detailed map (a stability diagram) showing exactly which shapes and spin speeds lead to disaster.

• What it means: It turns out that if you have a spinning, self-gravitating fluid body (like a star or a planet) with internal currents, it is very fragile. Even a tiny bit of friction can trigger a chain reaction that makes the shape collapse or change dramatically.
• The Takeaway: The author has built a universal toolkit that is faster and more accurate than previous methods. It allows scientists to predict the fate of these cosmic fluid balls with much greater precision, showing that the combination of spinning, internal sloshing, and even tiny amounts of friction creates a recipe for instability.

In short: The paper provides a new, faster way to calculate how spinning, squishy balls in space behave. It reveals that these balls are naturally unstable, and surprisingly, adding a little bit of "stickiness" (friction) can sometimes make them fall apart even faster, rather than holding them together.

Technical Summary

Technical Summary: On the Stability of Viscous Riemann Ellipsoids

Problem Statement
This study addresses the linear stability of Riemann ellipsoids, which are equilibrium configurations of self-gravitating, homogeneous, incompressible fluid masses characterized by a combination of solid-body rotation and uniform internal vorticity (specifically the S-type family). While the stability of Maclaurin spheroids and Jacobi ellipsoids has been extensively analyzed, a comprehensive linear stability theory for Riemann ellipsoids valid for arbitrary harmonic perturbations and including weak viscosity has been lacking. Previous inviscid analyses relied on the virial tensor method, which becomes computationally prohibitive at high harmonic degrees, or short-wavelength (WKB) approximations that lack uniform validity across all wavelengths. Furthermore, the specific mechanisms of viscosity-driven instabilities in these non-axisymmetric, internally sheared flows had not been systematically treated.

Methodology
The author develops a unified linear stability framework applicable to both the inviscid limit and the presence of weak viscosity (small Ekman number, $Ek \ll 1$).

1. Inviscid Formulation:

   • The study derives a generalized Poincaré equation governing small-amplitude oscillations about the equilibrium state.
   • A family of polynomial solutions is constructed for the hydrodynamic potential. This family extends the classical results of Cartan (1922) for strainless flows to include flows with a uniform strain field.
   • Unlike the virial tensor method, this approach allows for the direct computation of the full mode structure for ellipsoidal harmonics of arbitrary order ($n$).
   • An analytic dispersion relation is derived by expanding the disturbance fields in surface ellipsoidal harmonics and applying boundary conditions.

2. Viscous Formulation:

   • For weak viscosity, the analysis employs a boundary-layer approach based on Prandtl's theory.
   • The study focuses on the case $n=2$, where the inviscid velocity perturbations remain harmonic and satisfy the Navier-Stokes equations exactly; for $n>2$, bulk corrections would be required, which are not treated here.
   • Viscous effects are incorporated by solving coupled Prandtl-type equations within a thin boundary layer adjacent to the free surface.
   • This yields first-order viscous corrections to the inviscid dispersion relation, accounting for the stress-free boundary conditions and the interaction between the inviscid bulk flow and the viscous boundary layer.

Key Contributions

• Generalized Poincaré Equation: The derivation of a generalized Poincaré equation for rotating fluids with internal strain, accompanied by a specific family of polynomial solutions that generalizes Cartan's classical results.
• Analytic Dispersion Relations: The construction of explicit analytic dispersion relations for inviscid Riemann ellipsoids that remain computationally efficient at arbitrary harmonic degrees, contrasting with the limitations of the virial tensor method.
• Viscous Stability Theory: The development of a systematic procedure to calculate first-order viscous corrections via boundary-layer analysis, providing an explicit criterion for secular instability in Riemann ellipsoids.
• Stability Diagrams: The presentation of stability diagrams over the parameter space of admissible Riemann ellipsoids, illustrating the interplay between rotation, internal strain, and diffusion.

Results

• Inviscid Stability: The study confirms that almost all admissible S-type Riemann ellipsoids are unstable to second-order harmonic disturbances. The analysis recovers the known bifurcation points for Jacobi ellipsoids (where $f=0$) and identifies additional instabilities for $f \neq 0$ arising from the crossing of distinct frequency branches (elliptic instability).
• Viscous Instability: The inclusion of viscosity reveals the existence of secular instabilities driven by dissipation. Specifically, viscosity induces instabilities in Riemann ellipsoids outside the inviscid interval of instability.
• Comparison with Maclaurin Spheroids: While viscosity destabilizes Maclaurin spheroids at the bifurcation point to the Jacobi sequence (a symmetry-breaking event), the instability in Riemann ellipsoids occurs in intrinsically non-axisymmetric configurations. The mechanism involves the destruction of gyroscopic stabilization by diffusion, leading to exponential growth in modes that are stable in the inviscid limit.
• Avoided Crossings: In the viscous regime, the dispersion curves exhibit avoided crossings near the inviscid Hamilton–Hopf bifurcation points, a feature consistent with the behavior of dissipative systems.

Significance
The paper claims to fill a critical gap in the literature by providing a tractable, unified linear stability theory for Riemann ellipsoids that is valid for arbitrary harmonic perturbations and includes viscous effects. The methodology offers a computationally efficient alternative to the virial tensor method and avoids the restrictions of WKB approximations. The findings provide a systematic description of dissipation-induced instabilities in rotating, self-gravitating fluids with internal shear, offering a framework relevant to the secular evolution of rotating planets, stars, and fluid interiors in geophysical and astrophysical contexts. The work also sets the stage for future investigations into the Chandrasekhar–Friedman–Schutz (CFS) instability in differentially rotating bodies, though the current analysis remains within a Newtonian framework.