BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christian Y. Cardall

Christian Y. Cardall appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-phastro-ph.HEastro-ph.IMphysics.comp-phastro-ph.SRgr-qchep-phnucl-thphysics.flu-dynSoftware Engineering

Catalog footprint

What is connected

18works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

18 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fluid Dynamics in Curvilinear Coordinates without Fictitious Forces

Use of curvilinear coordinates is sometimes indicated by the inherent geometry of a fluid dynamics problem, but this introduces fictitious forces into the momentum equations that spoil strict conservative form. If one is willing to work in three dimensions, these fictitious forces can be eliminated by solving for rectangular (Cartesian) momentum components on a curvilinear mesh. A thoroughly geometric approach to fluid dynamics on spacetime demonstrates this transparently, while also giving insight into a greater unity of the relativistic and nonrelativistic cases than is usually appreciated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accelerating Our Understanding of Supernova Explosion Mechanism via Simulations and Visualizations with GenASiS

Core-collapse supernovae are among the most powerful explosions in the Universe, releasing about $10^{53}~\mbox{erg}$ of energy on timescales of a few tens of seconds. These explosion events are also responsible for the production and dissemination of most of the heavy elements, making life as we know it possible. Yet exactly how they work is still unresolved. One reason for this is the sheer complexity and cost of a self-consistent, multi-physics, and multi-dimensional core-collapse supernova simulation, which is impractical, and often impossible, even on the largest supercomputers we have available today. To advance our understanding we instead must often use simplified models, teasing out the most important ingredients for successful explosions, while helping us to interpret results from higher fidelity multi-physics models. In this paper we investigate the role of instabilities in the core-collapse supernova environment. We present here simulation and visualization results produced by our code GenASiS.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

GenASiS: General Astrophysical Simulation System. I. Refinable Mesh and Nonrelativistic Hydrodynamics

GenASiS (General Astrophysical Simulation System) is a new code being developed initially and primarily, though by no means exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper---the first in a series---demonstrates a centrally refined coordinate patch suitable for gravitational collapse and documents methods for compressible nonrelativistic hydrodynamics. We benchmark the hydrodynamics capabilities of GenASiS against many standard test problems; the results illustrate the basic competence of our implementation, demonstrate the strengths and limitations of the HLLC relative to the HLL Riemann solver in a number of interesting cases, and provide preliminary indications of the code's ability to scale and to function with cell-by-cell fixed-mesh refinement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

GenASiS: General Astrophysical Simulation System. II. Nonrelativistic Hydrodynamics

In this paper, the second in a series, we document the algorithms and solvers for compressible nonrelativistic hydrodynamics implemented in GenASiS (General Astrophysical Simulation System)---a new code being developed initially and primarily, though by no means exclusively, for the simulation of core-collapse supernovae. In the Mathematics division of GenASiS we introduce Solvers, which includes finite-volume updates for generic hyperbolic BalanceEquations and ordinary differential equation integration Steps. We also introduce the Physics division of GenASiS; this extends the Manifolds division of Mathematics into physical Spaces, defines StressEnergies, and combines these into Universes. We benchmark the hydrodynamics capabilities of GenASiS against many standard test problems; the results illustrate the basic competence of our implementation, demonstrate the manifest superiority of the HLLC over the HLL Riemann solver in a number of interesting cases, and provide preliminary indications of the code's ability to scale and to function with cell-by-cell fixed-mesh refinement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stochasticity and efficiency of simplified core-collapse supernova explosions

We present an initial report on 160 simulations of a highly simplified model of the post-bounce core-collapse supernova environment in three spatial dimensions (3D). We set different values of a parameter characterizing the impact of nuclear dissociation at the stalled shock in order to regulate the post-shock fluid velocity, thereby determining the relative importance of convection and the stationary accretion shock instability (SASI). While our convection-dominated runs comport with the paradigmatic notion of a `critical neutrino luminosity' for explosion at a given mass accretion rate (albeit with a nontrivial spread in explosion times just above threshold), the outcomes of our SASI-dominated runs are much more stochastic: a sharp threshold critical luminosity is `smeared out' into a rising probability of explosion over a $\sim 20\%$ range of luminosity. We also find that the SASI-dominated models are able to explode with 3 to 4 times less efficient neutrino heating, indicating that progenitor properties, and fluid and neutrino microphysics, conducive to the SASI would make the neutrino-driven explosion mechanism more robust.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Conservative 3+1 General Relativistic Boltzmann Equation

We present a new derivation of the conservative form of the general relativistic Boltzmann equation and specialize it to the 3+1 metric. The resulting transport equation is intended for use in simulations involving numerical relativity, particularly in the absence of spherical symmetry. The independent variables are lab frame coordinate basis spacetime position components and comoving frame curvilinear momentum space coordinates. With an eye towards astrophysical applications---such as core-collapse supernovae and compact object mergers---in which the fluid includes nuclei and/or nuclear matter at finite temperature, and in which the transported particles are neutrinos, we examine the relationship between lepton number and four-momentum exchange between neutrinos and the fluid.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Conservative 3+1 General Relativistic Variable Eddington Tensor Radiation Transport Equations

