Sessile nanofluid droplet drying

Nanofluid droplet evaporation has gained much audience nowadays due to its wide applications in painting, coating, surface patterning, particle deposition, etc. This paper reviews the drying progress and deposition formation from the evaporative sessile droplets with the suspended insoluble solutes, especially nanoparticles. The main content covers the evaporation fundamental, the particle self-assembly, and deposition patterns in sessile nanofluid droplet. Both experimental and theoretical studies are presented. The effects of the type, concentration and size of nanoparticles on the spreading and evaporative dynamics are elucidated at first, serving the basis for the understanding of particle motion and deposition process which are introduced afterward. Stressing on particle assembly and production of desirable residue patterns, we express abundant experimental interventions, various types of deposits, and the effects on nanoparticle deposition. The review ends with the introduction of theoretical investigations, including the Navier–Stokes equations in terms of solutions, the Diffusion Limited Aggregation approach, the Kinetic Monte Carlo method, and the Dynamical Density Functional Theory. Nanoparticles have shown great influences in spreading, evaporation rate, evaporation regime, fluid flow and pattern formation of sessile droplets. Under different experimental conditions, various deposition patterns can be formed. The existing theoretical approaches are able to predict fluid dynamics, particle motion and deposition patterns in the particular cases. On the basis of further understanding of the effects of fluid dynamics and particle motion, the desirable patterns can be obtained with appropriate experimental regulations.     

Recent progress in experiments for sessile droplet wetting on structured surfaces

Wetting of a sessile droplet on structured or patterned surface can be found in a broad range of applications. The researchers have been promoted to keep working on the topic. The review is on the basis of the recent experimental advances on the sessile droplet wetting on the hydrophobic, hydrophilic, or combined hydrophobic and hydrophilic surfaces under isothermal conditions, and on heating or cooling substrates having nonisothermal conditions. More attention has been paid on the wetting configuration between the sessile droplet and the structured substrate; the research gap has been discussed on identifying the three-phase line shape. Further, the three-dimensional measurement for the sessile droplets on the patterned surfaces with focusing more on the contact line of sessile droplets might reveal new physical insights. This review targets at building a holistic overview on the sessile droplet wetting behaviors on the structured substrate in the past 2 years.     

Simultaneous spreading and evaporation: Recent developments

The recent progress in theoretical and experimental studies of simultaneous spreading and evaporation of liquid droplets on solid substrates is discussed for pure liquids including nanodroplets, nanosuspensions of inorganic particles (nanofluids) and surfactant solutions. Evaporation of both complete wetting and partial wetting liquids into a nonsaturated vapour atmosphere are considered. However, the main attention is paid to the case of partial wetting when the hysteresis of static contact angle takes place. In the case of complete wetting the spreading/evaporation process proceeds in two stages. A theory was suggested for this case and a good agreement with available experimental data was achieved. In the case of partial wetting the spreading/evaporation of a sessile droplet of pure liquid goes through four subsequent stages: (i) the initial stage, spreading, is relatively short (1–2 min) and therefore evaporation can be neglected during this stage; during the initial stage the contact angle reaches the value of advancing contact angle and the radius of the droplet base reaches its maximum value, (ii) the first stage of evaporation is characterised by the constant value of the radius of the droplet base; the value of the contact angle during the first stage decreases from static advancing to static receding contact angle; (iii) during the second stage of evaporation the contact angle remains constant and equal to its receding value, while the radius of the droplet base decreases; and (iv) at the third stage of evaporation both the contact angle and the radius of the droplet base decrease until the drop completely disappears. It has been shown theoretically and confirmed experimentally that during the first and second stages of evaporation the volume of droplet to power 2/3 decreases linearly with time. The universal dependence of the contact angle during the first stage and of the radius of the droplet base during the second stage on the reduced time has been derived theoretically and confirmed experimentally. The theory developed for pure liquids is applicable also to nanofluids, where a good agreement with the available experimental data has been found. However, in the case of evaporation of surfactant solutions the process deviates from the theoretical predictions for pure liquids at concentration below critical wetting concentration and is in agreement with the theoretical predictions at concentrations above it.     

Liquid nanodroplets spreading on chemically patterned surfaces

