Droplet Morphology and Mobility on Lubricant-Impregnated Surfaces: A Molecular Dynamics Study

Study on the Effects of Wettability and Pressure in Shale Matrix Nanopore Imbibition during Shut-in Process by Molecular Dynamics Simulations

Shut-in after fracturing is generally adopted for wells in shale oil reservoirs, and imbibition occurring in matrix nanopores has been proven as an effective way to improve recovery. In this research, a molecular dynamics (MD) simulation was used to investigate the effects of wettability and pressure on nanopore imbibition during shut-in for a typical shale reservoir, Jimsar. The results indicate that the microscopic advancement mechanism of the imbibition front is the competitive adsorption between "interfacial water molecules" at the imbibition front and "adsorbed oil molecules" on the pore wall. The essence of spontaneous imbibition involves the adsorption and aggregation of water molecules onto the hydroxyl groups on the pore wall. The flow characteristics of shale oil suggest that the overall push of the injected water to the oil phase is the main reason for the displacement of adsorbed oil molecules. Thus, shale oil, especially the heavy hydrocarbon component in the adsorbed layer, tends to slip on the walls. However, the weak slip ability of heavy components on the wall surface is an important reason that restricts the displacement efficiency of shale oil during spontaneous imbibition. The effectiveness of spontaneous imbibition is strongly dependent on the hydrophilicity of the matrix pore's wall. The better hydrophilicity of the matrix pore wall facilitates higher levels of adsorption and accumulation of water molecules on the pore wall and requires less time for "interfacial water molecules" to compete with adsorbed oil molecules. During the forced imbibition process, the pressure difference acts on both the bulk oil and the boundary adsorption oil, but mainly on the bulk oil, which leads to the occurrence of wetting hysteresis. Meanwhile, shale oil still existing in the pore always maintains a good, stratified adsorption structure. Because of the wetting hysteresis phenomenon, as the pressure difference increases, the imbibition effect gradually increases, but the actual capillary pressure gradually decreases and there is a loss in the imbibition velocity relative to the theoretical value. Simultaneously, the decline in hydrophilicity further weakens the synergistic effect on the imbibition of the pressure difference because of the more pronounced wetting hysteresis. Thus, selecting an appropriate well pressure enables cost savings and maximizes the utilization of the formation's natural power for enhanced oil recovery (EOR).

Theoretical and Three-Dimensional Molecular Dynamics Study of Droplet Wettability and Mobility on Lubricant-Infused Porous Surfaces

Profiting from their slippery nature, lubricant-infused porous surfaces endow with droplets excellent mobility and consequently promise remarkable heat transfer improvement for dropwise condensation. To be a four-phase wetting system, the droplet wettability configurations and the corresponding dynamic characteristics on lubricant-infused porous surfaces are closely related to many factors, such as multiple interfacial interactions, surface features, and lubricant thickness, which keeps a long-standing challenge to promulgate the underlying physics. In this work, thermodynamically theoretical analysis and three-dimensional molecular dynamics simulations with the coarse-grained water and hexane models are carried out to explore droplet wettability and mobility on lubricant-infused porous surfaces. Combined with accessible theoretical criteria, phase diagrams of droplet configurations are constructed with a comprehensive consideration of interfacial interactions, surface structures, and lubricant thickness. Subsequently, droplet sliding and coalescence dynamics are quantitatively defined under different configurations. Finally, in terms of the promotion of dropwise condensation, a non-cloaking configuration with the encapsulated state underneath the droplet is recommended to achieve high droplet mobility owing to the low viscous drag of the lubricant and the eliminated pinning effect of the contact line. On the basis of the low oil-water and water-solid interactions, a stable lubricant layer with a relatively low thickness is suggested to construct slippery surfaces.

Interdependence of Surface Roughness on Icephobic Performance: A Review

Ice protection techniques have attracted significant interest, notably in aerospace and wind energy applications. However, the current solutions are mostly costly and inconvenient due to energy-intensive and environmental concerns. One of the appealing strategies is the use of passive icephobicity, in the form of coatings, which is induced by means of several material strategies, such as hydrophobicity, surface texturing, surface elasticity, and the physical infusion of ice-depressing liquids, etc. In this review, surface-roughness-related icephobicity is critically discussed to understand the challenges and the role of roughness, especially on superhydrophobic surfaces. Surface roughness as an intrinsic, independent surface property for anti-icing and de-icing performance is also debated, and their interdependence is explained using the related physical mechanisms and thermodynamics of ice nucleation. Furthermore, the role of surface roughness in the case of elastomeric or low-modulus polymeric coatings, which typically instigate an easy release of ice, is examined. In addition to material-centric approaches, the influence of surface roughness in de-icing evaluation is also explored, and a comparative assessment is conducted to understand the testing sensitivity to various surface characteristics. This review exemplifies that surface roughness plays a crucial role in incorporating and maintaining icephobic performance and is intrinsically interlinked with other surface-induced icephobicity strategies, including superhydrophobicity and elastomeric surfaces. Furthermore, the de-icing evaluation methods also appear to be roughness sensitive in a certain range, indicating a dominant role of mechanically interlocked ice.

