The coupling of surface charge and boundary slip at the solid-liquid interface and their combined effect on fluid drag: A review

Dependence of the formation kinetics of carbon dioxide hydrate on clay aging for solid carbon dioxide storage

Clay-based marine sediments have great potential for safe and effective carbon dioxide (CO2) encapsulation by storing enormous amounts of CO2 in solid gas hydrate form. However, the aging of clay with time changes the surface properties of clay and complicates the CO2 hydrate formation behaviors in sediments. Due to the long clay aging period, it is difficult to identify the role of clay aging in the formation of CO2 hydrate in marine sediments. Here, we used ultrasonication and plasma treatment to simulate the breakage and oxidation of clay nanoflakes in aging and investigated the influence of clay aging on CO2 hydrate formation kinetics. We found that the breakage and oxidation of clay nanoflakes would disrupt the siloxane rings and graft hydroxyl on the clay nanoflakes. This decreased the negative charge density of clay nanoflakes and weakened the interfacial interaction of clay nanoflakes with the surrounding water. Therefore, the small clay nanoflakes enriched in hydroxyl would disrupt the surrounding tetrahedral water structure analogous to the CO2 hydrate, resulting in the prolongation of CO2 hydrate nucleation. These results revealed the influence of the structure-function relationship of clay nanoflakes with CO2 hydrate formation and are favorable for the development of hydrate-based CO2 storage.

Dual mechanisms based on synergistic effects of evaporation potential and streaming potential for natural water evaporation

Electricity generation by natural water evaporation generators (NWEGs) using porous materials shows great potential for energy harvesting, but mechanistic investigations of NWEGs have mostly been limited to streaming potential studies. In this study, we propose the coexistence of an evaporation potential and streaming potential in a NWEG using ZSM-5 as the generation material. The iron probe method, salt concentration regulation, solution regulation, and side evaporation area regulation were used to analyze the NWEG mechanism. Our findings revealed that a streaming potential formed as water flowed inside the ZSM-5 nanochannels, driven by electrodynamic effects that increased from the bottom to the top of the generator. In addition, an evaporation potential existed at the surface interface between ZSM-5 and water, which decreased from the bottom to the top as the evaporation height of the generator increased. The resulting open-circuit voltage (Voc) depended on the balance between the evaporation and streaming potentials, both of which were influenced by the evaporation enthalpy (Ee) or vapor pressure. Generally, a higher Ee or lower vapor pressure led to a lower evaporation potential and subsequently a lower Voc. A dual mechanism involving synergistic evaporation potential and streaming potential is proposed to explain the mechanism of NWEGs.

Effect of the Surface Charge-Dependent Boundary Slip on the Electrophoresis of a Hydrophobic Polarizable Rigid Colloid

The electrophoresis of a hydrophobic charged rigid colloid is studied by considering the lateral movement of the adsorbed surface charge. The slip velocity condition at the hydrophobic surface is modified to take into account the impact of the frictional and electric forces created by the adsorbed laterally mobile surface charge. Though the dependency of the surface charge on the slip velocity in the context of electrophoresis has been addressed before, the effect of the laterally mobile adsorbed surface charge on the electrophoresis of hydrophobic colloids has not been studied. The dielectric colloid is considered to polarize and create an induced immobile surface charge when subjected to an imposed electric field. The impact of the mobile surface charge along with the immobile induced surface charge on electrophoresis of a hydrophobic colloid is elucidated by numerically solving the governing electrokinetic equations in their full form. We have also developed a simplified model under a weak applied field consideration, which can be further reduced to a closed-form analytic expression for the mobility under the Debye-Hückel approximation. This analytic model for mobility is in excellent agreement with the exact numerical solution for an entire range of the Debye length when the ζ-potential is in the order of the thermal potential. One of the notable features of this closed-form mobility expression is that it accounts for the mobile adsorbed surface charge on the hydrodynamic slip condition and the dielectric polarization of the particle. We find that the mobility of the surface charge decreases the electrophoretic mobility of the hydrophobic dielectric colloid. However, the mobile surface charge enhances the mobility of a conducting hydrophobic colloid.

