2019/07/27 by Benjamin K. Tapley, Tapley, Benjamin K., Helge I. Andersson +5
Computer Science · Mathematics · Physics and Astronomy · #Computational Physics (physics.comp-ph) #FOS: Mathematics #FOS: Physical sciences #Fluid Dynamics (physics.flu-dyn) #Numerical Analysis (math.NA) #cs.NA #math.NA #physics.comp-ph #physics.flu-dyn
paper · pdf · doi:10.48550/arxiv.1907.11936
arxiv created 2021/02/25 · arxiv updated 2021/02/26
A geometric numerical method for simulating suspensions of spherical and non-spherical particles with Stokes drag is proposed. The method combines divergence-free matrix-valued radial basis function interpolation of the fluid velocity field with a splitting method integrator that preserves the sum of the Lyapunov spectrum while mimicking the centrifuge effect of the exact solution. We discuss how breaking the divergence-free condition in the interpolation step can erroneously affect how the volume of the particulate phase evolves under numerical methods. The methods are tested on suspensions of 104 particles evolving in discrete cellular flow field. The results are that the proposed geometric methods generate more accurate and cost-effective particle distributions compared to conventional methods.