Directional Drop Rebound on Adhesive-Gradient Micro-Nanostructured Surfaces Formed by a Femtosecond Laser

The dynamic behavior of droplets hitting a solid surface has received extensive attention due to its broad application prospects. Additionally, controlling the rebound behavior of impacting droplets is an important research topic. Current methods for investigating this behavior focus on the construction of a differentiated wettability surface, which is characterized by contact angle measurements, or a differentiated topography surface, which is represented by geometric height. This information allows one to obtain the nonuniform kinetic energy distribution of rebounding droplets and to realize control of rebounding droplet behavior. In this paper, femtosecond laser processing is proposed for the fabrication of an anisotropic surface with differences in adhesion, which allows for the control of impacting droplet rebound behavior. The experimental results show that the micro-nanostructure of the surface affects its adhesion. By changing the micro-nanostructure of the solid surface, the difference in surface adhesion can be controlled, thereby realizing precise control of impacting droplet rebound behavior. This study demonstrates that the micro-nanostructured surface formed by a femtosecond laser can be used to control a droplet rebound direction and landing site, which is of great significance to the development of liquid transport, microfluidic devices, and other fields.

Programmable Transport of C60 by Straining Graphene Substrate

Controlling the maneuverability of nanocars and molecular machines on the surface is essential for the targeted transportation of materials and energy at the nanoscale. Here, we evaluate the motion of fullerene, as the most popular candidate for use as a nanocar wheel, on the graphene nanoribbons with strain gradients based on molecular dynamics (MD), and theoretical approaches. The strain of the examined substrates linearly decreases by 20%, 16%, 12%, 8%, 4%, and 2%. MD calculations were performed with the open source LAMMPS solver. The essential physics of the interactions is captured by Lennard-Jones and Tersoff potentials. The motion of C60 on the graphene nanoribbon is simulated in canonical ensemble, which is implanted by using a Nose-Hoover thermostat. Since the potential energy of C60 is lower on the unstrained end of nanoribbons, this region is energetically more favorable for the molecule. As the strain gradient of the surface increases, the trajectories of the motion and the C60 velocity indicate more directed movements along the gradient of strain on the substrate. Based on the theoretical relations, it was shown that the driving force and diffusion coefficient of the C60 motion respectively find linear and quadratic growth with the increase of strain gradient, which is confirmed by MD simulations. To understand the effect of temperature, at each strain gradient of substrate, the simulations are repeated at the temperatures of 100, 200, 300, and 400 K. The large ratio of longitudinal speed to the transverse speed of fullerene at 100 and 200 K refers to the rectilinear motion of molecule at low temperatures. Using successive strain gradients on the graphene in perpendicular directions, we steered the motion of C60 to the desired target locations. The programmable transportation of nanomaterials on the surface has a significant role in different processes at the nanoscale, such as bottom-up assembly.

The Effect of Liquid-Solid Interactions upon Nucleate Boiling on Rough Surfaces: Insights from Molecular Dynamics

Despite the fact that engineered surface enabling remarkable phase change heat transfer have elicited increasing attention due to their ubiquitous applications in thermal management, the underlying mechanisms of intrinsic rough structures as well as the surface wettability on bubble dynamics remain to be explored. Therefore, a modified molecular dynamics simulation of nanoscale boiling was conducted in the present work to investigate bubble nucleation on rough nanostructured substrates with different liquid-solid interactions. Specifically, the initial stage of nucleate boiling was mainly investigated and the bubble dynamic behaviors were quantitively studied under different energy coefficients. Results shows that as the contact angle decreases, the nucleation rate increases, because liquid obtains more thermal energy there compared with that on less wetting surfaces. The rough profiles of the substrate can provide nanogrooves, which can enhance initial nucleate embryos, thereby improving thermal energy transfer efficiency. Moreover, atomic energies are calculated and adopted to explain how bubble nuclei are formed on various wetting substrates. The simulation results are expected to provide guidance towards surface design in state-of-the art thermal management systems, such as the surface wettability and the nanoscale surface patterns.

