We describe the CRASH (Center for Radiative Shock Hydrodynamics) code, a block adaptive mesh code for multi-material radiation hydrodynamics. The implementation solves the radiation diffusion model with the gray or multigroup method and uses a flux limited diffusion approximation to recover the free-streaming limit. The electrons and ions are allowed to have different temperatures and we include a flux limited electron heat conduction. The radiation hydrodynamic equations are solved in the Eulerian frame by means of a conservative finite volume discretization in either one, two, or three-dimensional slab geometry or in two-dimensional cylindrical symmetry. An operator split method is used to solve these equations in three substeps: (1) solve the hydrodynamic equations with shock-capturing schemes, (2) a linear advection of the radiation in frequency-logarithm space, and (3) an implicit solve of the stiff radiation diffusion, heat conduction, and energy exchange. We present a suite of verification test problems to demonstrate the accuracy and performance of the algorithms. The CRASH code is an extension of the Block-Adaptive Tree Solarwind Roe Upwind Scheme (BATS-R-US) code with this new radiation transfer and heat conduction library and equation-of-state and multigroup opacity solvers. Both CRASH and BATS-R-US are part of the publicly available Space Weather Modeling Framework (SWMF).

It is now generally agreed that multidimensional, multigroup, radiation hydrodynamics is an indispensable element of any realistic model of stellar-core collapse, core-collapse supernovae, and protoneutron star instabilities. We have developed a new, two-dimensional, multigroup algorithm that can model neutrino-radiation-hydrodynamic flows in core-collapse supernovae. Our algorithm uses an approach that is similar to the ZEUS family of algorithms, originally developed by Stone and Norman. However, we extend that previous work in three significant ways: First, we incorporate multispecies, multigroup, radiation hydrodynamics in a flux-limited-diffusion approximation. Our approach is capable of modeling pair-coupled neutrino-radiation hydrodynamics, and includes effects of Pauli blocking in the collision integrals. Blocking gives rise to nonlinearities in the discretized radiation-transport equations, which we evolve implicitly in time. We employ parallelized Newton-Krylov methods to obtain a solution of these nonlinear, implicit equations. Our second major extension to the ZEUS algorithm is inclusion of an electron conservation equation, which describes evolution of electron-number density in the hydrodynamic flow. This permits following the effects of deleptonization in a stellar core. In our third extension, we have modified the hydrodynamics algorithm to accommodate realistic, complex equations of state, including those having non-convex behavior. In this paper, >... <abstract continues>

Understanding the explosion mechanism of core collapse supernovae is a problem that has plagued nuclear astrophysicists since the first computational models of this phenomenon were carried out in the 1960s. Our current theories of this violent phenomenon center around multi-dimensional effects involving radiation-hydrodynamic flows of hot, dense matter and neutrinos. Modeling these multi-dimensional radiative flows presents a computational challenge that will continue to stress high-performance computing beyond the teraflop to the petaflop level. In this paper we describe a few of the scientific discoveries that we have made via terascale computational simulations of supernovae under the auspices of the SciDAC-funded Terascale Supernova Initiative.

We discuss recent advances in the radiative-hydrodynamic modeling of core collapse supernovae in multi-dimensions. A number of earlier attempts at fully radiation-hydrodynamic models utilized either the grey approximation to describe the neutrino distribution or utilized more sophisticated multigroup transport methods restricted to radial rays. In both cases these models have also neglected the O(v/c) terms that couple the radiation and matter strongly in the optically thick regions of the collapsed core. In this paper we present some recent advances that resolve some shortcomings of earlier models.