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

Guillaume Laibe

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.IMastro-ph.GAastro-ph.COastro-ph.HEphysics.ao-phphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

29works
9topics
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

29 published item(s)

preprint2026arXiv

Rescaling Transforms for Local Models of Spherical Flows

Previously we developed a local model for a spherically contracting/expanding gas cloud that can be used to study turbulence and small scale instabilities in such flows. In this work we generalise the super-comoving variables used in studies of cosmological structure formation to our local spherical flow model, which make it significantly easier to derive analytical solutions and analyse the interactions of more complex flows with the background. We show that a wide class of solutions to the local spherical flow model can be obtained via a mapping from the corresponding solutions in regular Cartesian flows. The rescaling of time in the transformation results in a modification of the linear instabilities that can occur in spherical flows, causing them to have a time dependent growth rate in the physical time coordinate, and can prevent slower instabilities from operating. Finally, we show that the small scale flows in isotropic contraction/expansion can be mapped directly to Cartesian, inviscid, incompressible hydrodynamics, meaning that one expects a form of rescaled Kolmogorov-turbulence at the small scale of isotropically contracting/expanding flows.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A solution to the overdamping problem when simulating dust-gas mixtures with smoothed particle hydrodynamics

We present a fix to the overdamping problem found by Laibe & Price (2012) when simulating strongly coupled dust-gas mixtures using two different sets of particles using smoothed particle hydrodynamics. Our solution is to compute the drag at the barycentre between gas and dust particle pairs when computing the drag force by reconstructing the velocity field, similar to the procedure in Godunov-type solvers. This fixes the overdamping problem at negligible computational cost, but with additional memory required to store velocity derivatives. We employ slope limiters to avoid spurious oscillations at shocks, finding the van Leer Monotonized Central limiter most effective.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Channels for streaming instability in dusty discs

Streaming instability is a privileged channel to bridge the gap between collisional growth of dust grains and planetesimal formation triggered by gravity. This instability is thought to develop through its secular mode, which is long-time growing and may not develop easily in real discs. We address this point by revisiting its perturbation analysis. A third-order expansion with respect to the Stokes number reveals important features over-looked so far. The secular mode can be stable. Epicycles can be unstable, more resistant to viscosity and are identified by Green's function analysis as promising channels for planetesimals formation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the settling of small grains in dusty discs: analysis and formulas

Instruments achieve sharper and finer observations of micron-in-size dust grains in the top layers of young stellar discs. To provide accurate models, we revisit the theory of dust settling for small grains, when gas stratification, dust inertia and finite correlation times for the turbulence should be handled simultaneously. We start from a balance of forces and derive distributions at steady-state. Asymptotic expansions require caution since limits do not commute. In particular, non-physical bumpy distributions appear when turbulence is purely diffusive. This excludes very short correlation times for real discs, as predicted by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Protostellar collapse: the conditions to form dust rich protoplanetary disks

