Researcher profile

Colin P. McNally

Colin P. McNally contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

6 published item(s)

preprint2026arXiv

A constrained-transport embedded boundary method for compressible resistive magnetohydrodynamics

Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of compressible resistive magnetohydrodynamics. The method is built around a finite volume formulation of the equations in which a Riemann solver is used to compute fluxes on the faces between grid cells, and a face-centered constrained transport formulation of the induction equation. The small time step problem associated with the cut cells is avoided by always computing fluxes on the faces and edges of the Cartesian mesh. We extend the method to model a moving interface between two materials with different properties using a ghost-fluid approach, and show some preliminary results including shock-wave-driven and magnetically-driven dynamical compressions of magnetohydrostatic equilibria. We present a thorough verification of the method and show that it converges at second order in the absence of discontinuities, and at first order with a discontinuity in material properties.

preprint2026arXiv

Ion Temperature Inference from Neutron Counting in Maxwellian Deuterium Plasmas

A method is presented for inferring the deuterium fuel ion temperature from neutron counts measured with fast liquid scintillators in conditions where the ion velocity distribution is Maxwellian. Local neutron count rates at each scintillator position are combined to estimate total neutron yield from the plasma, where absolute detection efficiency is determined via MCNP neutron scattering simulation based on a 3D model of the experiment structure. This method is particularly advantageous for Magnetized Target Fusion applications as it yields a time-resolved diagnostic and does not require direct line-of-sight to the plasma or collimation of the neutrons. The instrumentation configuration, pulse-shape discrimination and pile-up correction algorithms, detector calibration, and ion temperature calculation method with uncertainty characterization are discussed. An application of the method to General Fusion's Plasma Injector~3 (PI3) spherical tokamak device is demonstrated and the results are compared to an Ion Doppler spectroscopy ion temperature diagnostic.

preprint2021arXiv

Polydisperse Streaming Instability II. Methods for solving the linear stability problem

Occurring in protoplanetary discs composed of dust and gas, streaming instabilities are a favoured mechanism to drive the formation of planetesimals. The Polydispserse Streaming Instability is a generalisation of the Streaming Instability to a continuum of dust sizes. This second paper in the series provides a more in-depth derivation of the governing equations and presents novel numerical methods for solving the associated linear stability problem. In addition to the direct discretisation of the eigenproblem at second order introduced in the previous paper, a new technique based on numerically reducing the system of integral equations to a complex polynomial combined with root finding is found to yield accurate results at much lower computational cost. A related method for counting roots of the dispersion relation inside a contour without locating those roots is also demonstrated. Applications of these methods show they can reproduce and exceed the accuracy of previous results in the literature, and new benchmark results are provided. Implementations of the methods described are made available in an accompanying Python package psitools.

preprint2021arXiv

Polydisperse Streaming Instability III. Dust evolution encourages fast instability

Planet formation via core accretion requires the production of km-sized planetesimals from cosmic dust. This process must overcome barriers to simple collisional growth, for which the Streaming Instability (SI) is often invoked. Dust evolution is still required to create particles large enough to undergo vigorous instability. The SI has been studied primarily with single size dust, and the role of the full evolved dust distribution is largely unexplored. We survey the Polydispserse Streaming Instability (PSI) with physical parameters corresponding to plausible conditions in protoplanetary discs. We consider a full range of particle stopping times, generalized dust size distributions, and the effect of turbulence. We find that, while the PSI grows in many cases more slowly with a interstellar power-law dust distribution than with a single size, reasonable collisional dust evolution, producing an enhancement of the largest dust sizes, produces instability behaviour similar to the monodisperse case. Considering turbulent diffusion the trend is similar. We conclude that if fast linear growth of PSI is required for planet formation, then dust evolution producing a distribution with peak stopping times on the order of 0.1 orbits and an enhancement of the largest dust significantly above the single power-law distribution produced by a fragmentation cascade is sufficient, along with local enhancement of the dust to gas volume mass density ratio to order unity.

preprint2020arXiv

Hydrodynamical turbulence in eccentric circumbinary discs and its impact on the in situ formation of circumbinary planets

Eccentric gaseous discs are unstable to a parametric instability involving the resonant interaction between inertial-gravity waves and the eccentric mode in the disc. We present 3D global hydrodynamical simulations of inviscid circumbinary discs that form an inner cavity and become eccentric through interaction with the central binary. The parametric instability grows and generates turbulence that transports angular momentum with stress parameter $α\sim 5 \times 10^{-3}$ at distances $\lesssim 7 \;a_{bin} $, where $a_{bin}$ is the binary semi-major axis. Vertical turbulent diffusion occurs at a rate corresponding to $α_{diff}\sim 1-2\times 10^{-3}$. We examine the impact of turbulent diffusion on the vertical settling of pebbles, and on the rate of pebble accretion by embedded planets. In steady state, dust particles with Stokes numbers ${\it St} \lesssim 0.1$ form a layer of finite thickness $H_d \gtrsim 0.1 H$, where $H$ is the gas scale height. Pebble accretion efficiency is then reduced by a factor $r_{acc}/H_d$, where $r_{acc}$ is the accretion radius, compared to the rate in a laminar disc. For accreting core masses with $m_p \lesssim 0.1\; M_\oplus$, pebble accretion for particles with ${\it St} \gtrsim 0.5$ is also reduced because of velocity kicks induced by the turbulence. These effects combine to make the time needed by a Ceres-mass object to grow to the pebble isolation mass, when significant gas accretion can occur, longer than typical disc lifetimes. Hence, the origins of circumbinary planets orbiting close to their central binary systems, as discovered by the Kepler mission, are difficult to explain using an in situ model that invokes a combination of the streaming instability and pebble accretion.

preprint2020arXiv

Low-mass planet migration in three dimensional wind-driven inviscid discs: a negative corotation torque

We present simulations of low-mass planet--disc interactions in inviscid three-dimensional discs. We show that a wind-driven laminar accretion flow through the surface layers of the disc does not significantly modify the migration torque experienced by embedded planets. More importantly, we find that 3D effects lead to a dramatic change in the behaviour of the dynamical corotation torque compared to earlier 2D theory and simulations. Although it was previously shown that the dynamical corotation torque could act to slow and essentially stall the inward migration of a low-mass planet, our results in 3D show that the dynamical corotation torque has the complete opposite effect and speeds up inward migration. Our numerical experiments implicate buoyancy resonances as the cause. These have two effects: (i) they exert a direct torque on the planet, whose magnitude relative to the Lindblad torque is measured in our simulations to be small; (ii) they torque the gas librating on horseshoe orbits in the corotation region and drive evolution of its vortensity, leading to the negative dynamical corotation torque. This indicates that at low turbulent viscosity, the detailed vertical thermal structure of the protoplanetary disc plays an important role in determining the migration behaviour of embedded planets. If this result holds up under a more refined treatment of disc thermal evolution, then it has important implications for understanding the formation and early evolution of planetary systems.