Molecular diffusion and slip boundary conditions at smooth surfaces with periodic and random nanoscale textures

Insights from molecular simulations on liquid slip over nanostructured surfaces

The current study focuses on non-equilibrium molecular dynamics (NEMD) simulations to investigate the slip properties of water flowing over different nanostructured surfaces. A simulation protocol is developed that applies constant shear stress throughout the fluid before measuring the slip length. Using pseudo-data, the reliability of this protocol in terms of both accuracy and noise of the results for high-slip and multiphase systems is demonstrated. In contrast to the NEMD techniques available in the literature, the protocol also enables a convenient way to compare the slip lengths of different surface coatings. The fluid slip lengths of surface coatings comprising carbon nanotubes on platinum are predicted using the proposed protocol with nitrogen gas trapped in the interstitial gaps. The role of these gas pockets in determining surface slip properties is investigated. The NEMD results from the proposed model compare well with a macroscopic theoretical model for nano-patterned surfaces. Finally, it is concluded that entrapped gas within nanostructures may offer significant drag reduction only if the gas surface coverage is above 95%.

A review on slip boundary conditions at the nanoscale: recent development and applications

The slip boundary condition for nanoflows is a key component of nanohydrodynamics theory, and can play a significant role in the design and fabrication of nanofluidic devices. In this review, focused on the slip boundary conditions for nanoconfined liquid flows, we firstly summarize some basic concepts about slip length including its definition and categories. Then, the effects of different interfacial properties on slip length are analyzed. On strong hydrophilic surfaces, a negative slip length exists and varies with the external driving force. In addition, depending on whether there is a true slip length, the amplitude of surface roughness has different influences on the effective slip length. The composition of surface textures, including isotropic and anisotropic textures, can also affect the effective slip length. Finally, potential applications of nanofluidics with a tunable slip length are discussed and future directions related to slip boundary conditions for nanoscale flow systems are addressed.

Influence of Surface Texture on the Variation of Electrokinetic Streaming Potentials

The electrokinetic streaming potential (V_s) obtained through electrolyte flow in a microchannel is shown to be related to the underlying surface pattern. Pillar, mesh, and groove patterns were studied for comparing the relative magnitudes of the V_s with air-/liquid-filled surfaces. A record value of the related figure of merit, in terms of the developed V_s per-unit applied pressure, of ∼0.127 mV/Pa, was observed in a mesh texture liquid-filled surface (LFS) impregnated with an electrolyte-immiscible oil. The study indicated that increasing the solid fraction of the pattern surface decreases the effective slip length while enhancing the overall channel ζ potential. Consequently, maximizing the obtained V_s implies a balancing of the slip with the surface potential, with plausibly more significance of the latter. The work has implications for higher-efficiency electrical voltage sources.

Molecular Dynamics Simulations of the Rotational and Translational Diffusion of a Janus Rod-Shaped Nanoparticle

The diffusion of a Janus rod-shaped nanoparticle in a dense Lennard-Jones fluid is studied using molecular dynamics (MD) simulations. The Janus particle is modeled as a rigid cylinder whose atoms at each end have different interaction energies with fluid molecules, thus comprising wetting and nonwetting surfaces. We found that both rotational and translational diffusion coefficients are larger for Janus particles with lower average wettability, and these values are bound between the two limiting cases of uniformly wetting and nonwetting particles. It was also shown that values of the diffusion coefficients for displacements parallel and perpendicular to the major axis of a uniformly wetting particle agree well with analytical predictions despite a finite slip at the particle surface present in MD simulations. It was further demonstrated that diffusion of Janus particles is markedly different from that of uniform particles; namely, Janus particles preferentially rotate and orient their nonwetting ends along the displacement vector to reduce drag. This correlation between translation and rotation is consistent with the previous results on diffusive dynamics of a spherical Janus particle with two hemispheres of different wettability.

Characterizing Structural Complexity in Disordered Carbons: From the Slit Pore to Atomistic Models

