SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

Dynamic contact angle of water-based titanium oxide nanofluid

Milad Radiom1, Chun Yang2* and Weng Kong Chan2

Author Affiliations

1 Department of Chemical Engineering, Virginia Tech, Blacksburg, VA 24060, USA

2 School of Mechanical and Aerospace Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798, Singapore

For all author emails, please log on.

Nanoscale Research Letters 2013, 8:282  doi:10.1186/1556-276X-8-282

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/8/1/282

Received:27 February 2013
Accepted:15 April 2013
Published:11 June 2013

© 2013 Radiom et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


This paper presents an investigation into spreading dynamics and dynamic contact angle of TiO2-deionized water nanofluids. Two mechanisms of energy dissipation, (1) contact line friction and (2) wedge film viscosity, govern the dynamics of contact line motion. The primary stage of spreading has the contact line friction as the dominant dissipative mechanism. At the secondary stage of spreading, the wedge film viscosity is the dominant dissipative mechanism. A theoretical model based on combination of molecular kinetic theory and hydrodynamic theory which incorporates non-Newtonian viscosity of solutions is used. The model agreement with experimental data is reasonable. Complex interparticle interactions, local pinning of the contact line, and variations in solid–liquid interfacial tension are attributed to errors.

Dynamic contact angle; Hydrodynamic theory; Molecular kinetic theory; Nanofluids; Nanoparticles; Non-Newtonian fluid; 68.08.Bc


Industrial operations such as spin coating, painting, and lubrication are based on spreading of fluids over solid surfaces. The fluid may be simple [1-3] or particulate such as paint, ink, or dye [4]. For many years, capillary flow of simple fluids has received considerable attention, and physics of capillary action is known for a long time [5-9]. In addition, capillary flow of micellar surfactant solutions which contain monodisperse and naturally stabilized nanoparticles has been studied [10-14]. However, the same study on liquids laden with metallic and oxide nanoparticles such as silver, copper, zinc oxide, and titanium oxide is scarce. These fluid suspensions are termed as nanofluids after the seminal work by Choi and Eastman [15]. The application of nanofluids is coined with enhanced heat transfer performance compared with their base fluids. They are proposed for applications in cooling of electronic devices, ventilation and air conditioning, and biomedical applications [14,16-24].

It is known that out-of-equilibrium interfacial energy (σ(cos θ0 − cos θ)) provides free energy of capillary flow where σ is the liquid-air surface tension and θ0 and θ are the equilibrium and dynamic contact angles, respectively. During capillary flow, the free energy is dissipated by two mechanisms [5]: (1) contact line friction (T ∑ l) which occurs in proximity of three-phase contact line (solid–liquid–air). The friction at the three-phase contact line is due to intermolecular interactions between solid molecules and liquid molecules. (2) Wedge film viscosity (W) which occurs in the wedge film region behind the three-phase contact line. Lubricating and rolling flow patterns in the wedge film region result in the dissipation of the free energy. For each mechanism of energy dissipation, a theory is developed: (1) molecular kinetic theory (MKT) [25,26] models the contact line friction, and (2) hydrodynamic theory (HDT) [27,28] models the wedge film viscosity. For partial wetting systems (θ0 > 10°), it is assumed that both dissipative mechanisms coexist and models that combine MKT and HDT are developed by Petrov [29] and De Ruijter [30]. In Petrov's model, it is assumed that the equilibrium contact angle θ0 is not constant and its change is described by MKT. In De Ruijter's model, it is assumed that θ0 is constant and the dissipation functions are added to form the total dissipation function, tot = T ∑ l + W. These models are developed for Newtonian fluids and show generally good agreement with experimental data [31].

This paper presents an investigation into spreading dynamics and dynamic contact angle of TiO2-deionized (DI) water nanofluids. Metal oxide TiO2 nanoparticle was chosen for its ease of access and popularity in enhanced heat removal applications. Various nanoparticle volume concentrations ranging from 0.05% to 2% were used. The denser solutions exhibit non-Newtonian viscosity at shear rate ranges that are common to capillary flow. To model experimental data a theoretical model based on combination of MKT and HDT similar to De Ruijter's model is used. The non-Newtonian viscosity of the solutions is incorporated in the model.


Preparation of nanofluids

