Rapidly rotating internally heated convection: bounds on long-time averages

This paper derives an asymptotically reduced model for rapidly rotating, internally heated convection at infinite Prandtl number to establish rigorous, optimized bounds on mean temperature and vertical heat transport that reveal two distinct scaling behaviors with respect to the Rayleigh and Ekman numbers.

The Gist

Imagine you are trying to understand how heat moves inside a giant, spinning pot of soup. This isn't just any soup; it's the kind found deep inside planets like Earth or Jupiter, or even inside icy moons like Europa. These places are spinning incredibly fast, and they are being heated from the inside (like a radioactive core), not just from the bottom like a pot on a stove.

This paper is a mathematical detective story. The authors are trying to figure out the "rules of the game" for how hot this soup gets and how fast the heat escapes, without having to run a computer simulation that would take a billion years to finish.

Here is the breakdown of their work using simple analogies:

1. The Problem: The "Too Fast, Too Hot" Dilemma

In the real world, the Earth's core spins so fast and is so hot that the fluid inside behaves in a very strange way.

• The Spin: Because the planet spins so fast, the fluid wants to move in tall, thin columns (like spinning dancers holding hands), rather than swirling in big, messy blobs.
• The Heat: The heat comes from everywhere inside the fluid, not just the bottom.
• The Computer Problem: If you try to simulate this on a computer, the math gets so complicated and the numbers get so huge that even the world's fastest supercomputers would crash. It's like trying to count every single grain of sand on a beach while the beach is moving.

2. The Solution: The "Shadow Puppet" Trick

Since they can't simulate the whole messy reality, the authors created a simplified "shadow puppet" version of the physics.

• The Analogy: Imagine a complex 3D sculpture. Instead of trying to measure every curve of the sculpture, you shine a light on it and measure the 2D shadow. The shadow is simpler, but it still tells you the most important things about the shape.
• The Math: They used a technique called "asymptotic reduction." They stripped away the tiny, fast wiggles of the fluid and focused only on the slow, big movements that matter most when the planet is spinning rapidly. This gave them a much simpler set of equations to work with.

3. The Goal: Setting the "Speed Limits"

The authors wanted to find the speed limits (mathematical bounds) for two things:

1. How hot does the soup get on average? (Mean Temperature)
2. How much heat is carried away by the swirling currents? (Convective Heat Flux)

In normal cooking (Rayleigh-Bénard convection), if you know how hot the bottom is, you can guess how hot the top will be. But in this "internally heated" spinning soup, the rules are different. The heat is generated everywhere, so the top and bottom can have very different temperatures depending on how the spin affects the flow.

4. The Method: The "Background Field" Strategy

To find these speed limits, they used a clever mathematical trick called the Auxiliary Functional Method (or the "Background Field" method).

• The Analogy: Imagine you are trying to prove that a car cannot go faster than 100 mph, but you don't know the exact speed of the car at every second. Instead, you imagine a "ghost car" driving alongside it. You design this ghost car to have a specific, steady speed profile.
• The Logic: You then prove that no matter how the real car tries to wiggle or speed up, it can never break the rules set by your ghost car without violating the laws of physics (like conservation of energy).
• The Result: By tweaking the shape of this "ghost car" (the background field), they found the tightest possible speed limits for the real fluid.

5. The Findings: Two Different Worlds

They discovered that the fluid behaves differently depending on how strong the heating is compared to the spinning:

• The "Just-Right" Spin (Moderate Heating): When the heating is moderate, the spin dominates. The heat transport and temperature follow one specific mathematical rule (scaling with the 6/7 power). It's like the fluid is dancing in a very organized, rhythmic line.
• The "Overheated" Spin (Intense Heating): When the heating gets very intense, the fluid starts to break its organized dance. The rules change, and the heat transport follows a different rule (scaling with the 3/5 power). It's like the dancers start to break formation and run in different directions.

6. Why Does This Matter?

This paper provides the first rigorous "rulebook" for this specific type of fluid motion.

