Modeling the dynamics of partial wetting

Abstract

The behavior of interfaces and contact lines arises from intermolecular interactions like Van der Waals forces. To consider this multi–phase behavior on the continuum scale, appropriate physical descriptions must be formulated. While the Continuum Surface Force model is well–engineered for the description of interfaces, there is still a lack of treatment of contact lines, which are represented by the intersection of a fluid–fluid interface and a solid boundary surface. In our approach we use the “non compensated Young force” to model contact line dynamics and therefore use an extension to the Navier–Stokes equations in analogy to the extension of a two–phase interface in the CSF model. Because particle–based descriptions are well–suited for changing and moving interfaces we use Smoothed Particle Hydrodynamics. In this way we are not only able to calculate the equilibrium state of a two–phase interface with a static contact angle, but also for instance able to simulate droplet shapes and their dynamical evolution with corresponding contact angles towards the equilibrium state, as well as different pore wetting behavior. Together with the capability to model density differences, this approach has a high potential to model recent challenges of two–phase transport in porous media. Especially with respect to moving contact lines this is a novelty and indispensable for problems, where the dynamic contact angle dominates the system behavior.