A phase-field moving contact line model is presented for a two-phase system with soluble surfactants. With the introduction of some scalar auxiliary variables, the original free energy functional is transformed into an equivalent form, and then a new governing system is obtained. The resulting model consists of two Cahn-Hilliard-type equations and incompressible Navier-Stokes equation with variable densities, together with the generalized Navier boundary condition for the moving contact line. We prove that the proposed model satisfies the total energy dissipation with time. To numerically solve such a complex system, we develop a nonlinearly coupled scheme with unconditional energy stability. A splitting method based on pressure stabilization is used to solve the Navier-Stokes equation. Some subtle implicit-explicit treatments are adopted to discretize convection and stress terms. A stabilization term is artificially added to balance the explicit nonlinear term associated with the surface energy at the fluid-solid interface. We rigorously prove that the proposed scheme can preserve the discrete energy dissipation. An efficient finite difference method on staggered grids is used for the spatial discretization. Numerical results in both two and three dimensions demonstrate the accuracy and energy stability of the proposed scheme. Using our model and numerical scheme, we investigate the wetting behavior of droplets on a solid wall. Numerical results indicate that surfactants can affect the wetting properties of droplet by altering the value of contact angles.
|Original language||English (US)|
|Journal||Journal of Computational Physics|
|State||Published - Dec 4 2019|
Bibliographical noteKAUST Repository Item: Exported on 2020-10-01
Acknowledged KAUST grant number(s): BAS/1/1351-01-01, REP/1/2879-01, URF/1/2993-01
Acknowledgements: Jun Yao and Guangpu Zhu acknowledge that this work is supported by the National Science and Technology Major Project (2016ZX05011-001), the Natural Science Foundation of China (51804325, 5167428041 and 51774317) and Shandong Provincial Natural Science Foundation (ZR2019JQ21). The work of Shuyu Sun and Jisheng Kou is supported by the KAUST research fund awarded to the Computational Transport Phenomena Laboratory at KAUST through the Grant BAS/1/1351-01-01, URF/1/2993-01, and REP/1/2879-01.