• For Planetary Scientists: It helps them understand how Earth's magnetic field is generated or how heat moves inside icy moons, even when we can't simulate those environments directly.
• For Mathematicians: It proves that even in chaotic, spinning systems, there are hard limits to how wild things can get. It's a guarantee that the universe has a "speed limit" even in the most turbulent conditions.

Summary

Think of this paper as creating a traffic law for a planet's core. The authors couldn't watch every car (fluid particle) drive around, so they built a simplified model, invented a "ghost car" to set the rules, and proved exactly how fast the heat can move and how hot the planet can get, depending on how fast it spins. They found that there are two distinct "speed zones" with different rules, ensuring that even in the most chaotic spinning fluids, physics still keeps things in check.

Technical Summary

1. Problem Statement

The paper addresses the long-time behavior of rapidly rotating internally heated convection (IHC), a flow regime relevant to geophysical and astrophysical systems such as planetary cores (driven by radiogenic heating) and subsurface oceans on icy moons (driven by tidal heating).

• The Challenge: In rapidly rotating systems, the Coriolis force dominates buoyancy, organizing flow into columnar structures (Taylor columns). Direct Numerical Simulations (DNS) and laboratory experiments struggle to reach the extreme parameter regimes required (very low Ekman numbers $E \sim 10^{-15}$ and high Rayleigh numbers $R \sim 10^{20}$) to observe geostrophic turbulence.
• The Gap: Standard energy identities fail to capture rotational effects because the Coriolis force does no work. Previous rigorous bounds for IHC existed for non-rotating cases or rotating cases with no-slip boundaries, but the scaling behavior for stress-free boundaries in the rapid rotation limit remained an open question.
• Objective: To derive rigorous, asymptotically reduced bounds on the long-time averages of mean temperature (mixing efficiency) and mean convective heat flux (asymmetry in heat transport) in the limit of infinite Prandtl number ($Pr \to \infty$).

2. Methodology

A. Asymptotically Reduced Model (NH-QG)

Instead of solving the full Navier-Stokes equations, the authors employ the Non-Hydrostatic Quasi-Geostrophic (NH-QG) model.

• Scaling: They utilize an anisotropic non-dimensionalization where the horizontal length scale $L$ scales as $E^{1/3}H$ (where $H$ is the vertical height). This captures the columnar nature of the flow.
• Limit: The analysis is conducted in the limit of infinite Prandtl number ($Pr \to \infty$), reducing the momentum equation to a forced Stokes flow.
• Governing Equations: The system is reduced to a set of coupled PDEs for the vertical velocity $w$, vertical vorticity $\zeta$, horizontal stream function $\psi$, horizontally averaged temperature $\Theta$, and temperature fluctuation $\theta$.

B. Auxiliary Functional Method (Background Field Method)

The authors use the background field method (a variant of the auxiliary functional method) to derive rigorous bounds.

• Decomposition: The flow fields are decomposed into a steady background field (chosen by the authors) and fluctuations.
• Variational Problem: The problem is transformed into an optimization problem. To find a lower bound for mean temperature or an upper bound for heat flux, one must maximize/minimize a functional subject to a spectral constraint.
• Spectral Constraint: This is a condition requiring a specific quadratic form (involving the background field and flow fluctuations) to be non-negative for all admissible wavenumbers.
• Ansatz Construction:
  • For Mean Temperature: A background field $\phi(Z)$ with quadratic boundary layers and a constant interior is constructed.
  • For Mean Heat Flux: A background field with a linear bulk and a quadratic top boundary layer is used.
• Optimization: The authors derive sufficient conditions for the spectral constraint to hold by analyzing large and small wavenumbers separately, reconciling them to find optimal scaling laws for the background field parameters ($A$ and $\delta$) in terms of $R$ and $E$.

3. Key Contributions

1. First Rigorous Bounds for Rapidly Rotating IHC: This is the first study to provide rigorous bounds on mean temperature and convective heat flux for internally heated convection under stress-free boundary conditions in the rapid rotation regime.
2. Derivation of Scaling Laws: The authors identify two distinct scaling regimes based on the reduced Rayleigh number $\tilde{R} = R E^{4/3}$:
   • Regime 1 (Near Onset): $d_0 E^{-4/3} \le R < d_1 E^{-4/3}$.
   • Regime 2 (Strong Rotation): $R \ge d_1 E^{-4/3}$.
