Fast method to simulate dynamics of two-phase medium with intense interaction between phases by smoothed particle hydrodynamics: gas-dust mixture with polydisperse particles, linear drag, one-dimensional tests

Abstract

To simulate the dynamics of fluid with polydisperse particles on macroscale level, one has to solve hydrodynamic equations with several relaxation terms, representing momentum transfer from fluid to particles and vice versa. For small particles, velocity relaxation time (stopping time) can be much shorter than dynamical time of fluid \textcolor{red}{that makes} this problem stiff and thus computationally expensive. We present a new fast method for computing several stiff drag terms in two-phase polydisperse medium with Smoothed Particle Hydrodynamics (SPH). In our implementation, fluid and every fraction of dispersed phase are simulated with different sets of particles. The method is based on (1) linear interpolation of velocity values in drag terms, (2) implicit approximation of drag terms that conserves momentum with machine precision\textcolor{red}{,} and (3) solution of system of $N$ linear algebraic equations with $O(N^2)$ arithmetic operation instead of $O(N^3)$. We studied the properties of the proposed method on one-dimensional problems with known solutions. We found that we can obtain acceptable accuracy of the results with numerical resolution independent of short stopping time values. All simulation results discussed in the paper are obtained with \textcolor{red}{open source} software.

25 pages, 10 figure

References in corpus (25)

• SPLASH: An interactive visualisation tool for Smoothed Particle Hydrodynamics simulations
• Concentrating small particles in protoplanetary disks through the streaming instability
• 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
• Including Dust Coagulation in Hydrodynamic Models of Protoplanetary Disks: Dust Evolution in the Vicinity of a Jupiter-mass Planet
• The early evolution of viscous and self-gravitating circumstellar disks with a dust component
• Dust and gas mixtures with multiple grain species - a one-fluid approach
• Dusty gas with SPH - II. Implicit timestepping and astrophysical drag regimes
• Asymptotic-Preserving methods and multiscale models for plasma physics
• DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures
• Grand challenges in protoplanetary disc modelling
• Two fluid dust and gas mixtures in SPH: A semi-implicit approach
• MULTIGRAIN: A smoothed particle hydrodynamics algorithm for multiple small dust grains and gas
• A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows VI: Dilute Disperse Gas-Particle Multiphase System
• Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction
• Smoothed particle hydrodynamics simulations of gas and dust mixtures
• A Hybrid Scheme for Gas-Dust Systems Stiffly Coupled via Viscous Drag
• Analysis of Numerical Algorithms for Computing Rapid Momentum Transfers between the Gas and Dust in Simulations of Circumstellar Disks
• Two-fluid dust and gas mixtures in smoothed particle hydrodynamics II: an improved semi-implicit approach
• Simulations of Dynamical Gas-Dust Circumstellar Disks: Going Beyond the Epstein Regime
• Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag
• A solution to the overdamping problem when simulating dust-gas mixtures with smoothed particle hydrodynamics
• Monodisperse gas-solid mixtures with intense interphase interaction in two-fluid Smoothed Particle Hydrodynamics
• Development and application of fast methods for computing momentum transfer between gas and dust in supercomputer simulation of planet formation