A numerical model of superspreading surfactants on hydrophobic surface

Many contributions significantly on experimental and mathematical studies are made to understand the mechanism of superspreading. Only few numerical methods have been proposed which solve the system of equations with soluble and insoluble surfactants. Among them, we propose a computational fluid dynamics model, based on the volume of fluid technique, with the piecewise linear interface calculation method. Interface reconstruction is applied to simulate the time evolution of the dynamics of drop spreading of surfactants on a thin water layer. We have allowed the occurrence of both the regimes relating to a series of trisiloxane (M(D′EnOH)M), sodium dodecyl sulphate, and Tergitol NP10 surfactants drop on a thin water layer with the influence of Marangoni stress. The numerical data seem consistent with those experimental for both regimes. It validates predictions for the spreading exponent in which the law of the radius of the circular area covered by the surfactant grows as tα, where 0 < α < 1. The comparison of the numerical and experimental predictions by Lee et al. [“Spreading of trisiloxanes over thin aqueous layers,” Colloid J. 71, 365–369 (2009)] is well represented in both regimes. The numerical study confirms that the spreading rates during the first stage increase as the solubility increases. This finding suggests that the model is adequate for describing the spreading of surfactants on thin fluid layers.