We present conservative 3+1 general relativistic variable Eddington tensor radiation transport equations, including greater elaboration of the momentum space divergence (that is, the energy derivative term) than in previous work. These equations are intended for use in simulations involving numerical relativity, particularly in the absence of spherical symmetry. The independent variables are the lab frame coordinate basis spacetime position coordinates and the particle energy measured in the comoving frame. With an eye towards astrophysical applications---such as core-collapse supernovae and compact object mergers---in which the fluid includes nuclei and/or nuclear matter at finite temperature, and in which the transported particles are neutrinos, we pay special attention to the consistency of four-momentum and lepton number exchange between neutrinos and the fluid, showing the term-by-term cancellations that must occur for this consistency to be achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Conservative Moment Equations for Neutrino Radiation Transport with Limited Relativity

We derive conservative, multidimensional, energy-dependent moment equations for neutrino transport in core-collapse supernovae and related astrophysical systems, with particular attention to the consistency of conservative four-momentum and lepton number transport equations. After taking angular moments of conservative formulations of the general relativistic Boltzmann equation, we specialize to a conformally flat spacetime, which also serves as the basis for four further limits. Two of these---the multidimensional special relativistic case, and a conformally flat formulation of the spherically symmetric general relativistic case---are given in appendices for the sake of comparison with extant literature. The third limit is a weak-field, `pseudo-Newtonian' approach \citep{kim_etal_2009,kim_etal_2012} in which the source of the gravitational potential includes the trace of the stress-energy tensor (rather than just the mass density), and all orders in fluid velocity $v$ are retained. Our primary interest here is in the fourth limit: `$\mathcal{O}(v)$' moment equations for use in conjunction with Newtonian self-gravitating hydrodynamics. We show that the concept of `$\mathcal{O}(v)$' transport requires care when dealing with both conservative four-momentum and conservative lepton number transport, and present two self-consistent options: `$\mathcal{O}(v)$-plus' transport, in which an $\mathcal{O}(v^2)$ energy equation combines with an $\mathcal{O}(v)$ momentum equation to give an $\mathcal{O}(v^2)$ number equation; and `$\mathcal{O}(v)$-minus' transport, in which an $\mathcal{O}(v)$ energy equation combines with an $\mathcal{O}(1)$ momentum equation to give an $\mathcal{O}(v)$ number equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Turbulent Magnetic Field Amplification from Spiral SASI Modes: Implications for Core-Collapse Supernovae and Proto-Neutron Star Magnetization

We extend our investigation of magnetic field evolution in three-dimensional flows driven by the stationary accretion shock instability (SASI) with a suite of higher-resolution idealized models of the post-bounce core-collapse supernova environment. Our magnetohydrodynamic simulations vary in initial magnetic field strength, rotation rate, and grid resolution. Vigorous SASI-driven turbulence inside the shock amplifies magnetic fields exponentially; but while the amplified fields reduce the kinetic energy of small-scale flows, they do not seem to affect the global shock dynamics. The growth rate and final magnitude of the magnetic energy are very sensitive to grid resolution, and both are underestimated by the simulations. Nevertheless our simulations suggest that neutron star magnetic fields exceeding $10^{14}$ G can result from dynamics driven by the SASI, \emph{even for non-rotating progenitors}.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Parallel FFT-based Poisson Solver for Isolated Three-dimensional Systems

We describe an implementation to solve Poisson's equation for an isolated system on a unigrid mesh using FFTs. The method solves the equation globally on mesh blocks distributed across multiple processes on a distributed-memory parallel computer. Test results to demonstrate the convergence and scaling properties of the implementation are presented. The solver is offered to interested users as the library PSPFFT.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Generation of Magnetic Fields by the Stationary Accretion Shock Instability

We begin an exploration of the capacity of the stationary accretion shock instability (SASI) to generate magnetic fields by adding a weak, stationary, and radial (but bipolar) magnetic field, and in some cases rotation, to an initially spherically symmetric fluid configuration that models a stalled shock in the post-bounce supernova environment. In axisymmetric simulations we find that cycles of latitudinal flows into and radial flows out of the polar regions amplify the field parallel to the symmetry axis, typically increasing the total magnetic energy by about two orders of magnitude. Nonaxisymmetric calculations result in fundamentally different flows and a larger magnetic energy increase: shearing associated with the SASI spiral mode contributes to a widespread and turbulent field amplification mechanism, boosting the magnetic energy by almost four orders of magnitude (a result which remains very sensitive to the spatial resolution of the numerical simulations). While the SASI may contribute to neutron star magnetization, these simulations do not show qualitatively new features in the global evolution of the shock as a result of SASI-induced magnetic field amplification.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Towards neutrino transport with flavor mixing in supernovae: the Liouville operator