The solutions were prepared by dispersing 15 nm TiO2 nanoparticles (anatase, 99%, Nanostructured and Amorphous Materials Inc., Houston, TX, USA) in DI water. Oleic acid is reported to stabilize TiO2 nanoparticles in DI water [20] and was added to the mixture at 0.01vol.% concentration. The solution was stirred for 8 h followed by 100 min sonication (Sonicator 3000, 20 kHz and 80 kW, MISONIX, Farmingdale, NY, USA). Temperature of the solution was maintained at 25°C during the sonication process. Clustering and morphology of nanoparticles are important factors in nanofluid spreading capability. We used transmission electron microscopy (TEM) nanographs of TiO2 nanoparticles to examine these factors. An aliquot of dilute solution was dropped and dried on a carbon-coated copper grid. TEM images were then taken immediately. Figure 1 shows that the solution contains irregular particle clusters in addition to monodispersed particles. The sizes of the single particles were found to be close to 15 nm as specified by the supplier. The morphology of the monodispersed particles is spherical. Sonication of the nanofluid solution and addition of surfactant molecules is critical to break down the particle agglomerations and stabilize particle dispersion. The effective nanoparticle size was 260 nm measured with a particle size analyzer (Brookhaven Instruments Corporation, Holtsville, NY, USA). Adsorption of oleic acid surfactant molecules to the surface of TiO2 particles and dissociation of proton from carboxylic acid head groups result in net negative charges on the surface of particles and thus formation of electric double layer around them. Thick electric double layers cause the deviation of particle-particle interactions from hard-sphere interactions. The (Debye) length in nanometer of an electric double layer of 1:1 electrolyte in water at 25°C can be approximated by <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M1">View MathML</a> (where M is the molar concentration). For 0.01 vol.% concentration of oleic acid in water (which is 3.15 × 10-4 molar), the Debye length is estimated to be about 16.9 nm. Such a small increase in the effective diameter of particles allows for an assumption of hard-sphere interactions between particles in the solution which is an important assumption in using Krieger's formula [32]. All other experimental measurements were carried out at 25°C.

thumbnailFigure 1. TEM nanographs of 15 nm TiO2nanoparticles.

Measurement of viscosity

Viscosity of the solutions was measured using a controllable low shear rate concentric cylinders rheometer (Contraves, Low Shear 40, Zurich, Switzerland). The viscosity was measured at shear rates ranging from 0 to 50 s−1. This range corresponds to the shear rates that are common to capillary flow.

Measurement of surface tension

Surface tension of the solutions was measured by pendant droplet method using FTA200 system (First Ten Angstroms, Inc., Portsmouth, VA, USA). To form the pendant droplets, the solutions were pumped out of a syringe system at a very low rate, namely 1 μl/s, to minimize inertia effects. To minimize errors due to evaporation, surface tension was measured right after the pendant droplet reached its maximum volume, namely 10 μl for the dense solutions.

Measurement of dynamic contact angle

Dynamic contact angle of the solutions was measured using the FTA200 system. A droplet of solution was generated at a very low rate (1 μl/s) and detached from the syringe needle tip as soon as it touched the borosilicate glass slide. A video was captured while the droplet was spreading over the glass slide from initial contact to equilibrium position (see Figure 2: the time frame elapses between (a) to (b) and (b) to (c) are 5 and 100 s, respectively). The consecutive photographs were used to measure the contact angles. The spatial resolution was estimated to be about 50 μm on the basis of the focused area and camera pixel size. The standard deviation for contact angle measurements was less than 1°. The temporal resolution was estimated based on the frame speed of the CCD camera as 30 fps. For each concentration, three experiments were performed and average was taken.

thumbnailFigure 2. Consecutive photographs of spreading droplet detached from syringe needle tip.


Empirical analysis of viscosity

From Figure 3, it is obvious that 0.5%, 1%, and 2% solutions exhibit shear thinning viscosity at shear rates below 20 s−1. At higher shear rates, Newtonian behavior was observed for all solutions. For dilute solutions, 0.1 vol.% and 0.05 vol.%, a weak shear thinning behavior was also observed at very low shear rates [19].

thumbnailFigure 3. Viscosity of TiO2-DI water solutions.

A power-law equation is used to model the shear rate and nanoparticle concentration dependent viscosity:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M2">View MathML</a>


