Surface roughness-hydrophobicity coupling in microchannel and nanochannel flows

Analysis of Electroviscous Effect for Flow of Micropolar Fluids in a Nanochannel with Overlapping Electric Double Layers at High Zeta Potential

Flows in nanofluidic channels under the influence of pressure gradient often lead to the overlapping of the electrical double layers (EDLs) augmenting the electroviscous effect. In this work, we analyze the electroviscous effect for the flow of a micropolar fluid in a parallel plate nanochannel, considering EDL overlapping and interfacial slip. Closed-form expressions of EDL potential, velocity, microrotation, yielded streaming potential, and volumetric flow rate are derived semi-analytically for high zeta potential without accounting for the Boltzmann distribution. We observe a significant change in streaming potential with inverse ionic Peclet number (R) for its lower values, and the range of streaming potential is more for thicker EDL. The micropolarity parameter does not have any influence on the streaming potential for weaker EDL overlapping at a lower R value. For thicker EDL, velocity decreases with the micropolarity parameter, while it increases drastically with interfacial slip. However, velocity shows a non-monotonic behavior with interfacial slip and micropolarity parameter for thinner EDL. The microrotation remains invariant with interfacial slip for thicker EDL, whereas its magnitude decreases with slip length for thinner EDL. Although the sensitivity of the flow rate on slip (Qs*) increases with R for thicker EDL, the behavior is non-monotonic for thinner EDL. Furthermore, Qs* varies non-monotonically with the micropolarity parameter at higher couple stress parameter values and vice versa. In fact, the volumetric flow rate is highly sensitive to slip for thicker EDL overlapping.

Bio-inspired microfluidics: A review

Biomicrofluidics, a subdomain of microfluidics, has been inspired by several ideas from nature. However, while the basic inspiration for the same may be drawn from the living world, the translation of all relevant essential functionalities to an artificially engineered framework does not remain trivial. Here, we review the recent progress in bio-inspired microfluidic systems via harnessing the integration of experimental and simulation tools delving into the interface of engineering and biology. Development of "on-chip" technologies as well as their multifarious applications is subsequently discussed, accompanying the relevant advancements in materials and fabrication technology. Pointers toward new directions in research, including an amalgamated fusion of data-driven modeling (such as artificial intelligence and machine learning) and physics-based paradigm, to come up with a human physiological replica on a synthetic bio-chip with due accounting of personalized features, are suggested. These are likely to facilitate physiologically replicating disease modeling on an artificially engineered biochip as well as advance drug development and screening in an expedited route with the minimization of animal and human trials.

Breakthrough pressure of oil displacement by water through the ultra-narrow kerogen pore throat from the Young-Laplace equation and molecular dynamic simulations

As one common unconventional reservoir, shale plays a pivotal role to compensate the depletion of conventional oil resources. There are numerous nanoscale pores and ultra-narrow pore throats (sub 2 nm) in shale media. To displace oil through ultra-narrow pore throats by water, one needs to overcome excessively-high capillary pressure. Understanding the water-oil two-phase displacement process through pore throats is critical to numerical simulation on tight/shale oil exploitation and ultimate oil recovery estimation. In this work, we use molecular dynamics simulations to investigate oil (represented by n-octane) displacement by water through a ~2 nm kerogen (represented by Type II-C kerogen) pore throat. Besides, the applicability of the Young-Laplace equation to the ultra-narrow kerogen pore throat has been assessed. We find that although the Type II-C kerogen is generally oil-wet, water has an excellent displacement efficiency without the oil film on the substrate, thanks to the hydrogen bonding formed between water and heteroatoms (such as O, N, and S) on the kerogen surface. Unlike previous studies, the capillary pressure obtained from the widely used Young-Laplace equation shows a good agreement with the breakthrough pressure obtained from MD simulations for the ∼2 nm kerogen pore throat. Our work indicates that explicitly considering intermolecular interactions as well as atomistic and molecular level characteristics is imperative to study the two-phase displacement process through ultra-narrow pore throats.

Multiscale analysis of the particles on demand kinetic model

We present a thorough investigation of the particles on demand kinetic model. After a brief introduction of the method, an appropriate multiscale analysis is carried out to derive the hydrodynamic limit of the model. In this analysis, the effect of the time-space dependent comoving reference frames are taken into account. This could be regarded as a generalization of the conventional Chapman-Enskog analysis applied to the lattice Boltzmann models which feature global constant reference frames. Further simulations of target benchmarks provide numerical evidence confirming the theoretical predictions.

Simulation of Water Flow in a Nanochannel with a Sudden Contraction or Expansion

