Pinning, retraction, and terracing of evaporating droplets containing nanoparticles

MPs and PFOS single and combined exposure significantly alter genetic expressions of growth hormone and insulin growth factor-related biomarkers during zebrafish embryonic development

Microplastics (MPs) and perfluorooctane sulfonate (PFOS) are emerging pollutants that are ubiquitously present in the environment and can cause series of ecotoxicological effects on aquatic animals. This study examined how the expression of genes related to insulin growth factor (igf1, igf2a, igf2b, igfra, and igfrb) and growth hormone (ghrh, gh1, ghra, and ghrb) changes during the development of zebrafish embryos exposed to 8 μm polyethylene microplastics (PE-MPs) and perfluorooctane sulfonate (PFOS) individually and in combination for 72 h. Our findings revealed that both low-concentrations of MP (50 μg/L) and PFOS (0.02 μg/L) treatments could significantly activate gene expression within a short period. High concentrations of MPs (500 μg/L) and PFOS (0.1 μg/L) not only rapidly activated gene expression but also sustained high expression levels for a longer duration. During combined exposures, peak gene expression in the low concentration groups (50 μg/L MPs and 0.02 μg/L PFOS; 50 μg/L MPs and 0.1 μg/L PFOS) primarily occurred within 12 h after treatment. In the high concentration groups (500 μg/L MPs and 0.02 μg/L PFOS), peak expression was also observed within 12 h. Notably, the combined exposure groups exhibited more pronounced effects on gene expression than the individual exposure groups. The activation of gene expression was both more significant and longer-lasting in the combined exposure, indicating a synergistic regulatory effect of MPs and PFOS. Overall, our study suggests that zebrafish embryo development can be significantly impacted by exposure to MPs, PFOS, and their combination, with combined exposures having a more lasting and profound effect on gene regulation compared to single exposures.

Droplet evaporation on soft solid substrates

Droplet evaporation on soft solid substrates is relevant to applications such as fabrication of microlenses and controlled particle deposition. Here, we develop a lubrication-theory-based model to advance fundamental understanding of the important limiting case of a planar droplet evaporating on a linear viscoelastic solid. A set of partial differential equations describing the time evolution of the liquid-air and liquid-solid interfaces is derived and solved with a finite-difference method. A disjoining-pressure/precursor-film approach is used to describe contact-line motion, and the one sided model is used to describe solvent evaporation. Parametric studies are conducted to investigate the effect of solid properties (thickness, viscosity, shear modulus, wettability) and evaporation rate on droplet dynamics. Our results indicate that softer substrates speed up droplet evaporation due to prolonged pinning of the contact line. Results from our model are able to qualitatively reproduce some key trends observed in experiments. Due to its systematic formulation, our model can readily be extended to more complex situations of interest such as evaporation of particle-laden droplets on soft solid substrates.

Evaporation of Binary-Mixture Liquid Droplets: The Formation of Picoliter Pancakelike Shapes

Small multicomponent droplets are of increasing importance in a plethora of technological applications ranging from the fabrication of self-assembled hierarchical patterns to the design of autonomous fluidic systems. While often far away from equilibrium, involving complex and even chaotic flow fields, it is commonly assumed that in these systems with small drops surface tension keeps the shapes spherical. Here, studying picoliter volatile binary-mixture droplets of isopropanol and 2-butanol, we show that the dominance of surface tension forces at small scales can play a dual role: Minute variations in surface tension along the interface can create Marangoni flows that are strong enough to significantly deform the drop, forming micron-thick pancakelike shapes that are otherwise typical of large puddles. We identify the conditions under which these flattened shapes form and explain why, universally, they relax back to a spherical-cap shape toward the end of drop lifetime. We further show that the formation of pancakelike droplets suppresses the "coffee-ring" effect and leads to uniform deposition of suspended particles. The quantitative agreement between theory and experiment provides a predictive capability to modulate the shape of tiny droplets with implications in a range of technologies from fabrication of miniature optical lenses to coating, printing, and pattern deposition.

