A Thermodynamically Consistent Model and Its Conservative Numerical Approximation for Moving Contact Lines with Soluble Surfactants
arXiv:2104.02626 · doi:10.1016/j.cma.2021.114033
Abstract
We derive a continuum sharp-interface model for moving contact lines with soluble surfactants in a thermodynamically consistent framework. The model consists of the isothermal two-phase incompressible Navier-Stokes equations for the fluid dynamic and the bulk\slash surface convection-diffusion equations for the surfactant transportation. The interface condition, the slip boundary condition, the dynamic contact angle condition, and the adsorption\slash desorption condition are derived based on the principle of the total free energy dissipation. In particular, we recover classical adsorption isotherms from different forms of the surface free energy. The model is then numerically solved in two spatial dimensions. We present an Eulerian weak formulation for the Navier-Stokes equations together with an arbitrary Lagrangian-Eulerian weak formulation for the surfactant transport equations. Finite element approximations are proposed to discretize the two weak formulations on the moving mesh. The resulting numerical method is shown to conserve the total mass of the surfactants exactly. By using the proposed model and its numerical method, we investigate the droplet spreading and migration in the presence of surfactants and study their dependencies on various dimensionless adsorption parameters.
References in corpus (9)
- Molecular scale contact line hydrodynamics of immiscible flows
- A variational approach to the moving contact line hydrodynamics
- Hydrodynamic theory of de-wetting
- Gradient dynamics models for liquid films with soluble surfactant
- Field Theory of Reaction-Diffusion: Mass Action with an Energetic Variational Approach
- On the distinguished limits of the Navier slip model of the moving contact line problem
- On the Stable Numerical Approximation of Two-Phase Flow with Insoluble Surfactant
- A Finite Element Method for Electrowetting on Dielectric
- An Energy-stable Finite Element Method for the Simulation of Moving Contact Lines in Two-phase Flows