Water flow in a nanoscale channel is known to be affected by strong water-wall interactions; as a result, the flow considerably deviates from the conventional continuum flow. Nanochannel with a sudden contraction or expansion is the most fundamental morphological nanostructure in many nanoporous systems such as shale matrix, mudrock, membrane, etc. However, the nanoconfinement effects of water flow in nanoporous systems and its effect on macroscopic flow behavior are still evolving research topics. In this work, our recently developed pore-scale lattice Boltzmann method (LBM) considering the nanoscale effects is extended to directly simulate water flow in nanoporous systems. The results show that the flow rate is dramatically decreased in hydrophobic nanopores because of additional flow resistances at the contraction and expansion junctions. This indicates that the bundle of capillary models or the permeability averaging method overestimates the water flow rate in nanoporous media if the contraction/expansion effects between different nanopores are ignored. This work highlights the importance of wettability of the nanochannel in the determination of dynamic water flow behaviors in the contraction/expansion nanosystem. Other important aspects, including velocity distribution, flow patterns, and vortex characteristics as well as pressure variation along the flow direction, are for the first time revealed and quantified. Large differences can be found comparing gas or larger-scale water flows through the same system. A new type of pressure variation curve along the axis of flow direction is found in the hydrophobic nanochannel with a sudden contraction/expansion. This work provides the fundamental understanding of water transport through the nanoscale system with contraction and expansion effects, giving implications to a wide range of industry applications.

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.

Impact of the prequench state of binary fluid mixtures on surface-directed spinodal decomposition

Using lattice Boltzmann simulations we investigate the impact of the amplitude of concentration fluctuations in binary fluid mixtures prior to demixing when in contact with a surface that is preferentially wet by one of the components. We find a bicontinuous structure near the surface for an initial, prequench state of the mixture close to the critical point where the amplitude of concentration fluctuations is large. In contrast, if the initial state of the mixture is not near the critical point and concentration fluctuations are relatively weak, then the morphology is not bicontinuous but remains layered until the very late stages of coarsening. In both cases, it is the morphology of a depletion layer rich in the nonpreferred component that dictates the growth exponent of the thickness of the fluid layer that is in direct contact with the substrate. In the early stages of demixing, we find a growth exponent consistent with a value of 1/4 for a prequench state away from the critical point, which is different from the usual diffusive scaling exponent of 1/3 that we recover for a prequench state close to the critical point. We attribute this to the structure of a depletion layer that is penetrated by tubes of the preferred fluid, connecting the wetting layer to the bulk fluid even in the early stages if the initial state is characterized by concentration fluctuations that are large in amplitude. Furthermore, we find that in the late stages of demixing the flow through these tubes results in significant in-plane concentration variations near the substrate, leading to dropletlike structures with a concentration lower than the average concentration in the wetting layer. This causes a deceleration in the growth of the wetting layer in the very late stages of the demixing. Irrespective of the prequench state of the mixture, the late stages of the demixing process produce the same growth law for the layer thickness, with a scaling exponent of unity usually associated with the impact of hydrodynamic flow fields.

Kinetic Simulations of Compressible Non-Ideal Fluids: From Supercritical Flows to Phase-Change and Exotic Behavior

We investigate a kinetic model for compressible non-ideal fluids. The model imposes the local thermodynamic pressure through a rescaling of the particle’s velocities, which accounts for both long- and short-range effects and hence full thermodynamic consistency. The model is fully Galilean invariant and treats mass, momentum, and energy as local conservation laws. The analysis and derivation of the hydrodynamic limit is followed by the assessment of accuracy and robustness through benchmark simulations ranging from the Joule–Thompson effect to a phase-change and high-speed flows. In particular, we show the direct simulation of the inversion line of a van der Waals gas followed by simulations of phase-change such as the one-dimensional evaporation of a saturated liquid, nucleate, and film boiling and eventually, we investigate the stability of a perturbed strong shock front in two different fluid mediums. In all of the cases, we find excellent agreement with the corresponding theoretical analysis and experimental correlations. We show that our model can operate in the entire phase diagram, including super- as well as sub-critical regimes and inherently captures phase-change phenomena.

Mesoscopic method to study water flow in nanochannels with different wettability

Molecular dynamics (MD) simulations is currently the most popular and credible tool to model water flow in nanoscale where the conventional continuum equations break down due to the dominance of fluid-surface interactions. However, current MD simulations are computationally challenging for the water flow in complex tube geometries or a network of nanopores, e.g., membrane, shale matrix, and aquaporins. We present a novel mesoscopic lattice Boltzmann method (LBM) for capturing fluctuated density distribution and a nonparabolic velocity profile of water flow through nanochannels. We incorporated molecular interactions between water and the solid inner wall into LBM formulations. Details of the molecular interactions were translated into true and apparent slippage, which were both correlated to the surface wettability, e.g., contact angle. Our proposed LBM was tested against 47 published cases of water flow through infinite-length nanochannels made of different materials and dimensions-flow rates as high as seven orders of magnitude when compared with predictions of the classical no-slip Hagen-Poiseuille (HP) flow. Using the developed LBM model, we also studied water flow through finite-length nanochannels with tube entrance and exit effects. Results were found to be in good agreement with 44 published finite-length cases in the literature. The proposed LBM model is nearly as accurate as MD simulations for a nanochannel, while being computationally efficient enough to allow implications for much larger and more complex geometrical nanostructures.

Hydrodynamic behavior of the pseudopotential lattice Boltzmann method for interfacial flows

