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.     

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.     

Estimation of surface free energies and hamaker constants for fibrous solids by wetting force measurements

Wetting force at three-phase line was measured by the Wilhelmy technique using fibrous solids/liquid/liquid systems. Advancing and receding contact angles were calculated from the wetting forces during fiber immersion and emersion. The obtained results showed that contact angle hysteresis was due to the heterogeneity of the fiber surfaces. The dispersive and polar components of surface free energies of the fibers were determined from the advancing and receding contact angles, respectively. The Hamaker constants of the fibers were estimated from the dispersive components of their surface free energies.     

Wetting: Inverse Dynamic Problem and Equations for Microscopic Parameters

Movement of a liquid meniscus in a low-diameter capillary while it is being filled or emptied is considered. The liquid is nonvolatile. Assuming low Reynolds number and low capillary number, the liquid-gas interface shape is studied. Angles of inclination of this boundary to the solid near the contact line are small. Consideration is given to the inverse problem in wetting dynamics: to establish an analytic expression for the universal constant that influences the dynamics of a three-phase contact line. Inverse relations for microscopic parameters in terms of macroscopic measured values obtained in experiments with a meniscus moving through a capillary are derived. The inverse relations are substantiated independently. To do so, numerical experiments for a van der Waals liquid have been carried out, using the de Gennes model of partial wetting. General formulas for microparameters agree well with numerical experiments. The article provides the similarity criterion which influences the wetting in the case of a van der Waals liquid meniscus. The inverse dynamic problem for both an advancing and a receding meniscus is solved. A relation for the critical speed of meniscus recession is proposed. Two contact angles for a meniscus are discussed. Behavior of dynamic contact angles in the vicinity of the critical speed is studied. One of the angles is shown to vanish at less than the critical speed, and the other one, exactly at the critical speed. In the case of an advancing meniscus the equations for microparameters are valid for both partial and complete wetting. The proposed inverse expression for complete wetting allows determination of the maximum precursor film thickness and its dependence on the motion speed (also determination of the Hamaker constant in the case of a van der Waals liquid). Copyright 2000 Academic Press.     

Evaporative characteristics of sessile nanofluid droplet on micro‐structured heated surface

Micro‐structure patterned substrates attract our attention due to the special and programmable wettabilities. The interaction between the liquid and micro/nano structures gives rise to controllable spreading and thus evaporation. For exploration of the application versatility, the introduction of nanoparticles in liquid droplet results in interaction among particles, liquid and microstructures. In addition, temperature of the substrates strongly affects the spreading of the contact line and the evaporative property. The evaporation of sessile droplets of nanofluids on a micro‐grooved solid surface is investigated in terms of liquid and surface properties. The patterned nickel surface used in the experiments is designed and fabricated with circular and rectangular shaped pillars whose size ratios between interval and pillars is fixed at 5. The behavior is firstly compared between nanofluid and pure liquid on substrates at room temperature. For pure water droplet, the drying time is relatively longer due to the receding of contact line which slows down the liquid evaporation. Higher concentrations of nanoparticles tend to increase the total evaporation time. With varying concentrations of graphite at nano scale from 0.02% to 0.18% with an interval at 0.04% in water droplets and the heating temperature from 22 to 85°C, the wetting and evaporation of the sessile droplets are systematically studied with discussion on the impact parameters and the resulted liquid dynamics as well as the stain. The interaction among the phases together with the heating strongly affects the internal circulation inside the droplet, the evaporative rate and the pattern of particles deposition.     

Modeling contact angle hysteresis of a liquid droplet sitting on a cosine wave-like pattern surface

A liquid droplet sitting on a hydrophobic surface with a cosine wave-like square-array pattern in the Wenzel state is simulated by using the Surface Evolver to determine the contact angle. For a fixed drop volume, multiple metastable states are obtained at two different surface roughnesses. Unusual and non-circular shape of the three-phase contact line of a liquid droplet sitting on the model surface is observed due to corrugation and distortion of the contact line by structure of the roughness. The contact angle varies along the contact line for each metastable state. The maximum and minimum contact angles among the multiple metastable states at a fixed viewing angle correspond to the advancing and the receding contact angles, respectively. It is interesting to observe that the advancing/receding contact angles (and contact angle hysteresis) are a function of viewing angle. In addition, the receding (or advancing) contact angles at different viewing angles are determined at different metastable states. The contact angle of minimum energy among the multiple metastable states is defined as the most stable (equilibrium) contact angle. The Wenzel model is not able to describe the contact angle along the three-phase contact line. The contact angle hysteresis at different drop volumes is determined. The number of the metastable states increases with increasing drop volume. Drop volume effect on the contact angles is also discussed.     

