Turbulent Heating between 0.2 and 1 au: A Numerical Study

arXiv:2603.01276 · physics.space-ph, astro-ph.SR, physics.plasm-ph · Submitted 2026-03-01 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Today's paper: "Turbulent Heating between 0.2 and 1 au: A Numerical Study".

Jocelyn: This numerical study investigates whether magnetohydrodynamic (MHD) turbulence can account for the observed radial decrease in proton temperature in the solar wind between 0.2 and 1 AU,

Vera: First, who's behind it and why it matters.

Title and authors: Vera: So Jocelyn, this paper by Montagud-Camps et al., "Turbulent Heating between zero point two and one au: A Numerical Study," is really diving into how MHD turbulence might explain that temperature drop we see in the solar wind between zero point two and one AU, which isn't matching the standard adiabatic cooling predictions.

Jocelyn: Exactly, Vera; it tackles that discrepancy head-on by using simulations to see if turbulent dissipation can actually produce the required heating profile, specifically targeting that observed one-over-R dependence.

Subrahmanyan: From a theoretical standpoint, this research is significant because it tries to connect the energy cascade rate within turbulent flows directly to the temperature evolution of the plasma as it moves outward from the Sun.

Vera: That makes sense; essentially, they're testing if turbulence can generate that specific heating profile we're seeing in the data.

Jocelyn: And what they are using is a model called the Expanding Box Model, which modifies standard MHD equations to account for that radial expansion due to the mean wind.

Subrahmanyan: That expansion term is crucial because it allows them to track how spatial derivatives change as the plasma volume stretches over time and distance, which isn't captured in simpler models.

Vera: I was looking at the methodology section, and they define their dissipative terms, like viscosity and resistivity, to decrease as time or distance increases according to a specific relation.

Jocelyn: It’s interesting because it ties the dissipation rates directly into the expansion dynamics of the simulation itself, which is quite complex setup for modeling solar wind conditions.

Subrahmanyan: And then they derive this critical heating requirement, showing that for a one/R temperature profile, you need a specific amount of heating defined by Q c = Q one = (one/two)TU zero/R.

Vera: So they establish this theoretical link between the required heating and the geometric expansion of the wind, which is a solid step in their argument.

Jocelyn: They then compare this critical heating to what we call a generalized Kolmogorov cascade rate, QK41, finding that for cold winds, the ratio of visco-resistive dissipation to this cascade rate stays around zero point one.

Subrahmanyan: That low ratio of approximately zero point one is consistent with what we see in some observational solar wind data, which lends some empirical support to their theoretical framework connecting turbulence to heating in that regime.

Title and authors: Vera: That connection between the simulation parameter and the actual observed dissipation ratio seems like a very important piece of evidence for this paper.

Jocelyn: They run these numerical simulations starting at zero point two AU, using parameters like an initial aspect ratio a x = five and they systematically vary things like the Mach number M and expansion parameter epsilon.

Subrahmanyan: Varying the spectral properties like the initial spectral slope m and extent k max allows them to probe different turbulent regimes, which is necessary because solar wind turbulence isn't uniform.

Vera: The simulations show that when they consider a "strong reduction of the initial spectral inertial range," they can actually get a temperature profile that is close to the one/R law we're trying to explain.

Jocelyn: That suggests that suppressing the small-scale energy content in those turbulent fluctuations is what helps them match these specific observational profiles, which is a key finding for this paper.

Subrahmanyan: The conclusion they draw is that the origin of this extra heating might come from developing a turbulent regime that generates substantial heating, and their simulations confirm radial temperature profiles close to one/R on average.

Vera: It sounds like they’re saying that the combination of adiabatic decrease and this specific type of turbulent dissipation under certain initial conditions can produce the observed temperature structure.

Jocelyn: They do have a point about distinguishing between two phases in energy evolution: an early phase with rapid and strong dissipation, followed by a longer phase where residual decay leads to a profile similar to the average one/R zero point nine profile measured by Totten et al. (one thousand nine hundred ninety-five).

Subrahmanyan: That distinction between the early, strong dissipation phase and the later, more gradual decay phase is important for understanding how energy is ultimately converted into heat in these turbulent solar wind environments.

Vera: So, to summarize this paper on "Turbulent Heating between zero point two and one au: A Numerical Study," it essentially shows that MHD turbulence can produce a one/R temperature decrease if the initial conditions are set correctly, particularly by suppressing small-scale energy.