3. Optimization of Bounds: The derived bounds are optimized within the chosen methodology, showing that the analytical bounds are within $\sim 2.4\%$ (for temperature) and $\sim 10\%$ (for heat flux) of the numerically computed maximum/minimum of the variational problem.
4. Comparison with Non-Rotating Cases: The paper explicitly contrasts the scaling of these bounds with non-rotating IHC, highlighting how rotation alters the dependence on the Rayleigh number.

4. Key Results

The paper proves two main theorems for the long-time volume averages $\langle \cdot \rangle_{V,t}$:

Theorem 1: Lower Bound on Mean Temperature ($\langle \Theta \rangle_{V,t}$)

The mean temperature (representing thermal dissipation/mixing) is bounded from below.

• Scaling: In the rapid rotation limit ($R \gtrsim E^{-4/3}$), the bound scales as:
  $$\langle \Theta \rangle_{V,t} \gtrsim E^{-1} R^{-3/4}$$
  (Specifically, $c_2 E^{-1} R^{-3/4} - c_3 E^{-2} R^{-3/2}$).
• Significance: Unlike non-rotating IHC where the bound scales as $R^{-5/17}$, the rotating bound decays faster with $R$ but remains higher in absolute magnitude for moderate $R$. This implies that rotation suppresses mixing efficiency less effectively than non-rotating convection at certain scales, or rather, the thermal state remains "hotter" (less mixed) due to the rotational constraint.

Theorem 2: Upper Bound on Mean Convective Heat Flux ($\langle w\theta \rangle_{V,t}$)

The mean convective heat flux (representing the asymmetry between top and bottom heat loss) is bounded from above.

• Scaling: In the rapid rotation limit, the bound scales as:
  $$\langle w\theta \rangle_{V,t} \lesssim E^{4/5} R^{3/5}$$
  (Specifically, $c_4 E^{4/5} R^{3/5} - \dots$).
• Physical Interpretation: As rotation increases (lower $E$), the heat flux asymmetry is constrained. The bound suggests that for very large $R$, the flow approaches a state where heat exits more uniformly, though the bound eventually exceeds the theoretical maximum of $1/2$ (indicating the breakdown of the reduced model or the need for the minimum principle constraint).

Constants and Thresholds

• Linear Onset: The linear stability analysis shows convection begins at $\tilde{R}_c \approx 71.75$.
• Nonlinear Onset: The bounds are valid for $R$ above a threshold $d_0 E^{-4/3}$ (where $d_0 \approx 5.8$ for temperature and $3.5$ for flux), which serves as a lower bound for nonlinear instability.

5. Significance and Implications

• Theoretical Rigor: The work provides a mathematically rigorous framework for understanding rotating IHC, filling a gap left by the inability of current simulations to reach planetary parameter regimes.
• Geophysical Relevance: The derived scaling laws ($R^{-3/4}E^{-1}$ for temperature and $R^{3/5}E^{4/5}$ for flux) offer testable predictions for the thermal structure and heat transport in planetary cores and icy moons, where rotation is dominant.
• Methodological Advancement: The paper demonstrates that the NH-QG model, combined with the background field method and integral estimates (avoiding Green's functions), is a powerful tool for deriving bounds in rotating flows. It successfully handles the spectral constraint without relying on the specific Green's function solutions used in previous Rayleigh-Bénard studies.
• Limitations and Future Work: The authors note that the bounds are currently restricted to infinite Prandtl number. Extending these results to finite $Pr$ (relevant for liquid metals in cores) requires perturbative approaches. Additionally, the validity of the reduced model breaks down when buoyancy dominates rotation (very large $R$), a regime where the bounds may exceed physical limits (e.g., heat flux $> 1/2$).

In summary, this paper establishes the first rigorous constraints on the long-term statistical behavior of rapidly rotating internally heated convection, providing a foundational theoretical tool for interpreting the dynamics of geophysical and astrophysical fluid systems.