COCONUT: A coronal model with an energy decomposition strategy
Fuente:
arXiv
Saved in:
| Main Authors: | , , , , , , , , , , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866911365153685504 |
|---|---|
| author | Wang, Haopeng Poedts, Stefaan Lani, Andrea Dhib, Rayan Linan, Luis Baratashvili, Tinatin Jeong, Hyun-Jin Zhou, Yuhao Li, Yucong Najafi-Ziyazi, Mahdi Wang, Juan Schmieder, Brigitte Wang, Wensi Husidic, Edin |
| author_facet | Wang, Haopeng Poedts, Stefaan Lani, Andrea Dhib, Rayan Linan, Luis Baratashvili, Tinatin Jeong, Hyun-Jin Zhou, Yuhao Li, Yucong Najafi-Ziyazi, Mahdi Wang, Juan Schmieder, Brigitte Wang, Wensi Husidic, Edin |
| contents | In this paper, we propose an energy decomposition method combined with an HLL Riemann solver that includes an additional dissipation term in the energy equation to improve the numerical stability of the fully implicit, time-evolving coronal model COCONUT and extend its applicability to solar-maximum phases. In MHD simulations that evolve conservative variables in time, the thermal pressure is typically computed by subtracting the magnetic and kinetic energies from the total energy. In low-beta (the ratio of thermal to magnetic pressure; $< 10^{-3}$) regions, discretization errors of magnetic energy can be comparable to the thermal pressure, potentially leading to negative thermal pressure and causing the simulation to crash. Therefore, we update the decomposed energy, excluding the magnetic energy, at each time step. It avoids subtracting a large magnetic energy from the total energy to obtain a very small thermal pressure in low-$β$ regions, thereby improving the numerical stability of MHD models. We validate the algorithm using a time-evolving solar-maximum Carrington rotation simulation in 2025, which the previous code failed to run to completion. We also perform quasi-steady-state coronal simulations and 2D benchmark tests to further assess the algorithm's performance. The simulation results show that the algorithm produces results nearly identical to those obtained using the traditional full energy equation during solar minimum, while significantly improving COCONUT's ability to simulate coronal evolution under strong magnetic fields, even including fields exceeding 100 Gauss with $β<10^{-3}$. This method provides a promising approach for performing quasi-realistic coronal simulations during solar maxima. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2508_20423 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | COCONUT: A coronal model with an energy decomposition strategy Wang, Haopeng Poedts, Stefaan Lani, Andrea Dhib, Rayan Linan, Luis Baratashvili, Tinatin Jeong, Hyun-Jin Zhou, Yuhao Li, Yucong Najafi-Ziyazi, Mahdi Wang, Juan Schmieder, Brigitte Wang, Wensi Husidic, Edin Solar and Stellar Astrophysics In this paper, we propose an energy decomposition method combined with an HLL Riemann solver that includes an additional dissipation term in the energy equation to improve the numerical stability of the fully implicit, time-evolving coronal model COCONUT and extend its applicability to solar-maximum phases. In MHD simulations that evolve conservative variables in time, the thermal pressure is typically computed by subtracting the magnetic and kinetic energies from the total energy. In low-beta (the ratio of thermal to magnetic pressure; $< 10^{-3}$) regions, discretization errors of magnetic energy can be comparable to the thermal pressure, potentially leading to negative thermal pressure and causing the simulation to crash. Therefore, we update the decomposed energy, excluding the magnetic energy, at each time step. It avoids subtracting a large magnetic energy from the total energy to obtain a very small thermal pressure in low-$β$ regions, thereby improving the numerical stability of MHD models. We validate the algorithm using a time-evolving solar-maximum Carrington rotation simulation in 2025, which the previous code failed to run to completion. We also perform quasi-steady-state coronal simulations and 2D benchmark tests to further assess the algorithm's performance. The simulation results show that the algorithm produces results nearly identical to those obtained using the traditional full energy equation during solar minimum, while significantly improving COCONUT's ability to simulate coronal evolution under strong magnetic fields, even including fields exceeding 100 Gauss with $β<10^{-3}$. This method provides a promising approach for performing quasi-realistic coronal simulations during solar maxima. |
| title | COCONUT: A coronal model with an energy decomposition strategy |
| topic | Solar and Stellar Astrophysics |
| url | https://arxiv.org/abs/2508.20423 |