Finite Element Framework for Describing Dynamic Wetting Phenomena

TL;DR

This paper develops a robust finite element framework for simulating dynamic wetting phenomena, providing guidelines for mesh resolution, implementation consistency, and potential extensions to include additional physical effects.

Contribution

The paper introduces a new FEM framework that simplifies boundary condition implementation at the contact line and enables comprehensive numerical analysis of wetting flows.

Findings

Framework removes complex rotations of momentum equations

Guidelines for mesh resolution based on non-dimensional parameters

Framework can be extended to include surface tension gradients

Abstract

The finite element simulation of dynamic wetting phenomena, requiring the computation of flow in a domain confined by intersecting a liquid-fluid free surface and a liquid-solid interface, with the three-phase contact line moving across the solid, is considered. For this class of flows, different finite element method (FEM) implementations have been proposed in the literature and these are seen to produce apparently contradictory results. The purpose of this paper is to develop a robust framework for the FEM simulation of these flows and then, by performing numerical experiments, provide guidelines for future investigations. In the new framework, the boundary conditions on the liquid-solid interface are implemented in a methodologically similar way to those on the free surface so that the equations at the contact line, where the interfaces meet, are applied without any ad-hoc alterations. The new implementation removes the need for complex rotations of the momentum equations, usually required to apply boundary conditions normal and tangent to the solid surface. The developed code allows the convergence of the solution to be studied as the spatial resolution of the computational mesh is varied over many orders of magnitude. This makes it possible to provide practical recommendations on the spatial resolution required by a numerical scheme for a given set of non-dimensional similarity parameters. Furthermore, one can examine various implementations used in the literature and evaluate their performance. Finally, it is shown how the framework may be generalized to account for additional physical effects, such as gradients in surface tensions. A user-friendly step-by-step guide specifying the entire implementation and allowing the reader to easily reproduce all presented results is provided in the Appendix.