The calculation of neutrino decoupling from nuclear matter requires a transport formalism capable of handling both collisions and flavor mixing. The first steps towards such a formalism are the construction of neutrino and antineutrino "distribution matrices," and a determination of the Liouville equations they satisfy in the noninteracting case. These steps are accomplished through study of a Wigner-transformed "density function," the mean value of paired neutrino quantum field operators.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Conservative Formulations of General Relativistic Radiative Transfer

Accurate accounting of particle number and 4-momentum in radiative transfer may be facilitated by the use of transport equations that allow transparent conversion between volume and surface integrals in both spacetime and momentum space. Such conservative formulations of general relativistic radiative transfer in multiple spatial dimensions are presented, and their relevance to core-collapse supernova simulations described.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

General Relativity and Neutrino-driven Supernova Winds

To date, detailed core-collapse supernova simulations have not succeeded in elucidating the explosion mechanism. But whatever the mechanism turns out to be, there will be a neutrino-heated outflow between the hot, newly-born neutron star and the outgoing supernova shock wave whose neutron richness, low density, and high temperature could provide promising conditions for the synthesis of heavy nuclei via rapid neutron capture (the r-process). General relativistic effects improve the prospects for this outflow as an r-process site. Among relativistic effects, the enhanced "gravitational potential" is more important than the gravitational redshift or trajectory bending of neutrinos.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

On Steady State Neutrino Heated Ultra-Relativistic Winds from Compact Objects

We study steady state winds from compact objects in the regime where the wind velocity at infinity is ultra-relativistic. This may have relevance to some models of Gamma-Ray-Bursts (GRB's). Particular attention is paid to the case where neutrinos provide the heating. Unless the neutrino luminosity is very large, L>10^{54} erg/s, the only allowed steady state solutions are those where energy deposition is dominated by neutrino-antineutrino annihilation at the sonic point. In this case, the matter temperature near the neutron star surface is low, less than 1MeV for typical neutrino luminosities. This is in contrast to the case for sub-relativistic winds discussed in the context of supernovae where the matter temperature near the neutron star approximates the temperature characterizing the neutrinos. We also investigate the setting of the neutron to proton ratio (n/p) in these winds and find that only for large (>10 MeV) electron neutrino or electron anti-neutrino temperatures is n/p entirely determined by neutrino capture. Otherwise, n/p retains an imprint of conditions in the neutron star.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

General relativistic effects in the neutrino-driven wind and r-process nucleosynthesis

We discuss general relativistic effects in the steady-state neutrino-driven ``wind'' which may arise from nascent neutron stars. In particular, we generalize previous analytic estimates of the entropy per baryon $S$, the mass outflow rate $\dot M$, and the dynamical expansion time scale $τ_{dyn}$. We show that $S$ increases and $τ_{dyn}$ decreases with increasing values of the mass-to-radius ratio describing the supernova core. Both of these trends indicate that a more compact core will lead to a higher number of neutrons per iron peak seed nucleus. Such an enhancement in the neutron-to-seed ratio may be required for successful r-process nucleosynthesis in neutrino-heated supernova ejecta.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Can a ``natural'' three-generation neutrino mixing scheme satisfy everything?

We examine the potential for a ``natural'' three-neutrino mixing scheme to satisfy available data and astrophysical arguments. By ``natural'' we mean no sterile neutrinos, and a neutrino mass hierarchy similar to that of the charged leptons. We seek to satisfy (or solve): 1. Accelerator and reactor neutrino oscillation constraints, including LSND; 2. The atmospheric muon neutrino deficit problem; 3. The solar neutrino problem; 4. Supernova $r$-process nucleosynthesis in neutrino-heated supernova ejecta; 5. Cold+hot dark matter models. We argue that putative supernova $r$-process nucleosynthesis bounds on two-neutrino flavor mixing can be applied directly to three-neutrino mixing in the case where one vacuum neutrino mass eigenvalue difference dominates the others. We show that in this ``one mass scale dominance'' limit, a natural three-neutrino oscillation solution meeting all the above constraints exists only if the atmospheric neutrino data {\em and} the LSND data can be explained with one neutrino mass difference. In this model, an explanation for the solar neutrino data can be effected by employing the {\em other} independent neutrino mass difference. Such a solution is only marginally allowed by the current data, and proposed long-baseline neutrino oscillation experiments can definitively rule it out. If it were ruled out, the simultaneous solution of the above constraints by neutrino oscillations would then require sterile neutrinos and/or a neutrino mass hierarchy of a different nature than that of the charged leptons.

0 citations0 reviewsNo code linkedOpen access