A self-regularising DVM–PSE method for the modelling of diffusion in particle methods

The modelling of diffusive terms in particle methods is a delicate matter and several models have been proposed in the literature. The Diffusion Velocity Method (DVM) consists in rewriting these terms in an advective way, thus defining a so-called diffusion velocity. In addition to the actual velocity, it is used to compute the particles displacement. On the other hand, the well-known and commonly used Particle Strength Exchange method (PSE) uses an approximation of the Laplacian operator in order to model diffusion. This approximation is based on an exchange of particles strength. Although DVM is particularly well suited to particle methods since it preserves their Lagrangian aspect, its major drawback stems in the fact that it suffers from severe singular behaviours. This paper intends to give insights and ideas for coping with these issues, based on an exact decomposition of the diffusion coefficient allowing a hybrid DVM–PSE treatment of diffusive terms