Modelling thermocapillary migration of a microfluidic droplet on a solid surface

Haihu Liu, Yonghao Zhang

Research output: Contribution to journalArticle

19 Citations (Scopus)

Abstract

A multiphase lattice Boltzmann model is developed to simulate immiscible thermocapillary flows with the presence of fluid–surface interactions. In this model, interfacial tension force and Marangoni stress are included by introducing a body force term based on the concept of continuum surface force, and phase segregation is achieved using the recolouring algorithm proposed by Latva-Kokko and Rothman. At a solid surface, fluid–surface interactions are modelled by a partial wetting boundary condition that uses a geometric formulation to specify the contact angle, and a colour-conserving boundary closure scheme to improve the numerical accuracy and suppress spurious velocities at the contact line. An additional convection–diffusion equation is solved by the passive scalar approach to obtain the temperature field, which is coupled to the hydrodynamic equations through an equation of state. This model is first validated by simulations of static contact angle and dynamic capillary intrusion process when a constant interfacial tension is considered. It is then used to simulate the thermocapillary migration of a microfluidic droplet on a horizontal solid surface subject to a uniform temperature gradient. We for the first time demonstrate numerically that the droplet motion undergoes two different states depending on the surface wettability: the droplet migrates towards the cooler regions on hydrophilic surfaces but reverses on hydrophobic surfaces. Decreasing the viscosity ratio can enhance the intensity of thermocapillary vortices, leading to an increase in migration velocity. The contact angle hysteresis, i.e., the difference between the advancing and receding contact angles, is always positive regardless of the contact angle and viscosity ratio. The contact angle hysteresis and the migration velocity both first decrease and then increase with the contact angle, and their minimum values occur at the contact angle of 90 degrees.
LanguageEnglish
Pages37-53
Number of pages17
JournalJournal of Computational Physics
Volume280
Early online date28 Sep 2014
DOIs
Publication statusPublished - 1 Jan 2015

Fingerprint

thermocapillary migration
Microfluidics
solid surfaces
Contact angle
surface reactions
Surface tension
Hysteresis
Wetting
interfacial tension
hysteresis
Viscosity
viscosity
convection-diffusion equation
Fluids
hydrodynamic equations
fluids
wettability
coolers
Equations of state
intrusion

Keywords

  • thermocapillary migration
  • lattice Boltzmann method
  • multiphase flow
  • wetting boundary condition
  • microfluidics
  • droplet dynamics

Cite this

@article{13dccad5c77e4a83b006ee211b08583a,
title = "Modelling thermocapillary migration of a microfluidic droplet on a solid surface",
abstract = "A multiphase lattice Boltzmann model is developed to simulate immiscible thermocapillary flows with the presence of fluid–surface interactions. In this model, interfacial tension force and Marangoni stress are included by introducing a body force term based on the concept of continuum surface force, and phase segregation is achieved using the recolouring algorithm proposed by Latva-Kokko and Rothman. At a solid surface, fluid–surface interactions are modelled by a partial wetting boundary condition that uses a geometric formulation to specify the contact angle, and a colour-conserving boundary closure scheme to improve the numerical accuracy and suppress spurious velocities at the contact line. An additional convection–diffusion equation is solved by the passive scalar approach to obtain the temperature field, which is coupled to the hydrodynamic equations through an equation of state. This model is first validated by simulations of static contact angle and dynamic capillary intrusion process when a constant interfacial tension is considered. It is then used to simulate the thermocapillary migration of a microfluidic droplet on a horizontal solid surface subject to a uniform temperature gradient. We for the first time demonstrate numerically that the droplet motion undergoes two different states depending on the surface wettability: the droplet migrates towards the cooler regions on hydrophilic surfaces but reverses on hydrophobic surfaces. Decreasing the viscosity ratio can enhance the intensity of thermocapillary vortices, leading to an increase in migration velocity. The contact angle hysteresis, i.e., the difference between the advancing and receding contact angles, is always positive regardless of the contact angle and viscosity ratio. The contact angle hysteresis and the migration velocity both first decrease and then increase with the contact angle, and their minimum values occur at the contact angle of 90 degrees.",
keywords = "thermocapillary migration, lattice Boltzmann method, multiphase flow, wetting boundary condition, microfluidics, droplet dynamics",
author = "Haihu Liu and Yonghao Zhang",
year = "2015",
month = "1",
day = "1",
doi = "10.1016/j.jcp.2014.09.015",
language = "English",
volume = "280",
pages = "37--53",
journal = "Journal of Computational Physics",
issn = "0021-9991",

}