The lattice Boltzmann method (LBM) is routinely employed in the simulation of complex multiphase flows comprising bulk phases separated by nonideal interfaces. The LBM is intrinsically mesoscale with a hydrodynamic equivalence popularly set by the Chapman-Enskog analysis, requiring that fields slowly vary in space and time. The latter assumptions become questionable close to interfaces where the method is also known to be affected by spurious nonhydrodynamical contributions. This calls for quantitative hydrodynamical checks. In this paper, we analyze the hydrodynamic behavior of the LBM pseudopotential models for the problem of the breakup of a liquid ligament triggered by the Plateau-Rayleigh instability. Simulations are performed at fixed interface thickness, while increasing the ligament radius, i.e., in the "sharp interface" limit. The influence of different LBM collision operators is also assessed. We find that different distributions of spurious currents along the interface may change the outcome of the pseudopotential model simulations quite sensibly, which suggests that a proper fine-tuning of pseudopotential models in time-dependent problems is needed before the utilization in concrete applications. Taken all together, we argue that the results of the proposed paper provide a valuable insight for engineering pseudopotential model applications involving the hydrodynamics of liquid jets.

Molecular dynamics simulation of nanofluidics

This review reports the progress on the recent development of molecular dynamics simulation of nanofluidics. Molecular dynamics simulations of nanofluidics in nanochannel structure, surface roughness of nanochannel, carbon nanotubes, electrically charged, thermal transport in nanochannels and gases in nanochannels are illustrated and discussed. This paper will provide an expedient and valuable reference to designers who intend to research molecular dynamics simulation of nanofluidic devices.

Particles on Demand for Kinetic Theory

A novel formulation of fluid dynamics as a kinetic theory with tailored, on-demand constructed particles removes restrictions on flow speed and temperature as compared to its predecessors, the lattice Boltzmann methods and their modifications. In the new kinetic theory, discrete particles are determined by a rigorous limit process which avoids ad hoc assumptions about their velocities. Classical benchmarks for incompressible and compressible flows demonstrate that the proposed discrete-particles kinetic theory opens up an unprecedented wide domain of applications for computational fluid dynamics.

Flagellated microswimmers: Hydrodynamics in thin liquid films

The hydrodynamics of a flagellated microswimmer moving in thin films is discussed. The fully resolved hydrodynamics is exploited by solving the Stokes equations for the actual geometry of the swimmer. Two different interfaces are used to confine the swimmer: a bottom solid wall and a top air-liquid interface, as appropriate for a thin film. The swimmer follows curved clockwise trajectories that can converge towards an asymptotically stable circular path or can result in a collision with one of the two interfaces. A bias towards the air-liquid interface emerges. Slight changes in the swimmer geometry and film thickness strongly affect the resulting dynamics suggesting that a very reach phenomenology occurs in the presence of confinement. Under specific conditions, the swimmer follows a "crown-like" path. Implications for the motion of bacteria close to an air bubble moving in a microchannel are discussed.

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.

Magnetic-field-driven alteration in capillary filling dynamics in a narrow fluidic channel

We investigated pressure-driven transport of an immiscible binary system, constituted by two electrically conducting liquids, in a narrow fluidic channel under the influence of an externally applied magnetic field. The surface wettability was taken into account in the analysis considering that the walls of the channel are chemically treated to obtain various predefined contact angles as required for the study. Alterations in the capillary filling and wetting dynamics in the channel stemming from a complex interplay among different forces acting over the interface were investigated. It was shown that an alteration in the strength of the magnetic field leads to an alteration in the dynamics of the interface, which in turn, alters the filling and wetting dynamics nontrivially upon interaction with the surface tension force due to the wetted walls of the channel. It is emphasized that a contrast in properties of constituents of the binary system gives rise to an alteration in the forces being applied across the interface, leading to an intricate control over the filling and wetting dynamics for a given flow configuration and an applied field strength. We believe that the results obtained from this analysis may aid the design of microfluidic devices used for multiphase transport.

Transport and adsorption under liquid flow: the role of pore geometry

We study here the interplay between transport and adsorption in porous systems with complex geometries under fluid flow. Using a lattice Boltzmann scheme extended to take into account the adsorption at solid/fluid interfaces, we investigate the influence of pore geometry and internal surface roughness on the efficiency of fluid flow and the adsorption of molecular species inside the pore space. We show how the occurrence of roughness on pore walls acts effectively as a modification of the solid/fluid boundary conditions, introducing slippage at the interface. We then compare three common pore geometries, namely honeycomb pores, inverse opal, and materials produced by spinodal decomposition. Finally, we quantify the influence of those three geometries on fluid transport and tracer adsorption. This opens perspectives for the optimization of materials' geometries for applications in dynamic adsorption under fluid flow.

Lattice Boltzmann simulations of droplet formation in confined channels with thermocapillary flows

Based on mesoscale lattice Boltzmann simulations with the "Shan-Chen" model, we explore the influence of thermocapillarity on the breakup properties of fluid threads in a microfluidic T-junction, where a dispersed phase is injected perpendicularly into a main channel containing a continuous phase, and the latter induces periodic breakup of droplets due to the cross-flowing. Temperature effects are investigated by switching on-off both positive-negative temperature gradients along the main channel direction, thus promoting a different thread dynamics with anticipated-delayed breakup. Numerical simulations are performed at changing the flow rates of both the continuous and dispersed phases, as well as the relative importance of viscous forces, surface tension forces, and thermocapillary stresses. The range of parameters is broad enough to characterize the effects of thermocapillarity on different mechanisms of breakup in the confined T-junction, including the so-called "squeezing" and "dripping" regimes, previously identified in the literature. Some simple scaling arguments are proposed to rationalize the observed behavior, and to provide quantitative guidelines on how to predict the droplet size after breakup.

