A splitting integration scheme for the SPH simulation of concentrated particle suspensions
A splitting integration scheme for the SPH simulation of concentrated particle suspensions
复制标题
DOI：
10.1016/j.cpc.2013.08.015
复制
发表时间：
2014
期刊：
影响因子：
--
通讯作者：
X. Bian;M. Ellero
中科院分区：
文献类型：
--
作者：
X. Bian;M. Ellero
由于短程润滑力的发散行为，模拟悬浮液中几乎接触的固体颗粒是一项具有挑战性的任务，这对显式积分方案造成了严重的时间步长限制。这种普遍的困难严重限制了浓缩悬浮液模拟的总持续时间。受[S. Litvinov，M. Ellero，X.Y.胡，NA亚当斯，J. 229（2010）5457-5464]对于高耗散流体的模拟，我们在这项工作中提出了用于直接模拟悬浮在牛顿液体中的固体颗粒的分裂积分方案。该方案分离了作用在固体颗粒上的不同力的贡献。特别是，中期和长期的多体水动力，这是计算从离散的Navier-Stokes方程使用光滑粒子流体动力学（SPH）方法，考虑到使用显式积分;对于短程润滑力，成对相互作用的固体颗粒的速度通过迭代地扫过所有相邻对而隐含地更新，直到获得解的收敛。通过使用分裂积分，模拟可以稳定和有效地运行到非常大的固体颗粒浓度。此外，所提出的计划是不限于SPH方法在这里提出的，但可以很容易地应用到其他模拟技术用于颗粒悬浮液。
Simulating nearly contacting solid particles in suspension is a challenging task due to the diverging behavior of short-range lubrication forces, which pose a serious time-step limitation for explicit integration schemes. This general difficulty limits severely the total duration of simulations of concentrated suspensions. Inspired by the ideas developed in [S. Litvinov, M. Ellero, X.Y. Hu, N.A. Adams, J. Comput. Phys. 229 (2010) 5457–5464] for the simulation of highly dissipative fluids, we propose in this work a splitting integration scheme for the direct simulation of solid particles suspended in a Newtonian liquid. The scheme separates the contributions of different forces acting on the solid particles. In particular, intermediate- and long-range multi-body hydrodynamic forces, which are computed from the discretization of the Navier–Stokes equations using the smoothed particle hydrodynamics (SPH) method, are taken into account using an explicit integration; for short-range lubrication forces, velocities of pairwise interacting solid particles are updated implicitly by sweeping over all the neighboring pairs iteratively, until convergence in the solution is obtained. By using the splitting integration, simulations can be run stably and efficiently up to very large solid particle concentrations. Moreover, the proposed scheme is not limited to the SPH method presented here, but can be easily applied to other simulation techniques employed for particulate suspensions.