Evaporation of Sessile Droplets of Polyelectrolyte/Surfactant Mixtures on Silicon Wafers

The wetting and evaporation behavior of droplets of aqueous solutions of mixtures of poly(diallyldimethylammonium chloride) solution, PDADMAC, with two different anionic surfactants, sodium laureth sulfate, SLES, and sodium N-lauroyl N-methyl taurate, SLMT, were studied in terms of the changes of the contact angle θ and contact length L of sessile droplets of the mixtures on silicon wafers at a temperature of 25 °C and different relative humidities in the range of 30–90%. The advancing contact angle θa was found to depend on the surfactant concentration, independent of the relative humidity, with the mixtures containing SLES presenting improved wetting behaviors. Furthermore, a constant droplet contact angle was not observed during evaporation due to pinning of the droplet at the coffee-ring that was formed. The kinetics for the first evaporation stage of the mixture were independent of the relative humidity, with the evaporation behavior being well described in terms of the universal law for evaporation.

Drop Behavior Influenced by the Correlation Length on Noisy Surfaces

We investigate numerically the role of the correlation length in drop behavior on noisy surfaces. To this aim, a phase field tool has been used. Theoretical results are confirmed by experiments of distilled water drops sitting on stainless steel and silicon surfaces textured by laser-induced periodic self-organized structures: an increase of the noise amplitude results in an amplification of the original behavior (i.e., hydrophobic is getting more hydrophobic, hydrophilic is getting more hydrophilic). Furthermore, computer simulations in two and three spatial dimensions allow for predictions of drop behavior on noisy sloped substrates under a gravitational force, a problem of large interest in controlled motion in micro- and nanofluidics.

A review on suppression and utilization of the coffee-ring effect

Evaporation of sessile droplets containing non-volatile solutes dispersed in a volatile solvent leaves behind ring-like solid stains. As the volatile species evaporates, pinning of the contact line gives rise to capillary flows that transport non-volatile solutes to the contact line. This phenomenon, called the coffee-ring effect, compromises the overall performance of industrially relevant manufacturing processes involving evaporation such as printing, biochemical analysis, manufacturing of nano-structured materials through colloidal and macromolecular patterning. Various approaches have been developed to suppress this phenomenon, which is otherwise difficult to avoid. The coffee-ring effect has also been leveraged to prepare new materials through convection induced assembly. This review underlines not only the strategies developed to suppress the coffee-ring effect but also sheds light on approaches to arrive at novel processes and materials. Working principles and applicability of these strategies are discussed together with a critical comparison.

Multiscale assembly of solution-processed organic electronics: the critical roles of confinement, fluid flow, and interfaces

Organic semiconducting small molecules and polymers provide a rich phase space for investigating the fundamentals of molecular and hierarchical assembly. Stemming from weak intermolecular interactions, their assembly sensitively depends on processing conditions, which in turn drastically modulate their electronic properties. Much work has gone into molecular design strategies that maximize intermolecular interactions and encourage close packing. Less understood, however, is the non-equilibrium assembly that occurs during the fabrication process (especially solution coating and printing) which is critical to determining thin film morphology across length scales. This encompasses polymorphism and molecular packing at molecular scale, assembly of π-bonding aggregates at the tens of nanometers scale, and the formation of domains at the micron-millimeter device scale. Here, we discuss three phenomena ubiquitous in solution processing of organic electronic thin films: the confinement effect, fluid flows, and interfacial assembly and the role they play in directing assembly. This review focuses on the mechanistic understanding of how assembly outcomes couple closely to the solution processing environment, supported by salient examples from the recent literature.

Drying of Droplets of Colloidal Suspensions on Rough Substrates

