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

TL;DR

This paper introduces a fast, implicit SPH-based method for simulating two-phase fluid dynamics with polydisperse particles, efficiently handling stiff drag interactions and short stopping times in one-dimensional tests.

Contribution

A novel implicit SPH method that efficiently computes stiff drag terms in polydisperse two-phase flows, reducing computational complexity and maintaining accuracy.

Findings

Method achieves accurate results with short stopping times.

Simulation resolution is independent of particle stopping times.

Open source implementation available.

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.