Thermal and Flow Analysis of Fully Developed Electroosmotic Flow in Parallel-Plate Micro- and Nanochannels with Surface Charge-Dependent Slip

This study analytically investigates the coupled effects of surface charge and boundary slip on the fully developed electroosmotic flow and thermal transfer in parallel plate micro and nanochannels under the high zeta potential. The electric potential, velocity, temperature, flow rate, and Nusselt number are obtained analytically. The main results are that the velocity of bulk flow is significantly reduced in the presence of the surface charge-dependent slip. Moreover, the maximum velocity at ζ = -125 mV is approximately twice as large as that at ζ = -25 mV. The velocity and dimensionless temperature increase as the zeta potential increases. The dimensionless temperature of the surface charge-dependent slip flow is larger than that of the surface charge-independent slip flow. For the surface charge-dependent slip flow, the maximum temperature at ζ = -125 mV is approximately four times larger than that at ζ = -25 mV. The Nusselt number decreases with Joule heating and increases with a positive heat transfer coefficient. The Nusselt number decreases as the electric field and the magnitude of the zeta potential increase. In the surface charge-dependent slip flows, the Nusselt number is smaller than that in the surface charge-independent slip flows.

Passive solute separation in AC electroosmosis including surface charge-coupled hydrodynamic slip effects

Separation of electrically neutral, mutually noninteracting passive solutes via AC electroosmotic slit channel flows is investigated for general asymmetric wall surface zeta potentials and apparent hydrodynamic slip lengths. We consider the nontrivial coupling between the surface potentials (or charge densities) and the apparent slip lengths, and focus our attention on the occurrence of a so called "crossover phenomenon" for separating out the slow diffusers when both slow and fast diffusers are present. Results show that regardless of the potential-slip coupling, wider bandwidths become available for crossover phenomenon to occur when the electroosmotic velocity gradient (magnitude) is greater. Contrarily, plug-like velocity profiles inhibit crossover phenomenon, and the potential-slip parametric combinations leading to such profiles can be easily identified by the conditions for minimal transport enhancement reported in recent literature. When separating out the slow diffuser or crossover phenomenon is desired, we recommend incorporating significant asymmetry in the surface potential and apparent slip boundary conditions such that the operating frequency and flow oscillation amplitude may be lowered to more practical values. Our results also agree with and strengthen the physical picture for explaining crossover phenomenon in macroscopic pressure-driven oscillatory flows.

AFM Slip Length Measurements for Water at Selected Phyllosilicate Surfaces

Most reported slip length measurements have been made at the surfaces of synthetic materials and modified synthetic materials. In contrast, few slip length measurements at the surface of unmodified natural mineral surfaces have been reported. In this regard, flow at the silica face surfaces of the phyllosilicate minerals, talc and mica, was considered. A slip boundary condition was expected at the nonpolar hydrophobic silica surface of talc leading to enhanced flow, and a no-slip boundary condition was expected at the hydrophilic silica surface of mica. Atomic force microscopy (AFM) slip length measurements were made at the talc and mica surfaces. The slip length results for the hydrophobic silica surface of talc were contrasted to the results for the hydrophilic silica surface of mica (no-slip flow). The results are discussed based on molecular dynamics simulations (MDS), as reported in the literature, and AFM images of surface nanobubbles. For nonpolar hydrophobic surfaces (such as talc), it is doubtful that the MDS interfacial water structure and the water exclusion zone (3.2 Å) account for the AFM slip flow with slip lengths as great as 95 nm. Rather, a better explanation for the AFM slip flow condition is based on reduced interfacial viscosity due to the presence of dissolved gas and the accommodation of pancake nanobubbles at the talc surface having a height dimension of magnitude similar to the slip length.

Effect of pH on Effective Slip Length and Surface Charge at Solid-Oil Interfaces of Roughness-Induced Surfaces

Featured Application

Drag reduction in many applications of micro-/nano-fluidic channel.

Abstract