In many technological applications, excess solvent must be removed from liquid droplets to deposit solutes onto substrates. Often, the substrates on which the droplets rest may possess some roughness, either intended or unintended. Motivated by these observations, we present a lubrication-theory-based model to study the drying of droplets of colloidal suspensions on a substrate containing a topographical defect. The model consists of a system of one-dimensional partial differential equations accounting for the shape of the droplet and depth-averaged concentration of colloidal particles. A precursor film and disjoining pressure are used to describe the contact-line region, and evaporation is included using the well-known one-sided model. Finite-difference solutions reveal that when colloidal particles are absent, the droplet contact line can pin to a defect for a significant portion of the drying time due to a balance between capillary-pressure gradients and disjoining-pressure gradients. The time-evolution of the droplet radius and contact angle exhibits the constant-radius and constant-contact-angle stages that have been observed in prior experiments. When colloidal particles are present and the defect is absent, the model predicts that particles will be deposited near the center of the droplet in a cone-like pattern. However, when a defect is present, pinning of the contact-line accelerates droplet solidification, leading to particle deposition near the droplet edge in a coffee-ring pattern. These predictions are consistent with prior experimental observations, and illustrate the critical role contact-line pinning plays in controlling the dynamics of drying droplets.

Evaporation of Sessile Droplets Laden with Particles and Insoluble Surfactants

We consider the flow dynamics of a thin evaporating droplet in the presence of an insoluble surfactant and noninteracting particles in the bulk. On the basis of lubrication theory, we derive a set of evolution equations for the film height, the interfacial surfactant, and bulk particle concentrations, taking into account the dependence of liquid viscosity on the local particle concentration. An important ingredient of our model is that it takes into account the fact that the surfactant adsorbed at the interface hinders evaporation. We perform a parametric study to investigate how the presence of surfactants affects the evaporation process as well as the flow dynamics with and without the presence of particles in the bulk. Our numerical calculations show that the droplet lifetime is affected significantly by the balance between the ability of the surfactant to enhance spreading, suppressing the effect of thermal Marangoni stresses-induced motion, and to hinder the evaporation flux through the reduction of the effective interfacial area of evaporation, which tend to accelerate and decelerate the evaporation process, respectively. For particle-laden droplets and in the case of dilute solutions, the droplet lifetime is found to be weakly dependent on the initial particle concentration. We also show that the particle deposition patterns are influenced strongly by the direct effect of the surfactant on the evaporative flux; in certain cases, the "coffee-stain" effect is enhanced significantly. A discussion of the delicate interplay between the effects of capillary pressure and solutal and thermal Marangoni stresses, which drive the liquid flow inside of the evaporating droplet giving rise to the observed results, is provided herein.

Numerical simulation of dip-coating in the evaporative regime

A hydrodynamic model is used for numerical simulations of a polymer solution in a dip-coating-like experiment. We focus on the regime of small capillary numbers where the liquid flow is driven by evaporation, in contrast to the well-known Landau-Levich regime dominated by viscous forces. Lubrication approximation is used to describe the flow in the liquid phase. Evaporation in stagnant air is considered (diffusion-limited evaporation), which results in a coupling between liquid and gas phases. Self-patterning due to the solutal Marangoni effect is observed for some ranges of the control parameters. We first investigate the effect of evaporation rate on the deposit morphology. Then the role of the spatial variations in the evaporative flux on the wavelength and mean thickness of the dried deposit is ascertained, by comparing the 2D and 1D diffusion models for the gas phase. Finally, for the very low substrate velocities, we discuss the relative importance of diffusive and advective components of the polymer flux, and consequences on the choice of the boundary conditions.

Modification of the Young-Laplace equation and prediction of bubble interface in the presence of nanoparticles

Bubbles are fundamental to our daily life and have wide applications such as in the chemical and petrochemical industry, pharmaceutical engineering, mineral processing and colloids engineering. This paper reviews the existing theoretical and experimental bubble studies, with a special focus on the dynamics of triple line and the influence of nanoparticles on the bubble growth and departure process. Nanoparticles are found to influence significantly the effective interfacial properties and the dynamics of triple line, whose effects are dependent on the particle morphology and their interaction with the substrate. While the Young-Laplace equation is widely applied to predict the bubble shape, its application is limited under highly non-equilibrium conditions. Using gold nanoparticle as an example, new experimental study is conducted to reveal the particle concentration influence on the behaviour of triple line and bubble dynamics. A new method is developed to predict the bubble shape when the interfacial equilibrium conditions cannot be met, such as during the oscillation period. The method is used to calculate the pressure difference between the gas and liquid phases, which is shown to oscillate across the liquid-gas interface and is responsible for the interface fluctuation. The comparison of the theoretical study with the experimental data shows a very good agreement, which suggests its potential application to predict bubble shape during non-equilibrium conditions.