The reliable characterization of nanoporous carbons is critical to the design and optimization of their numerous applications; however, the vast majority of carbons in industrial use are highly disordered, with complex structures whose understanding has long challenged researchers. The idealized slit pore model represents the most commonly used approximation to a carbon nanopore; nevertheless, it has been only partially successful in predicting adsorption isotherms and fails significantly in predicting transport properties because of its inability to capture structural disorder and its effect on fluid accessibility. Atomistic modeling of the structure has much potential for overcoming this limitation, and among such approaches, hybrid reverse Monte Carlo simulation has emerged as the most attractive. This method reconstructs the structure of a carbon based on the fitting of its experimentally measured pair distribution function and appropriate properties such as porosity while minimizing the energy. The method is shown to be best implemented using a multistage strategy, with the first stage used to attain a deep minimum of the energy and subsequent stages to refine the structure based on the fitting of specific properties. Methods to determine the accessibility of gases based on the atomistic structure are outlined, and it is shown that energy barriers are very sensitive to small differences in the sizes of constrictions and pore entries. The ability to accurately predict macroscopic transport coefficients of adsorbates in nanoporous carbons appears to be the greatest limitation of such models. Overcoming this will require the fitting of properties more sensitive to long-range disorder than the currently used pair distribution and the use of a suitable multiscaling strategy, which is suggested as a future direction for advancing atomistic models. The inclusion of heteroatoms in the structure is also an important area requiring further attention, particularly in the development of computationally efficient force fields incorporating their interactions.

The importance of chemical potential in the determination of water slip in nanochannels

We investigate the slip properties of water confined in graphite-like nanochannels by non-equilibrium molecular dynamics simulations, with the aim of identifying and analyze separately the influence of different physical quantities on the slip length. In a system under confinement but connected to a reservoir of fluid, the chemical potential is the natural control parameter: we show that two nanochannels characterized by the same macroscopic contact angle--but a different microscopic surface potential--do not exhibit the same slip length unless the chemical potential of water in the two channels is matched. Some methodological issues related to the preparation of samples for the comparative analysis in confined geometries are also discussed.

Diffusion of a Janus nanoparticle in an explicit solvent: A molecular dynamics simulation study

Molecular dynamics simulations are carried out to study the translational and rotational diffusion of a single Janus particle immersed in a dense Lennard-Jones fluid. We consider a spherical particle with two hemispheres of different wettabilities. The analysis of the particle dynamics is based on the time-dependent orientation tensor, particle displacement, as well as the translational and angular velocity autocorrelation functions. It was found that both translational and rotational diffusion coefficients increase with decreasing surface energy at the nonwetting hemisphere, provided that the wettability of the other hemisphere remains unchanged. We also observed that in contrast to homogeneous particles, the nonwetting hemisphere of the Janus particle tends to rotate in the direction of the displacement vector during the rotational relaxation time.

Flows and mixing in channels with misaligned superhydrophobic walls

Aligned superhydrophobic surfaces with the same texture orientation reduce drag in the channel and generate secondary flows transverse to the direction of the applied pressure gradient. Here we show that a transverse shear can be easily generated by using superhydrophobic channels with misaligned textured surfaces. We propose a general theoretical approach to quantify this transverse flow by introducing the concept of an effective shear tensor. To illustrate its use, we present approximate theoretical solutions and Dissipative Particle Dynamics simulations for striped superhydrophobic channels. Our results demonstrate that the transverse shear leads to complex flow patterns, which provide a new mechanism of a passive vertical mixing at the scale of a texture period. Depending on the value of Reynolds number two different scenarios occur. At relatively low Reynolds number the flow represents a transverse shear superimposed with two corotating vortices. For larger Reynolds number these vortices become isolated, by suppressing fluid transport in the transverse direction.

Interfacial slip on rough, patterned and soft surfaces: a review of experiments and simulations

Advancements in the fabrication of microfluidic and nanofluidic devices and the study of liquids in confined geometries rely on understanding the boundary conditions for the flow of liquids at solid surfaces. Over the past ten years, a large number of research groups have turned to investigating flow boundary conditions, and the occurrence of interfacial slip has become increasingly well-accepted and understood. While the dependence of slip on surface wettability is fairly well understood, the effect of other surface modifications that affect surface roughness, structure and compliance, on interfacial slip is still under intense investigation. In this paper we review investigations published in the past ten years on boundary conditions for flow on complex surfaces, by which we mean rough and structured surfaces, surfaces decorated with chemical patterns, grafted with polymer layers, with adsorbed nanobubbles, and superhydrophobic surfaces. The review is divided in two interconnected parts, the first dedicated to physical experiments and the second to computational experiments on interfacial slip of simple (Newtonian) liquids on these complex surfaces. Our work is intended as an entry-level review for researchers moving into the field of interfacial slip, and as an indication of outstanding problems that need to be addressed for the field to reach full maturity.