Dust plays a key role during star, disk and planet formation. Yet, its dynamics during the protostellar collapse remains a poorly investigated field. Recent studies seem to indicate that dust may decouple efficiently from the gas during these early stages. We aim to understand how much and in which regions dust grains concentrate during the early phases of the protostellar collapse, and see how it depends on the properties of the initial cloud and of the solid particles. We use the multiple species dust dynamics solver of the grid-based code RAMSES to perform various simulations of dusty collapses. We perform hydrodynamical and MHD simulations where we vary the maximum grain size, the thermal-to-gravitational energy ratio and the magnetic properties of the cloud. We simulate the simultaneous evolution of ten neutral dust grains species with grain sizes varying from a few nm to a few hundredth of microns. We obtain a significant decoupling between the gas and the dust for grains of typical sizes a few 10 microns. This decoupling strongly depends on the thermal-to-gravitational energy ratio, the grain sizes or the inclusion of a magnetic field. With a semi-analytic model calibrated on our results, we show that the dust ratio mostly varies exponentially with the initial Stokes number at a rate that depends on the local cloud properties. We find that larger grains tend to settle and drift efficiently in the first-core and in the newly formed disk. This can produce dust-to-gas ratios of several times the initial value. Dust concentrates in high density regions and is depleted in low density regions. The size at which grains decouple from the gas depends on the initial properties of the clouds. Since dust can not necessarily be used as a proxy for gas during the collapse, we emphasize on the necessity of including the treatment of its dynamics in collapse simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Acceleration of Superrotation in Simulated Hot Jupiter Atmospheres

Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady state can be reached. Aims. A consistent physical understanding of superrotation is needed for arbitrary drag and radiative timescales, and the relevance of considering linear steady states needs to be assessed. Methods. We obtain an analytical expression for the structure, frequency and decay rate of propagating waves in hot Jupiter atmospheres around a state at rest in the 2D shallow water beta plane limit. We solve this expression numerically and confirm the robustness of our results with a 3D linear wave algorithm. We then compare with 3D simulations of hot Jupiter atmospheres and study the non linear momentum fluxes. Results. We show that under strong day night heating the dynamics does not transit through a linear steady state when starting from an initial atmosphere in solid body rotation. We further show that non linear effects favour the initial spin up of superrotation and that the acceleration due to the vertical component of the eddy momentum flux is critical to the initial development of superrotation. Conclusions. Overall, we describe the initial phases of the acceleration of superrotation, including consideration of differing radiative and drag timescales, and conclude that eddy-momentum driven superrotating equatorial jets are robust, physical phenomena in simulations of hot Jupiter atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Plane-Parallel Wind Solution For Testing Numerical Simulations of Photoevaporation

Here we derive a Parker-wind like solution for a stratified, plane-parallel atmosphere undergoing photoionisation. The difference compared to the standard Parker solar wind is that the sonic point is crossed only at infinity. The simplicity of the analytic solution makes it a convenient test problem for numerical simulations of photoevaporation in protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Grand challenges in protoplanetary disc modelling

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On dust entrainment in photoevaporative winds

We investigate dust entrainment by photoevaporative winds in protoplanetary discs using dusty smoothed particle hydrodrodynamics (SPH). We use unequal-mass particles to resolve more than five orders of magnitude in disc/outflow density and a one-fluid formulation to efficiently simulate an equivalent magnitude range in drag stopping time. We find that only micron sized dust grains and smaller can be entrained in EUV driven winds. The maximum grain size is set by dust settling in the disc rather than aerodynamic drag in the wind. More generally, there is a linear relationship between the base flow density and the maximum entrainable grain size in the wind. A pileup of micron sized dust grains can occur in the upper atmosphere at critical radii in the disc as grains decouple from the low-density wind. Entrainment is a strong function of location in the disc, resulting in a size sorting of grains in the outflow---the largest grain being carried out between $10$--$20\,$AU. The peak dust density for each grain size occurs at the inner edge of its own entrainment region.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the maximum grain size entrained by photoevaporative winds

We model the behaviour of dust grains entrained by photoevaporation-driven winds from protoplanetary discs assuming a non-rotating, plane-parallel disc. We obtain an analytic expression for the maximum entrainable grain size in extreme-UV radiation-driven winds, which we demonstrate to be proportional to the mass loss rate of the disc. When compared with our hydrodynamic simulations, the model reproduces almost all of the wind properties for the gas and dust. In typical turbulent discs, the entrained grain sizes in the wind are smaller than the theoretical maximum everywhere but the inner disc due to dust settling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the origin of horseshoes in transitional discs

We investigate whether the rings, lopsided features and horseshoes observed at millimetre wavelengths in transitional discs can be explained by the dynamics of gas and dust at the edge of the cavity in circumbinary discs. We use 3D dusty smoothed particle hydrodynamics calculations to show that binaries with mass ratio $q \gtrsim 0.04$ drive eccentricity in the central cavity, naturally leading to a crescent-like feature in the gas density, which is accentuated in the mm dust grain population with intensity contrasts in mm-continuum emission of 10 or higher. We perform mock observations to demonstrate that these features closely match those observed by ALMA, suggesting that the origin of rings, dust horseshoes and other non-axisymmetric structures in transition discs can be explained by the presence of massive companions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Two mechanisms for dust gap opening in protoplanetary discs