Controlling the spatial distribution of liquid droplets on surfaces via surface energy patterning can be used to deliver material to specified regions via selective liquid/solid wetting. Although studies of the equilibrium shape of liquid droplets on heterogeneous substrates exist, much less is known about the corresponding wetting kinetics. Here we present large-scale atomistic simulations of liquid nanodroplets spreading on chemically patterned surfaces. Results are presented for lines of polymer liquid (droplets) on substrates consisting of alternating strips of wetting (equilibrium contact angle theta0 = 0 degrees) and nonwetting (theta0 approximately 90 degrees) material. Droplet spreading is compared for different wavelength lambda of the pattern and strength of surface interaction on the wetting strips. For small lambda, droplets partially spread on both the wetting and nonwetting regions of the substrate to attain a finite contact angle less than 90 degrees. In this case, the extent of spreading depends on the interaction strength in the wetting regions. A transition is observed such that, for large lambda, the droplet spreads only on the wetting region of the substrate by pulling material from nonwetting regions. In most cases, a precursor film spreads on the wetting portion of the substrate at a rate strongly dependent on the width of the wetting region.     

Evaporation of pure liquid sessile and spherical suspended drops: a review

A sessile drop is an isolated drop which has been deposited on a solid substrate where the wetted area is limited by a contact line and characterized by contact angle, contact radius and drop height. Diffusion-controlled evaporation of a sessile drop in an ambient gas is an important topic of interest because it plays a crucial role in many scientific applications such as controlling the deposition of particles on solid surfaces, in ink-jet printing, spraying of pesticides, micro/nano material fabrication, thin film coatings, biochemical assays, drop wise cooling, deposition of DNA/RNA micro-arrays, and manufacture of novel optical and electronic materials in the last decades. This paper presents a review of the published articles for a period of approximately 120 years related to the evaporation of both sessile drops and nearly spherical droplets suspended from thin fibers. After presenting a brief history of the subject, we discuss the basic theory comprising evaporation of micrometer and millimeter sized spherical drops, self cooling on the drop surface and evaporation rate of sessile drops on solids. The effects of drop cooling, resultant lateral evaporative flux and Marangoni flows on evaporation rate are also discussed. This review also has some special topics such as drop evaporation on superhydrophobic surfaces, determination of the receding contact angle from drop evaporation, substrate thermal conductivity effect on drop evaporation and the rate evaporation of water in liquid marbles.     

Evaporation and Marangoni driven convection in small heated water droplets

Evaporation dynamics of small sessile water droplets under microgravity conditions is investigated numerically. The water-air interface is free, and the surrounding air is assumed to be quasisteady. The droplet is described by Navier-Stokes and heat equations and its surrounding water/air gaseous phase with Laplace equation. In the thermodynamic conditions of the simulations presented herein, the evaporative mass flow is nonlinear. It shows a minimum that indicates the existence of qualitative changes in the evaporative regimes although the droplet is sessile. Due to temperature gradients on the free interface, Marangoni motion occurs and generates inside the droplet convection cells that furthermore exhibit small fluctuating motion as evaporation goes on.     

Evaporation of sessile droplets of dilute aqueous solutions containing sodium n-alkylates from polymer surfaces: influences of alkyl length and concentration of solute

The evaporation of sessile droplets placed on polymer surfaces was studied by microscopic observation of the changes in shape of aqueous solution droplets in which the alkyl lengths and the initial concentrations of sodium n-alkylates were varied. Although the initial contact angles of the droplets were not significantly different, the evaporation process varied significantly with the alkyl length of the sodium n-alkylate employed. For the sodium dodecanoate (C 12), showing the highest surface activity, the concentration was found to have a significant effect on the evaporation process of the droplets. In the evaporation of water droplets, variations in the three distinct stages were caused by the different concentration of solutes distributed near or at the air/water interface. It is revealed that the concentration of droplet solute near the air/water interface requires not only solvent evaporation but also some affinity of the solute for the interface. The initial C 12 concentration-dependence of the evaporation of C 12 solution droplets is discussed with particular emphasis on the sudden spreading or sudden contraction of the contact area near the end of evaporation. It is suggested that the cluster formation by C 12 molecules at the air/liquid interface during the evaporation causes Marangoni instability in an evaporating droplet, and the clusters are expected to move dynamically, depending on the droplet concentration of C 12, from the droplet center to the contact line and vice versa, showing Marangoni flow along the air/water interface.     