Thermal slip for liquids at rough solid surfaces

Molecular dynamics simulation is used to examine the thermal slip of liquids at rough solid surfaces as characterized by fractal Cantor structures. The temperature profiles, potential energy distributions, thermal slip, and interfacial thermal resistance are investigated and evaluated for a variety of surface topographies. In addition, the effects of liquid-solid interaction, surface stiffness, and boundary condition on thermal slip length are presented. Our results indicate that the presence of roughness expands the low potential energy regions in adjacent liquids, enhances the energy transfer at liquid-solid interface, and decreases the thermal slip. Interestingly, the thermal slip length and thermal resistance for liquids in contact with solid surfaces depends not only on the statistical roughness height, but also on the fractal dimension (i.e., topographical spectrum).

Green-Kubo relation for friction at liquid-solid interfaces

We have developed a Green-Kubo relation that enables accurate calculations of friction at solid-liquid interfaces directly from equilibrium molecular dynamics (MD) simulations and that provides a pathway to bypass the time-scale limitations of typical nonequilibrium MD simulations. The theory has been validated for a number of different interfaces and it is demonstrated that the liquid-solid slip is an intrinsic property of an interface. Because of the high numerical efficiency of our method, it can be used in the design of interfaces for applications in aqueous environments, such as nano- and microfluidics.

Effective slip-length tensor for a flow over weakly slipping stripes

We discuss the flow past a flat heterogeneous solid surface decorated by slipping stripes. The spatially varying slip length, b(y), is assumed to be small compared to the scale of the heterogeneities, L, but finite. For such weakly slipping surfaces, earlier analyses have predicted that the effective slip length is simply given by the surface-averaged slip length, which implies that the effective slip-length tensor becomes isotropic. Here we show that a different scenario is expected if the local slip length has steplike jumps at the edges of slipping heterogeneities. In this case, the next-to-leading term in an expansion of the effective slip-length tensor in powers of max[b(y)/L] becomes comparable to the leading-order term, but anisotropic, even at very small b(y)/L. This leads to an anisotropy of the effective slip and to its significant reduction compared to the surface-averaged value. The asymptotic formulas are tested by numerical solutions and are in agreement with results of dissipative particle dynamics simulations.

Tensorial slip theory for gas flows and comparison with molecular dynamics simulations using an anisotropic gas-wall collision mechanism

In this paper we examine the anisotropic slip theory for gas flows based on tangential accommodation coefficients and compare it with molecular dynamics (MD) results. A special gas-wall boundary condition is employed within MD simulations to mimic the anisotropic gas-wall collision mechanism. Results from MD simulations with different surface orientations show good agreement with the slip quantification proposed in this work.

Flow past superhydrophobic surfaces with cosine variation in local slip length

Anisotropic superhydrophobic surfaces have the potential to greatly reduce drag and enhance mixing phenomena in microfluidic devices. Recent work has focused mostly on cases of superhydrophobic stripes. Here, we analyze a relevant situation of cosine variation of the local slip length. We derive approximate formulas for maximal (longitudinal) and minimal (transverse) directional effective slip lengths that are in good agreement with the exact numerical solution and lattice-Boltzmann simulations. Compared to the case of superhydrophobic stripes, the cosine texture can provide a very large effective slip. However, the difference between eigenvalues of the slip-length tensor is smaller, indicating that the flow is less anisotropic.

Effects of surface morphology and anisotropy on the tangential-momentum accommodation coefficient between Pt(100) and Ar

In this paper, we study the influence of platinum (100) surface morphology on the tangential-momentum accommodation coefficient with argon using a molecular dynamics method. The coefficient is computed directly by beaming Ar atoms onto the surfaces and measuring the relative momentum changes. The wall is maintained at a constant temperature and its interaction with the gas atoms is governed by the Kulginov potential. To capture correctly the surface effect of the walls and the atoms' trajectories, the quantum Sutton-Chen multibody potential is employed between the Pt atoms. The effects of wall surface morphology, incident direction, and temperature are considered in this work and provide full information on the gas-wall interaction.