A new multifluid method for dusty astrophysical flows. Application to turbulent protostellar collapses
arXiv:2508.19127 · doi:10.1051/0004-6361/202554662
Abstract
Stars and planets form in collapsing clouds of gas and dust. The presence of dust grains and their local distribution play a significant role throughout the protostellar sequence, from the thermodynamics and the chemistry of molecular clouds to the opacity of collapsing protostellar cores and the coupling between the gas and the magnetic field and down to planet formation in young and evolved disks. We aim to simulate the dynamics of the dust, considering the whole range of grain sizes, from few nanometers to millimeters. We implemented a neutral pressureless multifluid that samples the dust size distribution in the RAMSES code. This multifluid is dynamically coupled to the gas via a drag source term and self-gravity, relying on the Eulerian approach. We designed a Riemann solver for the gas and dust mixture that prevents unphysical dust-to-gas ratio variations for well-coupled grains. We illustrated the capacities of the code by performing simulations of a protostellar collapse down to the formation of a first hydrostatic core, both for small and large dust grains. Grains over 100 microns significantly decouple from the gas. The spatial maps and the probability density functions indicate that dust enrichment within the first hydrostatic core and in some locations of the envelope increases as a function of the grain size and the level of initial turbulence. Thanks to the novel Riemann solver, we recovered the terminal velocity regime, even at low resolution. Moreover, we successfully extended it to regimes where the grain inertia matters. The multifluid module performs the coupling between the dust and the gas self-consistently all through the dynamical scales. The dust enrichment in the first hydrostatic core and the envelope have been revised here, assuming the initial turbulence and grain sizes. This enables us to probe new potential conditions for planet formation.
arXiv commentAccepted for publication in Astronomy & Astrophysics
References in corpus (42)
- Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
- Streaming Instabilities in Protoplanetary Disks
- The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
- Closed-form expressions for particle relative velocities induced by turbulence
- A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
- The global dust SED: Tracing the nature and evolution of dust with DustEM
- FARGO3D: A new GPU-oriented MHD code
- Millimeter-wave polarization of protoplanetary disks due to dust scattering
- Why do protoplanetary disks appear not massive enough to form the known exoplanet population?
- The initial conditions of star formation in the Ophiuchus main cloud: Kinematics of the protocluster condensations
- The Ubiquity of Micrometer-Sized Dust Grains in the Dense Interstellar Medium
- A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
- Can non-ideal magnetohydrodynamics solve the magnetic braking catastrophe?
- Radiation hydrodynamics with Adaptive Mesh Refinement and application to prestellar core collapse. I Methods
- Protostellar Disk Formation Enabled by Removal of Small Dust Grains
- Resonant Drag Instabilities in protoplanetary disks: the streaming instability and new, faster-growing instabilities
- Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D
- Dusty gas with SPH - I. Algorithm and test suite
- Dusty gas with one fluid
- MPI-AMRVAC 3.0: updates to an open-source simulation framework
- Streaming Instability for Particle-Size Distributions
- A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations
- Low dust emissivities and radial variations in the envelopes of Class 0 protostars: a signature of early grain growth?
- Simulations of protostellar collapse using multigroup radiation hydrodynamics. II. The second collapse
- Incorporating Ambipolar and Ohmic Diffusion in the AMR MHD code RAMSES
- DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures
- Protostellar collapse: the conditions to form dust rich protoplanetary disks
- Protostellar collapse simulations in spherical geometry with dust coagulation and fragmentation
- Indirect evidence of significant grain growth in young protostellar envelopes from polarized dust emission
- MULTIGRAIN: A smoothed particle hydrodynamics algorithm for multiple small dust grains and gas
- Small dust grain dynamics on adaptive mesh-refinement grids. I. Methods
- A Multi-Fluid Dust Module in Athena++: Algorithms and Numerical Tests
- THEMIS 2.0: A self-consistent model for dust extinction, emission, and polarisation
- Impact of the Hall effect in star formation and the issue of angular momentum conservation
- Multigroup radiation hydrodynamics with flux-limited diffusion and adaptive mesh refinement
- IDEFIX: a versatile performance-portable Godunov code for astrophysical flows
- Grain growth for astrophysics with Discontinuous Galerkin schemes
- Three-temperature radiation hydrodynamics with PLUTO: Tests and applications to protoplanetary disks
- The birth and early evolution of a low mass protostar
- FAUST IX. Multi-band, multi-scale dust study of L1527 IRS. Evidence for dust properties variations within the envelope of a Class 0/I YSO
- Dynamics of Dusty Vortices I: Extensions and limitations of the terminal velocity approximation
- Dust dynamics in RAMSES -- I. Methods and turbulent acceleration