Chimaera simulation of complex states of flowing matter

We discuss a unified mesoscale framework (chimaera) for the simulation of complex states of flowing matter across scales of motion. The chimaera framework can deal with each of the three macro-meso-micro levels through suitable 'mutations' of the basic mesoscale formulation. The idea is illustrated through selected simulations of complex micro- and nanoscale flows.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

Lattice Kinetic Theory in a Comoving Galilean Reference Frame

We prove that the fully discrete lattice Boltzmann method is invariant with respect to Galilean transformation. Based on this finding, a novel class of shifted lattices is proposed which dramatically increases the operating range of lattice Boltzmann simulations, in particular, for gas dynamics applications. A simulation of vortex-shock interaction is used to demonstrate the accuracy and efficiency of the proposed lattices. With one single algorithm it is now possible to simulate a broad range of applications, from low Mach number flows to transonic and supersonic flow regimes.

Effectiveness of the Young-Laplace equation at nanoscale

Using molecular dynamics (MD) simulations, a new approach based on the behavior of pressurized water out of a nanopore (1.3–2.7 nm) in a flat plate is developed to calculate the relationship between the water surface curvature and the pressure difference across water surface. It is found that the water surface curvature is inversely proportional to the pressure difference across surface at nanoscale, and this relationship will be effective for different pore size, temperature, and even for electrolyte solutions. Based on the present results, we cannot only effectively determine the surface tension of water and the effects of temperature or electrolyte ions on the surface tension, but also show that the Young-Laplace (Y-L) equation is valid at nanoscale. In addition, the contact angle of water with the hydrophilic material can be further calculated by the relationship between the critical instable pressure of water surface (burst pressure) and nanopore size. Combining with the infiltration behavior of water into hydrophobic microchannels, the contact angle of water at nanoscale can be more accurately determined by measuring the critical pressure causing the instability of water surface, based on which the uncertainty of measuring the contact angle of water at nanoscale is highly reduced.

Droplets Can Rebound toward Both Directions on Textured Surfaces with a Wettability Gradient

The impact of water droplets on superhydrophobic surfaces with a wettability gradient is studied using the lattice Boltzmann simulation. Droplets impacting such textured surfaces have been previously reported to rebound obliquely following the wettability gradient due to the unbalanced interfacial forces created by the heterogeneous architectures. Here we demonstrate that droplets can rebound toward both directions on textured surfaces with a wettability gradient. Our simulation results indicate that the rebound trajectory of droplets is determined by the competition between the lateral recoil of the liquid and the penetration and capillary emptying of the penetrated liquid from the textures in the vertical direction. When the time scale for the droplet penetration and capillary emptying process is smaller than the time for the lateral spreading, the droplet will rebound following the wettability gradient. By contrast, the droplet will display a bouncing against the wettability gradient direction because of the significant capillary penetration and emptying in the transverse direction. We believe that our study provides important insight for the design of micro/nanotextured surfaces for controlled droplet manipulation.

Effects of viscoelasticity on droplet dynamics and break-up in microfluidic T-Junctions: a lattice Boltzmann study

The effects of viscoelasticity on the dynamics and break-up of fluid threads in microfluidic T-junctions are investigated using numerical simulations of dilute polymer solutions at changing the Capillary number (Ca), i.e. at changing the balance between the viscous forces and the surface tension at the interface, up to Ca ≈ 3×10(-2). A Navier-Stokes (NS) description of the solvent based on the lattice Boltzmann models (LBM) is here coupled to constitutive equations for finite extensible non-linear elastic dumbbells with the closure proposed by Peterlin (FENE-P model). We present the results of three-dimensional simulations in a range of Ca which is broad enough to characterize all the three characteristic mechanisms of break-up in the confined T-junction, i.e. squeezing, dripping and jetting regimes. The various model parameters of the FENE-P constitutive equations, including the polymer relaxation time τP and the finite extensibility parameter L2, are changed to provide quantitative details on how the dynamics and break-up properties are affected by viscoelasticity. We will analyze cases with Droplet Viscoelasticity (DV), where viscoelastic properties are confined in the dispersed (d) phase, as well as cases with Matrix Viscoelasticity (MV), where viscoelastic properties are confined in the continuous (c) phase. Moderate flow-rate ratios Q ≈ O(1) of the two phases are considered in the present study. Overall, we find that the effects are more pronounced in the case with MV, as the flow driving the break-up process upstream of the emerging thread can be sensibly perturbed by the polymer stresses.

A lattice Boltzmann study of the effects of viscoelasticity on droplet formation in microfluidic cross-junctions