Effect of particle geometry on triple line motion of nano-fluid drops and deposit nano-structuring

We illustrate the importance of particle geometry on droplet contact line pinning, 'coffee-stain' formation and nano-structuring within the resulting rings. We present the fundamentals of pure liquid droplet evaporation and then discuss the effect of particles on the evaporation process. The resulting coffee-stain patterns and particle structuring within them are presented and discussed. In the second part, we turn our attention to the effect of particle geometry on the evaporation process. A wide range of particle shapes, categorised according to aspect ratio, from the simple shape of a sphere to the highly irregular shapes of platelets and tubes is discussed. Particle geometry effect on evaporation behaviour was quantified in terms of change in contact angle and contact radius for the stick-slip cases. Consequently the hysteretic energy barrier pinning the droplets was estimated, showing an increasing trend with particle aspect ratio. The three-phase contact line (TL) motion kinetics are complemented with analysis of the nano-structuring behaviour of each shape, leading to the identification of the two main parameters affecting nanoparticle self-assembly behaviour at the wedge. Flow velocity and wedge constraints were found to have antagonist effects on particle deposition, although these varied with particle shape. This description should help in understanding the drying behaviour of more complex fluids. Furthermore, knowing the fundamentals of this simple and inexpensive surface patterning technique should permit its tailoring to the needs of many potential applications.

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.

Sagging of evaporating droplets of colloidal suspensions on inclined substrates

A droplet of a colloidal suspension placed on an inclined substrate may sag under the action of gravity. Solvent evaporation raises the concentration of the colloidal particles, and the resulting viscosity changes may influence the sag of the droplet. To investigate this phenomenon, we have developed a mathematical model for perfectly wetting droplets based on lubrication theory and the rapid-vertical-diffusion approximation. Precursor films are assumed to be present, the colloidal particles are taken to be hard spheres, and particle and liquid dynamics are coupled through a concentration-dependent viscosity and diffusivity. Evaporation is assumed to be limited by how rapidly solvent molecules can transfer from the liquid to the vapor phase. The resulting one-dimensional system of nonlinear partial differential equations describing the evolution of the droplet height and particle concentration is solved numerically for a range of initial particle concentrations and substrate temperatures. The solutions reveal that the interaction between evaporation and non-Newtonian suspension rheology gives rise to several distinct regimes of droplet shapes and particle concentration distributions. The results provide insight into how evaporation and suspension rheology can be tuned to minimize sagging as well as the well-known coffee-ring effect, an outcome which is important for industrial coating processes.

Electrostatic suppression of the "coffee stain effect"

The dynamics of a slender, evaporating, particle-laden droplet under the effect of electric fields are examined. Lubrication theory is used to reduce the governing equations to a coupled system of evolution equations for the interfacial position and the local, depth-averaged particle concentration. The model incorporates the effects of capillarity, viscous stress, Marangoni stress, elecrostatically induced Maxwell stress, van der Waals forces, concentration-dependent rheology, and evaporation. Via a parametric numerical study, the one-dimensional model is shown to recover the expected inhomogeneous ring-like structures in appropriate parameter ranges due to a combination of enhanced evaporation close to the contact line, and resultant capillarity-induced flow. It is then demonstrated that this effect can be significantly suppressed via the use of carefully chosen electric fields. Finally, the three-dimensional behavior of the film and the particle concentration field is briefly examined.

Patterned deposition at moving contact lines