Jocelyn: And they found that the parameters regulating this heating rate are primarily determined by the combination of the rms Mach number and the expansion parameter epsilon, specifically M two/epsilon.

Title and authors: Subrahmanyan: That finding, with other parameters like plasma beta having only a minor effect up to now, is quite telling for how robust this turbulent mechanism is in shaping these temperature profiles.

Vera: It’s encouraging to see that the Mach number and expansion parameter are the main drivers here, rather than needing every single plasma parameter to be perfectly tuned.

Jocelyn: Speaking of tuning parameters, they did explore how changing the mean magnetic field amplitude B zero and even the initial spectral slope m impacts these results, which shows how sensitive this model is to those initial conditions.

Subrahmanyan: This sensitivity suggests that future observational constraints will be vital in narrowing down the range of physical parameters that govern solar wind heating processes.

Vera: So, looking ahead at what this implies for our understanding of solar wind dynamics, it seems we need models that can incorporate this turbulent dissipation mechanism to accurately predict temperature evolution across these distances.

Jocelyn: It opens up new avenues for AI systems to perform parameter inference, allowing them to predict the necessary dissipation rate based on observed temperature profiles and initial wind conditions.

Subrahmanyan: And for the broader cosmic picture, it reinforces the idea that turbulence is a primary driver in energy conversion processes within magnetized plasmas across various astrophysical environments.

Vera: It’s exciting to think about how this could feed into models of plasma physics in other stellar and planetary atmospheres, given the applicability to these solar wind conditions.

Jocelyn: I think the real impact here is on how we interpret future missions, as they will need to look for signatures that point toward these specific turbulent regimes identified in this numerical work.

Subrahmanyan: Ultimately, understanding this mechanism helps us map out how kinetic energy from large-scale flows gets distributed across different scales and converted into thermal energy in the heliosphere.

Vera: So, to wrap up on "Turbulent Heating between zero point two and one au: A Numerical Study," it’s a numerical study demonstrating that specific MHD turbulence can generate the one/R temperature profile through careful tuning of initial spectral conditions.

Jocelyn: And they point toward the need for more observational data that can constrain these initial parameters to fully validate this turbulent heating hypothesis.

Subrahmanyan: It's a solid piece of work that provides a computational pathway to link turbulence directly to macroscopic plasma properties in our solar system.

The paper's summary: Vera: So, to recap, this numerical study basically shows that if you have certain initial conditions for turbulence in the solar wind, it can generate that specific temperature drop we see between zero point two and one AU by having turbulent dissipation play a big role alongside the standard cooling.

Jocelyn: That makes sense from an observational standpoint; it’s taking those complex simulations and mapping them onto what we actually measure in the solar wind, showing a way to get closer to matching those temperature gradients than just relying on adiabatic models alone.

Subrahmanyan: From a theoretical astrophysics angle, this work suggests that the energy cascade within MHD turbulence isn't just about stirring up the plasma; it’s actively converting kinetic energy into thermal energy in a spatially dependent manner across the solar wind.

Vera: Exactly; they found that suppressing small-scale turbulent activity is actually what helps them arrive at those observed one/R profiles, which is a really specific detail.

Jocelyn: It’s exciting because it gives us a concrete physical mechanism to investigate when we look at temperature data from missions like ACE or Parker Solar Probe, suggesting that the "extra" heating isn't just noise; it might be a structured process.

Subrahmanyan: And the implication here is that we need to refine our models of how energy is distributed across different scales in magnetized flows, moving beyond simple fluid descriptions to something more detailed regarding turbulent dissipation rates.

Vera: I think the part about distinguishing between different phases of energy evolution—the early strong dissipation versus the later residual decay—is super important for understanding the entire lifecycle of these turbulent structures.

Jocelyn: And that phase distinction hints at how long-lasting these heating effects might be, which is something we need to look for in our next round of solar wind observations.

Subrahmanyan: Plus, their finding that the parameters regulating this heating are primarily M two/epsilon gives us a clear focus for future theoretical modeling; we know where to concentrate our efforts to see if that relationship holds across different environments.

Vera: It really puts things into perspective, showing how finely tuned those initial conditions have to be for the turbulence to produce the specific profile we are trying to understand.

