We present results of simulations to assess the feasibility of modeling outflows from massive stars using the Los Alamos 3-D radiation hydrodynamics code Cassio developed for inertial confinement fusion (ICF) applications. We find that a 1-D stellar envelope simulation relaxes into hydrostatic equilibrium using computing resources that would make the simulation tractable in 2-D. We summarize next steps to include more physics fidelity and model the response to a large and abrupt energy deposition at the base of the envelope.

Carefully accounting for neutrino transport is an essential component of many astrophysical studies. Solving the full transport equation is too expensive for most realistic applications, especially those involving multiple spatial dimensions. For such cases, resorting to approximations is often the only viable option for obtaining solutions. One such approximation, which recently became popular, is the M1 method. It utilizes the system of the lowest two moments of the transport equation and closes the system with an ad hoc closure relation. The accuracy of the M1 solution depends on the quality of the closure. Several closures have been proposed in the literature and have been used in various studies. We carry out an extensive study of these closures by comparing the results of M1 calculations with precise Monte Carlo calculations of the radiation field around spherically-symmetric protoneutron star models. We find that no closure performs consistently better or worse than others in all cases. The level of accuracy a given closure yields depends on the matter configuration, neutrino type, and neutrino energy. Given this limitation, the maximum entropy closure by Minerbo (1978) on average yields relatively accurate results in the broadest set of cases considered in this work.

The light curves of Type Ia Supernovae (SN Ia) are powered by gamma-rays emitted by the decay of radioactive elements such as $^{56}$Ni and its decay products. These gamma-rays are downscattered,absorbed, and eventually reprocessed into the optical emission which makes up the bulk of all supernova observations. Detection of the gamma-rays that escape the expanding star provide the only direct means to study this power source for SN Ia light curves. Unfortunately, disagreements between calculations for the gamma-ray lines have made it difficult to interpret any gamma-ray observations. Here we present a detailed comparison of the major gamma-ray line transport codes for a series of 1-dimensional Ia models. Discrepancies in past results were due to errors in the codes, and the corrected versions of the seven different codes yield very similar results. This convergence of the simulation results allows us to infer more reliable information from the current set of gamma-ray observations of SNe Ia. The observations of SNe 1986G, 1991T and 1998bu are consistent with explosion models based on their classification: sub-luminous, super-luminous and normally-luminous respectively.