Convective rolls and hydrothermal waves in evaporating sessile drops

Recent experiments on the evaporation of sessile droplets have revealed the spontaneous formation of various patterns including the presence of hydrothermal waves. These waves had previously been observed, in the absence of evaporation, in thin liquid layers subjected to an imposed, uniform temperature gradient. This is in contrast to the evaporating droplet case wherein these gradients arise naturally due to evaporation and are spatially and temporally varying. In the present paper, we present a theory of evaporating sessile droplets deposited on a heated surface and propose a candidate mechanism for the observed pattern formation using a linear stability analysis in the quasi-steady-state approximation. A qualitative agreement with experimental trends is observed.     

Influence of surface orientation on the organization of nanoparticles in drying nanofluid droplets

The influence of droplet orientation on the flow directed organization of nanoparticles in evaporating nanofluid droplets is important for the efficiency of foliar applied fertilizers and contamination adhesion to the exterior of buildings. The so called "coffee ring" deposit resulting from the evaporation of a sessile nanofluid drop on a hydrophilic surface has received much attention in the literature. Deposits forming on hydrophobic surfaces in the pendant drop position (i.e. hanging drop), which are of importance in foliar fertilizer and exterior building contamination, have received much less attention. In this study, the deposit patterns resulting from the evaporation of water droplets containing silica nanoparticles on hydrophobic surfaces orientated in the sessile or pendant configuration are compared. In the case of a sessile drop the well known coffee ring pattern surrounding a thin nanoparticle layer was formed. A deposit consisting of a thin coffee ring surrounding a bump was formed in the pendant position. A mechanism involving flow induced aggregation at the droplet waist, settling, orientation dependant accumulation within the drop and pinning of the contact line is suggested to explain the findings. Differences in the contact area and adhesion of deposits with surface orientation will affect the efficiency and rainfastness of foliar fertilizers and the cleanliness of building exteriors.     

Behavior of evaporating droplets at nonsoluble and soluble surfaces: Modeling with molecular resolution

Liquid droplets in equilibrium with vapor are simulated at solidlike surfaces using the cooperative motion algorithm (CMA). These droplets behave like real droplets, i.e., the densities of the coexistent liquid and vapor phases obey empirical relations such as rho l - rho v proportional, variant (1 - T/Tc)(1/3). Droplet evaporation was studied under various interaction conditions, i.e., nonsoluble and soluble substrates. In the last case, substrate particles migrate toward the liquid-vapor interface to minimize the droplet surface energy. This leads to the formation of a microwell surrounded by a ringlike deposit on the substrate surface. It is shown that the ring formation in the first stages of evaporation results in pinning of the droplet contact area.     

Facile patterning of upconversion NaYF4:Yb,Er nanoparticles

Different kinds of highly ordered patterns of NaYF(4):Yb,Er nanoparticles on gold substrates were fabricated using a simple method combining micro-contact printing and "breath figures" techniques. Ordered arrays of water droplets were first formed in the hydrophilic regions of patterned self-assembled monolayers (SAMs). This was subsequently submerged in a chloroform solution of NaYF(4):Yb,Er nanoparticles. The particles were spontaneously assembled at the interface of chloroform/water droplet surface, leading to different kinds of uniform patterns after solvent evaporation. The structures of NaYF(4):Yb,Er particles patterns depended on the dimension of the substrate, the concentration of the NaYF(4):Yb,Er nanoparticles and the water condensation process.     

Sessile droplet de-pinning: new life for gravimetric data

Using three different types of surfaces as exemplars, we report a gravimetric method as a viable tool for studying the de-pinning process. Namely, the de-pin time, tau(d) (the time required for a horizontal sessile droplet to de-pin at the triple phase line on a given substrate), is estimated without using a time consuming and expensive video imaging system. This is made possible by deciphering the non-linear portion of mass vs time data of an evaporating sessile droplet. Typical gravimetric glass-substrate evaporative mass loss vs time data has two regimes: a long, linear regime followed by a short, non-linear regime. Traditionally, researchers extract only the evaporation rate of a droplet from the linear regime but discard (by truncating the data) or ignore (thus deriving no information from) the non-linear regime. The origin of the linear to non-linear transition, found almost universally in gravimetric data, persists unremarked upon. By constructing three very different types of surfaces and comparing gravimetric data with video imaging data taken simultaneously, we report the transition is correlated to the onset of the de-pinning event in each case. This realization enables us to measure the de-pin time, tau(d), with gravimetric data only; i.e., without the video system, gathering more information from gravimetric data than previously considered. The method has application in estimating the de-pin time of a droplet deposited on a substrate that yields poor top-view contrast for videography, such as a water droplets on silicon wafers or glass substrates. Finally, gravimetric data is more accurate for evaporation modeling when substrate/droplet interaction areas are not circular.     