Cloaking Transition of Droplets on Lubricated Brushes

We study the equilibrium properties and the wetting behavior of a simple liquid on a polymer brush, with and without the presence of lubricant by multibody Dissipative Particle Dynamics simulations. The lubricant is modeled as a polymeric liquid consisting of short chains that are chemically identical with the brush polymers. We investigate the behavior of the brush in terms of the grafting density and the amount of lubricant present. Regarding the wetting behavior, we study a sessile droplet on top of the brush. The droplet consists of nonbonded particles that form a dense phase. Our model and choice of parameters result in the formation of a wetting ridge and in the cloaking of the droplet by the lubricant; i.e., the lubricant chains creep up onto the droplet and eventually cover its surface completely. Cloaking is a phenomenon that is observed experimentally and is of integral importance to the dynamics of sliding droplets. We quantify the cloaking in terms of its thickness, which increases with the amount of lubricant present. The analysis reveals a well-defined transition point where the cloaking sets in. We propose a thermodynamic theory to explain this behavior. In addition, we investigate the dependence of the contact angles on the size of the droplet and the possible effect of line tension. We quantify the variation of the contact angle with the curvature of the contact line on a lubricant free brush and find a negative value for the line tension. Finally we investigate the effect of cloaking/lubrication on the contact angles and the wetting ridge. We find that lubrication and cloaking reduce the contact angles by a couple of degrees. The effect on the wetting ridge is a reduction in the extension of the brush chains near the three phase contact line, an effect that was also observed in experiments of droplets on cross-linked gels.

Role of the microridges on cactus spines

Cactus spines have inspired a wide range of micro- and nano-structures that cause droplets to move spontaneously and directionally. The conical shape and the surface wettability gradient are two typical characteristics in such systems. The cross section of the existing conical fibers is usually assumed to be an ideal circle. In fact, microridges are observed on the spine surface of the cactus, and the function is not yet fully understood. The present work thus focuses on how microridges affect droplet self-transport. Structures mimicking microridges are first investigated by constructing pyramidal cross sections with concave or convex lateral faces. The dissipative particle dynamics method is then employed to numerically investigate and theoretically analyze the dynamic behaviors of droplets on these conical fibers with different cross sections. The results show that the microridges reduce the base radius and the contact area of the droplet, thereby increasing the driving force and reducing the friction force. Moreover, by mimicking the microridges structure, we propose a conical fiber with a triple concave cross section, which increases the droplet velocity and the distance traveled over the traditional circular fiber. This work reveals the role of the microridges in the droplet self-transport, which opens up new prospects for the manufacture of fiber systems for microfluidics and liquid manipulation.

Contact angle of Nepenthes slippery zone: results from measurement and model analysis

The slippery zone of Nepenthes alata depends on its highly evolved morphology and structure to show remarkable superhydrophobicity, which has gradually become a biomimetic prototype for developing superhydrophobic materials. However, the mechanism governing this phenomenon has not been fully revealed through a model analysis. In this paper, the superhydrophobicity of the slippery zone is studied by contact angle measurement, morphology/structure examination and model analysis. The slippery zone causes an ultrapure water droplet to produce a considerably high contact angle (155.11–158.30°) and has micro–nanoscale hierarchical structures consisting of lunate cells and wax coverings. According to the Cassie–Baxter equation and a self-defined infiltration coefficient, a model was established to analyse the effect of a structure characteristic on the contact angle. The analysis, result showed that the calculated contact angle (154.67–159.49°) was highly consistent with the measured contact angle, indicating that the established model can quantitatively characterise the relationship between the contact angle and the structure characteristic. The authors’ study provides some evidences to further reveal the superhydrophobic mechanism of the slippery zone of N. alata, as well as inspiring the biomimetic development of superhydrophobic surfaces.

Head-on Collision of Two Nanodroplets on a Solid Surface: A Molecular Dynamics Simulation Study