Wetting of real solid surfaces: new glance on well-known problems

The wetting of solid surfaces is treated. The Young and receding contact angles are experimentally unattainable values for the majority of solid surfaces. Actually, we always observe the apparent contact angle. This makes the characterization of wetting of real surfaces problematic. It is proposed in this paper to characterize wetting of real surfaces with the advancing contact angle and the minimum work of adhesion calculated according to the Dupre equation. The advancing contact angle, which depends slightly on the experimental technique used for its measurement, corresponds to the maximal solid/liquid surface tension and correspondingly to the minimal work of adhesion, calculated according to the Dupre equation.     

Comparison of contact angle hysteresis of different probe liquids on the same solid surface

Advancing and receding contact angles of water, formamide and diiodomethane were measured on 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) layers deposited on three different solid supports—glass, mica and poly(methyl methacrylate). Up to five statistical monolayers were deposited on the surfaces by spreading DPPC solution. It was found that even on five statistical DPPC monolayers, the hysteresis of a given liquid depends on the kind of solid support. Also on the same solid support the contact angle hysteresis is different for each probe liquid used. The AFM images show that the heights of roughness of the DPPC films cannot be the primary cause of the observed hysteresis because the heights are too small to cause the observed hystereses. It is believed that the hysteresis is due to the liquid film present right behind the three-phase solid surface/liquid drop/gas (vapour) contact line and the presence of Derjaguin pressure. The value of contact angle hysteresis depends on both the solid surface and liquid properties as well as on intermolecular interactions between them.     

Wetting on axially-patterned heterogeneous surfaces

Contact angle variability, leading to errors in interpretation, arises from various sources. Contact angle hysteresis (history-dependent wetting) and contact angle multiplicity (corrugation of three-phase contact line) are irrespectively the most frequent causes of this uncertainty. Secondary effects also derived from the distribution of chemical defects on solid surfaces, and so due to the existence of boundaries, are the known "stick/jump-slip" phenomena. Currently, the underlying mechanisms in contact angle hysteresis and their connection to "stick/jump-slip" effects and the prediction of thermodynamic contact angle are not fully understood. In this study, axial models of smooth heterogeneous surface were chosen in order to mitigate contact angle multiplicity. For each axial pattern, advancing, receding and equilibrium contact angles were predicted from the local minima location of the system free energy. A heuristic model, based on the local Young equation for spherical drops on patch-wise axial patterns, was fruitfully tested from the results of free-energy minimization. Despite the very simplistic surface model chosen in this study, it allowed clarifying concepts usually misleading in wetting phenomena.     

Truncated versus extended microfilms at a vapor-liquid contact line on a heated substrate

The microstructure of a contact line formed by a liquid and its pure vapor on a perfectly wetted superheated smooth substrate, with the disjoining pressure most often in the form of a positive inverse cubic law (nonpolar case), is routinely considered to end up in a microfilm extended over adjacent "dry" parts of the solid surface. Invoking the spreading coefficient as an additional independent parameter within this framework, we argue however that a regime with a truncated microfilm is chosen instead if the spreading coefficient is decreased below a positive (still perfect wetting) critical value dependent upon the superheat, in which case the extended-microfilm thickness is surpassed by that of the "pancake" introduced by de Gennes and co-workers. Conversely, for a given positive spreading coefficient, there is a critical superheat above which the microfilm gets truncated, whereas for a negative one (partial wetting) the truncated regime should be preferred at any superheat. A parametric study of the apparent contact angle (a nonlinear eigenvalue of the steady microstructure problem) versus the spreading coefficient is carried out. When the latter is negative, Young's law is asymptotically recovered. Microfilm fronts on a bare surface are shown to be advancing or receding in accordance with the selected regime. A slightly more general class of disjoining pressures is also touched upon. The analysis is based in part upon thermodynamic considerations and in part upon a standard one-sided model of an evaporating liquid layer in the lubrication approximation.     

Constitutive modeling of contact angle hysteresis

We introduce a phase field model of wetting of surfaces by sessile drops. The theory uses a two-dimensional non-conserved phase field variable to parametrize the Gibbs free energy of the three-dimensional system. Contact line tension and contact angle hysteresis arise from the gradient term in the free energy and the kinetic coefficient respectively. A significant advantage of this approach is in the constitutive specification of hysteresis. The advancing and receding angles of a surface, the liquid-vapor interfacial energy and three-phase line tension are the only required constitutive inputs to the model. We first simulate hysteresis on a smooth chemically homogeneous surface using this theory. Next we show that it is possible to study heterogeneous surfaces whose component surfaces are themselves hysteretic. We use this theory to examine the wetting of a surface containing a circular heterogeneous island. The contact angle for this case is found to be determined solely by the material properties at the contact line in accord with recent experimental data.     

