Moving Contact Line Hydrodynamics

Updated 12 June 2026

Moving contact line hydrodynamics is the study of fluid interfaces at solid boundaries, coupling molecular kinetics with macroscopic flows.
It employs regularization techniques like Navier-slip and phase-field models to resolve stress singularities observed in the classical no-slip assumption.
Advanced experiments and simulations validate numerical models that underpin practical applications in wetting, coating, boiling, and nucleation.

A moving contact line is the locus where two immiscible fluids meet a solid boundary and this intersection is dynamically displaced due to external or internal driving forces. The hydrodynamics of the moving contact line are inherently multiscale, coupling molecular kinetic processes at the nanometer scale with macroscopic interface deformation and bulk flow. The classical “contact-line paradox”—the prediction of a non-integrable viscous stress singularity in the Stokes solution for fluids meeting at a finite angle on a no-slip surface—has driven the development of diverse theoretical, computational, and experimental approaches to resolve this singularity, clarify dissipation mechanisms, and establish well-posed models of contact-line motion relevant for applications ranging from wetting, coating, and spreading to boiling and nucleation phenomena.

1. Fundamental Theoretical Frameworks

The starting point for contact-line hydrodynamics is the set of incompressible Navier–Stokes or Stokes equations for each fluid phase, supplemented by appropriate boundary and interfacial conditions. In the presence of a moving solid wall, the velocity field near the contact line, captured in the viscous-wedge or Huh–Scriven solution, exhibits a logarithmic divergence in stress as the distance to the contact line $r \to 0$ —a consequence of enforcing both no-slip at the solid wall and a finite contact-angle kinematic constraint at the interface (Gupta et al., 2023). The classical sharp-interface, zero-slip hydrodynamic models demonstrate this singularity, leading to non-integrable dissipation and ill-posedness (Lukyanov et al., 2017).

Resolution requires introducing a microscopic regularization that relaxes incompatible boundary conditions. The most established approaches include:

Navier-slip boundary condition: allows finite slip velocity with slip length $b$ at the wall, replacing no-slip by $u_t = b \partial_n u_t$ (Giri et al., 2022, Zhao et al., 2020).
Phase-field (diffuse-interface) models: replace the sharp liquid–gas interface with a continuous transition of width $\epsilon$ , eliminating the singular zone (Sibley et al., 2013, Li et al., 11 Mar 2025, Yu et al., 2017).
Interface formation models: introduce mass exchange between bulk and interfacial “surface phase”, relaxing the impermeability and tangential velocity constraints, which regularizes the singularity at the contact line (Lukyanov et al., 2017).
Contact angle–velocity laws: empirical or derived relations between dynamic contact angle $\theta_d$ and contact-line velocity $U_{CL}$ , such as the Cox–Voinov law $\theta_d^3 = \theta_0^3 + 9Ca\ln(L/\ell)$ with a microscopic cutoff $\ell$ (Giri et al., 2022, Yu et al., 2017).

For volatile fluids, nanoscale curvature-induced evaporation/condensation—the Kelvin effect—serves as a regularizing mechanism, making the cutoff length and apparent angle explicit functions of substrate temperature and volatility (Janecek et al., 2012).

2. Multiscale Nature and Singularities

The moving contact line is intrinsically multiscale, with key physical processes spanning orders of magnitude:

Molecular scale ( $\sim$ nm): Discrete fluid structure, molecular slip, and surface forces act within a few nanometers of the contact line. For atomically smooth surfaces, slip lengths extracted from molecular dynamics (MD) simulations are of order a fraction to a few nanometers and primarily determine the regularization scale for continuum models (Giri et al., 2022, Smith et al., 2018).
Intermediate/hydrodynamic scale (100 nm–100 µm): The continuum solution (e.g., Huh–Scriven wedge, Cox–Voinov framework) is valid except within a cutoff region around the contact line; beyond this, macroscopic dissipation dominates, and interface shapes reflect viscous bending captured by matched-asymptotic expansions (Gupta et al., 2023, Chan et al., 2013).
Macroscopic scale (>100 µm): Gravity, static meniscus, and bulk flow govern the far-field interface shape and flow.

Recent high-resolution experiments and simulations confirm that the apparent viscous bending of the interface predicted by theory is absent in the raw interface profiles when the outer meniscus is included in the comparison. Instead, the interface rapidly decelerates near the contact point, effectively resolving the stress singularity not by slip but by a natural reduction in interfacial flow (Gupta et al., 2023). In phase-field and diffuse-interface models, the stress and velocity remain finite and single-valued as the contact line is approached, provided the interface width $\epsilon$ is respected (Sibley et al., 2013, Yu et al., 2017, Li et al., 11 Mar 2025).