where ηb is the viscosity of DI water equal to 0.927 mPa s, F(ϕ) is a function of nanoparticle volume concentration (ϕ), <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M3">View MathML</a> is an indicator of shear thinning viscosity with K as the proportionality factor, and n as the power-law index. F(ϕ) is calculated using Krieger's formula [32]:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M4">View MathML</a>


where ϕmax is the fluidity limit that is empirically equal to 0.68 for hard spherical particles. In Equation 1, n and K are empirical constants which are obtained by fitting this equation to the experimental data shown in Figure 3. Table 1 shows the values of K and n for various nanoparticle volume concentrations. It is obvious that higher nanoparticle concentration results in a larger non-Newtonian behavior. Figure 3 also shows that the power-law Equation 1 is in good agreement with the experimental data.

Table 1. Power-law viscosity, surface tension, and equilibrium contact angle of TiO2-DI water solutions

Molecular kinetic theory

Schematic of a spreading droplet of radius r and contact angle θ that is inspired by De Gennes [5] and Blake [26] is depicted in Figure 4. Based on MKT [26], the rate of displacement of the three-phase contact line over adsorption sites on solid surface, U, is equal to the net frequency of molecular movements, KW (KW = K+ − K, where K+ is the frequency of forward motion and K is the frequency of backward motion), multiplied by average distance between the adsorption sites, λ:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M5">View MathML</a>


thumbnailFigure 4. Schematic of a spreading droplet.

The equilibrium frequency of the three-phase contact line motion (<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M6">View MathML</a>) is obtained from Eyring's theory of absolute reaction rates [26]:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M7">View MathML</a>


where <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M8">View MathML</a>,and T are the Boltzmann constant, Planck constant, Avogadro's number, and absolute temperature, respectively. In this equation, <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M9">View MathML</a> is the equilibrium free energy of capillary flow. An imbalance of the three interfacial tensions near the three-phase contact line, solid–liquid (σsl), solid-vapor (σsa), and liquid–vapor (σ), results in the out-of-equilibrium interfacial energy (σ(cos θ0 − cos θ)) which changes the total free energy of capillary flow. The frequency of the three-phase contact line motion in forward direction (+) and backward direction (−) is [26]:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M10">View MathML</a>


where n is the number of adsorption sites per unit area on solid surface. The net frequency of contact line motion is then as follows [26]:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M11">View MathML</a>


For small arguments of sinh, Equations 3 and 6 result in linear MKT [31]:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M12">View MathML</a>


where <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M13">View MathML</a> is in units of Pa s and is termed as the coefficient of friction at the three-phase contact line. It is noted that this equation is identical to equation twenty-two of [33] for U = 0 and σ cos(θ0) = σsa − σsl (Young's equation). Left hand side (LHS) of Equation 7 is the out-of-equilibrium interfacial energy which is the driving force of capillary flow. Right hand side (RHS) of Equation 7 only includes dissipation of the free energy due to the contact line friction. De Ruijter et al. [30] showed that the corresponding dissipation function (l) is:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M14">View MathML</a>


In the next section, the wedge film viscous dissipation is calculated and added to Equation 8 to form the total dissipation function from which the total drag force is calculated. The total drag force is then equated to the LHS of Equation 7 to form the complete equation of the three-phase contact line motion.

Hydrodynamic theory

To calculate the wedge film viscous dissipation (W), Navier–Stokes equation of motion is solved in the wedge film region. From Figure 4 for the film thickness (H) much smaller than the radial distance ρ (Hρ) and for capillary number Ca ≪ 1, lubrication theory is used:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M15">View MathML</a>


where p is the pressure and u is the velocity distribution at distance x inside the wedge film. For no stress boundary condition at the free fluid-air interface and no slip boundary condition at the solid surface, solution to Equation 9 gives:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M16">View MathML</a>


where ηn is replaced by its expression in Equation 1. The average cross-sectional fluid velocity in the wedge film (<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M17">View MathML</a>) is equal to the three-phase contact line velocity (<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M18">View MathML</a>). This results in:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M19">View MathML</a>


The viscous dissipation in the wedge film can be obtained as follows [5]:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M20">View MathML</a>


where τ is the shear stress (= ηn ∂ u/∂ z), and xm is the cutoff length similar to slip length in HDT [27,28]. Without consideration of xm, dissipation of energy at the wedge film grows infinitely close to the three-phase contact line. For a thin wedge film, solution to Equation 12 gives the wedge film viscous dissipation function:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M21">View MathML</a>


Dynamic contact angle

Combining Equations 8 and 13 gives the total dissipation function [30]:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M22">View MathML</a>


Taking derivative of the total dissipation function with respect to contact line velocity (∂ [l + W]/∂ U) results in the total drag force [5]:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M23">View MathML</a>


Finally, equating Equation 15 with LHS of Equation 7 gives:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M24">View MathML</a>


It is noted that for n = 1 (Newtonian fluid), the integral of Equation 12 results in logarithm ln(r/xm). In this case the final form of Equation 16 is similar to De Ruijter's model [30] (σ(cos θ0 − cos θ) = ζU + 6ηΦ(θ)U ln(r/a)) where Φ = sin 3θ/2 − 3 cos θ + cos 3θ and a is the cutoff length in De Ruijter's model).