Based on mesoscale lattice Boltzmann (LB) numerical simulations, we investigate the effects of viscoelasticity on the break-up of liquid threads in microfluidic cross-junctions, where droplets are formed by focusing a liquid thread of a dispersed (d) phase into another co-flowing continuous (c) immiscible phase. Working at small Capillary numbers, we investigate the effects of non-Newtonian phases in the transition from droplet formation at the cross-junction (DCJ) to droplet formation downstream of the cross-junction (DC) (Liu and Zhang, Phys. Fluids. 23, 082101 (2011)). We will analyze cases with Droplet Viscoelasticity (DV), where viscoelastic properties are confined in the dispersed phase, as well as cases with Matrix Viscoelasticity (MV), where viscoelastic properties are confined in the continuous phase. Moderate flow-rate ratios Q≈O(1) of the two phases are considered in the present study. Overall, we find that the effects are more pronounced with MV, where viscoelasticity is found to influence the break-up point of the threads, which moves closer to the cross-junction and stabilizes. This is attributed to an increase of the polymer feedback stress forming in the corner flows, where the side channels of the device meet the main channel. Quantitative predictions on the break-up point of the threads are provided as a function of the Deborah number, i.e., the dimensionless number measuring the importance of viscoelasticity with respect to Capillary forces.

Elucidating Nonwetting of Re-Entrant Surfaces with Impinging Droplets

Superomniphobic surfaces display both superoleophobic and superhydrophobic properties, having a contact angle greater than 150° for both water and oil droplets. In this work, lattice Boltzmann simulations on droplets impacting the surface textures of various topologies are performed to understand the mechanism of how the superomniphobic properties can be achieved by optimizing the geometry of re-entrant surfaces and the inherent hydrophobicity of substrates. Detailed kinetics for droplet impinging is analyzed for both liquid impalement and emptying, showing distinct dependences on geometrical details of re-entrant surfaces. The origins of the enhanced stability of Cassie states are ascribed to (i) the barrier of the Cassie-to-Wenzel transition for the impalement process, (ii) the driving force for liquid receding in the emptying process, and (iii) the contact line pinning from the entrance effect and the edge effect. Finally, we check the strategies proposed here by designing a new re-entrance structure that possesses an excellent property in maintaining the droplet Cassie state.

Lattice Boltzmann investigation of droplet inertial spreading on various porous surfaces

The spreading of liquid drops on solid surfaces is a wide-spread phenomenon of both fundamental and industrial interest. In many applications, surfaces are porous and spreading patterns are very complex with respect to the case on smooth surfaces. Focusing on the inertial spreading just before the Tanner-like viscous regime, this work investigates the spreading of a low-viscosity droplet on a porous surface using lattice Boltzmann numerical simulations. The case of a flat surface is first considered, and it reveals a dependence on the solid equilibrium contact angle θ(s)(eq), which is in good agreement with published experimental data. We conducted numerical experiments with various surfaces perforated by a regular pattern of holes of infinite length. The results show that the global spreading dynamics is independent of the porosity morphology. Through the assumption that, for wetting, the pores can be regarded as surface patches with a contact angle of θ(pore)(eq)=180°, we deduce an effective equilibrium contact angle θ(eff)(eq) on the porous surface from the Cassie-Baxter law. A spreading model is then proposed to describe both a prefactor and an exponent that are similar to a flat surface whose equilibrium contact angle is θ(eff)(eq). This model compares satisfactorily with a large number of numerical experiments under varying conditions.

Polymer melt flow through nanochannels: from theory and fabrication to application

This short review presents the theory, fabrication, and application of polymer melts through nanochannels.

Lattice Boltzmann modeling of droplet condensation on superhydrophobic nanoarrays

Droplet nucleation and growth on superhydrophobic nanoarrays is simulated by employing a multiphase, multicomponent lattice Boltzmann (LB) model. Three typical preferential nucleation modes of condensate droplets are observed through LB simulations with various geometrical parameters of nanoarrays, which are found to influence the wetting properties of nanostructured surfaces significantly. The droplets nucleated at the top of posts (top nucleation) or in the upside interpost space of nanoarrays (side nucleation) will generate a nonwetting Cassie state, while the ones nucleated at the bottom corners between the posts of nanoarrays (bottom nucleation) produce a wetting Wenzel state. The simulated time evolutions of droplet pressures at different locations are analyzed, which offers insight into the underlying physics governing the motion of droplets growing from different nucleation modes. It is demonstrated that the nanostructures with taller posts and a high ratio of post height to interpost space (H/S) are beneficial to produce the top- and side-nucleation modes. The simulated wetting states of condensate droplets on the nanostructures, having various geometrical configurations, compare reasonably well with experimental observations. The established relationship between the geometrical parameters of nanoarrays and the preferential nucleation modes of condensate droplets provides guidance for the design of nanoarrays with desirable anticondensation superhydrophobic properties.

Beyond Cassie equation: local structure of heterogeneous surfaces determines the contact angles of microdroplets

The application of Cassie equation to microscopic droplets is recently under intense debate because the microdroplet dimension is often of the same order of magnitude as the characteristic size of substrate heterogeneities, and the mechanism to describe the contact angle of microdroplets is not clear. By representing real surfaces statistically as an ensemble of patterned surfaces with randomly or regularly distributed heterogeneities (patches), lattice Boltzmann simulations here show that the contact angle of microdroplets has a wide distribution, either continuous or discrete, depending on the patch size. The origin of multiple contact angles observed is ascribed to the contact line pinning effect induced by substrate heterogeneities. We demonstrate that the local feature of substrate structure near the contact line determines the range of contact angles that can be stabilized, while the certain contact angle observed is closely related to the contact line width.

