Prandtl number dependence of rotating internally heated… — Plain-Language Explanation

✨

This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Cooking a Pot of Soup in a Spinning Pot

Imagine you have a giant, deep pot of soup. Instead of heating it from the bottom (like a stove), you have magic heating elements scattered throughout the entire soup, warming it from the inside out. This is what scientists call Internally Heated Convection (IHC).

In nature, this happens inside planets (like Earth's core) and stars. The soup wants to rise because it's hot, but it's also being spun around (like a planet rotating).

The researchers wanted to know: How does the "thickness" or "stickiness" of the fluid change the way heat moves?

In physics, this "stickiness" is called the Prandtl number ($Pr$).

Low $Pr$ (Runny Fluid): Think of water or liquid metal. It flows easily, and heat spreads out very fast.
High $Pr$ (Thick Fluid): Think of honey or molasses. It flows slowly, and heat gets stuck in place.

The team used powerful computers to simulate this soup pot, spinning it at different speeds and changing the fluid from runny water to thick honey. Here is what they found.

1. The Non-Rotating Case: The "Party" vs. The "Nap"

First, they looked at the pot without spinning it.

The Setup: Because the soup is heated from the inside, the top layer gets hot and wants to rise (unstable), while the bottom layer gets hot but is sitting on a cold plate, so it stays put (stable).
The Runny Fluid (Low $Pr$): Imagine a rowdy party. The turbulence from the top is so energetic that it stirs up the bottom layer, too. Even though the bottom should be calm, the chaotic motion from above "recovers" the symmetry, mixing everything up. The heat escapes efficiently through the bottom.
The Thick Fluid (High $Pr$): Imagine a library where everyone is whispering. The thick fluid is so sticky that the turbulence from the top can't reach the bottom. The bottom layer becomes a "Dead Zone." It sits there, perfectly calm and stratified, while all the action happens at the top. The heat struggles to get to the bottom.

The Surprise: Even though the bottom of the pot was either a "party" (runny) or a "dead zone" (thick), the average temperature of the whole pot stayed almost the same.

Analogy: It's like a classroom. Whether the back row is rowdy or silent, the average noise level of the whole room is mostly determined by the teacher at the front (the top boundary layer). The top layer controls the global temperature, regardless of what's happening at the bottom.

2. The Rotating Case: The "Spinning Top" Effect

Now, they started spinning the pot (simulating a rotating planet). This adds a force called the Coriolis force, which tries to organize the flow into vertical columns (like tornadoes).

The Runny Fluid (Low $Pr$): When they spin the runny soup, the heat spreads out so fast that the spinning force can't organize it into neat columns. The "magic" of rotation doesn't really help move heat faster. The efficiency stays the same.
The Thick Fluid (High $Pr$): When they spin the thick honey, something magical happens. The spinning creates Ekman pumping. Imagine the spinning creates a vacuum cleaner effect at the top and bottom, sucking fluid up and down in neat, organized columns.
- Result: This dramatically improves the efficiency of heat transport, but only if the fluid is thick enough ( $Pr \ge 1$ ). If the fluid is too runny, the heat diffuses away before the spinning can organize it.

The Takeaway: Rotation makes the heat transport better, but only if the fluid is "thick" enough to hold onto the organized structure.

3. The "Dead Zone" and Planetary Interiors

The most fascinating discovery is about the bottom of the pot.

In high-viscosity fluids (like the Earth's mantle or the deep interior of a gas giant), the bottom layer becomes a stagnant, dead zone. The fluid there is so thick and stable that the turbulence from above simply cannot penetrate it.

Real-world implication: If you are trying to model how Earth's core cools down, you can't just look at the top. You have to realize that the bottom might be "frozen" in a state of calm, effectively hiding from the mixing process. The "effective depth" of the active convection changes depending on how thick the fluid is.

Summary in a Nutshell

Top vs. Bottom: The top of the fluid layer controls the overall temperature, no matter how thick or thin the fluid is.
The Bottom Layer: The bottom layer is the sensitive one.
- Thin fluid: The bottom gets stirred up and active.
- Thick fluid: The bottom goes to sleep (a "dead zone").
Spinning Helps (Sometimes): Spinning the pot helps move heat faster, but only if the fluid is thick enough to let the spinning organize the flow. If the fluid is too runny, spinning does nothing to improve efficiency.

Why does this matter?
This helps scientists understand how planets and stars cool down over billions of years. It tells us that the "stickiness" of the material inside a planet determines whether the deep interior is a churning mix or a quiet, stagnant layer.

1. Problem Statement

The study investigates the influence of the Prandtl number ($Pr$) on Internally Heated Convection (IHC) in both non-rotating and rotating regimes. Unlike the classic Rayleigh-Bénard (RB) problem where heat is driven by boundary temperature gradients, IHC involves volumetric heating within the fluid, bounded by isothermal plates. This configuration creates a unique flow structure: an unstably stratified thermal boundary layer at the top and a stably stratified layer at the bottom.

The authors aim to resolve how the ratio of viscous to thermal diffusion ($Pr$) affects:

Global heat transport efficiency and mean temperature.
The morphology of the flow, particularly the interaction between the stable bottom layer and turbulent bulk.
The interplay between rotation (Coriolis force) and $Pr$, specifically whether the "rotational enhancement" of heat transport observed in RB convection (via Ekman pumping) persists in IHC across the full range of $Pr$ (from $0.1$ to $100$).

2. Methodology

Numerical Approach: The authors performed three-dimensional Direct Numerical Simulations (DNS) using the open-source code AFiD.
Governing Equations: The system is modeled using the Boussinesq approximation for the Navier-Stokes equations in a rotating frame, coupled with an advection-diffusion equation for temperature with a uniform volumetric heat source.
Parameter Space:
- Prandtl Number ($Pr$): Varied from 0.1 to 100 (covering low-$Pr$ liquid metals to high-$Pr$ planetary mantles).
- Rayleigh Number ($Ra$): Ranged from $3.16 \times 10^5$ to $10^{10}$ .
- Ekman Number ( $E$ ): Varied from $10^{-6}$ (strong rotation) to $\infty$ (non-rotating).
- Domain: Periodic horizontal boundaries with aspect ratios ( $\Gamma$ ) of 1 and 2 to ensure domain size independence.
Resolution: Simulations were resolved to satisfy exact global balance relations for viscous and thermal dissipation, ensuring accuracy within 1%.

3. Key Contributions

First DNS of Rotating IHC across a Wide $Pr$ Range: This study extends previous work (Ostilla-M´onico & Arslan, 2025) to include low-$Pr$ fluids ($Pr < 1$) and high-$Pr$ fluids ($Pr > 10$) in rotating IHC, a regime previously unexplored computationally.
Decoupling of Global and Local Dynamics: The paper demonstrates that while global mean temperature is robust to $Pr$, local flow structures and heat flux asymmetry are highly sensitive to it.
Mechanism of Rotational Enhancement: The study clarifies that the rotational enhancement of global heat transport (analogous to the Nusselt number in RB) is conditional on $Pr$, occurring only when $Pr \geq 1$ .

4. Key Results

A. Non-Rotating Regime ( $E = \infty$ )

Global Mean Temperature ( $\langle T \rangle$ ): Remarkably insensitive to $Pr$. It follows a scaling law of $\langle T \rangle \sim Ra^{-0.2}$ across all $Pr$. This is because $\langle T \rangle$ is primarily controlled by the dynamics of the unstably stratified top boundary layer, which remains largely unaffected by $Pr$.
Flow Morphology & Symmetry Recovery:
- Low $Pr $($ 0.1 - 1$): High thermal diffusivity allows turbulent stirring from the bulk to penetrate the stable bottom layer. This creates a "symmetry recovery" where the bottom layer becomes active, increasing vertical convective heat flux ( $\langle wT \rangle$ ) and reducing the heat fraction leaving the bottom ( $F_B$ ).
- High $Pr $($ 10 - 100$): High viscosity and low thermal diffusivity suppress turbulence in the stable bottom layer. This leads to a "dead zone" where the fluid near the bottom plate is nearly quiescent. Heat transport becomes dominated by plumes originating from the top.
Dissipation: Viscous dissipation shifts from being distributed (low $Pr$) to being dominated by the bulk and top layer (high $Pr$).

B. Rotating Regime ( $E < \infty$ )

Vertical Convective Flux ( $\langle wT \rangle$ ): Rotation enhances $\langle wT \rangle$ $⟨ w T ⟩$ across all $Pr$ values.
- At low $Pr$, the enhancement is driven by increased asymmetry between the top and bottom layers due to thermal diffusion maintaining the stable bottom.
- At high $Pr$, the enhancement is driven by Ekman pumping, which organizes flow into columnar structures that transport heat more efficiently.
Global Heat Transport Efficiency ( $\langle T \rangle$ ):
- $Pr \geq 1$ : Rotation significantly improves global cooling efficiency (reduces $\langle T \rangle$ ). This is attributed to the emergence of Ekman pumping, which effectively transports heat.
- $Pr < 1$: Rotation does not improve global cooling efficiency. High thermal diffusivity allows heat to escape the boundary layers before it can be vertically transported by columnar flows, rendering Ekman pumping ineffective for global heat removal.
Flow Structure:
- At low $Pr$, rotation has a minimal effect on the turbulent bottom layer.
- At high $Pr$, rotation induces a transition to geostrophic turbulence with vertically aligned Taylor columns.
- A bulk temperature gradient emerges once the inverse Rossby number ($1/Ro$) exceeds 1, a feature not present in non-rotating IHC.

C. Dissipation and Boundary Layers

Thermal Dissipation: In the rotating regime, the dominant source of thermal dissipation shifts from the top boundary layer to the bulk as rotation increases ( $E < 10^{-5}$ ), indicating a fundamental change in the heat transport bottleneck.
Boundary Layers: The top thermal boundary layer thickness remains relatively constant until strong rotation ( $1/Ro \approx 1$ ) causes it to thicken. The bottom thermal boundary layer thickens with increasing $Pr$ and rotation.

5. Significance

Planetary and Stellar Interiors: The findings are critical for modeling convection in planetary mantles (high $Pr$) and liquid metal cores (low $Pr$). The study suggests that in high-$Pr$ environments (like Earth's mantle), the lower portion of the convective layer may be dynamically "masked" or inactive due to stable stratification, effectively reducing the mixing depth.
Distinction from Rayleigh-Bénard: The paper highlights that while IHC shares scaling similarities with RB convection at the top boundary, the internal stratification creates a unique sensitivity to $Pr$ that is absent in boundary-driven systems.
Rotational Enhancement Limits: The results establish a clear threshold ( $Pr \geq 1$ ) for the effectiveness of rotation in enhancing global heat transport in internally heated systems, challenging the assumption that rotational effects are universally beneficial for cooling efficiency.

In summary, the paper establishes that while the global mean temperature in IHC is robust and controlled by the top boundary, the Prandtl number dictates the internal flow structure, the activity of the stable bottom layer, and the efficacy of rotation in enhancing global heat transport.

Prandtl number dependence of rotating internally heated convection