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

Haihu Liu, Yan Ba, Lei Wu, Zhen Li, Guang Xi, Yonghao Zhang

Research output: Contribution to journalArticle

18 Citations (Scopus)

Abstract

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.
LanguageEnglish
Pages381-412
Number of pages32
JournalJournal of Fluid Mechanics
Volume837
Early online date21 Dec 2017
DOIs
Publication statusPublished - 25 Feb 2018

Fingerprint

Finite difference method
Surface active agents
surfactants
Surface tension
interfacial tension
Shear flow
shear flow
Reynolds number
Coalescence
Microfluidics
coalescing
convection-diffusion equation
Capillarity
two phase flow
Equations of state
Two phase flow
Contacts (fluid mechanics)
equations of state

Keywords

  • droplet dynamics
  • surfactants
  • Boltzmann equation

Cite this

@article{ee7da70fbe8447caa4cb7ce1bb57a989,
title = "A hybrid lattice Boltzmann and finite difference method for droplet dynamics with insoluble surfactants",
abstract = "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.",
keywords = "droplet dynamics, surfactants, Boltzmann equation",
author = "Haihu Liu and Yan Ba and Lei Wu and Zhen Li and Guang Xi and Yonghao Zhang",
year = "2018",
month = "2",
day = "25",
doi = "10.1017/jfm.2017.859",
language = "English",
volume = "837",
pages = "381--412",
journal = "Journal of Fluid Mechanics",
issn = "0022-1120",
publisher = "Cambridge University Press",

}

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.

In: Journal of Fluid Mechanics, Vol. 837, 25.02.2018, p. 381-412.

Research output: Contribution to journalArticle

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 -