3. Modern Multiscale and Molecular Approaches

Molecular dynamics (MD) simulations provide direct access to the structure and dynamics at the moving contact line and allow one to extract slip lengths, interfacial friction coefficients, and verify continuum predictions at scales down to $b$ 0, where $b$ 1 is the molecular length (Giri et al., 2022, Smith et al., 2018, Lācis et al., 2020). MD–continuum coupling (domain overlap or parameterized closure) enables seamless multiscale computation and supports findings such as:

For atomically smooth substrates, dissipation is dominated by bulk viscous flow, with wall slip and contact-line localized friction accounting for less than $b$ 2 of total dissipation (Giri et al., 2022).
The dynamic contact angle–velocity relation follows the Cox–Voinov law, with the slip length $b$ 3 setting the microscopic cutoff $b$ 4 (Giri et al., 2022, Smith et al., 2018).
Hydrogen bonding or surface roughness can reduce slip length to sub-nm, as observed for water on silica (Lācis et al., 2020).

Hydrodynamic density-functional theory (HDDFT) offers a continuum approach incorporating molecular layering and nonlocal thermodynamics. It predicts spatially distinct regions near the contact line: an oscillatory (structured) slip region near the wall (width $b$ 5– $b$ 6) and an extended compressive dissipation region on the vapor side of the interface whose length grows rapidly with temperature (Nold et al., 2024). This model unites nanoscale layering, mesoscale dissipation, and macroscopic wetting dynamics with no ad hoc parameters.

4. Volatility and Evaporation Effects

For volatile fluids, interfacial phase change supplies a distinct regularization mechanism. The Kelvin effect—the curvature dependence of equilibrium vapor pressure—regulates the temperature at the interface. The evaporative/condensation flux $b$ 7 is set by conduction through the thin liquid film and balances the viscous flux divergence. This leads to a matched-asymptotic structure:

The microscopic cutoff (the Voinov length $b$ 8) is set by a Kelvin length $b$ 9, depending only on fluid properties and substrate temperature, not on a molecular slip length.
The macroscopic meniscus and dynamic angle (the Voinov angle $u_t = b \partial_n u_t$ 0) are explicit functions of wall superheat $u_t = b \partial_n u_t$ 1; $u_t = b \partial_n u_t$ 2 exceeds the Young angle $u_t = b \partial_n u_t$ 3 for $u_t = b \partial_n u_t$ 4 (Janecek et al., 2012).

This Kelvin-regularized model is highly predictive for boiling and evaporation: the apparent contact angle and onset of liquid entrainment both depend on $u_t = b \partial_n u_t$ 5, with the Kelvin length and functional dependence precisely derived from the inner problem.

5. Numerical and Computational Modeling

Multiple high-fidelity numerical frameworks have been established:

Sharp-interface, fitted-mesh finite element approaches implement Navier–Stokes/Stokes equations with Navier-slip and a dynamic contact-angle law, often through a Robin-type boundary condition (imposed weakly) (Zhao et al., 2020). These methods demonstrate unconditional energy stability and robust convergence for large viscosity/density ratios (Zhao et al., 2020).
Diffuse-interface and phase-field models (Cahn–Hilliard or Allen–Cahn) fully regularize the contact-line singularity by encoding surface tension and contact-angle energetics within spatially resolved order parameters, allowing for dynamic or relaxation boundary conditions at the wall (Yu et al., 2017, Li et al., 11 Mar 2025). Analytical results establish global existence and decay under certain geometries and parameter regimes (Li et al., 11 Mar 2025).
Level-set and volume-of-fluid (VOF) methods with local Navier slip or generalized Navier boundary conditions translate the continuum picture to the interface-capturing context, with careful treatment of singular force terms to impose dynamic contact-angle laws (Zhao et al., 2021). Resolution of the contact line is linked to grid scale or localized slip "bells" in VOF (Lācis et al., 2020).
Particle methods (SPH) enforce static or dynamic contact angles via correction of local surface normals and color-gradient interpolations, reducing unphysical shear at the triple point and reproducing classical $u_t = b \partial_n u_t$ 6 behavior even in meshfree, Lagrangian settings (Farrokhpanah et al., 2016, Bao et al., 2018).

