Slip length in sheared liquid films subject to mixed boundary conditions

Nikolai V. Priezjev, Anton A. Darhuber and Sandra M. Troian, Phys. Rev. E 71, 041608 (2005)


Fluid flow in confined geometries can be significantly affected by slip at the liquid/solid interface. The measure of slip is the so-called slip length, which is defined as an extrapolated distance relative to the wall where the tangential velocity component vanishes, see picture below. Factors that affect slip length include weak wall fluid attraction, surface roughness, high shear rates and nucleation of nanobubbles at hydrophobic surfaces. Molecular dynamics simulations provide a powerful tool for studying the effect of these parameters on boundary conditions at the molecular level.




As discussed in our paper, we investigated the behavior of the slip length in simple fluids subject to planar shear bounded by substrates with mixed boundary conditions. This study was inspired by recent experiments which suggest a possible presence of flat nanobubbles or dissolved gas at the liquid/solid interface. A significant drag reduction could be achieved because of no-shear stress (infinite slip) regions at the liquid/gas interface. We mimic this situation by considering a shear flow past a surface which is patterned with an infinite array of stripes representing alternating regions of no-shear and finite slip, see picture below. The question we try to answer is how the slip length depends on the distribution of non-wetting regions?




In this work the fluid was subject to steady planar shear by translating the upper homogeneous wall at a constant velocity U, while the lower, patterned wall remained stationary, see picture below. We reduced the attractive part of Lennard-Jones potential for fluid-wall interactions to simulate non-wetting (white) stripes.




In this movie we consider shear flow past a surface with a stripe width of about 6s, where s is the molecular diameter characterizing the Lennard-Jones interaction. Note the enhanced layering near the wetting patches and reduced fluid density above the non-wetting regions. When the stripe width is of order of molecular diameter, the alternating wall potential represents effectively a rough surface. These molecular scale corrugations created by the composite wall potential strongly reduces the effective slip length in the transverse flow orientation, leading eventually to a no-slip boundary condition at approximately a~s.




The difference in the contact densities and relative velocities above the wetting and non-wetting patches produces significant oscillations in the averaged velocity profiles for large stripe widths. We found that effective slip length increases monotonically with stripe period to a saturation value for both parallel and transverse flow orientations. Detailed comparison between the results of the hydrodynamic calculations and molecular dynamics simulations shows excellent agreement when the length scale of the substrate pattern geometry is roughly larger than ten molecular diameters.

In a recent study, we have considered a more general situation when the mean flow direction is not aligned with the symmetry axis of a patterned substrate and the apparent slip velocity acquires a non-zero transverse component. It was found that the angular dependence of the effective slip length obtained from MD simulations is in good agreement with hydrodynamic predictions provided that the stripe width is larger than several molecular diameters. Perhaps most interestingly, we find that the directional diffusion coefficient of fluid molecules in contact with patterned substrate correlates well with the effective slip length as a function of the shear flow direction. These findings lend support for the microscopic justification of the tensor formulation of the effective slip boundary conditions for noninertial flows of Newtonian fluids over smooth surfaces with nanoscale anisotropic textures. (DFD talk).


Molecular dynamics simulations described here were conducted with LAMMPS numerical code.

Back to Nikolai Priezjev's page