Does macroscopic flow geometry influence wetting dynamic?

The macroscopic flow geometry has long been assumed to have little impact on dynamic wetting behavior of liquids on solid surfaces. This study experimentally studied both spontaneous spreading and forced wetting of several kinds of Newtonian and non-Newtonian fluids to study the effect of the macroscopic flow geometry on dynamic wetting. The relationship between the dynamic contact angle, θ(D), and the velocity of the moving contact line, U, indicates that the macroscopic flow geometry does not influence the advancing dynamic wetting behavior of Newtonian fluids, but does influence the advancing dynamic wetting behavior of non-Newtonian fluids, which had not been discovered before.     

Wetting phenomena at the CO2/water/glass interface

A novel high-pressure apparatus and technique were developed to measure CO2/water/solid contact angles (theta) in situ for pressures up to 204 bar. For two glass substrates with different hydrophilicities, theta increased significantly with CO2 pressure. As the pressure was increased, an increase in the cohesive energy density of CO2 caused the substrate/CO2 and water/CO2 interfacial tensions (gamma) to decrease, whereas the water/substrate gamma value increased. theta for the more hydrophobic substrate was predicted accurately from the experimental water/CO2 gamma value and an interfacial model that included only long-range forces. However, for the more hydrophilic substrate, short-range specific interactions due to capping of the silanol groups by physisorbed CO2 resulted in an unusually large increase in the water/substrate gamma value, which led to a much larger increase in theta than predicted by the model. A novel type of theta hysteresis was discovered in which larger theta values were observed during depressurization than during pressurization, even down to ambient pressure. Effective receding angles were observed upon pressurization, and effective advancing angles were observed upon depressurization on the basis of movement of the three-phase contact line. The greater degree of hysteresis for the more hydrophilic silica can be attributed in part to the capping of silanol groups with CO2. The large effects of CO2 on the various interfacial energies play a key role in the enhanced ability of CO2, relative to many organic solvents, to dry silica surfaces as reported previously on the basis of FTIR spectroscopy (Tripp, C. P.; Combes, J. R. Langmuir 1998, 14, 7348-7352).     

Determination of nanolayer thickness and effective thermal conductivity of nanofluids

Nanofluids having high thermal conductivity enhancement relative to conventional pure fluids are fluids engineered by suspending solid nanoparticles into base fluids. In the present study, calculating the Van der Waals interaction energy between a nanoparticle and an ordered liquid nanolayer around it, the nanolayer thickness was determined, the average velocity of the Brownian motion of nanoparticles in a fluid was estimated, and by taking into account both the aggregation of nanoparticles and the presence of a nanolayer a new thermal conductivity model for nanofluids was proposed. It has been shown that the nanolayer thickness in nanofluids is independent on the radius of nanoparticles when the radius of the nanoparticles is much greater than the nanolayer thickness and determines by the specific interaction of the given liquid and solid nanoparticle through the Hamaker constant, the surface tension and the wetting angle. It was proved that the frequency of heat exchange by fluid molecules is two orders of magnitude higher than the frequency of heat transfer by nanoparticles, so that the contribution due to the Brownian motion of nanoparticles in the thermal conductivity of nanofluids can be neglected. The predictions of the proposed model of thermal conductivity were compared with the experimental data and a good correlation was achieved.     

Drops down the hill: theoretical study of limiting contact angles and the hysteresis range on a tilted plate

The limiting inclination angle (slip angle), for which a two-dimensional water drop may be at equilibrium on a chemically heterogeneous surface, is exactly calculated for a variety of cases. The main conclusion is that, in the cases studied, the contact angles at the upper and lower contact line do not always simultaneously equal the receding and advancing contact angles, respectively. On a hydrophobic surface, the lowest contact angle (at the upper contact line) tends to be approximately equal to the receding contact angle, while the highest contact angle (at the lower contact line) may be much lower than the advancing contact angle. For hydrophilic surfaces, the opposite is true. These conclusions imply that the hysteresis range cannot in general be measured by analyzing the shape of a drop on an inclined plane. Also, the limiting inclination angle cannot in general be calculated from the classical equation based only on the advancing and receding contact angles.     

Partial wetting of chemically patterned surfaces: the effect of drop size