In Equation 16, the base radius (r) is in millimeter length scale while the cutoff length (xm) is in nanometer length scale. Thus, rxm, and consequently r1nxm1−n for n ranging from 0.04 to 0.92 (see Table 1). Also, for a sessile droplet of spherical geometry (see Figure 2), the base radius is geometrically related to the dynamic contact angle:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M25">View MathML</a>


where V is the volume of the droplet. Neglecting xm1 − n and substituting r with Equation 17 gives:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M26">View MathML</a>


Equation 18 shows the dynamic contact angle (θ) as a function of contact line velocity (U), solid–liquid molecular interactions (ζ), and non-Newtonian viscosity (n, K). Finally, substituting U with dr/dt = (dr/) × (/dt) the following equation can be obtained for the time evolution of the dynamic contact angle:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M27">View MathML</a>


in which the dynamic contact angle θ = π − α. To compare with experimental data θ is used. Equation 19 is an implicit ordinary differential equation, which cannot be solved analytically, and thus numerical solutions to this equation will be sought.

Results and discussion

The effective diameter of nanoparticles was equal to 260 nm at the lowest solution concentration of 0.05 vol.%. At higher particle concentrations, the increased interparticle interactions result in larger clusters. This increases the possibility of clusters to deposit on the surface of solid and form a new hydrophilic surface. Due to their larger size, these clusters are less possible to deposit on the three-phase contact line, and thus a heterogeneous surface will form: within the wedge film and away from the three-phase contact line, deposition of TiO2 clusters results in a hydrophilic surface with higher surface energy (approximately 2.2 J/m2[34]) than the three-phase contact line where the bare borosilicate glass is present (approximately 0.11 J/m2[35]). The higher surface energy inside the droplet shrinks the wetted area by increasing the equilibrium contact angle (denser solutions are more hydrophilic inside than outside). As a result, solid–liquid interfacial tension increases which on the other hand enhances the equilibrium contact angle [21]. Surface tension of these solutions decreases with particle concentration that is in accordance with Gibb's adsorption isotherm. The shear thinning viscosity of the solutions is due to strong interparticle interaction of the nanoparticle clusters [19,23,36]. Other nanofluids such as ethylene glycol-based ZnO nanofluid [23] and CuO nanofluid [37] also exhibited shear thinning viscosity at low shear rates.

Equation 19 suggests that the contact line friction dissipation (first term on the RHS of Equation 19) and the wedge film viscous dissipation (second term on the RHS of Equation 19) can occur at different time scales [38]. The time dependence of these dissipations has been shown by our experimental data: Figure 5 shows the experimental three-phase contact line velocity (U = dr/dt) plotted versus σ cos θ, where the base radius r is calculated from the experimental dynamic contact angle θ using Equation 17. Figure 5 shows a linear trend that is in accordance with the contact line friction dissipation and a nonlinear trend (see inset of Figure 5) that is in accordance with the wedge film viscous dissipation. This suggests that at the start of capillary flow, the contact line friction is the dominant dissipative mechanism. As capillary flow slows down, the wedge film viscous dissipation becomes more dominant. This corresponds to the solution's higher viscosity at lower shear rates (see Figure 3). Transition to wedge film viscous dominant regime occurs earlier in dilute solutions; for example, Figure 6 shows that for 0.05% concentration the viscous forces start to dominate at time scales around 4 to 8 s while for 2% concentration at time scales around 25 to 32 s.

