A multi-fluid approach for polydisperse pebble accretion: From particles to fluids, establishing the multifluid framework

Abstract

Pebble accretion offers an efficient pathway to form planets, driven by a constant supply of inward drifting mass and an accretion efficiency enhanced by gas drag. While most studies assume a single pebble size (monodisperse), real discs contain a range of sizes (polydisperse) that drift, interact, and accrete at different rates. We aim to model polydisperse pebble accretion with a fluid approach, validating the method and exploring how gas disc evolution, solid-to-gas back-reaction, and a polydisperse size distribution affect growth. We used FARGO3D, modified to allow pebble accretion, to run 2D hydrodynamic simulations in a global disc with multiple pebble species representing an underlying continuous pebble size distribution. With our multi-fluid approach, we find values for pebble accretion efficiency consistent with earlier studies for a static gas disc. This confirms that our approach gives an accurate representation of pebble accretion. Evolving the gas disc, we find lower efficiencies compared to an unperturbed gas disc for high Stokes numbers ($\gtrsim 0.3$) and higher efficiencies for smaller Stokes numbers ($\lesssim 0.3$). This effect increases for higher planet masses. The accretion rate is mostly dominated by the highest Stokes numbers in our parameter study ($\mathrm{St}\in [10^{-2},10^0]$). The ratio we find between the polydisperse and monodisperse pebble accretion rates is higher than previous estimations. We constructed a multi-fluid model framework capable of accurately simulating polydisperse pebble accretion consistent with previous studies. This framework offers advantages for simulating higher planet masses and for modelling multiple pebble species coupled to the gas. We find that the protoplanet's perturbation of the gas-disc lowers the accretion rate when assuming an MRN-distribution of solids.

Accepted for publication in Astronomy and Astrophysics. 14 pages, 12 figures

References in corpus (20)

• Array Programming with NumPy
• PLUTO: a Numerical Code for Computational Astrophysics
• Coagulation, fragmentation and radial motion of solid particles in protoplanetary disks
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Separating gas-giant and ice-giant planets by halting pebble accretion
• Dust flow in gas disks in the presence of embedded planets
• Global MHD simulations of stratified and turbulent protoplanetary discs. II. Dust settling
• Dust growth and evolution in protoplanetary disks
• Predicting the observational signature of migrating Neptune-sized planets in low-viscosity disks
• Eccentricity excitation and merging of planetary embryos heated by pebble accretion
• An Analytical Theory for the Growth from Planetesimals to Planets by Polydisperse Pebble Accretion
• Influences of three-dimensional gas flow induced by protoplanets on pebble accretion --$\mathrm{I}$. shear regime
• Torques felt by solid accreting planets
• Pebble-driven migration of low-mass planets in the 2D regime of pebble accretion
• Dynamics of dust grains in turbulent molecular clouds. Conditions for decoupling and limits of different numerical implementations
• Polydisperse Formation of Planetesimals: The dust size distribution in clumps
• On the evolution of pebble-accreting planets in evolving protoplanetary discs
• Accretion of aerodynamically large pebbles
• Resonant Drag Instabilities for Polydisperse Dust, I. The Acoustic Resonant Drag Instability
• Resonant drag instabilities for polydisperse dust. II. The streaming and settling instabilities