Turbulent Heating between 0.2 and 1 au: A Numerical Study

This numerical study demonstrates that MHD turbulent simulations using the expanding box model can reproduce the observed 1/R radial temperature profile of the slow solar wind between 0.2 and 1 AU, provided the turbulence begins with a Mach number of unity and a quasi-2D spectral anisotropy, albeit with limitations imposed by the modest Reynolds numbers achievable at such high Mach speeds.

The Gist

The Big Mystery: Why is the Solar Wind Still Warm?

Imagine the Sun blowing a giant, constant wind made of charged particles (plasma). As this wind travels away from the Sun, it spreads out into the vastness of space.

The Problem:
According to the laws of physics, if you have a gas expanding into a vacuum, it should cool down very quickly. It's like opening a can of compressed air; the air rushes out and gets cold. If the solar wind behaved like a normal gas, the temperature of its protons (the particles) should drop drastically as it travels from the Sun to Earth.

The Reality:
Scientists have measured the solar wind, and it's not cooling down as fast as it should. Between 0.2 and 1 Astronomical Unit (AU)—which is the distance from the Sun to Earth—the wind stays surprisingly warm. It's like the can of compressed air is somehow being reheated while it flies through space.

The Theory:
Scientists suspect that turbulence is the culprit. Think of the solar wind not as a smooth flow, but as a chaotic river with swirling eddies, whirlpools, and choppy waves. As these turbulent swirls crash into each other, they break down into smaller and smaller swirls. Eventually, this chaotic motion turns into heat, warming up the wind.

The Experiment: A Virtual Solar Wind Lab

The authors of this paper wanted to prove that this turbulence theory actually works. They didn't just guess; they built a super-computer simulation (a "numerical study") to act as a virtual laboratory.

They created a digital "box" of solar wind and let it expand from 0.2 AU to 1 AU, mimicking the journey from near the Sun to Earth. They used a special tool called the Expanding Box Model (EBM).

The Analogy:
Imagine you are filming a video of a balloon being inflated.

• Normal filming: The balloon gets bigger, and the camera has to zoom out to keep it in the frame.
• The EBM approach: Instead of zooming out, the camera moves with the balloon, but the film itself stretches sideways as the balloon expands. This allows the scientists to study the tiny, swirling details inside the wind without losing track of the big picture.

The Recipe: What Makes the Wind Stay Warm?

The scientists tried thousands of different "recipes" for the solar wind to see which one produced the right amount of heat. They tweaked four main ingredients:

1. The "Chaos Level" (Mach Number): How wild and fast the turbulence is compared to the speed of sound.
2. The "Stretching Speed" (Expansion Parameter): How fast the wind is expanding relative to how fast the turbulence is swirling.
3. The "Spectrum" (The Size of the Swirls): This is crucial. They had to decide how many different sizes of whirlpools to start with.
4. The "Plasma Beta": The balance between magnetic pressure and gas pressure.

The Big Discovery: It's All About the "Swirls"

The most interesting finding was about the size of the initial turbulence.

The Failed Recipe (Run A):
When they started with a huge variety of swirl sizes (from giant eddies to tiny ones), the simulation went wrong.

• What happened: The tiny, fast swirls dissipated their energy too quickly at the start. It was like pouring a bucket of hot water into a cold lake; the heat was released all at once, then the water cooled down too fast.
• The Result: The wind got too hot too early, then cooled down too much later. It didn't match the real solar wind.

The Winning Recipe (Run E):
They realized that in the real solar wind, the "inertial range" (the zone where energy cascades from big swirls to small ones) is limited.

• The Fix: They started the simulation with a limited range of swirl sizes. They didn't start with tiny, high-energy eddies.
• The Result: By limiting the initial "fuel," the turbulence burned slowly and steadily. This created a steady, gentle heating effect that perfectly matched the observed 1/R temperature profile (where temperature drops slowly as distance increases).

The Metaphor:
Think of the solar wind like a campfire.

• The Bad Way: If you throw a whole log of dry wood (too much small-scale energy) onto the fire at once, it flares up instantly and burns out, leaving you cold later.
• The Good Way: If you carefully feed the fire with just the right amount of kindling and small sticks (limited spectral extent), the fire burns steadily and keeps you warm for a long time.

