BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Paul C. Duffell

Paul C. Duffell contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.EPastro-ph.GAphysics.comp-phphysics.flu-dynastro-ph.COastro-ph.IMastro-ph.SRgr-qc

Trust snapshot

Quick read

Trust 21 - Emerging
19works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

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.

Published work

19 published item(s)

preprint2022arXiv

Numerical Investigation of Dynamical and Morphological Trends in Relativistic Jets

Active galactic nuclei (AGN) show a range of morphologies and dynamical properties, which are determined not only by parameters intrinsic to the central engine but also their interaction with the surrounding environment. We investigate the connection of kiloparsec scale AGN jet properties to their intrinsic parameters and surroundings. This is done using a suite of 40 relativistic hydrodynamic simulations spanning a wide range of engine luminosities and opening angles. We explore AGN jet propagation with different ambient density profiles, including $r^{-2}$ (self-similar solution) and $r^{-1}$, which is more relevant for AGN host environments. The Fanaroff-Riley (FR) morphological dichotomy arises naturally in our models. Jets with low energy density compared to the ambient medium produce a center-brightened emissivity distribution, while emissivity from relatively higher energy density jets is dominated by a terminal bright spot. We observe recollimation shocks in our simulations that can generate bright spots along the spine of the jet, providing a possible explanation for "knots" observed in AGN jets. We additionally find a scaling relation between the number of knots and the jet-head-to-surroundings energy density ratio. This scaling relation is generally consistent with the observations of the jets in M87 and Cygnus A. Our model also correctly predicts M87 as FR I and Cygnus A as FR II. Our model can be used to relate jet dynamical parameters such as jet head velocity, jet opening angle, and external pressure to jet power and ambient density estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Circumbinary Disks: Accretion and Torque as a Function of Mass Ratio and Disk Viscosity

Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around $0.05$. This means that for mass ratios below $0.05$, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from $α= 0.03$ to $α= 0.15$, and find that this outward torque is proportional to the viscous torque, simply proportional to viscosity in this range. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between $0.01$ and unity. For mass ratio $q > 0.6$, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a "smoking gun" prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Effect of Jet-Ejecta Interaction on the Viewing Angle Dependence of Kilonova Light Curves

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 multi-dimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at 1.5h after merger. We find that on this timescale, jet shock heating does not affect the kilonova emission. 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 ($\lesssim 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An empirically-derived formula for the shape of planet-induced gaps in protoplanetary disks

This study uses numerical hydrodynamics calculations and a novel method for densely sampling parameter space to measure the precise shape of a gap opened by a planet in a gaseous disk, as a function of planet-to-star mass ratio, disk Mach number, and disk viscosity. Formulas for gap depth and width are determined, which are combined to form a complete formula for surface density as a function of radius in the disk. This new analytical formula is compared with numerically-derived gaps opened by planets ranging from very low masses up to a few times Jupiter's mass, and excellent agreement is found over a wide range of parameter space relevant to planet-disk interactions. A simple-to-use code is presented to rapidly generate synthetic disk profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A One-Dimensional Model for Rayleigh-Taylor Instability in Supernova Remnants

This study presents a method for approximating the multidimensional effects of Rayleigh-Taylor instability as a modification of the one-dimensional hydro equations. This modification is similar to the Shakura-Sunyaev α prescription for modeling the coarse-grained effects of turbulence in astrophysical disks. The model introduces several dimensionless tunable parameters that are calibrated by comparing with high-resolution two-dimensional axisymmetric numerical calculations of Rayleigh-Taylor unstable flows. A complete description of the model is presented, along with a handful of test problems that demonstrate the extent to which the one-dimensional model is able to reproduce multidimensional effects.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

DISCO: A 3D Moving-Mesh Magnetohydrodynamics Code Designed for the Study of Astrophysical Disks

This work presents the publicly available moving-mesh magnetohydrodynamics code DISCO. DISCO is efficient and accurate at evolving orbital fluid motion in two and three dimensions, especially at high Mach number. DISCO employs a moving-mesh approach utilizing a dynamic cylindrical mesh that can shear azimuthally to follow the orbital motion of the gas. The moving mesh removes diffusive advection errors and allows for longer timesteps than a static grid. Magnetohydrodynamics is implemented in DISCO using an HLLD Riemann solver and a novel constrained transport scheme which is compatible with the mesh motion. DISCO is tested against a wide variety of problems, which are designed to test its stability, accuracy and scalability. In addition, several magnetohydrodynamics tests are performed which demonstrate the accuracy and stability of the new constrained transport approach, including two tests of the magneto-rotational instability (MRI); one testing the linear growth rate and the other following the instability into the fully turbulent regime.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Narrow Short-Duration GRB Jet from a Wide Central Engine

We use two-dimensional relativistic hydrodynamic numerical calculations to show that highly collimated relativistic jets can be produced in neutron star merger models of short-duration gamma ray bursts without the need for a highly directed engine or a large net magnetic flux. Even a hydrodynamic engine generating a very wide sustained outflow on small scales can in principle produce a highly collimated relativistic jet, facilitated by a dense surrounding medium which provides a cocoon surrounding the jet core. An oblate geometry to the surrounding gas significantly enhances the collimation process. Previous numerical simulations have shown that the merger of two neutron stars produces an oblate, expanding cloud of dynamical ejecta. We show that this gas can efficiently collimate the central engine power much as the surrounding star does in long-duration GRB models. For typical short-duration GRB central engine parameters, we find jets with opening angles of order 10 degrees in which a large fraction of the total outflow power of the central engine resides in highly relativistic material. These results predict large differences in the opening angles of outflows from binary neutron star mergers versus neutron star-black hole mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Simple Analytical Model for Gaps in Protoplanetary Disks

An analytical model is presented for calculating the surface density as a function of radius $Σ(r)$ in protoplanetary disks in which a planet has opened a gap. This model is also applicable to circumbinary disks with extreme binary mass ratios. The gap profile can be solved for algebraically, without performing any numerical integrals. In contrast with previous one-dimensional gap models, this model correctly predicts that low-mass (sub-Jupiter) planets can open gaps in sufficiently low-viscosity disks, and it correctly recovers the power-law dependence of gap depth on planet-to-star mass ratio $q$, disk aspect ratio $h/r$, and dimensionless viscosity $α$ found in previous numerical studies. Analytical gap profiles are compared with numerical calculations over a range of parameter space in $q$, $h/r$, and $α$, demonstrating accurate reproduction of the "partial gap" regime, and general agreement over a wide range of parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Eccentric Jupiters via Disk-Planet Interactions

Numerical hydrodynamics calculations are performed to determine conditions under which giant planet eccentricities can be excited by parent gas disks. Unlike in other studies, Jupiter-mass planets are found to have their eccentricities amplified --- provided their orbits start eccentric. We disentangle the web of co-rotation, co-orbital, and external resonances to show that this finite-amplitude instability is consistent with that predicted analytically. Ellipticities can grow until they reach of order the disk's aspect ratio, beyond which the external Lindblad resonances that excite eccentricity are weakened by the planet's increasingly supersonic epicyclic motion. Forcing the planet to still larger eccentricities causes catastrophic eccentricity damping as the planet collides into gap walls. For standard parameters, the range of eccentricities for instability is modest; the threshold eccentricity for growth ($\sim$$0.04$) is not much smaller than the final eccentricity to which orbits grow ($\sim$$0.07$). If this threshold eccentricity can be lowered (perhaps by non-barotropic effects), and if the eccentricity driving documented here survives in 3D, it may robustly explain the low-to-moderate eccentricities $\lesssim 0.1$ exhibited by many giant planets (including Jupiter and Saturn), especially those without planetary or stellar companions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

From Engine to Afterglow: Collapsars Naturally Produce Top-Heavy Jets and Early-Time Plateaus in Gamma Ray Burst Afterglows

We demonstrate that the steep decay and long plateau in the early phases of gamma ray burst (GRB) X-ray afterglows are naturally produced in the collapsar model, by a means ultimately related to the dynamics of relativistic jet propagation through a massive star. We present two-dimensional axisymmetric hydrodynamical simulations which start from a collapsar engine and evolve all the way through the late afterglow phase. The resultant outflow includes a jet core which is highly relativistic after breaking out of the star, but becomes baryon-loaded after colliding with a massive outer shell, corresponding to mass from the stellar atmosphere of the progenitor star which became trapped in front of the jet core at breakout. The prompt emission produced before or during this collision would then have the signature of a high Lorentz factor jet, but the afterglow is produced by the amalgamated post-collision ejecta which has more inertia than the original highly relativistic jet core and thus has a delayed deceleration. This naturally explains the early light curve behavior discovered by Swift, including a steep decay and a long plateau, without invoking late-time energy injection from the central engine. The numerical simulation is performed continuously from engine to afterglow, covering a dynamic range of over ten orders of magnitude in radius. Light curves calculated from the numerical output demonstrate that this mechanism reproduces basic features seen in early afterglow data. Initial steep decays are produced by internal shocks, and the plateau corresponds to the coasting phase of the outflow.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Halting Migration: Numerical Calculations of Corotation Torques in the Weakly Nonlinear Regime

Planets in their formative years can migrate due to the influence of gravitational torques in the protoplanetary disk they inhabit. For low-mass planets in an isothermal disk, it is known that there is a strong negative torque on the planet due to its linear perturbation to the disk, causing fast inward migration. The current investigation demonstrates that in these same isothermal disks, for intermediate-mass planets, there is a strong positive nonlinear corotation torque due to the effects of gas being pulled through a gap on horseshoe orbits. For intermediate-mass planets, this positive torque can partially or completely cancel the linear (Type I) torque, leading to slower or outward migration, even in an isothermal disk. The effect is most significant for Super-Earth and Sub-Jovian planets, during the transition from a low-mass linear perturber to a non-linear gap-opening planet, when the planet has opened a so-called 'partial gap'. In this study, numerical calculations of planet-disk interactions calculate these torques explicitly, and scalings are empirically constructed for migration rates in this weakly nonlinear regime. These results find outward migration is possible for planets with masses in the 20 - 100 Earth Mass range.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-Frequency Voronoi Noise Reduced by Smoothed Mesh Motion

We describe a technique for improving the performance of hydrodynamics codes which employ a moving Voronoi mesh. Currently, such codes are susceptible to high-frequency noise produced by rapid adjustments in the grid topology on the smallest scales. The treatment for this grid noise is simple; instead of moving the mesh-generating marker points with the local fluid velocity, this velocity field is smoothed on small scales, so that neighboring marker points generally have similar velocities. We demonstrate significant improvement gained by this adjustment in several code tests relevant to the physics which moving-mesh codes are designed to capture.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shallow Cavities in Multiple-Planet Systems

Large cavities are often observed in protoplanetary disks, which might suggest the presence of planets opening gaps in the disk. Multiple planets are necessary to produce a wide cavity in the gas. However, multiple planets may also be a burden to the carving out of very deep gaps. When additional planets are added to the system, the time-dependent perturbations from these additional satellites can stir up gas in the gap, suppressing cavity opening. In this study, we perform two-dimensional numerical hydro calculations of gap opening for single and multiple planets, showing the effect that additional planets have on the gap depths. We show that multiple planets produce much shallower cavities than single planets, so that more massive planets are needed in the multiple-planet case to produce an equivalent gap depth as in the single-planet case. To deplete a gap by a factor of 100 for the parameters chosen in this study, one only requires $M_p \approx 3.5M_J$ in the single-planet case, but much more massive planets, $M_p \approx 7M_J$ are required in the multiple-planet case. This requirement of high-mass planets implies that such planets may be detectable in the next generation of direct imaging projects, in gaps whose depths are constrained to be sufficiently deep by ALMA.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Migration of Gap-Opening Planets is not Locked to Viscous Disk Evolution

Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work, it is demonstrated using highly accurate numerical calculations that the actual migration rate is dependent on disk and planet parameters, and can be significantly larger or smaller than the viscous drift rate. In the limiting case of a disk much more massive than the secondary, the migration rate saturates to a constant which is sensitive to disk parameters and is not necessarily of order viscous rate. In the opposite limit of a low-mass disk, the migration rate decreases linearly with disk mass. Steady-state solutions in the low disk mass limit show no pile-up outside the secondary's orbit, and no corresponding drainage of the inner disk.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A "Boosted Fireball" Model for Structured Relativistic Jets

We present a model for relativistic jets which generates a particular angular distribution of Lorentz factor and energy per solid angle. We consider a fireball with specific internal energy E/M launched with bulk Lorentz factor γ_B. This "boosted fireball" model is motivated by the phenomenology of collapsar jets, but is applicable to a wide variety of relativistic flows. In its center-of-momentum frame the fireball expands isotropically, converting its internal energy into radially expanding flow with asymptotic Lorentz factor η_0 ~ E/M. In the lab frame the flow is beamed, expanding with Lorentz factor Γ= 2 η_0 γ_B in the direction of its initial bulk motion and with characteristic opening angle θ_0 ~ 1/γ_B. The flow is jet-like with Γθ_0 ~ 2 η_0 such that jets with Γ> 1/θ_0 are naturally produced. The choice η_0 ~ γ_B ~ 10 yields a jet with Γ~ 200 on-axis and angular structure characterized by opening angle θ_0 ~ 0.1 of relevance for cosmological GRBs, while γ_B >~ 1 may be relevant for low-luminosity GRBs. The model produces a family of outflows, of relevance for different relativistic phenomena with structures completely determined by η_0 and γ_B. We calculate the energy per unit solid angle for the model and use it to compute light curves for comparison with the widely used top-hat model. The jet break in the boosted fireball light curve is greatly subdued when compared to the top-hat model because the edge of the jet is smoother than for a top-hat. This may explain missing jet breaks in afterglow light curves.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gap Opening by Extremely Low Mass Planets in a Viscous Disk

By numerically integrating the compressible Navier-Stokes equations in two dimensions, we calculate the criterion for gap formation by a very low mass (q ~10^{-4}) protoplanet on a fixed orbit in a thin viscous disk. In contrast with some previously proposed gap-opening criteria, we find that a planet can open a gap even if the Hill radius is smaller than the disk scale height. Moreover, in the low-viscosity limit, we find no minimum mass necessary to open a gap for a planet held on a fixed orbit. In particular, a Neptune-mass planet will open a gap in a minimum mass solar nebula with suitably low viscosity (α<10^{-4}). We find that the mass threshold scales as the square root of viscosity in the low mass regime. This is because the gap width for critical planet masses in this regime is a fixed multiple of the scale height, not of the Hill radius of the planet.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rayleigh-Taylor Instability in a Relativistic Fireball on a Moving Computational Grid

We numerically calculate the growth and saturation of the Rayleigh-Taylor instability caused by the deceleration of relativistic outflows with Lorentz factor Γ = 10, 30, and 100. The instability generates turbulence whose scale exhibits strong dependence on Lorentz factor, as only modes with angular size smaller than 1/Γ can grow. We develop a simple diagnostic to measure the kinetic energy in turbulent fluctuations, and calculate a ratio of turbulent kinetic energy to thermal energy of .03 in the region affected by the instability. Although our numerical calculation does not include magnetic fields, we argue that small scale turbulent dynamo amplifies magnetic fields to nearly this same fraction, giving a ratio of magnetic to thermal energy of ~ .01, to within a factor of two. The instability completely disrupts the contact discontinuity between the ejecta and the swept up circumburst medium. The reverse shock is stable, but is impacted by the Rayleigh-Taylor instability, which strengthens the reverse shock and pushes it away from the forward shock. The forward shock front is unaffected by the instability, but Rayleigh-Taylor fingers can penetrate of order 10% of the way into the energetic region behind the shock during the two-shock phase of the explosion. We calculate afterglow emission from the explosion and find the reverse shock emission peaks at a later time due to its reduced Lorentz factor and modified density and pressure at the shock front. These calculations are performed using a novel numerical technique that includes a moving computational grid. The moving grid is essential as it maintains contact discontinuities to high precision and can easily evolve flows with extremely large Lorentz factors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Global Calculations of Density Waves and Gap Formation in Protoplanetary Disks using a Moving Mesh

We calculate the global quasi-steady state of a thin disk perturbed by a low-mass protoplanet orbiting at a fixed radius using extremely high-resolution numerical integrations of Euler&#39;s equations in two dimensions. The calculations are carried out using a moving computational domain, which greatly reduces advection errors and allows for much longer time-steps than a fixed grid. We calculate the angular momentum flux and the torque density as a function of radius and compare them with analytical predictions. We discuss the quasi-steady state after 100 orbits and the prospects for gap formation by low mass planets.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

TESS: A Relativistic Hydrodynamics Code on a Moving Voronoi Mesh

We have generalized a method for the numerical solution of hyperbolic systems of equations using a dynamic Voronoi tessellation of the computational domain. The Voronoi tessellation is used to generate moving computational meshes for the solution of multi-dimensional systems of conservation laws in finite-volume form. The mesh generating points are free to move with arbitrary velocity, with the choice of zero velocity resulting in an Eulerian formulation. Moving the points at the local fluid velocity makes the formulation effectively Lagrangian. We have written the TESS code to solve the equations of compressible hydrodynamics and magnetohydrodynamics for both relativistic and non-relativistic fluids on a dynamic Voronoi mesh. When run in Lagrangian mode, TESS is significantly less diffusive than fixed mesh codes and thus preserves contact discontinuities to high precision while also accurately capturing strong shock waves. TESS is written for Cartesian, spherical and cylindrical coordinates and is modular so that auxilliary physics solvers are readily integrated into the TESS framework and so that the TESS framework can be readily adapted to solve general systems of equations. We present results from a series of test problems to demonstrate the performance of TESS and to highlight some of the advantages of the dynamic tessellation method for solving challenging problems in astrophysical fluid dynamics.

0 citations0 reviewsNo code linkedOpen access