thumbnailFigure 5. Experimental three-phase contact line velocity (U = dr/dt) plotted versus σ cos θ.

thumbnailFigure 6. Dynamic contact angle of TiO2-DI water solutions.

Figure 6 shows the dynamic contact angle of TiO2-DI water nanofluids at various nanoparticle volume concentrations ranging from 0.05% to 2%. Due to limitation in camera frame per second speed (30 fps), the onset of pendant droplet touching the surface of solid cannot be determined accurately. Hence, the time axis in Figure 6 was shifted to where all of the captured images were readable to the FTA200 software. From Figure 6, it is obvious that for higher nanoparticle concentrations, the contact angles are higher. Figure 6 also shows that the spreading of these nanofluids starts from a primary region where the contact angle changes rapidly followed by a region where the contact angle changes more gradually (note that in a very short period of time (less than 300 ms), the contact angle evolves from 180° at point of contact to angles that are readable to our software and are plotted in Figure 6 at the shifted zero time). In the primary region, the contact line friction dissipation predominates the wedge film viscous dissipation causing fast reduction in the contact angle; then the wedge film viscous dissipation controls the droplet spreading [31].

Using Equation 19, ζ is obtained for the best fit of theory to experimental data that gives the least squared error. Figure 7 shows a reasonable comparison between experimental data and theory. The error is calculated from the following equation and is reported in Table 2:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M28">View MathML</a>


where N is the number of experimental points, <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M29">View MathML</a> is the fitted dynamic contact angle, <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M30">View MathML</a> is the experimental dynamic contact angle, and e = 1° is the experimental error associated with the standard deviation in the measurement of the contact angles. A close examination of table three in [31] and table four in [38] reveals that the agreement between experiment and theory in our case is reasonable considering the complexity of the solution.

thumbnailFigure 7. Dynamic contact angle of TiO2-DI water nanofluid, comparison of experiment and theory.

Table 2. Coefficient of contact line frictionζ, theoretical equilibrium contact angle<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M31">View MathML</a>, and error of comparison between theory and experiment

Table 2 shows values of ζ for various nanoparticle volume concentrations. From solution concentration of 0.05% to 0.5% ζ only changes by 5%; however, it drops rapidly for denser solutions. It is possible that the relative higher hydrophobicity at the three-phase contact line for denser solutions lowers the affinity of surface molecules to water molecules, thereby lowering the friction. At dense concentrations, the presence of large amount of nanoparticles in the wedge film varies the flow field structure. Without nanoparticles, it has been stated that there are two flow patterns in the wedge film: rolling and lubricating patterns [5]. Nanoparticles in the wedge film can change these flow patterns and result in more complex flow structures. As a result of these interparticle interactions, dissipation is more pronounced in the wedge film. Equation 19 gives better results at lower nanoparticle concentrations since complex interparticle interactions are less frequent in dilute solutions (see Table 2). Other sources of disagreement between experiment and theory can be local variations in the concentration of the nanoparticles in the nanofluid [21], pinning of the contact line, and variations in solid–liquid interfacial tension (σsl) [18,21]. It is not possible to model all these effects in theory, and only simple models which can accommodate some of these effects can be developed. Also shown in Table 2 are the theoretical equilibrium contact angles, <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M33">View MathML</a>, which are in reasonable agreement with the experimental equilibrium contact angles, <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/282/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/282/mathml/M34">View MathML</a> (see Table 1).


Due to a wide range of industrial applications, studying capillary flow of liquids laden with metallic and metal oxide nanoparticles is important. Metal oxide TiO2 nanoparticles are especially interesting in enhanced heat removal applications. Agglomeration of nanoparticles results in clusters that have larger effective diameter than the actual particle size. These clusters can deposit on the surface of solid substrates and form a heterogeneous surface condition inside the droplet away from the three-phase contact line that increases the equilibrium contact angle. Dynamic contact angle of metal oxide TiO2 nanoparticles dispersed in DI water revealed two stages of spreading: rapid reduction in contact angle coincides with contact line friction dissipation governed by MKT while gradual reduction in contact angle coincide with wedge film viscous dissipation governed by HDT. Non-Newtonian viscosity of the solution is incorporated in HDT model to give reasonable comparison with experimental data. Nanoparticles in the wedge film change lubricating and rolling flow patterns and result in complex flow field structures. Including all physical aspects of such complex flow in theory is not feasible at the current stage. Simple theoretical equations can only give reasonable comparisons with experiment.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MR, CY, and WKC contributed equally in carrying out the experimental and theoretical studies. All authors read and approved the final manuscript.


