Researcher profile

Donald E. Willcox

Donald E. Willcox contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.HE astro-ph.SR astro-ph.IM physics.comp-ph

Trust snapshot

Quick read

Trust 15 - Baseline

3works

0followers

4topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

thornado+Flash-X: A Hybrid DG-IMEX and Finite-Volume Framework for Neutrino-Radiation Hydrodynamics in Core-Collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement (AMR)-based multiphysics simulation framework Flash-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $O(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino-matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino-electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved with a hybrid DG-FV representation of the fluid variables and operator-split evolution in Flash-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within Flash-X. thornado's neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

preprint2019arXiv

Type Ia Supernova Explosions from Hybrid Carbon-Oxygen-Neon White Dwarf Progenitors That Have Mixed During Cooling

The creation of "hybrid" white dwarfs, made of a C-O core within a O-Ne shell has been proposed, and studies indicate that ignition in the C-rich central region makes these viable progenitors for thermonuclear (type Ia) supernovae. Recent work found that the C-O core is mixed with the surrounding O-Ne as the white dwarf cools prior to accretion, which results in lower central C fractions in the massive progenitor than previously assumed. To further investigate the efficacy of hybrid white dwarfs as progenitors of thermonuclear supernovae, we performed simulations of thermonuclear supernovae from a new series of hybrid progenitors that include the effects of mixing during cooling. The progenitor white dwarf model was constructed with the one-dimensional stellar evolution code MESA and represented a star evolved through the phase of unstable interior mixing followed by accretion until it reached conditions for the ignition of carbon burning. This MESA model was then mapped to a two-dimensional initial condition for explosions simulated with FLASH. For comparison, similar simulations were performed for a traditional C-O progenitor white dwarf. By comparing the yields of the explosions, we find that, as with earlier studies, the lower C abundance in the hybrid progenitor compared to the traditional C-O progenitor leads to a lower average yield of 56Ni. Although the unmixed hybrid WD showed a similar decrement also in total iron group yield, the mixed case does not and produces a smaller fraction of iron group elements in the form of 56Ni. We attribute this to the higher central density required for ignition and the location, center or off-center, of deflagration ignition.

preprint2016arXiv

Type Ia Supernova Explosions from Hybrid Carbon-Oxygen-Neon White Dwarf Progenitors

Motivated by recent results in stellar evolution that predict the existence of hybrid white dwarf (WD) stars with a C-O core inside an O-Ne shell, we simulate thermonuclear (Type Ia) supernovae from these hybrid progenitors. We use the FLASH code to perform multidimensional simulations in the deflagration to detonation transition (DDT) explosion paradigm. Our hybrid progenitor models were produced with the MESA stellar evolution code and include the effects of the Urca process, and we map the progenitor model to the FLASH grid. We performed a suite of DDT simulations over a range of ignition conditions consistent with the progenitor's thermal and convective structure assuming multiple ignition points. To compare the results from these hybrid WD stars to previous results from C-O white dwarfs, we construct a set of C-O WD models with similar properties and similarly simulate a suite of explosions. We find that despite significant variability within each suite, trends distinguishing the explosions are apparent in their $^{56}$Ni yields and the kinetic properties of the ejecta. We comment on the feasibility of these explosions as the source of some classes of observed subluminous events.