Phase-Field Models, Sharp Interface Limits, and Numerical Schemes for Contact Line Dynamics

We study phase-field and sharp-interface models for contact line dynamics of a liquid droplet on a solid substrate within a unified variational framework.

The motion of the contact line, where liquid, gas, and solid phases meet, poses a fundamental difficulty in continuum modeling due to the classical stress singularity of no-slip hydrodynamics.

Phase-field models regularize this singularity by introducing a thin transition layer of thickness and encoding interfacial effects through a Ginzburg-Landau free energy augmented by a wall energy on the substrate.

Starting from the total free energy $E = E_b + E_w$, we analyze two phase-field models: the Allen-Cahn equation and the Cahn-Hilliard equation.

Using matched asymptotic expansions as $\delta \to 0$, we recover their corresponding sharp interface limits.

In the Allen-Cahn case, the limit yields motion by mean curvature with a contact line law driven by deviations of the dynamic contact angle from Young's angle.

In the Cahn-Hilliard case, the limit leads to a Mullins-Sekerka problem with the same form of contact line dynamics.

A central result of this work is the identification of consistent gradient-flow structures across both models.

The Allen-Cahn dynamics correspond to an $L^2$-gradient flow, while the Cahn-Hilliard dynamics correspond to an $H^{-1}$-gradient flow, and both converge to sharp-interface evolutions that preserve the same energy-dissipation structure.

This provides a unified interpretation of contact line motion as a consequence of a single variational principle.

Finally, we develop energy-stable numerical schemes based on the minimizing movement principle and establish discrete energy dissipation and well-posedness of the fully discrete problem.

Numerical examples confirm that both schemes relax toward the same stationary sharp interface solution while their dynamics reflect the different dissipation mechanisms.