Modelling thermocapillary migration of a microfluidic droplet on a solid surface. / Liu, Haihu; Zhang, Yonghao.

In: Journal of Computational Physics, Vol. 280, 01.01.2015, p. 37-53.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Modelling thermocapillary migration of a microfluidic droplet on a solid surface

AU - Liu, Haihu

AU - Zhang, Yonghao

PY - 2015/1/1

Y1 - 2015/1/1

N2 - A multiphase lattice Boltzmann model is developed to simulate immiscible thermocapillary flows with the presence of fluid–surface interactions. In this model, interfacial tension force and Marangoni stress are included by introducing a body force term based on the concept of continuum surface force, and phase segregation is achieved using the recolouring algorithm proposed by Latva-Kokko and Rothman. At a solid surface, fluid–surface interactions are modelled by a partial wetting boundary condition that uses a geometric formulation to specify the contact angle, and a colour-conserving boundary closure scheme to improve the numerical accuracy and suppress spurious velocities at the contact line. An additional convection–diffusion equation is solved by the passive scalar approach to obtain the temperature field, which is coupled to the hydrodynamic equations through an equation of state. This model is first validated by simulations of static contact angle and dynamic capillary intrusion process when a constant interfacial tension is considered. It is then used to simulate the thermocapillary migration of a microfluidic droplet on a horizontal solid surface subject to a uniform temperature gradient. We for the first time demonstrate numerically that the droplet motion undergoes two different states depending on the surface wettability: the droplet migrates towards the cooler regions on hydrophilic surfaces but reverses on hydrophobic surfaces. Decreasing the viscosity ratio can enhance the intensity of thermocapillary vortices, leading to an increase in migration velocity. The contact angle hysteresis, i.e., the difference between the advancing and receding contact angles, is always positive regardless of the contact angle and viscosity ratio. The contact angle hysteresis and the migration velocity both first decrease and then increase with the contact angle, and their minimum values occur at the contact angle of 90 degrees.

AB - A multiphase lattice Boltzmann model is developed to simulate immiscible thermocapillary flows with the presence of fluid–surface interactions. In this model, interfacial tension force and Marangoni stress are included by introducing a body force term based on the concept of continuum surface force, and phase segregation is achieved using the recolouring algorithm proposed by Latva-Kokko and Rothman. At a solid surface, fluid–surface interactions are modelled by a partial wetting boundary condition that uses a geometric formulation to specify the contact angle, and a colour-conserving boundary closure scheme to improve the numerical accuracy and suppress spurious velocities at the contact line. An additional convection–diffusion equation is solved by the passive scalar approach to obtain the temperature field, which is coupled to the hydrodynamic equations through an equation of state. This model is first validated by simulations of static contact angle and dynamic capillary intrusion process when a constant interfacial tension is considered. It is then used to simulate the thermocapillary migration of a microfluidic droplet on a horizontal solid surface subject to a uniform temperature gradient. We for the first time demonstrate numerically that the droplet motion undergoes two different states depending on the surface wettability: the droplet migrates towards the cooler regions on hydrophilic surfaces but reverses on hydrophobic surfaces. Decreasing the viscosity ratio can enhance the intensity of thermocapillary vortices, leading to an increase in migration velocity. The contact angle hysteresis, i.e., the difference between the advancing and receding contact angles, is always positive regardless of the contact angle and viscosity ratio. The contact angle hysteresis and the migration velocity both first decrease and then increase with the contact angle, and their minimum values occur at the contact angle of 90 degrees.

KW - thermocapillary migration

KW - lattice Boltzmann method

KW - multiphase flow

KW - wetting boundary condition

KW - microfluidics

KW - droplet dynamics

UR - http://www.sciencedirect.com/science/article/pii/S0021999114006512

U2 - 10.1016/j.jcp.2014.09.015

DO - 10.1016/j.jcp.2014.09.015

M3 - Article

VL - 280

SP - 37

EP - 53

JO - Journal of Computational Physics

T2 - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -