H. Grimm-Strele, F. Kupka, B. Löw-Baselli, E. Mundprecht, F. Zaussinger, P. Schiansky
We have implemented open boundary conditions into the ANTARES code to increase the realism in our simulations of stellar surface convection. Even though we greatly benefit from the high accuracy due to our 5th-order numerical scheme (WENO5), the broader stencils needed for the numerical scheme complicate the implementation of boundary conditions. Moreover, one-dimensional starting models may not provide enough data to set up the simulation domain large enough in each direction, or the initial data may lead to numerical instabilities. We show how to overcome these difficulties and demonstrate to what extent numerical simulations of stellar surface convection are sensitive to the numerical setup and the boundary conditions. A sloppy choice of parameters contained in the boundary conditions can have a severe impact. Only if the initial model, the extent and position of the simulation box, and the parameters from the boundary conditions are chosen adequately, numerical simulations of stellar surface convection are (physically) meaningful and realistic throughout the entire simulation domain.
View original:
http://arxiv.org/abs/1305.0743
No comments:
Post a Comment