When a simple or complex liquid recedes from a smooth solid substrate it often leaves a homogeneous or structured deposit behind. In the case of a receding non-volatile pure liquid the deposit might be a liquid film or an arrangement of droplets depending on the receding speed of the meniscus and the wetting properties of the system. For complex liquids with volatile components as, e.g., polymer solutions and particle or surfactant suspensions, the deposit might be a homogeneous or structured layer of solute--with structures ranging from line patterns that can be orthogonal or parallel to the receding contact line via hexagonal or square arrangements of drops to complicated hierarchical structures. We review a number of recent experiments and modelling approaches with a particular focus on mesoscopic hydrodynamic long-wave models. The conclusion highlights open question and speculates about future developments.

Gradient dynamics description for films of mixtures and suspensions: dewetting triggered by coupled film height and concentration fluctuations

A thermodynamically consistent gradient dynamics model for the evolution of thin layers of liquid mixtures, solutions, and suspensions on solid substrates is presented which is based on a film-height- and mean-concentration-dependent free energy functional. It is able to describe a large variety of structuring processes, including coupled dewetting and decomposition processes. As an example, the model is employed to investigate the dewetting of thin films of liquid mixtures and suspensions under the influence of effective long-range van der Waals forces that depend on solute concentration. The occurring fluxes are discussed, and it is shown that spinodal dewetting may be triggered through the coupling of film height and concentration fluctuations. Fully nonlinear calculations provide the time evolution and resulting steady film height and concentration profiles.

Structural transitions in a ring stain created at the contact line of evaporating nanosuspension sessile drops

Monodisperse nanosuspension droplets, placed on a flat surface, evaporated following the stick-slip motion of the three-phase contact line. Unexpectedly, a disordered region formed at the exterior edge of a closely packed nanocolloidal crystalline structure during the "stick" period. In order to assess the role of particle velocity on particle structuring, we did experiments in a reduced pressure environment which allowed the enhancement of particle velocity. These experiments revealed the promotion of hexagonal packing at the very edge of the crystallite with increasing velocity. Quantification of particle velocity and comparison with measured deposit shape for each case allowed us to provide a tentative description of the underlying mechanisms that govern particle deposition of nanoparticles at the triple line of an evaporating droplet. Behavior is governed by an interplay between the fluid, and hence particle, flow velocity (main ordering parameter) and wedge constraints, and consequently disjoining pressure (main disordering parameter). Furthermore, the formation of a second disordered particle region at the interior edge of the deposit (towards bulk fluid) was found and attributed to the rapid motion of the triple line during the "slip" regime. Additionally, the magnitude of the pinning forces acting on the triple line of the same drops was calculated. These findings provide further insight into the mechanisms of the phenomenon and could facilitate its exploitation in various nanotechnological applications.

Stitching chemically converted graphene on solid surfaces by solvent evaporation

The suspension of chemically converted graphene (CCG) provides a cost-effective and facile approach to construct graphene-based materials. However, wrinkles and aggregates usually occur when transferring graphene from suspension to solid-state, which significantly alter the optical, electrical, and electrochemical properties of deposited graphene. Our effort is devoted to the control of the morphology of individual graphene sheet deposited on solid surfaces by the solvent evaporation. Here we have studied the effects of additional components (e.g., organic solvent and electrolyte) in the CCG suspension on the CCG morphology. It was found that the CCG sheets could be stitched together and the graphene monolayer could be flattened by the addition of appropriate additives to the CCG suspension.

Nonlinear dynamics of thin liquid films consisting of two miscible components

Recently, we systematically derived a system of two coupled conservation equations governing a thin liquid layer with a deformable surface composed of two completely miscible components [Phys. Fluids 22, 104102 (2010)]. One equation describes the location of the free surface and the second one the evolution of the mean concentration. This lubrication model was investigated previously in linearized form. The study is now extended to the fully nonlinear case of thin liquid films of a binary mixture (in one and two horizontal spatial dimensions) with and without heat transport. For an initially flat and motionless film heated from below, we analyze the component separation induced by the Soret effect. Nonlinear simulations show that the Soret effect can cause a multitude of interesting behaviors, such as oscillatory patterns and solitonlike structures (localized traveling drops or holes). A stronger component separation induced by stronger Soret effects favors faster-moving localized structures. For isothermal systems, we study the fusion and the mixing of two thin liquid films of different but perfectly miscible liquids. Marangoni-driven forces can cause delayed coalescence, ripple formation, and fingering patterns at the borderline between the two liquid layers. A systematic analysis for ripple pattern formation and finger instabilities at different diffusion constants shows that these phenomena appear more pronounced for lower diffusion in the system.

