Realistic Simulations of Stellar Surface Convection with ANTARES: I. Boundary Conditions and Model Relaxation

TL;DR

This paper introduces open boundary conditions into the ANTARES code to improve the realism of stellar surface convection simulations, highlighting the importance of initial and boundary conditions for accurate modeling.

Contribution

The paper presents a novel implementation of open boundary conditions in the ANTARES code and analyzes their impact on stellar surface convection simulations.

Findings

Boundary conditions significantly influence simulation outcomes.

Thermal stratification changes are limited by Kelvin-Helmholtz time scale.

Proper initial and boundary parameter choices are crucial for realistic results.

Abstract

We have implemented open boundary conditions into the ANTARES code to increase the realism of our simulations of stellar surface convection. Even though we greatly benefit from the high accuracy of our fifth order numerical scheme (WENO5), the broader stencils needed for the numerical scheme complicate the implementation of boundary conditions. We show that the effective temperature of a numerical simulation cannot be changed by corrections at the lower boundary since the thermal stratification does only change on the Kelvin-Helmholtz time scale. Except for very shallow models, this time scale cannot be covered by multidimensional simulations due to the enormous computational requirements. We demonstrate to what extent numerical simulations of stellar surface convection are sensitive to the initial conditions and the boundary conditions. An ill-conceived choice of parameters for the boundary conditions can have a severe impact. Numerical simulations of stellar surface convection will only be (physically) meaningful and realistic if the initial model, the extent and position of the simulation box, and the parameters from the boundary conditions are chosen adequately.