A Two-Phase SPH Model for Dynamic Contact Angles Including Fluid–Solid Interactions at the Contact Line

Abstract

The description of wetting phenomena on the continuum scale is a challenging problem, since intermolecular interactions, like van der Waals forces between liquid and solid, alter the flow field at the contact line. Recently, these effects were included in the smoothed particle hydrodynamics method by introducing a contact line force (CLF) on the continuum scale. This physically based contact line force model is employed here to simulate two-phase flow in a wide range of wetting dynamics parametrized by capillary number. In particular, dynamic contact angles at various capillary number values are calculated by CLF model and compared to measured values. We find that there is significant disagreement between simulated and measured results, specially at low wetting speeds. It is indeed expected that most of the driving force is dissipated to overcome strong liquid-solid interactions, which are not adequately accounted for in the existing CLF model. Therefore, we have extended that model to account for stick-slip (SSL) behavior of the contact line caused by solid-fluid interactions. The new SSL model results in dynamic contact angle values that are in good agreement with experimental data for the full range of wetting dynamics. Keywords Dynamic contact angle • Contact line force • Two-phase flow • Stick-slip behavior • Smoothed particle hydrodynamics (SPH) List of symbols A Arbitrary field variable B Transition function for continuous surface force at solid boundary c Color function B P. Kunz

Simulation of Electrolyte Imbibition in Gas Diffusion Electrodes

“…Details on the algorithm and boundary conditions are presented in , . The accuracy of the applied two‐phase flow algorithm was shown in a previous work .…”

Pore-network modeling of flow in shale nanopores: Network structure, flow principles, and computational algorithms

“…DNS methods such as computational fluid dynamics (CFD) and Lattice Boltzmann methods (LBM) solve pore-scale governing equations in actual pore space obtained from digitalized porous media without modification or idealization (Chen et al, 2015;Guo et al, 2018). In the DNS methods for two phases, modeling the interface by sharp interface models or diffuse interface models, has been an intractable problem (Kunz et al, 2018;Yin et al, 2019;Zhu et al, 2020;Zhu et al, 2019). The moving contact line, defined as the intersection of the fluid-fluid interface with the solid wall, is incompatible with the no-slip boundary condition on the solid wall surface.…”