Fast evaporation of spreading droplets of colloidal suspensions

When a coffee droplet dries on a countertop, a dark ring of coffee solute is left behind, a phenomenon often referred to as the coffee-ring effect. A closely related yet less-well-explored phenomenon is the formation of a layer of particles, or skin, at the surface of the droplet during drying. In this work, we explore the behavior of a mathematical model that can qualitatively describe both phenomena. We consider a thin axisymmetric droplet of a colloidal suspension on a horizontal substrate undergoing spreading and evaporation. In contrast to prior work, precursor films (rather than pinned contact lines) are present at the droplet edge, and evaporation is assumed to be limited by how quickly molecules can transfer out of the liquid phase (rather than by how quickly they can diffuse through the gas phase). The lubrication approximation is applied to simplify the mass and momentum conservation equations, and the colloidal particles are allowed to influence the droplet rheology through their effect on the viscosity. By describing the transport of the colloidal particles with the full convection-diffusion equation, we are able to capture depthwise gradients in particle concentration and thus describe skin formation, a feature neglected in prior models of droplet evaporation. The highly coupled model equations are solved for a range of problem parameters using a finite-difference scheme based on a moving overset grid. The presence of evaporation and a large particle Peclet number leads to the accumulation of particles at the liquid-air interface. Whereas capillarity creates a flow that drives particles to the droplet edge to produce a coffee ring, Marangoni flows can compete with this and promote skin formation. Increases in viscosity due to particle concentration slow down droplet dynamics and can lead to a reduction in the spreading rate.

Temperature distribution inside an evaporating two-dimensional droplet lying on curved or flat substrates

The temperature field inside an evaporating two-dimensional droplet resting on curved or flat isothermal substrates is studied under various evaporation conditions. An analytical solution for the temperature is derived, which can be directly used in the cases of pinned and depinned contact lines as well as in stick-slip evaporation modes. It is found that the temperature drop at the free surface of a droplet on a convex, hydrophobic substrate is far greater than that for flat or concave substrates of the same hydrophobicity. The analytical solution for the temperature field allows the direct estimate of the local temperature gradient and, hence, of the local surface tension gradient that gives rise to Marangoni flow directed from the contact lines to the top of the droplet.

Flow transition within an evaporating binary mixture sessile drop

The flow field along the base of an evaporating ethanol-water droplet and its evolution time was measured by particle image velocimetry. Three stages are revealed, a first stage dominated by multiple vortices, a second transition stage characterized by a remarkable spike in outward flow not previously identified, and a third stage dominated by outward flow identical to that found for pure water. Stage I is thought to be driven by surface tension gradients arising from local concentration variation. The spike in outward flow is explained in terms of a transition corresponding to almost total depletion of ethanol. An exponential decay in vorticity during the transition stage is explained in terms of ethanol diffusion from the bulk to the interface. We speculate on the existence of a zero-concentration wave propagating from the apex to the contact line corresponding to the final total depletion of ethanol.

Quantitatively probing the means of controlling nanoparticle assembly on surfaces