Jocelyn: So, it’s not just about finding *a* heating mechanism; it’s about figuring out *which* turbulent regime produces the exact shape of the temperature drop we see as we get further from the Sun.

Subrahmanyan: And that connects back to our background on pickup He+ tori; if turbulence is key here, it might be a common feature in other complex astrophysical plasmas where energy is being transferred across scales.

Vera: It sounds like this paper provides a vital computational tool for us to bridge the gap between microscopic turbulent physics and macroscopic solar wind observations.

The paper's improvements: Vera: So, to wrap up on that discussion, this paper outlines several ways future work can push this research further regarding turbulent heating in the solar wind between zero point two and one AU.

Jocelyn: That's true; they point out that while their current model is quite robust for certain conditions, there are still parameters like plasma beta that could be explored more deeply to see how much influence they really have on the final temperature profile shape.

Subrahmanyan: They also suggest extending the simulation domain further out, perhaps testing these turbulent heating mechanisms in regions beyond one AU to see if the physics scales up in a predictable way.

Vera: That makes sense; expanding the spatial extent of those simulations would give us a much better picture of how this turbulence behaves across different radial distances.

Jocelyn: And they mention focusing more on the initial spectral slope, m, to see if that’s truly as critical as they initially suggested, which is something we can test with more detailed observational constraints from future missions.

Subrahmanyan: From a theoretical standpoint, exploring those parameter spaces will help us build a more comprehensive framework for how kinetic energy gets converted into thermal energy in these expanding solar wind environments.

Vera: I think the implication here is that the next step is to move from finding *if* turbulence works to understanding *exactly* what combination of initial conditions produces which specific temperature profile, like the one we are trying to match.

Jocelyn: It gives us a clear roadmap for experimental and observational teams; they know exactly what parameters they need to target when looking at temperature profiles from instruments like STEP or future spacecraft.

Subrahmanyan: And if this approach proves successful, it could help us explain heating mechanisms in other astrophysical plasmas where the expansion dynamics are similarly complex.

Vera: It’s really exciting to think about how these detailed simulation results can inform our understanding of energy transport in various stellar and planetary atmospheres too.

Jocelyn: So, the main thing is that they're not just stopping at a single solution; they're providing a framework for exploring the whole parameter space of turbulent heating.

Subrahmanyan: And I think the future work on extending the domain is where we could see if this mechanism has broader relevance across different solar wind regimes.

Vera: That’s right; it shows that this area of research isn't settled, and there's a lot more ground to cover with these suggested improvements.

Conclusion: Vera: So, to wrap up on this discussion of "Turbulent Heating between zero point two and one au: A Numerical Study," this paper demonstrates that MHD turbulence can produce the observed temperature decrease by carefully tuning initial spectral conditions to suppress small-scale energy.

Jocelyn: It really shows how crucial those initial setup parameters are, and it gives us a tangible link between simulation inputs and the resulting temperature structure we see in space.

Subrahmanyan: The overall impact is that it strengthens the theoretical argument that turbulent dissipation is a significant contributor to heating profiles in solar wind environments, which has wide implications for understanding energy conversion across astrophysical scales.

Vera: I think the most exciting part is seeing how the authors successfully modeled those different phases of energy evolution, showing that both rapid early dissipation and later residual decay can contribute to the final profile.

Jocelyn: That phase distinction is something we need to keep in mind when we analyze future observational data from missions that are looking at temperature gradients over long distances.

Subrahmanyan: And for the cosmic picture, it reinforces the idea that these processes are fundamental ways energy gets distributed within magnetized plasmas as they expand into the heliosphere.

Vera: It’s a solid piece of work, and it gives us a powerful numerical tool to test our ideas about how kinetic energy is converted into heat in these specific solar wind conditions.

Jocelyn: I think we should definitely keep an eye out for how this specific turbulent heating mechanism relates to the observations from Solar Orbiter's STEP instrument.

Subrahmanyan: Indeed, and understanding this paper’s results will help us better interpret those complex structures we see in the solar wind environment.

LPP, Ecole Polytechnique, CNRS · Universita di Firenze, Dipartimento di Fisica e Astronomia

physics.space-ph, astro-ph.SR, physics.plasm-ph

Submitted: 2026-03-01

Updated: 2026-03-01

Journal ref: ApJ 853 153 (2018)

