Simulations of stellar convection with CO5BOLD

• Published in: Journal of computational physics 2012-02, Vol.231 (3), p.919-959
• Main Authors: Freytag, B., Steffen, M., Ludwig, H.-G., Wedemeyer-Böhm, S., Schaffenberger, W., Steiner, O.
• Format: Article
• Language: English
• Subjects: Numerical simulations; Radiation (magneto)hydrodynamics; Stellar surface convection
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2011.09.026

High-resolution images of the solar surface show a granulation pattern of hot rising and cooler downward-sinking material – the top of the deep-reaching solar convection zone. Convection plays a role for the thermal structure of the solar interior and the dynamo acting there, for the stratification of the photosphere, where most of the visible light is emitted, as well as for the energy budget of the spectacular processes in the chromosphere and corona. Convective stellar atmospheres can be modeled by numerically solving the coupled equations of (magneto)hydrodynamics and non-local radiation transport in the presence of a gravity field. The CO5BOLD code described in this article is designed for so-called “realistic” simulations that take into account the detailed microphysics under the conditions in solar or stellar surface layers (equation-of-state and optical properties of the matter). These simulations indeed deserve the label “realistic” because they reproduce the various observables very well – with only minor differences between different implementations. The agreement with observations has improved over time and the simulations are now well-established and have been performed for a number of stars. Still, severe challenges are encountered when it comes to extending these simulations to include ideally the entire star or substellar object: the strong stratification leads to completely different conditions in the interior, the photosphere, and the corona. Simulations have to cover spatial scales from the sub-granular level to the stellar diameter and time scales from photospheric wave travel times to stellar rotation or dynamo cycle periods. Various non-equilibrium processes have to be taken into account. Last but not least, realistic simulations are based on detailed microphysics and depend on the quality of the input data, which can be the actual accuracy limiter. This article provides an overview of the physical problem and the numerical solution and the capabilities of CO5BOLD, illustrated with a number of applications.