A conservative SPH method for surfactant dynamics

In this paper, a Lagrangian particle method is proposed for the simulation of multiphase flows with surfactant. The model is based on the multiphase smoothed particle hydrodynamics (SPH) framework of Hu and Adams (2006) [1]. Surface-active agents (surfactants) are incorporated into our method by a scalar quantity describing the local concentration of molecules in the bulk phase and on the interface. The surfactant dynamics are written in conservative form, thus global mass of surfactant is conserved exactly. The transport model of the surfactant accounts for advection and diffusion. Within our method, we can simulate insoluble surfactant on an arbitrary interface geometry as well as interfacial transport such as adsorption or desorption. The flow-field dynamics and the surfactant dynamics are coupled through a constitutive equation, which relates the local surfactant concentration to the local surface-tension coefficient. Hence, the surface-tension model includes capillary and Marangoni-forces. The present numerical method is validated by comparison with analytic solutions for diffusion and for surfactant dynamics. More complex simulations of an oscillating bubble, the bubble deformation in a shear flow, and of a Marangoni-force driven bubble show the capabilities of our method to simulate interfacial flows with surfactants.