We identify two distinct physical mechanisms for dust gap opening by embedded planets in protoplanetary discs based on the symmetry of the drag-induced motion around the planet: I) A mechanism where low mass planets, that do not disturb the gas, open gaps in dust by tidal torques assisted by drag in the inner disc, but resisted by drag in the outer disc; and II) The usual, drag assisted, mechanism where higher mass planets create pressure maxima in the gas disc which the drag torque then acts to evacuate further in the dust. The first mechanism produces gaps in dust but not gas, while the second produces partial or total gas gaps which are deeper in the dust phase. Dust gaps do not necessarily indicate gas gaps.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A fast and explicit algorithm for simulating the dynamics of small dust grains with smoothed particle hydrodynamics

We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves solving one extra diffusion equation for the dust fraction in addition to the usual equations of hydrodynamics. This "diffusion approximation for dust" is valid when the dust stopping time is smaller than the computational timestep. We present a numerical implementation using Smoothed Particle Hydrodynamics (SPH) that is conservative, accurate and fast. It does not require any implicit timestepping and can be straightforwardly ported into existing 3D codes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

ALMA images of discs: are all gaps carved by planets?

Protoplanetary discs are now routinely observed and exoplanets, after the numerous indirect discoveries, are starting to be directly imaged. To better understand the planet formation process, the next step is the detection of forming planets or of signposts of young planets still in their disc, such as gaps. A spectacular example is the ALMA science verification image of HL Tau showing numerous gaps and rings in its disc. To study the observability of planet gaps, we ran 3D hydrodynamical simulations of a gas and dust disc containing a 5 M J gap-opening planet and characterised the spatial distribution of migrating, growing and fragmenting dust grains. We then computed the corresponding synthetic images for ALMA. For a value of the dust fragmentation threshold of 15 m s --1 for the collisional velocity, we identify for the first time a self-induced dust pile up in simulations taking fragmentation into account. This feature, in addition to the easily detected planet gap, causes a second apparent gap that could be mistaken for the signature of a second planet. It is therefore essential to be cautious in the interpretation of gap detections.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On planet formation in HL Tau

We explain the axisymmetric gaps seen in recent long-baseline observations of the HL Tau protoplanetary disc with the Atacama Large Millimetre/Submillimetre Array (ALMA) as being due to the different response of gas and dust to embedded planets in protoplanetary discs. We perform global, three dimensional dusty smoothed particle hydrodynamics calculations of multiple planets embedded in dust/gas discs which successfully reproduce most of the structures seen in the ALMA image. We find a best match to the observations using three embedded planets with masses of 0.2, 0.27 and 0.55 $M_{\rm J}$ in the three main gaps observed by ALMA, though there remain uncertainties in the exact planet masses from the disc model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The marriage of gas and dust

Dust-gas mixtures are the simplest example of a two fluid mixture. We show that when simulating such mixtures with particles or with particles coupled to grids a problem arises due to the need to resolve a very small length scale when the coupling is strong. Since this is occurs in the limit when the fluids are well coupled, we show how the dust-gas equations can be reformulated to describe a single fluid mixture. The equations are similar to the usual fluid equations supplemented by a diffusion equation for the dust-to-gas ratio or alternatively the dust fraction. This solves a number of numerical problems as well as making the physics clear.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two phase mixtures in SPH - A new approach

We present a new approach to simulating mixtures of gas and dust in smoothed particle hydrodynamics (SPH). We show how the two-fluid equations can be rewritten to describe a single-fluid 'mixture' moving with the barycentric velocity, with each particle carrying a dust fraction. We show how this formulation can be implemented in SPH while preserving the conservation properties (i.e. conservation of mass of each phase, momentum and energy). We also show that the method solves two key issues with the two fluid approach: it avoids over-damping of the mixture when the drag is strong and prevents a problem with dust particles becoming trapped below the resolution of the gas. We also show how the general one-fluid formulation can be simplified in the limit of strong drag (i.e. small grains) to the usual SPH equations plus a diffusion equation for the evolution of the dust fraction that can be evolved explicitly and does not require any implicit timestepping. We present tests of the simplified formulation showing that it is accurate in the small grain/strong drag limit. We discuss some of the issues we have had to solve while developing this method and finally present a preliminary application to dust settling in protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Diversity in the outcome of dust radial drift in protoplanetary discs