In the development of micro/nano fluid control systems, fluid resistance has always been one of the key factors restricting its development. According to previous studies, it is found that the boundary slip effect of the solid-liquid interface can effectively reduce the resistance of the microfluid and improve the transport efficiency of the microfluid. The boundary slip length is mainly affected by surface wettability, roughness, and surface charge density. Among them, the influence mechanism of surface charge density on the boundary slip is the most complicated, and there is a lack of relevant research, and further investigation is needed. In this paper, we present research on quantification of effective slip length and surface charge density, where the roughness effect is considered. The electrostatic and hydrodynamic force data obtained from atomic force microscopy (AFM) measurements were fitted and processed for comparative analysis. We obtained the variation of surface charge density and effective slip length when different oleophobic surface samples were immersed in ethylene glycol with different pH values. The effect of pH on the surface charge density and effective slip length was investigated by their variations. The mechanism of the effect of pH on the surface charge density was discussed. The experimental results show that in the ethylene glycol solution, no matter whether the pH value of the solution increases or decreases, the charge density of the surface with the same properties decreases, and the effective boundary slip length also shows a downward trend. In deionized water, the surface charge density and effective boundary slip length decreases with the decrease of PH value.

Liquid-Solid Slip on Charged Walls: The Dramatic Impact of Charge Distribution

Nanofluidic systems show great promise for applications in energy conversion, where their performance can be enhanced by nanoscale liquid-solid slip. However, efficiency is also controlled by surface charge, which is known to reduce slip. Combining molecular dynamics simulations and analytical developments, we show the dramatic impact of surface charge distribution on the slip-charge coupling. Homogeneously charged graphene exhibits a very favorable slip-charge relation (rationalized with a new theoretical model correcting some weaknesses of the existing ones), leading to giant electrokinetic energy conversion. In contrast, slip is strongly affected on heterogeneously charged surfaces, due to the viscous drag induced by counterions trapped on the surface. In that case slip should depend on the detailed physical chemistry of the interface controlling the fraction of bound ions. Our numerical results and theoretical models provide new fundamental insight into the molecular mechanisms of liquid-solid slip, and practical guidelines for searching new functional interfaces with optimal energy conversion properties, e.g., for blue energy or waste heat harvesting.

Functionalized Surface-Charged SiO2 Nanoparticles Induce Pro-Inflammatory Responses, but Are Not Lethal to Caco-2 Cells

Nanoparticles (NPs) are widely used in food, and analysis of their potential gastrointestinal toxicity is necessary. The present study was designed to determine the effects of silica dioxide (SiO₂), titanium dioxide (TiO₂), and zinc oxide (ZnO) NPs on cultured THP-1 monocyte-derived macrophages and human epithelial colorectal adenocarcinoma (Caco-2) cells. Exposure to ZnO NPs for 24 h increased the production of redox response species (ROS) and reduced cell viability in a dose-dependent manner in THP-1 macrophages and Caco-2 cells. Although TiO₂ and SiO₂ NPs induced oxidative stress, they showed no apparent cytotoxicity against both cell types. The effects of functionalized SiO₂ NPs on undifferentiated and differentiated Caco-2 cells were investigated using fluorescently labeled SiO₂ NPs with neutral, positive, or negative surface charge. Exposure of both types of cells to the three kinds of SiO₂ NPs significantly increased their interaction in a dose-dependent manner. The largest interaction with both types of cells was noted with exposure to more negatively surface-charged SiO₂ NPs. Exposure to either positively or negatively, but not neutrally, surface-charged SiO₂ NPs increased NO levels in differentiated Caco-2 cells. Exposure of differentiated Caco-2 cells to positively or negatively surface-charged SiO₂ NPs also upregulated interleukin-8 expression. We conclude that functionalized surface-charged SiO₂ NPs can induce pro-inflammatory responses but are noncytotoxic.

Cross-Sectional Dimension Dependence of Electroosmotic Flow in Fractal Treelike Rectangular Microchannel Network