A key challenge remains the accurate representation of dynamic contact angle hysteresis and the rate-dependent cross-coupling between local interfacial viscous stresses, surface roughness, and the imposed slip length or interfacial friction (Bao et al., 2018).

6. Contact-Line Dynamics on Heterogeneous and Curved Substrates

Modeling contact-line movement over chemically or topographically inhomogeneous surfaces introduces spatial variation in the microscopic contact angle $u_t = b \partial_n u_t$ 7, leading to a dynamically disordered energy landscape for the contact line. Linear-response and modal analyses reduce the multi-dimensional lubrication problem to nonlocal 1+1D “rheological” laws featuring memory kernels and stochastic noise, facilitating the study of depinning, thermal activation, and line elasticity (Perrin et al., 2017).

For curved walls, perturbations of the classical Huh–Scriven wedge flow yield regular solutions where the leading $u_t = b \partial_n u_t$ 8 singularity is cancelled by geometric correction, enabling non-singular Dirichlet velocity boundary conditions for moving contact-point simulations on general substrates (Holmgren et al., 2019, Holmgren et al., 2017). This approach extends to multiscale computational coupling, where phase-field “microboxes” provide dynamic angle–velocity relations driving effective boundary conditions at macroscopic scale, resolving the singularity (Holmgren et al., 2017).

7. Experimental Validation and Physical Implications

Recent synchronized interface-shape and PIV measurements in plate-immersion experiments have validated the theoretical constructs across four decades in viscosity and capillary numbers. The main findings (Gupta et al., 2023) include:

Absence of predicted viscous bending in the raw interface profile; composite models including full meniscus geometry are essential for correct shape prediction.
Quantitative agreement between modulated-wedge solution and experimentally measured (u,v) fields everywhere but within ~100 μm of the contact line.
Deceleration of material points along the interface as they approach the contact line, directly removing the classical stress singularity.
The measured deceleration provides critical boundary condition data for numerical models, facilitating the development of new interface-formation or inner-region regularization theories.

This experimental validation substantiates the matched-asymptotic and multiscale paradigms and sets new benchmarks for theoretical and computational development in contact-line hydrodynamics.

References:

(Janecek et al., 2012) Janeček & Nikolayev, "Moving contact line of a volatile fluid", Phys. Rev. Lett. 111, 154501 (2013)
(Gupta et al., 2023) Gupta et al., "An experimental study of flow near an advancing contact line: a rigorous test of theoretical models"
(Giri et al., 2022) Bonn et al., "Resolving the microscopic hydrodynamics at the moving contact line"
(Li et al., 11 Mar 2025) Ding et al., "Navier-Stokes/Allen-Cahn system with moving contact line"
(Nold et al., 2024) Ledesma et al., "Hydrodynamic density-functional theory for the moving contact-line problem reveals fluid structure and emergence of a spatially distinct pattern"
(Lukyanov et al., 2017) Lukyanov & Pryer, "Hydrodynamics of moving contact lines: macroscopic versus microscopic"
(Sibley et al., 2013) Sibley et al., "The contact line behaviour of solid-liquid-gas diffuse-interface models"
(Yu et al., 2017) Yang et al., "Numerical Approximations for a Phase-Field Moving Contact Line Model with Variable Densities and Viscosities"
(Farrokhpanah et al., 2016) Farrokhpanah et al., "Applying Contact Angle to a 2D Multiphase Smoothed Particle Hydrodynamics Model"
(Chan et al., 2013) Chan et al., "Hydrodynamics of air entrainment by moving contact lines"
(Bao et al., 2018) Bao et al., "A Modified Smoothed Particle Hydrodynamics Approach for Modelling Dynamic Contact Angle Hysteresis"
(Perrin et al., 2017) Perrin et al., "A one dimensional modal approach for flows controlled by contact line motion"
(Smith et al., 2018) De Ruijter et al., "Moving Contact Lines: Linking Molecular Dynamics and Continuum-Scale Modeling"
(Holmgren et al., 2019) Holmgren & Kreiss, "A Hydrodynamic Model of Movement of a Contact Line Over a Curved Wall"
(Holmgren et al., 2017) Holmgren & Kreiss, "A Computational Multiscale Model for Contact Line Dynamics"
(Zhao et al., 2021) Zhao et al., "A Level Set Method for the Simulation of Moving Contact Lines in Three Dimensions"
(Zhao et al., 2020) Zhao & Ren, "An Energy-stable Finite Element Method for the Simulation of Moving Contact Lines in Two-phase Flows"