(Zingale, M., Woosley, S. E. Rendleman, C. A., Day, M. S., Bell, J. B. 2005, ApJ, submitted)
[astro-ph preprint]
in collaboration with the CCSE/LBL group.
We have extended our two-dimensional reactive Rayleigh-Taylor study to three-dimensions, focusing on the flamelet regime. We want to understand the nature of flame-generated turbulence in 3d, and to determine whether the less computationally expensive 2d calculations provide an accurate description of the flame. All of these calculations are fully resolved---no flame model is used. These are large calculations, encompassing over 170 million zones at the end.
This calculation is at a density of 1.5 x 107 g/cc, appropriate to the late stages of a Type Ia supernova. The energy generation is a single C12+C12 thermonuclear fusion rate.
Power spectrum of the 3d RT. The solid curve is the uncorrected velocity power spectrum. The dashed curve is after projecting out the expansion components, and the dotted straight line is a n=-5/3 power law. At late times, we get a Kolmogorov power spectrum for the kinetic energy.
| Image | Power Spectrum | |
| step = 0, time = 0 s | ||
|
|
|
| | ||
| step = 120, time = 6.62e-4 s | ||
|
|
|
| | ||
| step = 240, time = 9.43e-4 s | ||
|
|
|
| | ||
| step = 360, time = 1.16e-3 s | ||
|
|
|
An AVI movie of the carbon mass fraction.
These calculations were carried out on the NERSC IBM/SP seaborg and the NASA/Ames Columbia machine. This research is supported by DOE grant No. DE-FC02-01ER41176 to the Supernova Science Center/UCSC.