### Abstract

Our method is first validated for the surfactant-laden droplet deformation in a three- dimensional (3D) extensional flow and a 2D shear flow, and then applied to investigate the effect of surfactants on droplet dynamics in a 3D shear flow. Numerical results show that, at low capillary numbers, surfactants increase droplet deformation, due to reduced interfacial tension by the average surfactant concentration, and non-uniform effects from non-uniform capillary pressure and Marangoni stresses. The role of surfactants on critical capillary number (Ca

_{cr}) of droplet breakup is investigated for various confinements (defined as the ratio of droplet diameter to wall separation) and Reynolds numbers. For clean droplets, Ca

_{cr}first decreases and then increases with confinement, and the minimum value of Ca

_{cr}is reached at the confinement of 0.5; for surfactant-laden droplets, Ca

_{cr}exhibits the same variation in trend for the confinements lower than 0.7, but for higher confinements, Ca

_{cr}is almost a constant. The presence of surfactants decreases Ca

_{cr}for each confinement, and the decrease is also attributed to the reduction in average interfacial tension and non-uniform effects, which are found to prevent droplet breakup at low confinements but promote breakup at high confinements. In either clean or surfactant-laden cases, Ca

_{cr}first keeps almost unchanged and then decreases with increasing Reynolds number, and a higher confinement or Reynolds number favors ternary breakup. Finally, we study the collision of two equal-sized droplets in a shear flow in both surfactant-free and surfactant-contaminated systems with the same effective capillary numbers. It is identified that the non-uniform effects in near-contact interfacial region immobilize the interfaces when two droplets are approaching each other and thus inhibit their coalescence.

Language | English |
---|---|

Pages | 381-412 |

Number of pages | 32 |

Journal | Journal of Fluid Mechanics |

Volume | 837 |

Early online date | 21 Dec 2017 |

DOIs | |

Publication status | Published - 25 Feb 2018 |

### Fingerprint

### Keywords

- droplet dynamics
- surfactants
- Boltzmann equation

### Cite this

*Journal of Fluid Mechanics*,

*837*, 381-412. https://doi.org/10.1017/jfm.2017.859

}

*Journal of Fluid Mechanics*, vol. 837, pp. 381-412. https://doi.org/10.1017/jfm.2017.859

**A hybrid lattice Boltzmann and finite difference method for droplet dynamics with insoluble surfactants.** / Liu, Haihu; Ba, Yan; Wu, Lei; Li, Zhen; Xi, Guang; Zhang, Yonghao.

Research output: Contribution to journal › Article

TY - JOUR

T1 - A hybrid lattice Boltzmann and finite difference method for droplet dynamics with insoluble surfactants

AU - Liu, Haihu

AU - Ba, Yan

AU - Wu, Lei

AU - Li, Zhen

AU - Xi, Guang

AU - Zhang, Yonghao

PY - 2018/2/25

Y1 - 2018/2/25