The Key Takeaways

1. Turbulence is the Heater: The paper confirms that MHD (magnetohydrodynamic) turbulence is indeed capable of heating the solar wind to the levels we observe.
2. The "Goldilocks" Zone: To get the right heating, the turbulence needs to start with a specific "recipe." It needs a high level of chaos (Mach number near 1) but a limited range of swirl sizes. If you have too many tiny, energetic swirls at the start, the system breaks down.
3. The Expansion Factor: The way the wind expands as it travels is just as important as the turbulence itself. The balance between how fast the wind stretches and how fast the turbulence decays determines the final temperature.
4. No Magic Needed: They didn't need to invent new physics. By using standard equations and tweaking the initial conditions to match reality, they successfully reproduced the solar wind's behavior.

Conclusion

In simple terms, this paper solved a puzzle by running a virtual experiment. They found that the solar wind stays warm because its internal "storms" (turbulence) are carefully balanced. If the storms are too chaotic with too many tiny eddies, the wind cools too fast. But if the storms are structured just right, they provide a steady, gentle heat that keeps the solar wind warm all the way from the Sun to Earth.

It's a bit like tuning a radio: if you have too much static (too much small-scale energy), you can't hear the music. But if you tune it just right, you get a clear signal that matches the real world perfectly.

Technical Summary

1. Problem Statement

The radial evolution of proton temperature ($T_p$) in the solar wind remains a significant unsolved problem in heliophysics. Observations indicate that between 0.2 and 1 AU, the temperature decreases much slower than the adiabatic prediction for a spherically expanding flow ($T \propto R^{-4/3}$). Instead, the observed profile follows a slower power law, approximately $T \propto R^{-\xi}$ where $\xi \in [0.8, 1]$.

The prevailing hypothesis is that this "extra" heating is caused by the dissipation of magnetohydrodynamic (MHD) turbulence. However, previous models often relied on simplified nonlinear couplings or free parameters to fit observations. The specific challenge addressed in this paper is to determine, using direct numerical simulations (DNS) with full nonlinear MHD couplings, under what specific initial conditions the turbulent dissipation rate naturally balances the expansion cooling to produce the observed $1/R$ temperature profile.

2. Methodology

Numerical Framework: Expanding Box Model (EBM)

The authors utilize the Expanding Box Model (EBM) to simulate the solar wind evolution from 0.2 AU to 1 AU.

• Coordinate System: The simulation uses a comoving frame that expands with the radial wind. The domain is initially elongated radially (aspect ratio $a_x = 5$) and stretches in the transverse directions as the plasma moves outward, eventually becoming a cube at 1 AU.
• Governing Equations: The study solves the ideal MHD equations modified for expansion, including:
  • Continuity, momentum, induction, and energy equations with expansion terms (e.g., $-2\rho(\epsilon/a)$).
  • Dissipation: Viscous, resistive, and thermal conduction terms are included. Crucially, the transport coefficients ($\mu, \eta, \kappa$) are scaled as $1/a(t)$ to moderate the decrease in the Reynolds number as the turbulent amplitude decays.
• Resolution: Simulations are performed on a $512^3$ grid.

Initial Conditions and Parameters

The study focuses on the slow solar wind, characterized by quasi-2D spectral anisotropy (turbulence primarily perpendicular to the mean magnetic field). The initial state at 0.2 AU is defined by:

• Mach Number ($M$): Ranging from 0.3 to 1.0.
• Expansion Parameter ($\epsilon$): Ranging from 0.12 to 0.4.
• Plasma Beta ($\beta$): Varying from 0.29 to 2.75.
• Spectral Properties: The initial energy spectrum slope ($m$) and extent ($k_{max}$) are varied.
• Key Dimensionless Group: The authors derive a critical parameter $M^2/\epsilon$ based on the balance between turbulent cascade heating and adiabatic cooling. Theoretical analysis suggests a critical value of $M^2/\epsilon \approx 4.4$ is required to sustain a $T \propto 1/R$ profile.

3. Key Contributions and Theoretical Framework