Effects of the hierarchical structure of rough solid surfaces on the wetting of microdroplets

We used the lattice Boltzmann method to investigate how the hierarchical structure of a rough solid surface, which in this work is modeled as the microstructure (micropillars) covered with nanostructures (nanopillars), affects the contact angle of microdroplets atop of the solid surface and the wetting transition between the Wenzel and Cassie states. Our simulation results show that the Wenzel-to-Cassie state transition can be achieved by decreasing the fluid-solid attraction, increasing the micropillar spacing, or coating the microstructures with nanostructures. For the effect of the hierarchical structure on the contact angle, we find that the micropillars show a negligible effect on the contact angle, but they may affect the sliding angle. In contrast, it is the nanostructure that determines the contact angle. The contact angle increases with the nanopillar length until reaching a maximal value, but its dependence on the nanopillar spacing becomes more complicated. The contact angle may first increase with the nanopillar spacing and then decreases, or decreases monotonously, depending on whether the liquid enters the nanostructure or not. In this work, we also demonstrate in the presence of contact line pinning, that the pinning effect affects the apparent contact angle.

Nonlinear amplification in electrokinetic pumping in nanochannels in the presence of hydrophobic interactions

We discover a nonlinear coupling between the hydrophobicity of a charged substrate and electrokinetic pumping in narrow fluidic confinements. Our analyses demonstrate that the effective electrokinetic transport in nanochannels may get massively amplified over a regime of bare surface potentials and may subsequently get attenuated beyond a threshold surface charging condition because of a complex interplay between reduced hydrodynamic resistance on account of the spontaneous inception of a less dense interfacial phase and ionic transport within the electrical double layer. We also show that the essential physics delineated by our mesoscopic model, when expressed in terms of a simple mathematical formula, agrees remarkably with that portrayed by molecular dynamics simulations. The nontrivial characteristics of the initial increment followed by a decrement of the effective zeta potential with a bare surface potential may open up the realm of hitherto-unexplored operating regimes of electrohydrodynamically actuated nanofluidic devices.

Thin liquid film flow over substrates with two topographical features

A multicomponent lattice Boltzmann scheme is used to investigate the surface coating of substrates with two topographical features by a gravity-driven thin liquid film. The considered topographies are U- and V-shaped grooves and mounds. For the case of substrates with two grooves, our results indicate that for each of the grooves there is a critical width such that if the groove width is larger than the critical width, the groove can be coated successfully. The critical width of each groove depends on the capillary number, the contact angle, the geometry, and the depth of that groove. The second groove critical width depends on, in addition, the geometry and the depth of the first groove; for two grooves with the same geometries and depths, it is at least equal to that of the first groove. If the second groove width lies between the critical widths, the second groove still can be coated successfully on the condition that the distance between the grooves is considered larger than a critical distance. For considered contact angles and capillary numbers our results indicate that the critical distance is a convex function of the capillary number and the contact angle. Our study also reveals similar results for the case of substrates with a mound and a groove.

Mesoscopic model for microscale hydrodynamics and interfacial phenomena: slip, films, and contact-angle hysteresis

We present a model based on the lattice Boltzmann equation that is suitable for the simulation of dynamic wetting. The model is capable of exhibiting fundamental interfacial phenomena such as weak adsorption of fluid on the solid substrate and the presence of a thin surface film within which a disjoining pressure acts. Dynamics in this surface film, tightly coupled with hydrodynamics in the fluid bulk, determine macroscopic properties of primary interest: the hydrodynamic slip; the equilibrium contact angle; and the static and dynamic hysteresis of the contact angles. The pseudo-potentials employed for fluid-solid interactions are composed of a repulsive core and an attractive tail that can be independently adjusted. This enables effective modification of the functional form of the disjoining pressure so that one can vary the static and dynamic hysteresis on surfaces that exhibit the same equilibrium contact angle. The modeled fluid-solid interface is diffuse, represented by a wall probability function that ultimately controls the momentum exchange between solid and fluid phases. This approach allows us to effectively vary the slip length for a given wettability (i.e., a given static contact angle) of the solid substrate.

High-order thermal lattice Boltzmann models derived by means of Gauss quadrature in the spherical coordinate system

We use the spherical coordinate system in the momentum space and an appropriate discretization procedure to derive a hierarchy of lattice Boltzmann (LB) models with variable temperature. The separation of the integrals in the momentum space into angular and radial parts allows us to compute the moments of the equilibrium distribution function by means of Gauss-Legendre and Gauss-Laguerre quadratures, as well as to find the elements of the discrete momentum set for each LB model in the hierarchy. The capability of the high-order models in this hierarchy to capture specific effects in microfluidics is investigated through a computer simulation of Couette flow by using the Shakhov collision term to get the right value of the Prandtl number.

Macroscopic phase-field model of partial wetting: bubbles in a capillary tube