The authors gratefully acknowledge the financial support of the research grant (MOE2009-T2-2-102) from the Ministry of Education of Singapore to CY and the Singapore A*STAR scholarship to MR.


  1. Sikalo S, Tropea C, Ganic EN: Dynamic wetting angle of a spreading droplet.

    Experimental Thermal and Fluid Science 2005, 29:795-802. Publisher Full Text OpenURL

  2. Carre A, Woehl P: Spreading of silicone oils on glass in two geometries.

    Langmuir 2006, 22:134-139. PubMed Abstract | Publisher Full Text OpenURL

  3. Wang MJ, Lin FH, Hung YL, Lin SY: Dynamic behaviors of droplet impact and spreading: water on five different substrates.

    Langmuir 2009, 25:6772-6780. PubMed Abstract | Publisher Full Text OpenURL

  4. Smith JT, Viglianti BL, Reichert WM: Spreading diagrams for the optimization of quill pin printed microarray density.

    Langmuir 2002, 18:6289-6293. Publisher Full Text OpenURL

  5. De Gennes PG: Wetting - statics and dynamics.

    Rev Mod Phys 1985, 57:827-863. Publisher Full Text OpenURL

  6. Marmur A: Equilibrium and spreading of liquids on solid-surfaces.

    Adv Colloid Interface Sci 1983, 19:75-102. Publisher Full Text OpenURL

  7. Fraaije J, Cazabat AM: Dynamics of spreading on a liquid substrate.

    J Colloid Interface Sci 1989, 133:452-460. Publisher Full Text OpenURL

  8. Chen JD, Wada N: Edge profiles and dynamic contact angles of a spreading drop.

    J Colloid Interface Sci 1992, 148:207-222. Publisher Full Text OpenURL

  9. Sikalo S, Wilhelm HD, Roisman IV, Jakirlic S, Tropea C: Dynamic contact angle of spreading droplets: experiments and simulations.

    Phys Fluids 2005, 17:062103. Publisher Full Text OpenURL

  10. Kolev VL, Kochijashky II, Danov KD, Kralchevsky PA, Broze G, Mehreteab A: Spontaneous detachment of oil drops from solid substrates: governing factors.

    J Colloid Interface Sci 2003, 257:357-363. PubMed Abstract | Publisher Full Text OpenURL

  11. Kralchevsky PA, Danov KD, Kolev VL, Gurkov TD, Temelska MI, Brenn G: Detachment of oil drops from solid surfaces in surfactant solutions: molecular mechanisms at a moving contact line.

    Ind Eng Chem Res 2005, 44:1309-1321. Publisher Full Text OpenURL

  12. Nikolov A, Kondiparty K, Wasan D: Nanoparticle self-structuring in a nanofluid film spreading on a solid surface.

    Langmuir 2010, 26:7665-7670. PubMed Abstract | Publisher Full Text OpenURL

  13. Wasan DT, Nikolov AD: Spreading of nanofluids on solids.

    Nature 2003, 423:156-159. PubMed Abstract | Publisher Full Text OpenURL

  14. Matar OK, Craster RV, Sefiane K: Dynamic spreading of droplets containing nanoparticles.

    Physical Review E 2007, 76:056315. OpenURL

  15. Choi SUS, Eastman JA: Enhancing thermal conductivity of fluids with nanoparticles. San Francisco, CA; 1995. [ASME International Mechanical Engineering Congress and Exposition] OpenURL

  16. Duursma G, Sefiane K, Kennedy A: Experimental studies of nanofluid droplets in spray cooling.

    Heat Transfer Engineering 2009, 30:1108-1120. Publisher Full Text OpenURL

  17. Sefiane K, Bennacer R: Nanofluids droplets evaporation kinetics and wetting dynamics on rough heated substrates.

    Adv Colloid Interface Sci 2009, 147–148:263-271. OpenURL

  18. Sefiane K, Skilling J, MacGillivray J: Contact line motion and dynamic wetting of nanofluid solutions.

    Adv Colloid Interface Sci 2008, 138:101-120. PubMed Abstract | Publisher Full Text OpenURL

  19. He Y, Jin Y, Chen H, Ding Y, Cang D, Lu H: Heat transfer and flow behaviour of aqueous suspensions of TiO2 nanoparticles (nanofluids) flowing upward through a vertical pipe.

    Int J Heat Mass Transf 2007, 50:2272-2281. Publisher Full Text OpenURL

  20. Murshed SMS, Leong KC, Yang C: Enhanced thermal conductivity of TiO2-water based nanofluids.

    Int J Therm Sci 2005, 44:367-373. Publisher Full Text OpenURL

  21. Vafaei S, Borca-Tasciuc T, Podowski MZ, Purkayastha A, Ramanath G, Ajayan PM: Effect of nanoparticles on sessile droplet contact angle.

    Nanotechnology 2006, 17:2523. PubMed Abstract | Publisher Full Text OpenURL

  22. Vafaei S, Purkayastha A, Jain A, Ramanath G, Borca-Tasciuc T: The effect of nanoparticles on the liquid–gas surface tension of Bi 2 Te 3 nanofluids.

    Nanotechnology 2009, 20:185702. PubMed Abstract | Publisher Full Text OpenURL

  23. Yu W, Xie H, Chen L, Li Y: Investigation of thermal conductivity and viscosity of ethylene glycol based ZnO nanofluid.

    Thermochim Acta 2009, 491:92-96. Publisher Full Text OpenURL

  24. Wong KV, De Leon O: Applications of nanofluids: current and future.

    Advances in Mechanical Engineering 2010, 2010:519659. OpenURL

  25. Blake TD, Haynes JM: Kinetics of liquid/liquid displacement.

    J Colloid Interface Sci 1969, 30:421-423. Publisher Full Text OpenURL

  26. Blake TD: The physics of moving wetting lines.

    J Colloid Interface Sci 2006, 299:1-13. PubMed Abstract | Publisher Full Text OpenURL

  27. Voinov OV: Hydrodynamics of wetting.

    Fluid Dynamics 1976, 11:714-721. OpenURL

  28. Cox RG: The dynamics of the spreading of liquids on a solid surface. Part 1. Viscous flow.

    J Fluid Mech 1986, 168:169-194. Publisher Full Text OpenURL

  29. Petrov P, Petrov I: A combined molecular-hydrodynamic approach to wetting kinetics.

    Langmuir 1992, 8:1762-1767. Publisher Full Text OpenURL

  30. De Ruijter MJ, De Coninck J, Oshanin G: Droplet spreading: partial wetting regime revisited.

    Langmuir 1999, 15:2209-2216. Publisher Full Text OpenURL

  31. Seveno D, Vaillant A, Rioboo R, Adão H, Conti J, De Coninck J: Dynamics of wetting revisited.

    Langmuir 2009, 25:13034-13044. PubMed Abstract | Publisher Full Text OpenURL

  32. Phillips RJ, Armstrong RC, Brown RA, Graham AL, Abbott JR: A constitutive equation for concentrated suspensions that accounts for shear-induced particle migration.

    Physics of Fluids A: Fluid Dynamics 1992, 4:30-40. Publisher Full Text OpenURL

  33. Starov VM: Equilibrium and hysteresis contact angles.

    Adv Colloid Interface Sci 1992, 39:147-173. OpenURL

  34. Naicker PK, Cummings PT, Zhang HZ, Banfield JF: Characterization of titanium dioxide nanoparticles using molecular dynamics simulations.

    J Phys Chem B 2005, 109:15243-15249. PubMed Abstract | Publisher Full Text OpenURL

  35. Rhee SK: Surface energies of silicate-glasses calculated from their wettability data.

    J Mater Sci 1977, 12:823-824. Publisher Full Text OpenURL

  36. Chen H, Ding Y, Tan C: Rheological behaviour of nanofluids.

    New J Phys 2007, 9:367. Publisher Full Text OpenURL

  37. Kwak K, Kim C: Viscosity and thermal conductivity of copper oxide nanofluid dispersed in ethylene glycol.

    Korea-Australia Rheology Journal 2005, 17:35-40. OpenURL

  38. De Ruijter MJ, Charlot M, Voué M, De Coninck J: Experimental evidence of several time scales in drop spreading.

    Langmuir 2000, 16:2363-2368. Publisher Full Text OpenURL