Most researchers focus on the collision of a single droplet with a solid surface, while it is common for a droplet to collide with a sessile droplet on a solid surface in reality. This study performed the head-on collision of two nanodroplets on a solid surface using the molecular dynamics simulation method. The effects of impact velocity, interaction intensity between solid and liquid atoms, and the solid fraction of the surface on the collision process are studied with independent simulation cases. The maximum spreading factor and the dimensionless maximum spreading time are recorded and calculated to describe the collision process quantitatively. The simulation results indicate that the maximum spreading factor depends more on the solid fraction than the interaction intensity since it does not fundamentally change the wetting state of the droplet at its maximum spreading state. Because of two different effects, the maximum dimensionless spreading time decreases first and then increases with the interaction intensity, and both effects weaken with the increase of impact velocity. As the solid fraction increases, the maximum spreading factor increases significantly at high impact velocity, and the maximum dimensionless spreading time first decreases and then increases because the wetting state of the coalescent droplet at the maximum spreading moment gradually changes from the Wenzel state to the Cassie state. In general, the initial wetting state of the sessile droplet and the wetting state of the coalescent droplet at the maximum spreading moment have important effects on the maximum spreading factor and the maximum spreading time. We establish a theoretical prediction model for the maximum spreading factor on a smooth surface based on energy conservation with quite good accuracy. This research has improved our understanding of the head-on collision process of two nanodroplets on a solid surface.

Dropwise condensation: From fundamentals of wetting, nucleation, and droplet mobility to performance improvement by advanced functional surfaces

As a ubiquitous vapor-liquid phase-change process, dropwise condensation has attracted tremendous research attention owing to its remarkable efficiency of energy transfer and transformative industrial potential. In recent years, advanced functional surfaces, profiting from great progress in modifying micro/nanoscale features and surface chemistry on surfaces, have led to exciting advances in both heat transfer enhancement and fundamental understanding of dropwise condensation. In this review, we discuss the development of some key components for achieving performance improvement of dropwise condensation, including surface wettability, nucleation, droplet mobility, and growth, and discuss how they can be elaborately controlled as desired using surface design. We also present an overview of dropwise condensation heat transfer enhancement on advanced functional surfaces along with the underlying mechanisms, such as jumping condensation on nanostructured superhydrophobic surfaces, and new condensation characteristics (e.g., Laplace pressure-driven droplet motion, hierarchical condensation, and sucking flow condensation) on hierarchically structured surfaces. Finally, the durability, cost, and scalability of specific functional surfaces are focused on for future industrial applications. The existing challenges, alternative strategies, as well as future perspectives, are essential in the fundamental and applied aspects for the practical implementation of dropwise condensation.

Anisotropic Ice Adhesion of Micro-Nano-Structured Metal Surface by a Femtosecond Laser

The icephobic materials induced using micro-nano-structured surfaces have aroused great attention for promising applications. Previously, the characterization of ice adhesion of icephobic materials by shear force is usually performed without direction discrimination along the surface whatever the surface is anisotropic or not. In this work, we studied the direction-dependent ice adhesion strength on groove-shaped micro-nano-structured aluminum alloy surfaces formed using a femtosecond laser. It is found that the ice adhesion strength on the surfaces exhibits anisotropy, which corresponds to a smaller ice adhesion strength in the direction parallel to the groove than that orthogonal to the groove. Furthermore, it is found that the ice adhesion strength decreases with the increase in groove width in the orthogonal direction, while it does not change much in the parallel direction. The anisotropic ice adhesion strength is attributed to the change of wettability and morphology in the two directions. The findings in this work suggest that anisotropic ice adhesion should be fully considered when designing an icephobic micro-nano-structured metal structure, which is of great significance to the characterization and application of icephobic materials.

Droplet Mobility on Slippery Lubricant Impregnated and Superhydrophobic Surfaces under the Effect of Air Shear Flow

The focus of this study is to investigate and compare the behavior of a droplet on superhydrophobic (SHS) and slippery lubricant impregnated (SLIPS) surfaces under the effect of air shear flow. In this regard, both experimental and numerical analyses have been conducted to compare their performance on droplet mobility under different air speeds. Two different lubricants have been utilized to scrutinize their effect on droplet movement. The numerical simulations have been performed based on the volume of fluid method coupled with the large eddy simulation turbulent model in conjunction with the dynamic contact angle method in addition to a model that can represent the effect of lubricants on slippery surfaces. The numerical simulations are compared with the experimental study in order to shed light on the underlying mechanisms. The results showed that under the same conditions, the critical velocity for droplet movement on the superhydrophobic surfaces is lower than that on the slippery lubricant impregnated surfaces due to the smaller droplet base diameter and the larger contact angle. The hydrodynamics of droplet mobility on superhydrophobic surfaces exhibits a rolling behavior while for the slippery lubricant impregnated surfaces a combination of rolling and sliding is observed. Beyond the critical airflow speed, a complete droplet shedding on all surfaces occurs. The wetting length and position of the droplet on superhydrophobic and slippery surfaces have been measured. On slippery surfaces, the speed of droplets is greatly affected by the lubricant properties while similar behavior in the wetting lengths is observed.

