A multicomponent pseudopotential lattice Boltzmann model for liquid–liquid systems with soluble surfactants

Journal Article (2025)
Author(s)

Mohammad Pourtousi (University of Limerick)

Arman Safdari (University of Limerick)

Orest Shardt (University of Limerick)

H.E.A. van den Akker (University of Limerick, TU Delft - ChemE/Transport Phenomena)

Research Group
ChemE/Transport Phenomena
DOI related publication
https://doi.org/10.1016/j.ijmultiphaseflow.2025.105255
More Info
expand_more
Publication Year
2025
Language
English
Research Group
ChemE/Transport Phenomena
Volume number
189
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

The intermolecular interactions in the pseudo-potential lattice Boltzmann (PPLB) method can readily be extended to more than two components. We report about a three-component PPLB approach to explore whether the effect of a surfactant could be included in describing droplet behaviour in (liquid–liquid) emulsions. The two main liquid components are taken to follow the Carnahan-Starling equation of state (EoS), while the surfactant obeys an ideal EoS. We investigate the nature of the phases present at equilibrium and the dependence of the interfacial tension between the two liquid phases on the amount of surfactant. The response of a droplet subjected to simple shear is investigated in the absence and the presence of a surfactant. Our exploratory simulations show how during droplet deformation the surfactant re-distributes itself due to the action of the shear and flows towards the far ends of the deformed droplet, up to the moment the droplet breaks up. This inhomogeneous surfactant distribution along the interface increases the shear rate that is needed for droplet breakup such that the critical capillary number for breakup increases and the breakup process is delayed. The simulations also reveal the detailed flow fields inside and outside the deforming droplet.