Partial wetting of chemically heterogeneous substrates is simulated. Three-dimensional sessile drops in equilibrium with smooth surfaces supporting ordered chemical patterns are considered. Significant features are observed as a result of changing the drop volume. The number of equilibrated drops is found either to remain constant or to increase with growing drop volume. The shape of larger drops appears to approach that of a spherical cap and their three-phase contact line seems, on a larger scale, more circular in shape than that of smaller drops. In addition, as the volume is increased, the average contact angle of drops whose free energy is lowest among all equilibrium-shaped drops of the same volume appears to approach the angle predicted by Cassie. Finally, contrary to results obtained with two-dimensional drops, contact angle hysteresis observed in this system is shown to exhibit a degree of volume dependence in the advancing and receding angles. Qualitative differences in the wetting behavior associated with the two different chemical patterns considered here, as well as differences between results obtained with two-dimensional and three-dimensional drops, can possibly be attributed to variations in the level of constraint imposed on the drop by the different patterns and by the dimensionality of the system.     

Drops sitting on a tilted plate: receding and advancing pinning

The wetting behavior of a liquid drop sitting on an inclined plane is investigated experimentally and theoretically. Using Surface Evolver, the numerical simulations are performed based on the liquid-induced defect model, in which only two thermodynamic parameters (solid-liquid interfacial tensions before and after wetting) are required. A drop with contact angle (CA) equal to θ is first placed on a horizontal plate, and then the plate is tilted. Two cases are studied: (i) θ is adjusted to the advancing CA (θ(a)) before tilting, and (ii) θ is adjusted to the receding CA (θ(r)) before tilting. In the first case, the uphill CA declines and the downhill CA remains unchanged upon inclination. When the tilted drop stays at rest, the pinning of the receding part of the contact line (receding pinning) and the depinning of the advancing part of the contact line (advancing depinning) are observed. The free energy analysis reveals that upon inclination, the reduction of the solid-liquid free energy dominates over the increment of the liquid-gas free energy associated with shape deformation. In the second case, the downhill CA grows and the uphill CA remains the same upon inclination. Advancing pinning and receding depinning are noted for the tilted drop at rest. The free energy analysis indicates that upon inclination, the decrease of the liquid-gas free energy compensates the increment of the solid-liquid free energy. The experimental results are in good agreement with those of simulations.     

Spreading of triple line and dynamics of bubble growth inside nanoparticle dispersions on top of a substrate plate

This work investigates the feasibility of engineering surface wettability by using different nanoparticles. As an illustration, detailed formation of gas bubbles on top of a stainless steel substrate plate in a quiescent pool of aqueous gold and alumina nanofluids is studied. The presence of nanoparticles is shown to be able to modify the dynamics of triple line and bubble growth significantly. An early pinning of the bubble triple line is observed and a larger bubble contact angle is found for bubbles growing in a gold nanofluid, whereas an opposite phenomenon is observed for bubbles growing in an alumina nanofluid compared to those of pure water. Other bubble parameters such as departure volume, bubble frequency, and waiting time of bubble formation are also affected by the presence of nanoparticles. The variation of solid surface tensions due to the existence of nanoparticles and the resultant force at the triple line should be responsible for such differences. Such results illustrate the big potential of nanoparticle in engineering surface wettability of a solid-liquid-gas system.     

Contact angle hysteresis of microbead suspensions

Microbead suspensions are often used in microfluidic devices for transporting biomolecules. An experimental investigation on the wettability of microbead suspension is presented in this study. The variation in the surface tension and the equilibrium contact angle with the change in the volume fraction of the microbead is presented here. The surface tension of the microbead suspension is measured with the pendant drop technique, whereas the dynamic contact angle measurements, i.e., advancing and receding contact angles, are measured with the sessile drop technique. An equilibrium contact angle of a suspension with particular volume fraction is determined by computing an average over the measured advancing and receding contact angles. It is observed that the surface tension and the equilibrium contact angle determined from advancing and receding contact angles vary with the magnitude of the microbeads volume fraction in the suspension. A decrease in the surface tension with an increase in the volume fraction of the microbead suspension is observed. The advancement and the recession in contact line for dynamic contact angle measurements are achieved with the motorized dosing mechanism. For microbead suspensions, the advancement of the contact line is faster as compared to the recession of the contact line for the same flow rate. The presence of microbeads assists in the advancement and the recession of the contact line of the suspension. A decrease in the equilibrium contact angles with an increase in the microbead suspension volume fraction is observed. Inclusion of microbeads in the suspension increases the wetting capability for the considered combination of the microbead suspension and substrate. Finally, empirical correlations for the surface tension and the contact angle of the suspension as a function of microbead volume fraction are proposed. Such correlations can readily be used to develop mechanistic models for the capillary transport of microbead suspensions related to LOC applications.     

Contact angle hysteresis explained

A view of contact angle hysteresis from the perspectives of the three-phase contact line and of the kinetics of contact line motion is given. Arguments are made that advancing and receding are discrete events that have different activation energies. That hysteresis can be quantified as an activation energy by the changes in interfacial area is argued. That this is an appropriate way of viewing hysteresis is demonstrated with examples.