Drops and bubbles are nonspreading, local, compactly supported features. They are also equilibrium configurations in partial wetting phenomena. Yet, current macroscopic theories of capillary-dominated flow are unable to describe these systems. We propose a framework to model multiphase flow in porous media with nonspreading equilibrium configurations. We illustrate our approach with a one-dimensional model of two-phase flow in a capillary tube. Our model allows for the presence of compactons: nonspreading steady-state solutions in the absence of external forces. We show that local rate dependency is not needed to explain globally rate-dependent displacement patterns, and we interpret dynamic wetting transitions as the route from equilibrium, capillary-dominated configurations, towards viscous-dominated flow. Mathematically, these transitions are possible due to nonclassical shock solutions and the role of bistability and higher-order terms in our model.

Consistent description of electrohydrodynamics in narrow fluidic confinements in the presence of hydrophobic interactions

Electrohydrodynamics in the presence of hydrophobic interactions in narrow confinements is traditionally represented from a continuum viewpoint by a Navier slip-based conceptual paradigm, in which the slip length carries the sole burden of incorporating the effects of substrate wettability on interfacial electromechanics, precluding any explicit dependence of the interfacial potential distribution on the substrate wettability. Here we show that this traditional way of treating electrokinetics-wettability coupling may lead to serious discrepancies in predicting the resultant transport characteristics as manifested through an effective zeta potential. We suggest that an alternative consistent description of the underlying physics through a free-energy-based formalism, in conjunction with considerations of hydrodynamic and electrical property variations consistent with the pertinent phase-field description, may represent the underlying consequences in a more rational manner, as compared to the traditional slip-based model coupled with a two-layer description. Our studies further reveal that the above discrepancies may not occur solely due to the slip-based route of representing the interfacial wettability, but may be additionally attributed to the act of "discretizing" the interfacial phase fraction distribution through an artificial two-layer route.

Wall-mass effects on hydrodynamic boundary slip

This paper investigates the combined effects of surface stiffness κ and wall particles' mass m(w) on the slip length. It aims to enhance our understanding of the momentum and energy transfer across solid-liquid interfaces. Elastic spring potentials are employed to simulate the thermal solid walls and model the surface stiffness κ. The thermal oscillation amplitude is primarily dictated by values of stiffness, whereas the oscillating frequency is proportional to √(κ/m(w)). It is shown that for cases with variable wall mass the relation of slip length and thermal oscillating frequencies can be approximated by a "master" curve according to which the length initially increases, then approaches a peak value, and afterwards is reduced toward an asymptotic value.

Enhancement of surface wettability via the modification of microtextured titanium implant surfaces with polyelectrolytes

Micrometer- and submicrometer-scale surface roughness enhances osteoblast differentiation on titanium (Ti) substrates and increases bone-to-implant contact in vivo. However, the low surface wettability induced by surface roughness can retard initial interactions with the physiological environment. We examined chemical modifications of Ti surfaces [pretreated (PT), R(a) ≤ 0.3 μm; sand blasted/acid etched (SLA), R(a) ≥ 3.0 μm] in order to modify surface hydrophilicity. We designed coating layers of polyelectrolytes that did not alter the surface microstructure but increased surface ionic character, including chitosan (CHI), poly(L-glutamic acid) (PGA), and poly(L-lysine) (PLL). Ti disks were cleaned and sterilized. Surface chemical composition, roughness, wettability, and morphology of surfaces before and after polyelectrolyte coating were examined by X-ray photoelectron spectroscopy (XPS), contact mode profilometry, contact angle measurement, and scanning electron microscopy (SEM). High-resolution XPS spectra data validated the formation of polyelectrolyte layers on top of the Ti surface. The surface coverage of the polyelectrolyte adsorbed on Ti surfaces was evaluated with the pertinent SEM images and XPS peak intensity as a function of polyelectrolyte adsorption time on the Ti surface. PLL was coated in a uniform thin layer on the PT surface. CHI and PGA were coated evenly on PT, albeit in an incomplete monolayer. CHI, PGA, and PLL were coated on the SLA surface with complete coverage. The selected polyelectrolytes enhanced surface wettability without modifying surface roughness. These chemically modified surfaces on implant devices can contribute to the enhancement of osteoblast differentiation.

Molecular transport and flow past hard and soft surfaces: computer simulation of model systems

The equilibrium and flow of polymer films and drops past a surface are characterized by the interface and surface tensions, viscosity, slip length and hydrodynamic boundary position. These parameters of the continuum description are extracted from molecular simulations of coarse-grained models. Hard, corrugated substrates are modelled by a Lennard-Jones solid while polymer brushes are studied as prototypes of soft, deformable surfaces. Four observations are discussed. (i) If the surface becomes strongly attractive or is coated with a brush, the Navier boundary condition fails to describe the effect of the surface independently of the strength and type of the flow. This failure stems from the formation of a boundary layer with an effective, higher viscosity. (ii) In the case of brush-coated surfaces, flow induces a cyclic, tumbling motion of the tethered chain molecules. Their collective motion gives rise to an inversion of the flow in the vicinity of the grafting surfaces and leads to strong, non-Gaussian fluctuations of the molecular orientations. The flow past a polymer brush cannot be described by Brinkman's equation. (iii) The hydrodynamic boundary condition is an important parameter for predicting the motion of polymer droplets on a surface under the influence of an external force. Their steady-state velocity is dictated by a balance between the power that is provided by the external force and the dissipation. If there is slippage at the liquid-solid interface, the friction at the solid-liquid interface and the viscous dissipation of the flow inside the drop will be the dominant dissipation mechanisms; dissipation at the three-phase contact line appears to be less important on a hard surface. (iv) On a soft, deformable substrate like a polymer brush, we observe a lifting-up of the three-phase contact line. Controlling the grafting density and the incompatibility between the brush and the polymer liquid we can independently tune the softness of the surface and the contact angle and thereby identify the parameters for maximizing the deformation at the three-phase contact.