The growth of dust particles into planet embryos needs to circumvent the radial-drift barrier, i.e. the accretion of dust particles onto the central star by radial migration. The outcome of the dust radial migration is governed by simple criteria between the dust-to-gas ratio and the exponents p and q of the surface density and temperature power laws. The transfer of radiation provides an additional constraint between these quantities because the disc thermal structure is fixed by the dust spatial distribution. To assess which discs are primarily affected by the radial-drift barrier, we used the radiative transfer code MCFOST to compute the temperature structure of a wide range of disc models, stressing the particular effects of grain size distributions and vertical settling. We find that the outcome of the dust migration process is very sensitive to the physical conditions within the disc. For high dust-to-gas ratios (> 0.01) or flattened disc structures (H/R < 0.05), growing dust grains can efficiently decouple from the gas, leading to a high concentration of grains at a critical radius of a few AU. Decoupling of grains can occur at a large fraction (> 0.1) of the initial radius, for a dust-to-gas ratio greater than ~ 0.05. The exact value of the required dust-to-gas ratio for dust to stop its migration is strongly dependent on the disc temperature structure. Non growing dust grains are accreted for discs with flat surface density profiles (p<0.7) while they always remain in the disc if the surface density is steep enough (p>1.2). Both the presence of large grains and vertical settling tend to favour the accretion of non growing dust grains onto the central object, but it slows down the migration of growing dust grains. All the disc configurations are found to have favourable temperature profiles over most of the disc to retain their planetesimals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dust and gas mixtures with multiple grain species - a one-fluid approach

We derive the single-fluid evolution equations describing a mixture made of a gas phase and an arbitrary number of dust phases, generalising the approach developed in Laibe & Price (2014a). A generalisation for continuous dust distributions as well as analytic approximations for strong drag regimes are also provided. This formalism lays the foundation for numeri- cal simulations of dust populations in a wide range of astrophysical systems while avoiding limitations associated with a multiple-fluid treatment. The usefulness of the formalism is illustrated on a series of analytical problems, namely the dustybox, dustyshock and dustywave problems as well as the radial drift of grains and the streaming instability in protoplanetary discs. We find physical effects specific to the presence of several dust phases and multiple drag timescales, including non-monotonic evolution of the differential velocity between phases and increased efficiency of the linear growth of the streaming instability. Interestingly, it is found that under certain conditions, large grains can migrate outwards in protoplanetary discs. This may explain the presence of small pebbles at several hundreds of astronomical units from their central star.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dusty gas with one fluid

In this paper, we show how the two-fluid equations describing the evolution of a dust and gas mixture can be reformulated to describe a single fluid moving with the barycentric velocity of the mixture. This leads to evolution equations for the total density, momentum, the differential velocity between the dust and the gas phases and either the dust-to-gas ratio or the dust fraction. The equations are similar to the usual equations of gas dynamics, providing a convenient way to extend existing codes to simulate two-fluid mixtures without modifying the code architecture. Our approach avoids the inherent difficulties related to the standard approach where the two phases are separate and coupled via a drag term. In particular, the requirements of infinite spatial and temporal resolution as the stopping time tends to zero are no longer necessary. This means that both small and large grains can be straightforwardly treated with the same method, with no need for complicated implicit schemes. Since there is only one resolution scale the method also avoids the problem of unphysical trapping of one fluid (e.g. dust) below the resolution of the other. We also derive a simplified set of equations applicable to the case of strong drag/small grains, consisting of the standard fluid equations with a modified sound speed, plus an advection-diffusion equation for the dust-to-gas ratio. This provides a simple and fast way to evolve the mixture when the stopping time is smaller than the Courant timestep. We present a Smoothed Particle Hydrodynamics implementation in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dusty gas with one fluid in smoothed particle hydrodynamics