The present work theoretically and numerically studies the electroosmotic flow (EOF) within a fractal treelike rectangular microchannel network with uniform channel height. To obtain minimum EOF fluidic resistance, the microchannel cross-sectional dimensions of the fractal network are optimized. It is found that the cross-sectional dimension dependence of EOF fluidic resistance within a symmetric fractal network is only dependent on the channel width when the total channel volume is constant, and the optimal microchannel widths to reach the minimum EOF fluidic resistance satisfy the scaling law of κ = N^-1 (where κ is the width ratio of the rectangular channels at two successive branching levels, N is the branching number); however, for the symmetric fractal network with constant total surface area , the optimal cross-sectional dimensions should simultaneously satisfy κ = N^-1 and (where H is the channel height, S is the total channel surface area, l₀ is the channel length at the original branching level, γ is the channel length ratio at two successive branching levels and m is the total branching level) to obtain the minimum EOF fluidic resistance. The optimal scaling laws established in present work can be used for the optimization design of the fractal rectangular microchannel network for EOF to reach maximum transport efficiency.

Giant Thermoelectric Response of Nanofluidic Systems Driven by Water Excess Enthalpy

Nanofluidic systems could in principle be used to produce electricity from waste heat, but current theoretical descriptions predict a rather poor performance as compared to thermoelectric solid materials. Here we investigate the thermoelectric response of NaCl and NaI solutions confined between charged walls, using molecular dynamics simulations. We compute a giant thermoelectric response, 2 orders of magnitude larger than the predictions of standard models. We show that water excess enthalpy-neglected in the standard picture-plays a dominant role in combination with the electro-osmotic mobility of the liquid-solid interface. Accordingly, the thermoelectric response can be boosted using surfaces with large hydrodynamic slip. Overall, the heat harvesting performance of the model systems considered here is comparable to that of the best thermoelectric materials, and the fundamental insight provided by molecular dynamics suggests guidelines to further optimize the performance, opening the way to recycle waste heat using nanofluidic devices.

Effects of surface charge and boundary slip on time-periodic pressure-driven flow and electrokinetic energy conversion in a nanotube

Time-periodic pressure-driven slip flow and electrokinetic energy conversion efficiency in a nanotube are studied analytically. The slip length depends on the surface charge density. Electric potential, velocity and streaming electric field are obtained analytically under the Debye-Hückel approximation. The electrokinetic energy conversion efficiency is computed using these results. The effects of surface charge-dependent slip and electroviscous effect on velocity and electrokinetic energy conversion efficiency are discussed. The main results show that the velocity amplitude and the electrokinetic energy conversion efficiency of the surface charge-dependent slip flow are reduced compared with those of the surface charge-independent slip flow.

Establishment of a Standard Method for Boundary Slip Measurement on Smooth Surfaces Based on AFM

Typically, it is difficult to analyze and design a micro/nanofluid system, and the design process cannot follow the traditional law of hydrodynamics. The boundary condition is very important in the applications of a micro/nanofluid system. The existence of boundary slip can reduce the hydrodynamic resistance and enhance fluid flow. How to accurately determine the dynamic boundary conditions is increasingly concerned by researchers. Atomic force microscope (AFM) is proven to be the most advanced experimental instrument for studying the characteristics of the surface and the interaction interface. Most studies on the application of atomic force microscopy to the measurement of the boundary slip do not describe a systematic standard process, leading to many differences in the measurement results. In this paper, a standard process of measuring slip on smooth and flat surfaces is developed, including the data processing methods that minimize the interference factors in the original data as well as simplify the data expression. Thus, the boundary slip can be obtained more easily and accurately.

Role of Solid Wall Properties in the Interface Slip of Liquid in Nanochannels

A two-dimensional molecular dynamics model of the liquid flow inside rough nanochannels is developed to investigate the effect of a solid wall on the interface slip of liquid in nanochannels with a surface roughness constructed by rectangular protrusions. The liquid structure, velocity profile, and confined scale on the boundary slip in a rough nanochannel are investigated, and the comparison of those with a smooth nanochannel are presented. The influence of solid wall properties, including the solid wall density, wall-fluid coupling strength, roughness height and spacing, on the interfacial velocity slip are all analyzed and discussed. It is indicated that the rough surface induces a smaller magnitude of the density oscillations and extra energy losses compared with the smooth solid surface, which reduce the interfacial slip of liquid in a nanochannel. In addition, once the roughness spacing is very small, the near-surface liquid flow dominates the momentum transfer at the interface between liquid and solid wall, causing the role of both the corrugation of wall potential and wall-fluid coupling strength to be less obvious. In particular, the slip length increases with increasing confined scales and shows no dependence on the confined scale once the confined scale reaches a critical value. The critical confined scale for the rough channel is larger than that of the smooth scale.