N2 - Droplet dynamics in microfluidic applications is significantly influenced by surfactants. It remains a research challenge to model and simulate droplet behavior including deformation, breakup and coalescence, especially in the confined microfluidic environment. Here we propose a hybrid method to simulate interfacial flows with insoluble surfactants. The immiscible two-phase flow is solved by an improved lattice Boltzmann color-gradient model that incorporates a Marangoni stress resulting from non-uniform interfacial tension, while the convection-diffusion equation which describes the evolution of surfactant concentration in the entire fluid domain is solved by a finite difference method. The lattice Boltzmann and finite difference simulations are coupled through an equation of state, which describes how surfactant concentration influences interfacial tension.Our method is first validated for the surfactant-laden droplet deformation in a three- dimensional (3D) extensional flow and a 2D shear flow, and then applied to investigate the effect of surfactants on droplet dynamics in a 3D shear flow. Numerical results show that, at low capillary numbers, surfactants increase droplet deformation, due to reduced interfacial tension by the average surfactant concentration, and non-uniform effects from non-uniform capillary pressure and Marangoni stresses. The role of surfactants on critical capillary number (Cacr) of droplet breakup is investigated for various confinements (defined as the ratio of droplet diameter to wall separation) and Reynolds numbers. For clean droplets, Cacr first decreases and then increases with confinement, and the minimum value of Cacr is reached at the confinement of 0.5; for surfactant-laden droplets, Cacr exhibits the same variation in trend for the confinements lower than 0.7, but for higher confinements, Cacr is almost a constant. The presence of surfactants decreases Cacr for each confinement, and the decrease is also attributed to the reduction in average interfacial tension and non-uniform effects, which are found to prevent droplet breakup at low confinements but promote breakup at high confinements. In either clean or surfactant-laden cases, Cacr first keeps almost unchanged and then decreases with increasing Reynolds number, and a higher confinement or Reynolds number favors ternary breakup. Finally, we study the collision of two equal-sized droplets in a shear flow in both surfactant-free and surfactant-contaminated systems with the same effective capillary numbers. It is identified that the non-uniform effects in near-contact interfacial region immobilize the interfaces when two droplets are approaching each other and thus inhibit their coalescence.

AB - Droplet dynamics in microfluidic applications is significantly influenced by surfactants. It remains a research challenge to model and simulate droplet behavior including deformation, breakup and coalescence, especially in the confined microfluidic environment. Here we propose a hybrid method to simulate interfacial flows with insoluble surfactants. The immiscible two-phase flow is solved by an improved lattice Boltzmann color-gradient model that incorporates a Marangoni stress resulting from non-uniform interfacial tension, while the convection-diffusion equation which describes the evolution of surfactant concentration in the entire fluid domain is solved by a finite difference method. The lattice Boltzmann and finite difference simulations are coupled through an equation of state, which describes how surfactant concentration influences interfacial tension.Our method is first validated for the surfactant-laden droplet deformation in a three- dimensional (3D) extensional flow and a 2D shear flow, and then applied to investigate the effect of surfactants on droplet dynamics in a 3D shear flow. Numerical results show that, at low capillary numbers, surfactants increase droplet deformation, due to reduced interfacial tension by the average surfactant concentration, and non-uniform effects from non-uniform capillary pressure and Marangoni stresses. The role of surfactants on critical capillary number (Cacr) of droplet breakup is investigated for various confinements (defined as the ratio of droplet diameter to wall separation) and Reynolds numbers. For clean droplets, Cacr first decreases and then increases with confinement, and the minimum value of Cacr is reached at the confinement of 0.5; for surfactant-laden droplets, Cacr exhibits the same variation in trend for the confinements lower than 0.7, but for higher confinements, Cacr is almost a constant. The presence of surfactants decreases Cacr for each confinement, and the decrease is also attributed to the reduction in average interfacial tension and non-uniform effects, which are found to prevent droplet breakup at low confinements but promote breakup at high confinements. In either clean or surfactant-laden cases, Cacr first keeps almost unchanged and then decreases with increasing Reynolds number, and a higher confinement or Reynolds number favors ternary breakup. Finally, we study the collision of two equal-sized droplets in a shear flow in both surfactant-free and surfactant-contaminated systems with the same effective capillary numbers. It is identified that the non-uniform effects in near-contact interfacial region immobilize the interfaces when two droplets are approaching each other and thus inhibit their coalescence.

KW - droplet dynamics

KW - surfactants

KW - Boltzmann equation

UR - https://www.cambridge.org/core/journals/journal-of-fluid-mechanics

U2 - 10.1017/jfm.2017.859

DO - 10.1017/jfm.2017.859

M3 - Article

VL - 837

SP - 381

EP - 412

JO - Journal of Fluid Mechanics

T2 - Journal of Fluid Mechanics

JF - Journal of Fluid Mechanics

SN - 0022-1120

ER -