Power-Law Slip Profile of the Moving Contact Line in Two-Phase Immiscible Flows

Large scale molecular dynamics (MD) simulations on two-phase immiscible flows show that associated with the moving contact line, there is a very large $1/x$ partial-slip region where $x$ denotes the distance from the contact line. This power-law partial-slip region is verified in large-scale adaptive continuum simulations based on a local, continuum hydrodynamic formulation, which has proved successful in reproducing MD results at the nanoscale. Both MD and continuum simulations indicate the existence of a universal slip profile in the Stokes-flow regime, well described by $v^{slip}(x)/V_w=1/(1+{x}/{al_s})$, where $v^{slip}$ is the slip velocity, $V_w$ the speed of moving wall, $l_s$ the slip length, and $a$ is a numerical constant. Implications for the contact-line dissipation are discussed.Comment: 13 pages, 3 figure

An Energetic Variational Approach for ion transport

The transport and distribution of charged particles are crucial in the study of many physical and biological problems. In this paper, we employ an Energy Variational Approach to derive the coupled Poisson-Nernst-Planck-Navier-Stokes system. All physics is included in the choices of corresponding energy law and kinematic transport of particles. The variational derivations give the coupled force balance equations in a unique and deterministic fashion. We also discuss the situations with different types of boundary conditions. Finally, we show that the Onsager's relation holds for the electrokinetics, near the initial time of a step function applied field

Molecular hydrodynamics of the moving contact line in two-phase immiscible flows

The ``no-slip'' boundary condition, i.e., zero fluid velocity relative to the solid at the fluid-solid interface, has been very successful in describing many macroscopic flows. A problem of principle arises when the no-slip boundary condition is used to model the hydrodynamics of immiscible-fluid displacement in the vicinity of the moving contact line, where the interface separating two immiscible fluids intersects the solid wall. Decades ago it was already known that the moving contact line is incompatible with the no-slip boundary condition, since the latter would imply infinite dissipation due to a non-integrable singularity in the stress near the contact line. In this paper we first present an introductory review of the problem. We then present a detailed review of our recent results on the contact-line motion in immiscible two-phase flow, from MD simulations to continuum hydrodynamics calculations. Through extensive MD studies and detailed analysis, we have uncovered the slip boundary condition governing the moving contact line, denoted the generalized Navier boundary condition. We have used this discovery to formulate a continuum hydrodynamic model whose predictions are in remarkable quantitative agreement with the MD simulation results at the molecular level. These results serve to affirm the validity of the generalized Navier boundary condition, as well as to open up the possibility of continuum hydrodynamic calculations of immiscible flows that are physically meaningful at the molecular level.Comment: 36 pages with 33 figure

A variational approach to the moving contact line hydrodynamics

In immiscible two-phase flows, contact line denotes the intersection of the fluid-fluid interface with the solid wall. When one fluid displaces the other, the contact line moves along the wall. A classical problem in continuum hydrodynamics is the incompatibility between the moving contact line and the no-slip boundary condition, as the latter leads to a non-integrable singularity. The recently discovered generalized Navier boundary condition (GNBC) offers an alternative to the no-slip boundary condition which can resolve the moving contact line conundrum. We present a variational derivation of the GNBC through the principle of minimum energy dissipation (entropy production), as formulated by Onsager for small perturbations away from the equilibrium. Through numerical implementation of a continuum hydrodynamic model, it is demonstrated that the GNBC can quantitatively reproduce the moving contact line slip velocity profiles obtained from molecular dynamics simulations. In particular, the transition from complete slip at the moving contact line to near-zero slip far away is shown to be governed by a power-law partial slip regime, extending to mesoscopic length scales. The sharp (fluid-fluid) interface limit of the hydrodynamic model, together with some general implications of slip versus no-slip, are discussed.Comment: 44 pages, 8 figure