Effective Boundary Slip Induced by Surface Roughness and Their Coupled Effect on Convective Heat Transfer of Liquid Flow

As a significant interfacial property for micro/nano fluidic system, the effective boundary slip can be induced by the surface roughness. However, the effect of surface roughness on the effective slip is still not clear, both increased and decreased effective boundary slip were found with increased roughness. The present work develops a simplified model to study the effect of surface roughness on the effective boundary slip. In the created rough models, the reference position of the rough surfaces to determinate effective boundary slip was set based on ISO/ASME standard and the surface roughness parameters including Ra (arithmetical mean deviation of the assessed profile), Rsm (mean width of the assessed profile elements) and shape of the texture varied to form different surface roughness. Then, the effective boundary slip of fluid flow through the rough surface was analyzed by using COMSOL 5.3. The results show that the effective boundary slip induced by surface roughness of fully wetted rough surface keeps negative and further decreases with increasing Ra or decreasing Rsm. Different shape of roughness texture also results in different effective slip. A simplified corrected method for the measured effective boundary slip was developed and proved to be efficient when the Rsm is no larger than 200 nm. Another important finding in the present work is that the convective heat transfer firstly increases followed by an unobvious change with increasing Ra, while the effective boundary slip keeps decreasing. It is believed that the increasing Ra enlarges the area of solid-liquid interface for convective heat transfer, however, when Ra is large enough, the decreasing roughness-induced effective boundary slip counteracts the enhancement effect of roughness itself on the convective heat transfer.

Automated image segmentation-assisted flattening of atomic force microscopy images

Atomic force microscopy (AFM) images normally exhibit various artifacts. As a result, image flattening is required prior to image analysis. To obtain optimized flattening results, foreground features are generally manually excluded using rectangular masks in image flattening, which is time consuming and inaccurate. In this study, a two-step scheme was proposed to achieve optimized image flattening in an automated manner. In the first step, the convex and concave features in the foreground were automatically segmented with accurate boundary detection. The extracted foreground features were taken as exclusion masks. In the second step, data points in the background were fitted as polynomial curves/surfaces, which were then subtracted from raw images to get the flattened images. Moreover, sliding-window-based polynomial fitting was proposed to process images with complex background trends. The working principle of the two-step image flattening scheme were presented, followed by the investigation of the influence of a sliding-window size and polynomial fitting direction on the flattened images. Additionally, the role of image flattening on the morphological characterization and segmentation of AFM images were verified with the proposed method.

Interface conditions of roughness-induced superoleophilic and superoleophobic surfaces immersed in hexadecane and ethylene glycol

Interface conditions are an important property that can affect the drag of fluid flow. For surfaces with different oleophobicity, the boundary slip at the solid-oil interface is mostly larger than that at the solid-water interface. Roughness is a key factor for the wettability of superoleophilic/superoleophobic surfaces, and it has been found to affect the effective value of slip length in measurements. Moreover, there are no studies on the effect of roughness on slip at interfaces between oil and superoleophilic/superoleophobic surfaces. A theoretical description of the real surface roughness is yet to be found. Results show that the effective slip length is negative and decreases with an increasing root mean squared (RMS) roughness of surfaces, as the increasing roughness enhances the area with discontinuous slip at the solid-liquid interface. The underlying mechanisms are analyzed. The amplitude parameters of surface roughness could significantly inhibit the degree of boundary slip on both superoleophilic surfaces in Wenzel state and superoleophobic surfaces in Cassie state immersed in oil. The oleic systems were likely to enhance boundary slip and resulted in a corresponding reduction in drag with decreasing roughness on the solid-oil interfaces.

Effect of Surface Charge on the Nanofriction and Its Velocity Dependence in an Electrolyte Based on Lateral Force Microscopy