Simulations of slip flow on nanobubble-laden surfaces

On microstructured hydrophobic surfaces, geometrical patterns may lead to the appearance of a superhydrophobic state, where gas bubbles at the surface can have a strong impact on the fluid flow along such surfaces. In particular, they can strongly influence a detected slip at the surface. We present two-phase lattice Boltzmann simulations of a flow over structured surfaces with attached gas bubbles and demonstrate how the detected slip depends on the pattern geometry, the bulk pressure, or the shear rate. Since a large slip leads to reduced friction, our results give assistance in the optimization of microchannel flows for large throughput.

Random-roughness hydrodynamic boundary conditions

We report results of lattice Boltzmann simulations of a high-speed drainage of liquid films squeezed between a smooth sphere and a randomly rough plane. A significant decrease in the hydrodynamic resistance force as compared with that predicted for two smooth surfaces is observed. However, this force reduction does not represent slippage. The computed force is exactly the same as that between equivalent smooth surfaces obeying no-slip boundary conditions, but located at an intermediate position between peaks and valleys of asperities. The shift in hydrodynamic thickness is shown to depend on the height and density of roughness elements. Our results do not support some previous experimental conclusions on a very large and shear-dependent boundary slip for similar systems.

Boundary slip dependency on surface stiffness

The paper investigates the effects of surface stiffness on the slip process aiming to obtain a better insight of the momentum transfer at nanoscale. The surface stiffness is modeled through the stiffness, κ, of spring potentials, which are employed to construct the thermal walls. It is shown that variations of stiffness, κ, influence the slip mechanism either toward slip or stick conditions. Increasing the values of κ alters the oscillation frequency and the mean displacement of the wall particles toward higher and lower values, respectively. Our results suggest that the amount of slip produced as a function of stiffness follows a common pattern that can be modeled through a fifth-order polynomial function.

Fully reversible transition from Wenzel to Cassie-Baxter states on corrugated superhydrophobic surfaces

Liquid drops on textured surfaces show different dynamical behaviors depending on their wetting states. They are extremely mobile when they are supported by composite solid-liquid-air interfaces (Cassie-Baxter state) and immobile when they fully wet the textured surfaces (Wenzel state). By reversibly switching between these two states, it will be possible to achieve control over the fluid dynamics. Unfortunately, these wetting transitions are usually prevented by surface energy barriers. We demonstrate here a new, simple design paradigm consisting of parallel grooves with an appropriate aspect ratio that allows for the controlled, barrierless, reversible switching of the wetting states upon application of electrowetting. We report a direct observation of the barrierless dynamical pathway for the reversible transitions between the Wenzel (collapsed) and Cassie-Baxter (suspended) states and present a theory that accounts for these transitions, including detailed lattice Boltzmann simulations.

Energy transfer through streaming effects in time-periodic pressure-driven nanochannel flows with interfacial slip

We analytically investigate the prospect of using electrokinetic phenomena to transfer hydrostatic energy to electrical power with high energy transfer efficiencies, by exploiting time periodic pressure-driven flows in narrow fluidic confinements. An expression for the energy transfer efficiency for such pulsating pressure-driven flows is derived by considering wall-slip effects due to hydrophobic interactions, strong electrical double layer interactions in the confined flow passages, possibilities of exploring the regimes of large wall potentials, and the adverse consequences of the finite conductance of the Stern layer. It is revealed from our studies that high-frequency pressure pulsations may be employed in practice to improve the concerned energy transfer efficiency to a considerable extent, instead of using a steady-state pressure field. Such favorable effects are found to be best exploited by utilizing "slipping" electro-hydrodynamics in thick electrical double layers in the presence of high surface potentials.

Capillary filling in microchannels with wall corrugations: a comparative study of the Concus-Finn criterion by continuum, kinetic, and atomistic approaches

We study the impact of wall corrugations in microchannels on the process of capillary filling by means of three broadly used methods: computational fluid dynamics (CFD), lattice Boltzmann equations (LBE), and molecular dynamics (MD). The numerical results of these approaches are compared and tested against the Concus-Finn (CF) criterion, which predicts pinning of the contact line at rectangular ridges perpendicular to flow for contact angles of theta > 45 degrees . Whereas for theta = 30, 40 (no flow), and 60 degrees (flow) all methods are found to produce data consistent with the CF criterion, at theta = 50 degrees the numerical experiments provide different results. Whereas the pinning of the liquid front is observed both in the LB and CFD simulations, MD simulations show that molecular fluctuations allow front propagation even above the critical value predicted by the deterministic CF criterion, thereby introducing a sensitivity to the obstacle height.