Effect of charge inversion on nanoconfined flow of multivalent ionic solutions.

A comprehensive understanding of fluid dynamics of dilute electrolyte solutions in nanoconfinement is essential to develop more efficient nanofluidic devices. In nanoconduits, the electrical double layer can occupy a considerable part of the channel cross-section, therefore the transport properties of a nanoconfined electrolyte solution can be altered by interfacial phenomena such as the charge inversion (CI). CI is an electrokinetic effect that has been associated with the presence of hydrated multivalent cations in nanoconfinement. Here, all-atom molecular dynamics simulations are employed to study the structure and dynamics of aqueous multivalent electrolyte solutions within slit-shaped silica channels. All simulations are conducted for more than 100 ns to capture the equilibrium ion distribution, the interfacial hydrodynamic properties, and to reveal the influence of CI on nanoconfined fluid transport. The electrolyte solutions consist of water as solvent, chloride as co-ion and different counter-ions, i.e., sodium, magnesium and aluminum. We find that the interfacial viscosity is related to the concentration and valence of the counter-ions in the solution. Our results suggest that higher CI is correlated to the presence of a layer of fluid with augmented viscosity adjacent to the channel wall. As the thickness of this interfacial high-viscosity fluid increases, lower flow rates are measured whereas higher interfacial viscosities and friction coefficients are computed.

Effect of an external electric field on capillary filling of water in hydrophilic silica nanochannels.

Development of functional nanofluidic devices requires understanding the fundamentals of capillary driven flow in nanochannels. In this context, we conduct molecular dynamics simulations of water capillary imbibition in silica nanoslits under externally applied electric (E) fields with strengths between 0 and 1 V nm-1. For increasing E-fields, we observe a systematic lowering in the meniscus contact angle and a decrease in the corresponding water filling rates. These results contrast markedly the classical Washburn-Bosanquet's equation which predicts an increase in filling rates for lower water contact angles. Our study provides evidence that the observed decrease in water filling rates can be attributed to the interplay between two underlying mechanisms, a reduced fluidity of interfacial water and a systematic alignment of the water molecules in the bulk as a response to the particular strength of the applied E-field. Therefore, during water capillary filling a constant E-field applied in the direction parallel to the water imbibition leads to a lower than expected filling rate caused by a viscosity increase in the bulk and an altered solid-liquid friction on the channel walls. These coupled mechanisms governing capillarity under the action of applied E-fields could be manipulated for controlling imbibition of polar liquid solutions in nanofluidic devices.

Water thermophoresis in carbon nanotubes: the interplay between thermophoretic and friction forces.

Thermophoresis is the phenomenon wherein particles experience a net drift induced by a thermal gradient. In this work, molecular dynamics simulations are conducted to study with atomistic detail the thermophoresis of water nanodroplets inside carbon nanotubes (CNTs) and its interplay with the retarding liquid-solid friction. Different applied temperatures, thermal gradients, and droplet sizes are used to reveal the dynamics of the two kinetic regimes of the thermophoretic motion in CNTs. The results indicate that during the droplet motion, the thermophoretic force is independent of the velocity of the droplet, whereas the magnitude of the retarding friction force exhibits a linear dependence. In fact, in the initial regime the magnitude of the friction force increases linearly with the droplet velocity, until the thermophoretic force is balanced by the friction force as the droplet reaches its terminal velocity in the final regime. In addition, an increase in the magnitude of the thermophoretic force is found for longer water droplets. These findings provide a deeper understanding of liquid transport driven by temperature gradients in nanoconfined geometries where liquid-solid interfaces govern fluidics.

Slip divergence of water flow in graphene nanochannels: the role of chirality.

Graphene has attracted considerable attention due to its characteristics as a 2D material and its fascinating properties, providing a potential building block for nanofabrication. In nanochannels the solid-liquid interface plays a non-negligible role in determining the fluid dynamics. Therefore, for an optimal design of nanofluidic devices, a comprehensive understanding of the slippage in a water flow confined between graphene walls is important. In nanoconfinement, experimental and computational studies have found the slip length to increase nonlinearly when the shear rate is larger than a critical value. Here, by conducting molecular dynamics simulations, we study the influence of the graphene crystallographic orientation on the slip boundary conditions inside a nanoslit channel. The flow in channels with heights of 2.0, 2.4 and 2.8 nm is driven parallel to the zig-zag and arm-chair crystallographic directions. We extract flow rates, velocity profiles, slip velocities and slip lengths. The slip velocity displays a linear relationship to the shear stress up to a critical value, which is not size dependent. Moreover, the slip length is found to be shear stress dependent above a critical shear stress value of 0.4 MPa. Furthermore, our results indicate that after this critical shear stress is reached, the flow rates are significantly influenced (up to 10%) by the particular orientation of the graphene topology.

Effect of the meniscus contact angle during early regimes of spontaneous imbibition in nanochannels.

Nanoscale capillarity has been extensively investigated; nevertheless, many fundamental questions remain open. In spontaneous imbibition, the classical Lucas-Washburn equation predicts a singularity as the fluid enters the channel consisting of an anomalous infinite velocity of the capillary meniscus. Bosanquet's equation overcomes this problem by taking into account fluid inertia predicting an initial imbibition regime with constant velocity. Nevertheless, the initial constant velocity as predicted by Bosanquet's equation is much greater than those observed experimentally. In the present study, large scale atomistic simulations are conducted to investigate capillary imbibition of water in slit silica nanochannels with heights between 4 and 18 nm. We find that the meniscus contact angle remains constant during the inertial regime and its value depends on the height of the channel. We also find that the meniscus velocity computed at the channel entrance is related to the particular value of the meniscus contact angle. Moreover, during the subsequent visco-inertial regime, as the influence of viscosity increases, the meniscus contact angle is found to be time dependent for all the channels under study. Furthermore, we propose an expression for the time evolution of the dynamic contact angle in nanochannels which, when incorporated into Bosanquet's equation, satisfactorily explains the initial capillary rise.

Early regimes of water capillary flow in slit silica nanochannels.

Molecular dynamics simulations are conducted to investigate the initial stages of spontaneous imbibition of water in slit silica nanochannels surrounded by air. An analysis is performed for the effects of nanoscopic confinement, initial conditions of liquid uptake and air pressurization on the dynamics of capillary filling. The results indicate that the nanoscale imbibition process is divided into three main flow regimes: an initial regime where the capillary force is balanced only by the inertial drag and characterized by a constant velocity and a plug flow profile. In this regime, the meniscus formation process plays a central role in the imbibition rate. Thereafter, a transitional regime takes place, in which, the force balance has significant contributions from both inertia and viscous friction. Subsequently, a regime wherein viscous forces dominate the capillary force balance is attained. Flow velocity profiles identify the passage from an inviscid flow to a developing Poiseuille flow. Gas density profiles ahead of the capillary front indicate a transient accumulation of air on the advancing meniscus. Furthermore, slower capillary filling rates computed for higher air pressures reveal a significant retarding effect of the gas displaced by the advancing meniscus.