In a companion paper we have shown how the equations describing gas and dust as two fluids coupled by a drag term can be reformulated to describe the system as a single fluid mixture. Here we present a numerical implementation of the one-fluid dusty gas algorithm using Smoothed Particle Hydrodynamics (SPH). The algorithm preserves the conservation properties of the SPH formalism. In particular, the total gas and dust mass, momentum, angular momentum and energy are all exactly conserved. Shock viscosity and conductivity terms are generalised to handle the two-phase mixture accordingly. The algorithm is benchmarked against a comprehensive suit of problems: dustybox, dustywave, dustyshock and dustyoscill, each of them addressing different properties of the method. We compare the performance of the one-fluid algorithm to the standard two-fluid approach. The one-fluid algorithm is found to solve both of the fundamental limitations of the two- fluid algorithm: it is no longer possible to concentrate dust below the resolution of the gas (they have the same resolution by definition), and the spatial resolution criterion h < csts, required in two-fluid codes to avoid over-damping of kinetic energy, is unnecessary. Implicit time stepping is straightforward. As a result, the algorithm is up to ten billion times more efficient for 3D simulations of small grains. Additional benefits include the use of half as many particles, a single kernel and fewer SPH interpolations. The only limitation is that it does not capture multi-streaming of dust in the limit of zero coupling, suggesting that in this case a hybrid approach may be required.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Growing dust grains in protoplanetary discs - I. Radial drift with toy growth models

In a series of papers, we present a comprehensive analytic study of the global motion of growing dust grains in protoplanetary discs, addressing both the radial drift and the vertical settling of the particles. Here we study how the radial drift of dust particles is affected by grain growth. In a first step, toy models in which grain growth can either be constant, accelerate or decelerate are introduced. The equations of motion are analytically integrable and therefore the grains dynamics is easy to understand. The radial motion of growing grains is governed by the relative efficiency of the growth and migration processes which is expressed by the dimensionless parameter Lambda, as well as the exponents for the gas surface density and temperature profiles, denoted p and q respectively. When Lambda is of order unity, growth and migration are strongly coupled, providing the most efficient radial drift. For the toy models considered, grains pile up when -p+q+1/2<0. Importantly, we show the existence of a second process which can help discs to retain their solid materials. For accelerating growth, grains end up their migration at a finite radius, thus avoiding being accreted onto the central star.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Growing dust grains in protoplanetary discs - II. The Radial drift barrier problem

We aim to study the migration of growing dust grains in protoplanetary discs, where growth and migration are tightly coupled. This includes the crucial issue of the radial-drift barrier for growing dust grains. We therefore extend the study performed in Paper I, considering models for grain growth and grain dynamics where both the migration and growth rate depend on the grain size and the location in the disc. The parameter space of disc profiles and growth models is exhaustively explored. In doing so, interpretations for the grain motion found in numerical simulations are also provided. We find that a large number of cases is required to characterise entirely the grains radial motion, providing a large number of possible outcomes. Some of them lead dust particles to be accreted onto the central star and some of them don't. We find then that q<1 is required for discs to retain their growing particles, where q is the exponent of the radial temperature profile T(R) proportional to R^-q. Additionally, the initial dust-to gas ratio has to exceed a critical value for grains to pile up efficiently, thus avoiding being accreted onto the central star. Discs are also found to retain efficiently small dust grains regenerated by fragmentation. We show how those results are sensitive to the turbulent model considered. Even though some physical processes have been neglected, this study allows to sketch a scenario in which grains can survive the radial-drift barrier in protoplanetary discs as they grow.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Growing dust grains in protoplanetary discs - III. Vertical settling

We aim to derive a simple analytic model to understand the essential properties of vertically settling growing dust grains in laminar protoplanetary discs. Separating the vertical dynamics from the motion in the disc midplane, we integrate the equations of motion for both a linear and an exponential grain growth rate. Numerical integrations are performed for more complex growth models. We find that the settling efficiency depends on the value of the dimensionless parameter gamma, which characterises the relative efficiency of grain growth with respect to the gas drag. Since gamma is expected to be of order as the initial dust-to-gas ratio in the disc (of order 10^-2), grain growth enhances the energy dissipation of the dust particles and improve the settling efficiency in protoplanetary discs. This behaviour is mostly independent of the growth model considered as well as of the radial drift of the particles.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dusty gas with SPH - I. Algorithm and test suite

