A Hybrid Scheme for Gas-Dust Systems Stiffly Coupled via Viscous Drag
arXiv:1001.4794 · doi:10.1016/j.jcp.2010.01.034
Abstract
We present a stable and convergent method for studying a system of gas and dust, coupled through viscous drag in both non-stiff and stiff regimes. To account for the effects of dust drag in the update of the fluid quantities, we employ a fluid description of the dust component and study the modified gas-dust hyperbolic system following the approach in Miniati & Colella (2007). In addition to two entropy waves for the gas and dust components, respectively, the extended system includes three waves driven partially by gas pressure and partially by dust drift, which, in the limit of vanishing coupling, tend to the two original acoustic waves and the unhindered dust streaming. Based on this analysis we formulate a predictor step providing first order accurate reconstruction of the time-averaged state variables at cell interfaces, whence a second order accurate estimate of the conservative fluxes can be obtained through a suitable linearized Riemann solver. The final source term update is carried out using a one-step, second order accurate, L-stable, predictor corrector asymptotic method (the alpha-QSS method suggested by Mott et. al. 2000). This procedure completely defines a two-fluid method for gas-dust system. Using the updated fluid solution allows us to then advance the individual particle solutions, including self-consistently the time evolution of the gas velocity in the estimate of the drag force. This is done with a suitable particle scheme also based on the alpha-QSS method. A set of benchmark problems shows that our method is stable and convergent. When dust is modeled as a fluid (two-fluid) second order accuracy is achieved in both stiff and non-stiff regimes, whereas when dust is modeled with particles (hybrid) second order is achieved in the non-stiff regime and first order otherwise.
41 pages, 3 figures, 14 tables, accepted to J. Comp. Phys.
References in corpus (4)
- Global MHD simulations of stratified and turbulent protoplanetary discs. II. Dust settling
- Block Structured Adaptive Mesh and Time Refinement for Hybrid, Hyperbolic + N-body Systems
- SPH simulations of grain growth in protoplanetary disks
- Three-Dimensional Simulations of Kelvin-Helmholtz Instability in Settled Dust Layers in Protoplanetary Disks
Cited by in corpus (28)
- Adding particle collisions to the formation of asteroids and Kuiper belt objects via streaming instabilities
- Particle-Gas Dynamics with Athena: Method and Convergence
- High-resolution simulations of planetesimal formation in turbulent protoplanetary discs
- 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
- Simulating galactic dust grain evolution on a moving mesh
- Resistive Magnetic Field Generation at Cosmic Dawn
- On the dynamics of planetesimals embedded in turbulent protoplanetary discs
- Streaming Instability with Multiple Dust Species: I. Favourable Conditions for the Linear Growth
- DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures
- Grand challenges in protoplanetary disc modelling
- Linear growth of streaming instability in pressure bumps
- MULTIGRAIN: A smoothed particle hydrodynamics algorithm for multiple small dust grains and gas
- Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction
- Constrained-Transport Magnetohydrodynamics with Adaptive-Mesh-Refinement in CHARM
- Analysis of Numerical Algorithms for Computing Rapid Momentum Transfers between the Gas and Dust in Simulations of Circumstellar Disks
- A Particle Module for the PLUTO Code: III -- Dust
- Dust entrainment in photoevaporative winds: The impact of X-rays
- Effect of dust on Kelvin-Helmholtz instabilities
- Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag
- Dusty disc-planet interaction with dust-free simulations
- Spontaneous concentrations of solids through two-way drag forces between gas and sedimenting particles
- A solution to the overdamping problem when simulating dust-gas mixtures with smoothed particle hydrodynamics
- 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
- Chemo-dynamical simulations of dwarf galaxy evolution
- Observation of aerodynamic instability in the flow of a particle stream in a dilute gas
- A Staggered Semi-Analytic Method for Simulating Dust Grains Subject to Gas Drag