Statistics – Computation
Scientific paper
May 2002
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2002aas...200.1401g&link_type=abstract
American Astronomical Society, 200th AAS Meeting, #14.01; Bulletin of the American Astronomical Society, Vol. 34, p.663
Statistics
Computation
5
Scientific paper
We consider a Type Ia supernova explosion originating as a deflagration in the center of a carbon-oxygen Chandrasekhar-mass white dwarf (WD) with initial composition 0.5C+0.5O, central density 2.0 x 109 g/cm3 and initial radius 2.1 x 108 cm. A three-dimensional (3D) numerical model is based on reactive Euler equations of fluid dynamics coupled with an equation of state for a degenerate matter and a simplified kinetics of energy release. The energy-release model provides the correct propagation velocity for a laminar flame and takes into account carbon burning, as well as nuclear statistical quasi-equilibrium and equilibrium relaxations. The model for the turbulent burning on scales that are not resolved in the simulations is based on the assumption that burning on small scales is driven by the gravity-induced Rayleigh-Taylor (RT) instability. We performed 3D calculations for the first 1.9 seconds of explosion using an adaptively refined structured mesh. For the highest-resolution case, the minimum cell size was 2.6 x 105 cm, and the mesh consisted of 108 computational cells by the end of the simulation. The flame, started as a sphere with the radius 3 x 106 cm, becomes very convoluted due to the RT and Kelvin-Helmholtz instabilities on resolved scales and develops multiple buoyant plumes. As the plumes grow, the unburnt material either sinks towards the center or expands more slowly than the burnt material inside the plumes. The material burns at all distances from the center even when the larger flame plumes reach the outer layers of the star. By 1.9 seconds, some of these plumes approach the surface of the expanding WD that extends to (5-6) x 108 cm from the center. About 50% of the material burns out releasing 1.3 x 1051 ergs of nuclear energy. The expansion velocity at the surface reaches 1.2 x 109 cm/s and continues to grow. A convergence study shows that at high resolutions, the results become practically independent on the computational cell size and insensitive to subgrid model parameters.
Gamezo Vadim N.
Khokhlov Alexei M.
Oran Elaine S.
No associations
LandOfFree
Three-Dimensional Numerical Simulations of Type Ia Supernovae: Numerical Convergence for Deflagration Stage does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Three-Dimensional Numerical Simulations of Type Ia Supernovae: Numerical Convergence for Deflagration Stage, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Three-Dimensional Numerical Simulations of Type Ia Supernovae: Numerical Convergence for Deflagration Stage will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1717443