Thermal infrared mapping of the freezing phase change activity of micro liquid droplet

This article is dedicated to develop an experimental approach for directly visualizing the global freezing phase change behavior of micro liquid droplets. The infrared (IR) thermograph was proposed to image the basic solidification phenomena of droplet and to acquire its temperature variations during the transient process. In particular, the volumetric recalescence event, regarded as initiation of freezing, was revealed by IR images for the first time. Preliminary results demonstrated that the involved temperature transition due to release of the latent heat can be accurately characterized by evident color break in IR images. Further, experiments were also performed simultaneously on three kinds of droplets made of pure water, dimethylsulfoxide (DMSO) and nano liquid to grasp more precise temporal and spatial temperature distribution. Types of the occurring solidification and the initial frozen volume produced from the recalescence were generally discussed. The IR monitoring method suggests a straightforward way for detecting the freezing phase change activity and its temperature evolution at micro scale.     

Effect of temperature on ring-shaped-deposit formated at evaporation of droplets of silver-nanoparticle dispersions

The effect of temperature is studied on the geometric parameters and conductivity of ring-shaped deposits formed at evaporation of droplets of dispersions of silver nanoparticles on hydrophilic (glass) and hydrophobic (copper) substrates. It has been shown that increasing temperature leads to substantial changes in the deposit profile. Therewith, the effects of temperature on droplet evaporation on glass and copper substrates are different. It has been found that the lateral conductivity of a ring-shaped deposit formed on a glass substrate increases stepwise similarly to a percolation transition at a droplet-evaporation temperature of 58°C. It has been suggested that the reason for the temperature effect is related to a change in the ratio between the rates of physicochemical processes occurring at different stages of droplet evaporation.     

Effect of nonionic surfactant on wetting behavior of an evaporating drop under a reduced pressure environment

The evaporation of sessile drops at reduced pressure is investigated. The evaporation of water droplets on aluminum and PTFE surfaces at reduced pressure was compared. It was found that water droplets on an aluminum surface exhibit a 'depinning jump' at subatmospheric pressures. This is when a pinned droplet suddenly depins, with an increase in contact angle and a simultaneous decrease in the base width. The evaporation of sessile water droplets with a nonionic surfactant (Triton X-100) added to an aluminum surface was then studied. The initial contact angle exhibited a minimum at 0.001 wt% Triton X-100. A maximum in the evaporation rate was also observed at the same concentration. Droplets with low surfactant concentrations are found to exhibit the 'depinning jump.' It is thought that the local concentration of the surfactant causes a gradient of surface tension. The balance at the contact angle is dictated by complex phenomena, including surfactant diffusion and adsorption processes at interfaces. Due to the strong evaporation near the triple line, an accumulation of the surfactant will lead to a surface tension gradient along the interface. The gradient of surface tension will influence the wetting behavior (Marangoni effect). At low surfactant concentrations the contact line depins under the strong effect of surface tension gradient that develops spontaneously over the droplet interface due to surfactant accumulation near the triple line. The maximum evaporation rate corresponds to a minimum contact angle for a pinned droplet.     

Modeling of wetting: a study of nanowetting at rough and heterogeneous surfaces

Molecular dynamics simulations were performed to study the behavior of nanoscale water droplets at solid surfaces. Simulations of droplets on heterogeneous patterned surfaces show that the relative sizes of the domains and the droplets play an important role as do the interactions between the solid and the liquid, particularly when the domain width is comparable to the droplet radius. For pillar surfaces, a transition is observed between the Wenzel and the Cassie and Baxter regimes with increasing pillar height. The effects of pillar width and the gap between the pillars were also examined. The simulations show clearly the importance of the detailed topography and composition of the solid surface.     

