Source author record

Ben A. Ayliffe

Ben A. Ayliffe appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP astro-ph.SR

Catalog footprint

What is connected

7works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2014arXiv

Ambipolar diffusion in smoothed particle magnetohydrodynamics

In partially ionised plasmas, the magnetic field can become decoupled from the neutral gas and diffuse through it in a process known as ambipolar diffusion. Although ambipolar diffusion has been implemented in several grid codes, we here provide an implementation in smoothed particle magnetohydrodynamics (SPMHD). We use the strong coupling approximation in which the ion density is negligible, allowing a single fluid approach. The equations are derived to conserve energy, and to provide a positive definite contribution to the entropy. We test the implementation in both a simple 1D SPMHD code and the fully 3D code PHANTOM. The wave damping test yields agreement within 0.03-2 per cent of the analytical result, depending on the value of the collisional coupling constant. The oblique C-shocks test yields results that typically agree within 4 per cent of the semi-analytical result. Our algorithm is therefore suitable for exploring the effect ambipolar diffusion has on physical processes, such as the formation of stars from molecular clouds.

preprint2012arXiv

On the accumulation of planetesimals near disc gaps created by protoplanets

We have performed three-dimensional two-fluid (gas-dust) hydrodynamical models of circumstellar discs with embedded protoplanets (3 - 333 M\oplu) and small solid bodies (radii 10cm to 10m). We find that high mass planets (\gtrsim Saturn mass) open sufficiently deep gaps in the gas disc such that the density maximum at the outer edge of the gap can very efficiently trap metre-sized solid bodies. This allows the accumulation of solids at the outer edge of the gap as solids from large radii spiral inwards to the trapping region. This process of accumulation occurs fastest for those bodies that spiral inwards most rapidly, typically metre-sized boulders, whilst smaller and larger objects will not migrate sufficiently rapidly in the discs lifetime to benefit from the process. Around a Jupiter mass planet we find that bound clumps of solid material, as large as several Earth masses, may form, potentially collapsing under self-gravity to form planets or planetesimals. These results are in agreement with Lyra et al. (2009), supporting their finding that the formation of a second generation of planetesimals or of terrestrial mass planets may be triggered by the presence of a high mass planet.

preprint2012arXiv

The growth and hydrodynamic collapse of a protoplanet envelope

We have conducted three-dimensional self-gravitating radiation hydrodynamical models of gas accretion onto high mass cores (15-33 Earth masses) over hundreds of orbits. Of these models, one case accretes more than a third of a Jupiter mass of gas, before eventually undergoing a hydrodynamic collapse. This collapse causes the density near the core to increase by more than an order of magnitude, and the outer envelope to evolve into a circumplanetary disc. A small reduction in the mass within the Hill radius (R_H) accompanies this collapse as a shock propagates outwards. This collapse leads to a new hydrostatic equilibrium for the protoplanetary envelope, at which point 97 per cent of the mass contained within the Hill radius is within the inner 0.03 R_H which had previously contained less than 40 per cent. Following this collapse the protoplanet resumes accretion at its prior rate. The net flow of mass towards this dense protoplanet is predominantly from high latitudes, whilst at the outer edge of the circumplanetary disc there is net outflow of gas along the midplane. We also find a turnover of gas deep within the bound envelope that may be caused by the establishment of convection cells.

preprint2011arXiv

Migration of protoplanets with surfaces through discs with steep temperature gradients

We perform three-dimensional self-gravitating radiative transfer simulations of protoplanet migration in circumstellar discs to explore the impact upon migration of the radial temperature profiles in these discs. We model protoplanets with masses ranging between 10-100 M\bigoplus, in discs with surface density profiles of Σ \varpropto r^-1/2, and temperature profiles of the form T \varpropto r^-β, where β ranges 0-2. We find that steep (β > 1) temperature profiles lead to outward migration of low mass protoplanets in interstellar grain opacity discs, but in more optically thin discs the migration is always inwards. The trend in migration rates with changing β obtained from our models shows good agreement with those obtained using recent analytic descriptions which include consideration of the co-orbital torques and their saturation. We find that switching between two models of the protoplanet, one in which accretion acts by evacuating gas and one in which gas piles up on a surface to form an atmosphere, leads to a small shift in the migration rates. If comparing these models in discs with conditions which lead to a marginally inward migration, the small shift can lead to outward migration. However, the direction and speed of migration is dominated by disc conditions rather than by the specific prescription used to model the flow near the protoplanet.

preprint2010arXiv

Planet migration: self-gravitating radiation hydrodynamical models of protoplanets with surfaces

We calculate radial migration rates of protoplanets in laminar minimum mass solar nebula discs using three-dimensional self-gravitating radiation hydrodynamical (RHD) models. The protoplanets are free to migrate, whereupon their migration rates are measured. For low mass protoplanets (10-50 M_\oplus) we find increases in the migration timescales of up to an order of magnitude between locally-isothermal and RHD models. In the high-mass regime the migration rates are changed very little. These results are arrived at by calculating migration rates in locally-isothermal models, before sequentially introducing self-gravity, and radiative transfer, allowing us to isolate the effects of the additional physics. We find that using a locally-isothermal equation of state, without self-gravity, we reproduce the migration rates obtained by previous analytic and numerical models. We explore the impact of different protoplanet models, and changes to their assumed radii, upon migration. The introduction of self-gravity gives a slight reduction of the migration rates, whilst the inertial mass problem, which has been proposed for high mass protoplanets with circumplanetary discs, is reproduced. Upon introducing radiative transfer to models of low mass protoplanets (\approx 10 M_\oplus), modelled as small radius accreting point masses, we find outward migration with a rate of approximately twice the analytic inward rate. However, when modelling such a protoplanet in a more realistic manner, with a surface which enables the formation of a deep envelope, this outward migration is not seen.

preprint2009arXiv

Circumplanetary disc properties obtained from radiation hydrodynamical simulations of gas accretion by protoplanets

We investigate the properties of circumplanetary discs formed in three-dimensional, self-gravitating radiation hydrodynamical models of gas accretion by protoplanets. We determine disc sizes, scaleheights, and density and temperature profiles for different protoplanet masses, in solar nebulae of differing grain opacities. We find that the analytical prediction of circumplanetary disc radii in an evacuated gap (R_Hill/3) from Quillen & Trilling (1998) yields a good estimate for discs formed by high mass protoplanets. The radial density profiles of the circumplanetary discs may be described by power-laws between r^-2 and r^-3/2. We find no evidence for the ring-like density enhancements that have been found in some previous models of circumplanetary discs. Temperature profiles follow a ~r^-7/10 power-law regardless of protoplanet mass or nebula grain opacity. The discs invariably have large scaleheights (H/r > 0.2), making them thick in comparison with their encompassing circumstellar discs, and they show no flaring.