Engineering topics

Duffell, Paul C.

Publications and source records attributed to Duffell, Paul C..

4 records

DOE OSTI2024

The Santa Barbara Binary-disk Code Comparison

We have performed numerical calculations of a binary interacting with a gas disk, using 11 different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, 0.2(a 2 Ω B /v) is the kinematic viscosity. The largest discrepancies between codes resulted from the dimensionality of the setup (3D vs. 2D disks). We find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

The effect of jet–ejecta interaction on the viewing angle dependence of kilonova light curves

ABSTRACT The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates of a suite of 2D relativistic hydrodynamic simulations of jet–ejecta interaction as initial conditions for multidimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at $1.5\, \mathrm{h}$ after merger. We find that on this time-scale, jet shock heating does not affect the kilonova emission for the jet parameters we survey. However, the jet disruption to the density structure of the ejecta does change the light curves. The jet carves a channel into the otherwise spheroidal ejecta, revealing the hot, inner regions. As seen from near (≲30°) the jet axis, the kilonova is brighter by a factor of a few and bluer. The strength of this effect depends on the jet parameters, since the light curves of more heavily disrupted ejecta are more strongly affected. The light curves and spectra are also more heavily modified in the ultraviolet than in the optical.

Klion, Hannah↗

DOE OSTI2018

Jet Dynamics in Compact Object Mergers: GW170817 Likely Had a Successful Jet

Here, we use relativistic hydrodynamic numerical calculations to study the interaction between a jet and a homologous outflow produced dynamically during binary neutron star mergers. We quantify how the thermal energy imparted by the jet and the ability of the jet to escape the ejecta depend on the parameters of the jet engine and the ejecta. Under our assumptions, a collimated jet initiated at early times compared to the engine duration, we show that successful breakout of the forward cocoon shock necessitates a jet that successfully escapes the ejecta. This is because the ejecta is expanding, and the forward shock from a failed jet stalls before it reaches the edge of the ejecta. This conclusion can be circumvented only for very energetic wide angle jets, with parameters that are uncomfortable given short-duration GRB observations. For successful jets, we find two regimes of jet breakout from the ejecta: early breakout on timescales shorter than the engine duration, and late breakout well after the engine shuts off. A late breakout can explain the observed delay between gravitational waves and gamma rays in GW170817. We show that for the entire parameter space of jet parameters surveyed here (covering energies ~10 48 –10 51 erg and opening angles θ j ~ 0.07–0.4) the thermal energy deposited by the jet is less than that produced by r-process heating on second timescales by at least an order of magnitude. Shock heating is thus energetically subdominant in setting the luminosity of thermally powered transients coincident with neutron star mergers (kilonovae).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Rayleigh–Taylor Instability in Interacting Supernovae: Implications for Synchrotron Magnetic Fields

Synchrotron emission from a supernova (SN) necessitates a magnetic field, but it is unknown how strong the relevant magnetic fields are, and what mechanism generates them. In this study, we perform high-resolution numerical gas dynamics calculations in axisymmetry to determine the growth of turbulence due to Rayleigh–Taylor (RT) instability, and the resulting kinetic energy in turbulent fluctuations, as a means of inferring how strong magnetic fields can become when amplified by this turbulence. Assuming rough equipartition between kinetic and magnetic energy in the turbulence, we find that RT instability may produce turbulent fluctuations strong enough to amplify magnetic fields to a few percent of equipartition with the thermal energy. This turbulence stays concentrated near the reverse shock, but averaging this magnetic energy throughout the shocked region (weighting by emissivity) sets the magnetic fields at a minimum of 0.3 percent of equipartition. This line of argument predicts a minimum effective magnetic field strength (ϵ{sub B}>0.003) that should be present in all interacting SNe. This provides a prediction for what should be found in highly resolved, three-dimensional magnetohydrodynamics calculations. The strength and spatial distribution of turbulently generated magnetic fields would have implications for the shape and luminosity of SN radio light curves.

79 ASTRONOMY AND ASTROPHYSICS↗