We present a new algorithm for simulating two-fluid gas and dust mixtures in Smoothed Particle Hydrodynamics (SPH), systematically addressing a number of key issues including the generalised SPH density estimate in multi-fluid systems, the consistent treatment of variable smoothing length terms, finite particle size, time step stability, thermal coupling terms and the choice of kernel and smoothing length used in the drag operator. We find that using double-hump shaped kernels improves the accuracy of the drag interpolation by a factor of several hundred compared to the use of standard SPH bell-shaped kernels, at no additional computational expense. In order to benchmark our algorithm, we have developed a comprehensive suite of standardised, simple test problems for gas and dust mixtures: dustybox, dustywave, dustyshock, dustysedov and dustydisc, the first three of which have known analytic solutions. We present the validation of our algorithm against all of these tests. In doing so, we show that the spatial resolution criterion Δ< cs ts is a necessary condition in all gas+dust codes that becomes critical at high drag (i.e. small stopping time ts) in order to correctly predict the dynamics. Implicit timestepping and the implementation of realistic astrophysical drag regimes are addressed in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dusty gas with SPH - II. Implicit timestepping and astrophysical drag regimes

In a companion paper (Laibe & Price 2011b), we have presented an algorithm for simulating two-fluid gas and dust mixtures in Smoothed Particle Hydrodynamics (SPH). In this paper, we develop an implicit timestepping method that preserves the exact conservation of the both linear and angular momentum in the underlying SPH algorithm, but unlike previous schemes, allows the iterations to converge to arbitrary accuracy and is suited to the treatment of non- linear drag regimes. The algorithm presented in Paper I is also extended to deal with realistic astrophysical drag regimes, including both linear and non-linear Epstein and Stokes drag. The scheme is benchmarked against the test suite presented in Paper I, including i) the analytic solutions of the dustybox problem and ii) solutions of the dustywave, dustyshock, dustysedov and dustydisc obtained with explicit timestepping. We find that the implicit method is 1- 10 times faster than the explicit temporal integration when the ratio r between the the timestep and the drag stopping time is 1 < r < 1000.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures

In this paper we present the analytic solutions for two test problems involving two-fluid mixtures of dust and gas in an astrophysical context. The solutions provide a means of benchmarking numerical codes designed to simulate the non-linear dynamics of dusty gas. The first problem, dustybox, consists of two interpenetrating homogeneous fluids moving with relative velocity difference. We provide exact solutions to the full non-linear problem for a range of drag formulations appropriate to astrophysical fluids (i.e., various prescriptions for Epstein and Stokes drag in different regimes). The second problem, dustywave consists of the propagation of linear acoustic waves in a two-fluid gas-dust mixture. We provide the analytic solution for the case when the two fluids are interacting via a linear drag term. Both test problems are simple to set up in any numerical code and can be run with periodic boundary conditions. The solutions we derive are completely general with respect to both the dust-to-gas ratio and the amplitude of the drag coefficient. A stability analysis of waves in a gas-dust system is also presented, showing that sound waves in an astrophysical dust-gas mixture are linearly stable.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Revisiting the "radial-drift barrier" of planet formation and its relevance in observed protoplanetary discs

Context. To form metre-sized pre-planetesimals in protoplanetary discs, growing grains have to decouple from the gas before they are accreted onto the central star during their phase of fast radial migration and thus overcome the so-called "radial-drift barrier" (often inaccurately referred to as the "metre-size barrier"). Aims. To predict the outcome of the radial motion of dust grains in protoplanetary discs whose surface density and temperature follow power-law profiles, with exponent p and q respectively. We investigate both the Epstein and the Stokes drag regimes which govern the motion of the dust. Methods. We analytically integrate the equations of motion obtained from perturbation analysis. We compare these results with those from direct numerical integration of the equations of motion. Then, using data from observed discs, we predict the fate of dust grains in real discs. Results. When a dust grain reaches the inner regions of the disc, the acceleration due to the increase of the pressure gradient is counterbalanced by the increase of the gas drag. We find that most grains in the Epstein (resp. the Stokes) regime survive their radial migration if-p+q+1/2 \leq0 (resp. if q\leq 2/3). The majority of observed discs satisfies both-p+q+ 1/2 \leq0 and q\leq 2/3: a large fraction of both their small and large grains remain in the disc, for them the radial drift barrier does not exist.

0 citations0 reviewsNo code linkedOpen access