Slip length in sheared liquid films subject to mixed
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.
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
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
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