Phase-space averaging for stellar convection
- 1Univ Toulouse, CNRS, CNES, IRAP
Abstract
Context. Surface convection in cool stars sets the entropy jump between the atmosphere and the deep convective envelope and therefore affects the radius and structure of one-dimensional stellar models. In standard local treatments, this entropy jump is controlled through an effective convective efficiency. Aims. We have developed a phase-space closure that relates the mean stratification to the local distribution of entropy and velocity among convective elements. Our aim is to describe how the surface entropy transition emerges from the statistics of these elements rather than prescribing it directly. Methods. We first constructed a maximum-entropy model for the convective bulk, assuming that the realised states are bounded by a deep entropy level S. This prediction was then confronted with three-dimensional radiation-hydrodynamics simulations from the M3DIS grid. The comparison motivated an extension to a two-population closure in which an optimal bulk component coexists with a radiatively cooled plume component. Results. The maximum-entropy construction predicts a one-sided exponential entropy distribution, and its low-entropy tail is recovered in the deep convective part of the simulations. Closer to the surface, the distribution is reorganised by radiative cooling and separates into two thermodynamic components. Using a radiative entropy branch, an adiabatically descending cooled branch, and a population fraction controlled by optical depth, the two-population closure reproduces the mean entropy stratification of the four reference simulations over most of the displayed range without adjusting the entropy profile itself. Conclusions. The surface transition can therefore be interpreted as a continuous population decomposition, rather than as a sharp boundary between radiative and convective layers. This provides a natural route towards one-dimensional stellar models in which the surface entropy jump is determined from the statistical organisation of the flow instead of being set by a calibrated mixing-length efficiency.
Versions
Related
- Phase-space averaging for stellar convection I. Liouvillian dynamics
- Phase-space averaging for stellar convection II. Maximum-entropy closures for mixed radiative-convective envelopes
- LACHESIS: Robust stellar parameters through Bayesian model averaging of isochrone grids
- Modelling stellar convective transport with plumes – II. Transport properties of locally and non-locally driven convection
- LACHESIS: Robust stellar parameters through Bayesian model averaging of isochrone grids
BibTeX
@article{p2026phasespace,
title = {{Phase-space averaging for stellar convection}},
author = {Houdayer P. S. and Rieutord M.},
journal = {Astronomy \& Astrophysics},
year = {2026},
month = oct,
doi = {10.1051/0004-6361/202661149},
url = {https://doi.org/10.1051/0004-6361/202661149},
note = {CC BY},
}
RIS
TY - JOUR
TI - Phase-space averaging for stellar convection
AU - Houdayer P. S.
AU - Rieutord M.
PY - 2026
DA - 2026/10/23
JO - Astronomy & Astrophysics
DO - 10.1051/0004-6361/202661149
UR - https://doi.org/10.1051/0004-6361/202661149
AB - Context. Surface convection in cool stars sets the entropy jump between the atmosphere and the deep convective envelope and therefore affects the radius and structure of one-dimensional stellar models. In standard local treatments, this entropy jump is controlled through an effective convective efficiency. Aims. We have developed a phase-space closure that relates the mean stratification to the local distribution of entropy and velocity among convective elements. Our aim is to describe how the surface entropy transition emerges from the statistics of these elements rather than prescribing it directly. Methods. We first constructed a maximum-entropy model for the convective bulk, assuming that the realised states are bounded by a deep entropy level S. This prediction was then confronted with three-dimensional radiation-hydrodynamics simulations from the M3DIS grid. The comparison motivated an extension to a two-population closure in which an optimal bulk component coexists with a radiatively cooled plume component. Results. The maximum-entropy construction predicts a one-sided exponential entropy distribution, and its low-entropy tail is recovered in the deep convective part of the simulations. Closer to the surface, the distribution is reorganised by radiative cooling and separates into two thermodynamic components. Using a radiative entropy branch, an adiabatically descending cooled branch, and a population fraction controlled by optical depth, the two-population closure reproduces the mean entropy stratification of the four reference simulations over most of the displayed range without adjusting the entropy profile itself. Conclusions. The surface transition can therefore be interpreted as a continuous population decomposition, rather than as a sharp boundary between radiative and convective layers. This provides a natural route towards one-dimensional stellar models in which the surface entropy jump is determined from the statistical organisation of the flow instead of being set by a calibrated mixing-length efficiency.
SN - 0004-6361
SN - 1432-0746
ER -