As a means of developing a simple, cost-effective, and reliable method for probing nanoparticle behavior, we have used atomic force microscopy to gain a quantitative 3D visual representation of the deposition patterns of citrate-capped Au nanoparticles on a substrate as a function of (a) sample preparation, (b) the choice of substrate, (c) the dispersion solvent, and (d) the number of loading steps. Specifically, we have found that all four parameters can be independently controlled and manipulated in order to alter the resulting pattern and quantity of as-deposited nanoparticles. From these data, the sample preparation technique appears to influence deposition patterns most broadly, and the dispersion solvent is the most convenient parameter to use in tuning the quantity of nanoparticles deposited onto the surface under spin-coating conditions. Indeed, we have quantitatively measured the effect of surface coverage for both mica and silicon substrates under preparation techniques associated with (i) evaporation under ambient air, (ii) heat treatment, and (iii) spin-coating preparation conditions. In addition, we have observed a decrease in nanoparticle adhesion to a substrate when the ethylene glycol content of the colloidal dispersion solvent is increased, which had the effect of decreasing interparticle-substrate interactions. Finally, we have shown that substrates prepared by these diverse techniques have potential applicability in surface-enhanced Raman spectroscopy.

Wetting and spreading of nanofluids on solid surfaces driven by the structural disjoining pressure: statics analysis and experiments

The wetting and spreading of nanofluids composed of liquid suspensions of nanoparticles have significant technological applications. Recent studies have revealed that, compared to the spreading of base liquids without nanoparticles, the spreading of wetting nanofluids on solid surfaces is enhanced by the structural disjoining pressure. Here, we present our experimental observations and the results of the statics analysis based on the augmented Laplace equation (which takes into account the contribution of the structural disjoining pressure) on the effects of the nanoparticle concentration, nanoparticle size, contact angle, and drop size (i.e., the capillary and hydrostatic pressure); we examined the effects on the displacement of the drop-meniscus profile and spontaneous spreading of a nanofluid as a film on a solid surface. Our analyses indicate that a suitable combination of the nanoparticle concentration, nanoparticle size, contact angle, and capillary pressure can result not only in the displacement of the three-phase contact line but also in the spontaneous spreading of the nanofluid as a film on a solid surface. We show here, for the first time, that the complete wetting and spontaneous spreading of the nanofluid as a film driven by the structural disjoining pressure gradient (arising due to the nanoparticle ordering in the confined wedge film) is possible by decreasing the nanoparticle size and the interfacial tension, even at a nonzero equilibrium contact angle. Experiments were conducted on the spreading of a nanofluid composed of 5, 10, 12.5, and 20 vol % silica suspensions of 20 nm (geometric diameter) particles. A drop of canola oil was placed underneath the glass surface surrounded by the nanofluid, and the spreading of the nanofluid was monitored using an advanced optical technique. The effect of an electrolyte, such as sodium chloride, on the nanofluid spreading phenomena was also explored. On the basis of the experimental results, we can conclude that a nanofluid with an effective particle size (including the electrical double layer) of about 40 nm, a low equilibrium contact angle (<3°), and a high effective volume concentration (>30 vol %) is desirable for the dynamic spreading of a nanofluid system with an interfacial tension of 0.5 mN/m. Our experimental observations also validate the major predications of our theoretical analysis.

Drying of a solution in a meniscus: a model coupling the liquid and the gas phases

A model simulating the drying of a solution in a meniscus in contact with a moving substrate is developed. It takes into account the hydrodynamics in the solution in the framework of the lubrication approximation, the vapor diffusion in the gas phase, and the variation of physical properties during drying. The free surface profile and spatial evaporation flux are not imposed a priori but result from the simulation of the mass transfer in the liquid/gas system (1.5-sided model). Several regimes are observed depending on the substrate velocity. For a large substrate velocity, the classical Landau-Levich regime is obtained. For smaller velocities, a drying front appears that is characterized by a strong concentration gradient and a peak in the evaporation flux. The coupling between the evaporation flux and the meniscus shape in this regime is analyzed. Another regime appears at a very low substrate velocity and seems to be driven by a competition between advection and diffusion. This macroscopic model simulates recent experimental results, namely, the dependence of the deposit thickness on the substrate velocity, which scales as 1/V in the regime dominated by evaporation.

Bubble formation in a quiescent pool of gold nanoparticle suspension

