Influence of numerical domain design in simulations of solar coronal heating with the MAGEC code
Abstract
Numerical simulations are essential for studying the solar atmosphere, especially in the context of the coronal heating problem. Understanding how numerical factors and boundary conditions (BC) influence the energetic structure of these models is crucial for unraveling coronal heating mechanisms. We investigate how the computational domain dimensions, upper BC, and vertical grid resolution influence the relaxed thermal structure, energetics, and long-term stability of 2D solar atmosphere simulations with a vertical magnetic field. Using the MAGEC code, we conducted parametric studies varying domain width and height, vertical grid stretching, and top BC, including hot-plate, closed, unrestricted temperature, background heating, and non-reflecting characteristic configurations. We evaluated average temperature profiles for parameter dependence and temporal stability, and analyzed the energetic budget as a function of height, isolating contributions from radiative losses, thermal conduction, pressure work, and internal energy flux. Models with larger vertical domains show higher average coronal temperatures and steeper transition region gradients. Among the less artificial BC, non-reflecting characteristics most effectively delay coronal cooling. We determined the threshold height at which a dynamic top BC can mimic fixed-temperature setups, and established the maximum allowable grid stretching that preserves physical accuracy while reducing computational costs. The energetic terms change significantly with the chosen BC, with specific mechanisms acting as cooling agents in different atmospheric layers. Our results demonstrate that domain height, vertical grid design, and BC critically dictate the relaxed thermal structure and energetic balance of the simulated solar atmosphere.
Versions
Related
- Understanding how Coronal Heating and Magnetic Field Scaling Affect Global Properties of the Low and Middle Corona
- Enhanced Coronal Heating Detected across the Prominence–Corona Boundary Using DKIST/Cryo-NIRSP Spectropolarimetric Observations
- Photospheric Kelvin–Helmholtz Vortices as Possible Drivers of Coronal Heating: Implications of the DKIST Observations
- Revisiting the Role of Global Magnetic Variables in Coronal Heating of Solar Active Regions
- Heating of Solar Coronal Loops by Converging Unipolar Magnetic Flux: Evidence from SDO Observations and a Radiative MHD Simulation
BibTeX
@misc{navarro2026influence,
title = {{Influence of numerical domain design in simulations of solar coronal heating with the MAGEC code}},
author = {Anamaría Navarro and E. Khomenko and T. Felipe and S. J. González Manrique},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.22012},
eprint = {2610.22012},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.22012},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Influence of numerical domain design in simulations of solar coronal heating with the MAGEC code
AU - Anamaría Navarro
AU - E. Khomenko
AU - T. Felipe
AU - S. J. González Manrique
PY - 2026
DA - 2026/10/15
PB - arXiv
DO - 10.48550/arxiv.2610.22012
UR - https://arxiv.org/abs/2610.22012
AB - Numerical simulations are essential for studying the solar atmosphere, especially in the context of the coronal heating problem. Understanding how numerical factors and boundary conditions (BC) influence the energetic structure of these models is crucial for unraveling coronal heating mechanisms. We investigate how the computational domain dimensions, upper BC, and vertical grid resolution influence the relaxed thermal structure, energetics, and long-term stability of 2D solar atmosphere simulations with a vertical magnetic field. Using the MAGEC code, we conducted parametric studies varying domain width and height, vertical grid stretching, and top BC, including hot-plate, closed, unrestricted temperature, background heating, and non-reflecting characteristic configurations. We evaluated average temperature profiles for parameter dependence and temporal stability, and analyzed the energetic budget as a function of height, isolating contributions from radiative losses, thermal conduction, pressure work, and internal energy flux. Models with larger vertical domains show higher average coronal temperatures and steeper transition region gradients. Among the less artificial BC, non-reflecting characteristics most effectively delay coronal cooling. We determined the threshold height at which a dynamic top BC can mimic fixed-temperature setups, and established the maximum allowable grid stretching that preserves physical accuracy while reducing computational costs. The energetic terms change significantly with the chosen BC, with specific mechanisms acting as cooling agents in different atmospheric layers. Our results demonstrate that domain height, vertical grid design, and BC critically dictate the relaxed thermal structure and energetic balance of the simulated solar atmosphere.
ER -