DOI: 10.3847/1538-4357/aaa1ea

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 62/100

The gist: This numerical study investigates whether magnetohydrodynamic (MHD) turbulence can account for the observed radial decrease in proton temperature in the solar wind between 0.2 and 1 AU, which

Key concepts

Expanding Box Model (EBM)
This is a set of modified MHD equations used to simulate plasma evolution in the solar wind. It accounts for the systematic expansion of the plasma volume due to its mean radial flow, allowing researchers to model how turbulence and dissipation affect temperature over distance.
Critical Heating ($Q_c$)
This is the minimum heating rate required in a turbulent system to produce a specific temperature decrease profile. The study found that for a 1/R temperature profile, this critical heating is directly proportional to the initial thermal energy and inversely proportional to the distance from the Sun.
Kolmogorov Cascade Rate ($Q_{K41}$)
This represents the rate at which energy flows down through different scales in a turbulent fluid, similar to how energy cascades in water turbulence. The ratio of dissipation to this cascade rate helps determine if turbulence is effectively heating the plasma according to observed solar wind data.

Terminology

Summary

This numerical study investigates whether magnetohydrodynamic (MHD) turbulence can account for the observed radial decrease in proton temperature in the solar wind between 0.2 and 1 AU, which deviates from adiabatic predictions. The research aims to determine if turbulent dissipation can generate the required heating profile, specifically a 1/R dependence, by simulating slow solar wind conditions characterized by quasi-2D spectral anisotropy.

Model and Governing Equations

The study utilizes the Expanding Box Model (EBM) equations, which are modifications of standard MHD equations incorporating expansion due to the mean radial wind. The plasma evolution is governed by these equations, which include terms accounting for the systematic velocity field perpendicular to the radial direction. Key aspects of this model include:

  1. The domain expansion is defined by a normalized heliospheric distance parameter, where the expansion rate is given by an expression involving the initial nonlinear time and turnover time:

  2. The spatial derivatives are modified to account for increasing lateral stretching of the plasma volume with time/distance, using coordinates comobile with this transverse expansion:

  3. Dissipative terms are defined as viscous and resistive terms, where viscosity, resistivity, and conductivity decrease as time/distance increases according to the relation:

  4. The turbulent heating rate is derived from these dissipative terms and is expressed by the term: Qν = µ(˜ω2 + 4/3 (∇ · ˜ u)2) + ηJ˜2.

Critical Heating and Temperature Profile Derivation

The paper establishes a theoretical link between the energy cascade rate and the temperature profile. It derives an equation expressing the critical heating required to produce a temperature decrease following a power law, specifically:

Qα = (1/2)TU0/R

For the specific case of a 1/R profile, which corresponds to setting α = 1, the critical heating is defined as Qc:

Qc = Q1 = (1/2)TU0/R

The paper then relates this critical heating to the turbulent cascade rate by defining a generalized Kolmogorov cascade rate, QK41. It is found that for cold winds, the ratio of visco-resistive dissipation to the Kolmogorov energy cascade rate (RV = Qν/QK41) holds at approximately 0.1, which is consistent with observations in solar wind data.

Numerical Simulations and Parameter Variation

The research employs a numerical box with a resolution of Nx = Ny = Nz = 512, starting simulations at 0.2 AU with an initial aspect ratio ax = 5 along the radial direction. The simulations are conducted across various initial conditions to test different regimes:

  1. The Mach number (M) and expansion parameter (ǫ) are varied, with typical values chosen as M=1 and ǫ=0.2 for representative runs.

  2. The simulation parameters are systematically explored by varying the Mach number, expansion parameter, plasma beta (β), mean magnetic field amplitude (B0), and the spectral properties such as the initial spectral slope (m) and extent (kmax).

  3. The results show that when considering a strong reduction of the initial spectral inertial range, a temperature profile close to a 1/R law is obtained, suggesting that suppressing small-scale energy content is crucial for matching observations.

Key Findings on Turbulence and Heating

The simulations provide several key insights into the physics of turbulent heating:

The origin of this extra heating may be attributed to the development of a turbulent regime in the wind which generates a substantial heating.

We find radial temperature profiles close to 1/R in average.