This paper begins with an extensive review of the formation of gas bubbles, with a particular focus on the dynamics of triple lines, in a pure liquid and progresses into an experimental study of bubble formation on a micrometer-sized nozzle immersed in a quiescent pool of aqueous gold nanofluid. Unlike previous studies of triple line dynamics in a nanofluid under evaporation or boiling conditions, which are mainly caused by the solid surface modification due to particle sedimentation, this work focuses on the roles of nanoparticles suspended in the liquid phase. The experiments are conducted under a wide range of flow rates and nanoparticle concentrations, and many interesting phenomena are revealed. It is observed that nanofluids prevent the spreading of the triple line during bubble formation, i.e. the triple line is pinned somewhere around the middle of the tube wall during the rapid bubble formation stage whereas it spreads to the outer edge of the tube for pure water. A unique 'stick-slip' movement of the triple line is also observed for bubbles forming in nanofluids. At a given bubble volume, the radius of the contact line is found to be smaller for higher particle concentrations, but a reverse trend is found for the dynamic bubble contact angle. With the increase of particle concentration, the bubble frequency is raised and the bubble departure volume is decreased. The bubble shape is found to be in a good agreement with the prediction from Young-Laplace equation for given flow rates. The influence of nanoparticles on other detailed characteristics related to bubble growth inside, including the variation of bubble volume expansion rate, the radius of the curvature at the apex, the bubble height and bubble volume, is revealed. It is suggested that the variation of surface tensions and the resultant force balance at the triple line might be responsible for the modified dynamics of the triple line.

Nanoparticle self-structuring in a nanofluid film spreading on a solid surface

Liquids containing nanoparticles (nanofluids) exhibit different spreading or thinning behaviors on solids than liquids without nanoparticles. Previous experiments and theoretical investigations have demonstrated that the spreading of nanofluids on solid surfaces is enhanced compared to the spreading of base fluids without nanoparticles. However, the mechanisms for the observed enhancement in the spreading of nanofluids on solid substrates are not well understood. The complex nature of the interactions between the particles in the nanofluid and with the solid substrate alters the spreading dynamics [Wasan, D. T.; Nikolov, A. D. Nature 2003, 423, 156]. Here, we report, for the first time, the results of an experimental observation of nanoparticles self-structuring in a nanofluid film formed between an oil drop and a solid surface. Using a silica-nanoparticle aqueous suspension (with a nominal diameter of 19 nm and 10 vol %) and reflected light interferometry, we show the nanoparticle layering (i.e., stratification) phenomenon during film thinning on a smooth hydrophilic glass surface. Our experiments revealed that the film thickness stability on a solid substrate depends on the film size (i.e., the drop size). A film formed from a small drop (with a high capillary pressure) is thicker and contains more particle layers than a film formed from a large drop (with a lower capillary pressure). The data for the film-meniscus contact angle verses film thickness (corresponding to the different number of particle layers) were obtained and used to calculate the film structural energy isotherm. These results may provide a better understanding of the complex phenomena involved in the enhanced spreading of nanofluids on solid surfaces.

Modelling approaches to the dewetting of evaporating thin films of nanoparticle suspensions

We review recent experiments on dewetting thin films of evaporating colloidal nanoparticle suspensions (nanofluids) and discuss several theoretical approaches to describe the ongoing processes including coupled transport and phase changes. These approaches range from microscopic discrete stochastic theories to mesoscopic continuous deterministic descriptions. In particular, we describe (i) a microscopic kinetic Monte Carlo model, (ii) a dynamical density functional theory and (iii) a hydrodynamic thin film model. Models (i) and (ii) are employed to discuss the formation of polygonal networks, spinodal and branched structures resulting from the dewetting of an ultrathin 'postcursor film' that remains behind a mesoscopic dewetting front. We highlight, in particular, the presence of a transverse instability in the evaporative dewetting front, which results in highly branched fingering structures. The subtle interplay of decomposition in the film and contact line motion is discussed. Finally, we discuss a simple thin film model (iii) of the hydrodynamics on the mesoscale. We employ coupled evolution equations for the film thickness profile and mean particle concentration. The model is used to discuss the self-pinning and depinning of a contact line related to the 'coffee-stain' effect. In the course of the review we discuss the advantages and limitations of the different theories, as well as possible future developments and extensions.