The nanofriction between a silicon nitride probe and both a silicon wafer and an octadecyltrichlorosilane (OTS)-coated surface is studied in saline solution by using lateral force microscopy (LFM). The effects of surface charge on the nanofriction in an electrolyte as well as its velocity dependence are studied, while the surface charge at the solid-liquid interface is adjusted by changing the pH value of the electrolyte. The results show that the nanofrictional behavior between the probe and the samples in an electrolyte depends strongly on the surface charge at the solid-liquid interface. When the probe and the sample in the electrolyte are charged with the same sign, a repulsive electrostatic interaction between the probe and the sample is produced, leading to a reduction in nanofriction. In contrast, when the two surfaces are charged with the opposite sign, nanofriction is enhanced by the attractive electrostatic interaction between the probe and the sample. The velocity dependence of nanofriction in an electrolyte is believed to be tied to charge regulation referring to a decreasing trend in surface charge densities for the two approaching charged surfaces in an electrolyte. When the probe slides on the sample at a low velocity, charge regulation occurs and weakens the electrostatic interaction between the probe and the sample. As a result, nanofriction is reduced for surfaces charged with the opposite sign, and it is enhanced for surfaces charged with the same sign. When the sliding velocity between the probe and the sample is high, there is insufficient time for charge regulation to occur. Thus, the friction pair shows a larger nanofriction when the surfaces are charged with the opposite sign and a smaller nanofriction when the surfaces are charged with the same sign when compared to the case of a lower sliding velocity.

Utilizing the Discrete Element Method for the Modeling of Viscosity in Concentrated Suspensions

The rheological behavior of concentrated suspensions is a complicated problem because it originates in the collective motion of particles and their interaction with the surrounding fluid. For this reason, it is difficult to accurately model the effect of various system parameters on the viscosity even for highly simplified systems. We model the viscosity of a hard-sphere suspension subjected to high shear rates using the dynamic discrete element method (DEM) in three spatial dimensions. The contact interaction between particles was described by the Hertz model of elastic spheres (soft-sphere model), and the interaction of particles with flow was accounted for by the two-way coupling approach. The hydrodynamic interaction between particles was described by the lubrication theory accounting for the slip on particle surfaces. The viscosity in a simple-shear model was evaluated from the force balance on the wall. The obtained results are in close agreement with literature data for systems with hard spheres. Namely, the viscosity is shown to be independent of shear rate and primary particle size for monodisperse suspensions. In accordance with theory and experimental data, the viscosity grows rapidly with particle volume fraction. We show that this rheological behavior is predominantly caused by the lubrication forces. A novel approach based on the slip of water on a particle surface was developed to overcome the divergent behavior of lubrication forces. This approach was qualitatively validated with literature data from AFM measurements using a colloidal probe. The model presented in this work represents a new, robust, and versatile approach to the modeling of viscosity in suspensions with the possibility to include various interaction models and study their effect on viscosity.

Electroviscous effect on fluid drag in a microchannel with large zeta potential

The electroviscous effect has been widely studied to investigate the effect of surface charge-induced electric double layers (EDL) on the pressure-driven flow in a micro/nano channel. EDL has been reported to reduce the velocity of fluid flow and increase the fluid drag. Nevertheless, the study on the combined effect of EDL with large zeta potential up to several hundred millivolts and surface charge depenedent-slip on the micro/nano flow is still needed. In this paper, the nonlinear Poisson-Boltzmann equation for electrical potential and ion distribution in non-overlapping EDL is first analytically solved. Then, the modified Navier-Stokes equation for the flow considering the effect of surface charge on the electrical conductivity of the electrolyte and slip length is analytically solved. This analysis is used to study the effect of non-overlapping EDL with large zeta potential on the pressure-driven flow in a microchannel with no-slip and charge-dependent slip conditions. The results show that the EDL leads to an increase in the fluid drag, but that slip can reduce the fluid drag. When the zeta potential is large enough, the electroviscous effect disappears for flow in the microchannel under a no-slip condition. However, the retardation of EDL on the flow and the enhancement of slip on the flow counteract each other under a slip condition. The underlying mechanisms of the effect of EDL with large zeta potential on fluid drag are the high net ionic concentration near the channel wall and the fast decay of electrical potential in the EDL when the zeta potential is large enough.