• Derivation of Critical Heating: The authors re-derive the condition for critical heating ($Q_c$) required to maintain a $T \propto 1/R$ profile. They establish a relationship linking the turbulent cascade rate, wind speed, and temperature, leading to the constraint $M^2/\epsilon \approx 4.4$ for cold/slow winds.
• Role of Spectral Extent: A major contribution is the identification that the initial spectral extent (the range of the inertial range) is a critical factor. High Mach numbers ($M \approx 1$) require high viscosity to prevent numerical instability (shocks), which artificially truncates the inertial range.
• Validation of Turbulent Scaling: The study validates that in the correct parameter regime, the ratio of actual dissipation to the Kolmogorov cascade rate ($Q_\nu / Q_{K41}$) converges to $\approx 0.1$, consistent with observational data from Vasquez et al. (2007).

4. Results

The Importance of Spectral Truncation

• Run A (Extended Spectrum): With a large initial spectral extent ($k_{max}=64$) and $M=1$, the simulation exhibits excessive initial heating. The turbulent energy is dissipated too rapidly in the early phase (0.2–0.5 AU), leading to a temperature profile that is too steep initially and then flattens out, failing to match the $1/R$ profile consistently.
• Runs B, C, E (Truncated Spectrum): By reducing the initial spectral extent ($k_{max}=4$) and adjusting the slope, the authors suppress the artificial early-phase dissipation.
  • Result: These runs produce a temperature profile that closely follows the $1/R$ law (specifically $T \propto R^{-0.9}$) across the 0.2–1 AU range.
  • Mechanism: Limiting the spectral extent ensures that the turbulent energy reservoir is not depleted too quickly, allowing for a sustained heating rate that balances adiabatic cooling.

Parameter Sensitivity

• $M^2/\epsilon$: The parameter $M^2/\epsilon$ acts as a primary control knob. Runs with values near the theoretical critical value (e.g., $M^2/\epsilon \approx 5$) yield the best temperature profiles.
• Mach Number: Lowering $M$ reduces the turbulent energy reservoir relative to internal energy, decreasing the heating rate.
• Expansion Parameter ($\epsilon$): A lower $\epsilon$ favors heating initially but can lead to faster energy depletion over the long transport distance due to the increased number of nonlinear turnover times.
• Plasma Beta ($\beta$) and Angle ($\theta$): Variations in $\beta$ (up to 2.75) and the initial magnetic field angle had negligible effects on the resulting temperature profile compared to $M$ and $\epsilon$.

Energy Conservation and Decay

• Decay Rate: The turbulent amplitude ($Z$) decays as $Z \propto R^{-0.6}$, which is slightly faster than the WKB prediction ($R^{-0.5}$) but consistent with the dominance of expansion damping over turbulent dissipation in these regimes.
• Residual Terms: The study analyzed the deviation from turbulent energy conservation ($Q_{NL}$). They found that the primary cause of non-conservation is compressible exchange between turbulent and internal energy, rather than the expansion itself.

5. Significance and Conclusion

This paper provides robust numerical evidence that MHD turbulence alone is sufficient to explain the observed slow cooling of the solar wind proton temperature between 0.2 and 1 AU, without invoking exotic heating mechanisms.

Key Takeaways:

1. Feasibility: Direct numerical simulations with full nonlinear couplings can reproduce the $1/R$ temperature profile.
2. Constraint on Initial Conditions: To achieve this, the turbulence must start with a specific balance of Mach number and expansion rate ($M^2/\epsilon \approx 4.4$) and, crucially, a limited inertial range.
3. Physical Insight: The "limit" on the inertial range is not just a numerical artifact but reflects the physical reality that high-Mach-number turbulence in the inner heliosphere may be constrained by viscosity and shock formation, preventing a full Kolmogorov cascade from developing immediately.
4. Validation: The simulations reproduce the observed ratio of dissipation to cascade rate ($Q_\nu/Q_{K41} \approx 0.1$), confirming that the numerical setup accurately captures the physics of solar wind turbulence.

The study concludes that the observed temperature profile is a natural consequence of the interplay between adiabatic expansion and turbulent dissipation, provided the initial turbulent state satisfies specific spectral and dynamic constraints. Future work will aim to lower the Mach number to achieve larger Reynolds numbers and better resolve the $5/3$ power-law scaling.