CCSN Simulation with Spectral Two-Moment Neutrino Transport Using FLASH

We are developing the toolkit for high-order neutrino-radiation hydrodynamics (thornado) to model neutrino transport in core-collapse supernova (CCSN) explosion simulations. thornado, which implements spectral two-moment neutrino transport using the discontinuous Galerkin method and implicit-explicit time stepping\footnote{\tiny Chu et al. 2019, JCP, {389}, 62; Laiu et al. 2020, J. Phys.: Conf. Ser. 1623, 012013; Laiu et al. 2020, ApJS (submitted)}, as well as WeakLib\footnote{\tiny github.com/starkiller-astro/weaklib}, a library providing tabulated microphysics, has been coupled with FLASH\footnote{\tiny Fryxell et al. 2000, ApJS, 131, 273} as an external library. With this enhanced FLASH code, we aim to simulate CCSN explosions in multiple spatial dimensions. Here we present (1) a detailed comparison between Boltztran\footnote{\tiny Mezzacappa and Bruenn 1993, ApJ, 405, 669} and thornado in the context of a fixed, spherically symmetric, post-bounce profile from a simulation that used the LS220 equation of state (EoS) and ``Bruenn~85'' neutrino opacities, and (2) the spherically symmetric gravitational collapse of a 15 solar mass progenitor using the SFHo EoS, Bruenn~85 opacities, FLASH's Newtonian hydrodynamics with self-gravity, and thornado's neutrino transport.