The analysis distinguishes between different phases of energy evolution:

  1. An early phase characterized by rapid and strong dissipation where turbulent dissipation dominates expansion decay.

  2. A longer-lasting phase where the residual decay is smaller, leading to a temperature decrease that is not far from the average 1/R0.9 profile measured by Totten et al. (1995).

Furthermore, the study concludes that the parameters regulating the heating rate are the rms Mach number and the expansion parameter ǫ, combined as M2/ǫ, while other parameters like plasma beta and mean field angle have a minor effect up to now. The results support MHD turbulence as a mechanism driving temperature profiles decreasing significantly more slowly than adiabatic predictions in the 0.2 < R < 1 AU range.

Conclusion and Future Work

In summary, the numerical results demonstrate that radial temperature profiles as 1/R can result from the combination of adiabatic decrease and turbulent dissipation under specific initial conditions (e.g., M=1, ǫ=0.2). The study suggests that the observed profile is achieved when considering a strong reduction of the initial spectral inertial range, which suppresses excessive heating from small scales.

Improvements for AI systems

As a fastidious and diligent researcher, I have analyzed this paper, TURBULENT HEATING BETWEEN 0.2 AND 1 AU: A NUMERICAL STUDY, which investigates how Magnetohydrodynamic (MHD) turbulence drives the observed radial temperature profile of the solar wind.

Here are specific improvements to AI systems that can be derived from this scientific research:

) Specific Improvements for AI Systems Based on This Paper:

  1. AI System capable of performing Turbulent Heating Parameter Inference for Solar Wind Models (e.g., Magnetohydrodynamic simulations).

  2. AI System capable of predicting the required turbulent dissipation rate to match observed temperature profiles, accounting for Mach number, expansion parameter, and spectral properties.

  3. AI System capable of distinguishing between different turbulent regimes (e.g., Run A vs. Runs B/C/E) based on the resulting temperature profile shape (power-law index) and energy budget dynamics (e.g., Qν/Qc ratio).

) What the Improved AI System Can Do:

  1. AI can take observational data of proton temperature profiles in the solar wind (from instruments like ACE or Parker Solar Probe) and, using the derived critical heating formula (Eq. 34), infer whether the observed profile is consistent with a specific turbulent regime, such as one characterized by a near-critical heating ratio of 1.0 (i.e., Run B/C/E).

  2. AI can analyze high-fidelity MHD simulation outputs and calculate the Kolmogorov rate components (Eq. 35) versus the visco-resistive dissipation rate (Eq. 27). It can then assess if the turbulence exhibits the expected low ratio of approximately 0.1 (Run B/C/E) that is characteristic of cold solar wind, effectively validating or rejecting the underlying turbulent physics model used in that simulation run.

  3. AI can serve as a rapid parameter optimizer for MHD simulations: Given desired output profiles (e.g., a specific 1/R power law index between 4/3 and 1), the AI can iteratively adjust initial conditions—specifically the spectral extent (kmax) and slope (m)—to find the optimal set of parameters that yields this result, bypassing lengthy manual exploration of parameter space.

  4. AI can predict how changes in initial physical parameters—such as increasing Mach number or changing plasma beta—will affect the resultant temperature profile shape, providing a quantitative prediction for solar wind properties under different environmental conditions (e.g., high-beta vs. low-beta environments).

Abstract

The heating of the solar wind is a key to understand its dynamics and acceleration process. The observed radial decrease of proton temperature in the solar wind is slow compared to the adiabatic prediction and it is thought to be caused by turbulent dissipation. To generate the observed 1/R decrease, the dissipation rate has to reach a specific level which varies in turn with temperature, wind speed, and heliocentric distance. We want to prove that MHD turbulent simulations can lead to the 1/R profile. We consider here the slow solar wind, characterized by a quasi-2D spectral anisotropy. We use the EBM (expanding box model) equations, which incorporate into 3D MHD equations the expansion due to the mean radial wind, allowing to follow the plasma evolution between 0.2 and 1 AU. We vary the initial parameters which are: Mach number, expansion parameter, plasma beta, and properties of the energy spectrum as the spectral range and slope. Assuming turbulence starts at 0.2 AU with a Mach number equal to unity, with a 3D spectrum mainly perpendicular to the mean field, we find radial temperature profiles close to 1/R in average. This is done at the price of limiting the initial spectral extent, corresponding to the small number of modes in the inertial range available, due to the modest Reynolds number reachable with high Mach numbers.

Related papers