Failure and Recovery of Droplet Nucleation and Growth on Damaged Nanostructures: A Molecular Dynamics Study

The condensate flooding during dropwise condensation causes serious deterioration in heat transfer performance. In this study, the three-dimensional large-scale molecular dynamics simulation is carried out to investigate the droplet state transition from local flooding mode to Wenzel or from Wenzel to Cassie due to the droplet coalescence under the effect of nanostructure size. In particular, the effect of nanostructure breakage on droplet nucleation and growth is discussed to reveal the mechanism of dropwise condensation heat transfer deterioration. As a potential solution, the lubricant-impregnated surface is proposed to recover the preferred Cassie state by regulating the dynamic wetting characteristics of droplets, and thus the detrimental effect of nanostructure breakage could be effectively avoided.

Heterogeneous nucleation of argon vapor on the nanostructure surface with molecular dynamics simulation

The study of vapor condensation on surface has important engineering significance. The condensation nucleation process of vapor on substrates with different pillar structures is studied through molecular dynamic simulation. The condensation droplets on the pillar with various heights and solid fractions exhibit Wenzel state, Cassie state and the transformation from Wenzel state to Cassie state. The results show that the condensation efficiencies are correlated with the state of droplet and it is explained from the perspective of heat transfer. For the Wenzel state, the droplet fills the gap of pillar and the form of the heat conduction change with the growth of cluster. In the initial of condensation droplets, the heat conduction is similar with various heights of pillar. As the condensation droplets grow, the efficiency of heat conduction enhances with the increasing of height of pillar. For the Cassie state, the form of heat conduction is perpetual during the condensation process with the thermal resistance of the droplet dominated due to the droplets suspended on the pillar. The efficiency of heat transfer is insensitive to the height of pillar. The form of heat conduction for the transformation state transforms from Wenzel state to Cassie state leading to the reduction of condensation rate. The droplet formed in the Wenzel state has the higher transfer efficiency than the Cassie one.

Dynamic Wettability on the Lubricant-Impregnated Surface: From Nucleation to Growth and Coalescence

The surface dynamic wettability during droplet nucleation and growth involved with phase change is different from the static wettability formed from a sessile drop. Revealing this dynamic wettability of the lubricant-impregnated surfaces (LISs) and identification of the consistency between the wettability during condensation and the static wettability are of significant importance. In this study, we investigated condensation of water droplets on LISs using molecular dynamics simulations. All possible morphologies on LISs were investigated considering the effects of interfacial tension and lubricant thickness. The exploration of droplet behaviors from nucleation to growth and coalescence revealed four nucleation mechanisms and six growth modes. The lubricant was observed to be beneficial for the formation of droplets and maintaining dropwise condensation mode. The present investigation also established that the consistency between the wettability during condensation and the static wettability was determined by the solid-water-oil interface and the lubricant thickness. A map was proposed which helps in deciding whether the wettability during condensation is the same as the static wettability on LIS.

Inner surface of Nepenthes slippery zone: ratchet effect of lunate cells causes anisotropic superhydrophobicity

Inner surface of Nepenthes slippery zone shows anisotropic superhydrophobic wettability. Here, we investigate what factors cause the anisotropy via sliding angle measurement, morphology/structure observation and model analysis. Static contact angle of ultrapure-water droplet exhibits the value of 154.80°-156.83°, and sliding angle towards pitcher bottom and up is 2.82 ± 0.45° and 5.22 ± 0.28°, respectively. The slippery zone under investigation is covered by plenty of lunate cells with both ends bending downward, and a dense layer of wax coverings without directional difference in morphology/structure. Results indicate that the slippery zone has a considerable anisotropy in superhydrophobic wettability that is most likely caused by the lunate cells. A model was proposed to quantitatively analyse how the structure characteristics of lunate cells affect the anisotropic superhydrophobicity, and found that the slope/precipice structure of lunate cells forms a ratchet effect to cause ultrapure-water droplet to roll towards pitcher bottom/up in different order of difficulty. Our investigation firstly reveals the mechanism of anisotropic superhydrophobic wettability of Nepenthes slippery zone, and inspires the bionic design of superhydrophobic surfaces with anisotropic properties.

Life and death of liquid-infused surfaces: a review on the choice, analysis and fate of the infused liquid layer

We review the rational choice, the analysis, the depletion and the properties imparted by the liquid layer in liquid-infused surfaces – a new class of low-adhesion surface.