Ke-Jung Chen, Alexander Heger, Ann Almgren
Mapping one-dimensional stellar profiles onto multidimensional grids as
initial conditions for hydrodynamics calculations can lead to numerical
artifacts, one of the most severe of which is the violation of conservation
laws for physical quantities such as energy and mass. Here we introduce a
numerical scheme for mapping one-dimensional spherically-symmetric data onto
multidimensional meshes so that these physical quantities are conserved. We
validate our scheme by porting a realistic 1D Lagrangian stellar profile to the
new multidimensional Eulerian hydro code CASTRO. Our results show that all
important features in the profiles are reproduced on the new grid and that
conservation laws are enforced at all resolutions after mapping.
View original:
http://arxiv.org/abs/1112.3033
No comments:
Post a Comment