Fast, high-throughput creation of size-tunable micro/nanoparticle clusters via evaporative self-assembly in picoliter-scale droplets of particle suspension

We report a fast, high-throughput method to create size-tunable micro/nanoparticle clusters via evaporative assembly in picoliter-scale droplets of particle suspension. Mediated by gravity force and surface tension force of a contacting surface, picoliter-scale droplets of the suspension are generated from a nanofabricated printing head. Rapid evaporative self-assembly of the particles on a hydrophobic surface leads to fast clustering of micro/nanoparticles and forms particle clusters of tunable sizes and controlled spacing. The evaporating behavior of the droplet is observed in real-time, and the clustering characteristics of the particles are understood based on the physics of evaporative-assembly. With this method, multiplex printing of various particle clusters with accurate positioning and alignment are demonstrated. Also, size-unifomity of the cluster arrays is thoroughly analyzed by examining the metallic nanoparticle cluster-arrays based on surface-enhanced Raman spectroscopy (SERS).     

Stick-slip of evaporating droplets: substrate hydrophobicity and nanoparticle concentration

The dynamics of the three-phase contact line for water and ethanol is experimentally investigated using substrates of various hydrophobicities. Different evolutions of the droplet profile (contact line, R, and contact angle, θ) are found to be dependent on the hydrophobicity of the substrate. A simple theoretical approach based on the unbalanced Young force is used to explain the depinning of the contact line on hydrophilic surfaces or the monotonic slip on hydrophobic substrates. The second part of the article involves the addition of different quantities of titanium oxide nanoparticles to water, and a comparison of the evaporative behavior of these novel fluids with the base liquid (water) on substrates varying in hydrophobicity (i.e., silicon, Cytop, and PTFE) is presented. The observed stick-slip behavior is found to be dependent on the nanoparticle concentration. The evaporation rate is closely related to the dynamics of the contact line. These findings may have an important impact when considering the evaporation of droplets on different substrates and/or those containing nanoparticles.     

Onset of Marangoni convection for evaporating sessile droplets

We have generated stability parameters using a linear stability analysis to predict the onset criteria for Marangoni convection in evaporating sessile droplets for two types of substrates, insulating and conducting. The stability problem was formulated with boundary conditions that allow for a temperature discontinuity at the liquid-vapour interface and the inclusion of an expression for the evaporation flux that considers this temperature discontinuity. We introduce no fitting coefficients; therefore, the stability parameters we generate contain only physical variables. The results indicate that spherical sessile droplets evaporating on insulating substrates are predicted to have a similar onset criteria with sessile droplets evaporating on conducting substrates. The onset prediction for sessile droplets evaporating on insulating substrates is found to be considerably different than the case of liquids evaporating from conical funnels constructed of insulating materials owing to the modification of the boundary condition from the geometrical shift and the corresponding retention of modes in the solution. A parametric analysis demonstrates how the input variables impact the stability of evaporating sessile droplets.     

Analysis of droplet evaporation on a superhydrophobic surface

The evaporation process for small, 1-2-mm-diameter droplets of water from patterned polymer surfaces is followed and characterized. The surfaces consist of circular pillars (5-15 microm diameter) of SU-8 photoresist arranged in square lattice patterns such that the center-to-center separation between pillars is 20-30 microm. These types of surface provide superhydrophobic systems with theoretical initial Cassie-Baxter contact angles for water droplets of up to 140-167 degrees, which are significantly larger than can be achieved by smooth hydrophobic surfaces. Experiments show that on these SU-8 textured surfaces water droplets initially evaporate in a pinned contact line mode, before the contact line recedes in a stepwise fashion jumping from pillar to pillar. Provided the droplets of water are deposited without too much pressure from the needle, the initial state appears to correspond to a Cassie-Baxter one with the droplet sitting upon the tops of the pillars. In some cases, but not all, a collapse of the droplet into the pillar structure occurs abruptly. For these collapsed droplets, further evaporation occurs with a completely pinned contact area consistent with a Wenzel-type state. It is shown that a simple quantitative analysis based on the diffusion of water vapor into the surrounding atmosphere can be performed, and estimates of the product of the diffusion coefficient and the concentration difference (saturation minus ambient) are obtained.