This study is a numerical study on the coagulation as well as the sedimentation effect of nanofluids using the Brownian dynamics method. Three cases are simulated, focusing on the effects of the sizes, volume fraction, and ζ potentials of nano-particles on the formation of coagulation and sedimentation of nanofluids. The rms fluctuation of the particle number concentration, as well as the flatness factor of it, is employed to study the formation and variation of the coagulation process. The results indicate a superposition of coagulation and sedimentation effect of small nano-particles. Moreover, it is stable of nanofluids with the volume fraction of particles below the limit of "resolution" of the fluids. In addition, the effect of ζ potentials is against the formation of coagulation and positive to the stability of nanofluids.
The nanofluid is characterized by the fluid with nanometer-sized solid particles dispersed in solution , which can increase the heat transfer coefficient [2-6], enhance the critical heat flux in boiling heat transfer [7-9], reduce the wall friction force , improve the optical characteristics , etc. Nano-sized particles are utilized because of its better stability than the suspension of micro-sized particles. For a badly stable suspension, sedimentation or coagulation (agglomeration) may occur. It compromises the above-mentioned advantages of the nano-suspension.
As is well known , the occurrences of coagulation and sedimentation are the two main factors for the instability of nanofluid. The phenomenon of coagulation is characterized by the formation of particle clusters, i.e., particles are in contact with each other and the cohesion takes place. Then, the clusters grow up. Many researchers investigated the coagulation effect of particles by the Brownian dynamics simulation, focusing on the formation of gelation , coagulation rates , particle network , etc. For example, Hütter  identified the characteristic coagulation time scales in colloidal suspensions, and measured their dependencies on the solid content and potential interaction parameters. He also deduced different cluster-cluster bonding mechanisms in the presence of an energy barrier, etc. Besides, the sedimentation always occurs after a big particle cluster is established, i.e., the particles within the cluster sediment flow downward because of the increased effect of the gravity of the cluster over the buoyancy force of it, and reduced the effect of Brownian motion to the big cluster. Many researches were devoted to the sedimentation [16-18] using the Brownian dynamics simulation too. For example, Soppe and Jannsen  studied the sediment formation of colloidal particle by a process of irreversible single-particle accretion. They used the algorithm of Ermak and McCammon, incorporating the inter-particle forces and hydrodynamic interaction on the two-particle level, and analyzed the effect of two-particle hydrodynamic interactions on the sediment structure, etc. They found that the process of sediment formation by colloidal particle is the result of a delicate balance of sediment field strength, DLVO interactions, and hydrodynamic interactions.
However, there is an important issue about which few researches have been concerned: the interaction between the coagulation and sedimentation for the instability of nanofluids. For example, the processing of coagulation causes the particle clusters to grow up, and then the large clusters are more prone to sedimentation than that of small clusters because of the intensive gravity effect. In other words, the coagulation effect is able to augment the sedimentation effect. Thus, this study is intended to carry out some research on this issue, exploring the complex interaction as well as the close relation between the coagulation and the sedimentation phenomena.
In this study, the Brownian dynamics technique is employed to investigate the motion of nanoparticles. The governing equation is the so-called Langevin equation, which is formulated as follows :
where the superscript 0 indicates that the variable is corresponding to the beginning of the time step Δt; ri is the ith component of the position vector of particle; Dij is the element of the diffusion tensor indexed by (i, j); Fj is the force experienced by the jth particle; kB is the Boltzmann constant; and T denotes the temperature. The displacement Ri(Δt) is a random displacement of a Gaussian distribution with a zero expectation and a ( ) covariance ( ). In this study, the Rotne-Prager tensor  is utilized as approximations to the hydrodynamic interaction:
where η is the viscosity, a is the particle radius, δij is the Kronecker delta, is the vector from the center of particle i to the center of particle j, and is the unit tensor.
Moreover, three forces, i.e., the attractive Van der Walls force, fv, the repulsive electrostatic force by the electric double layer, fe, and the gravity force fg are considered as the forces experienced by any particle, which are formulated as follows :
where A, d, ε, κ, ζ, ρf, ρ, and g are the Hamaker constant, the particle diameter, the electric permittivity of the fluid, the inverse of the double-layer thickness, the zeta potential of the suspension, the density of the fluid, the density of the particle, and the gravity acceleration, respectively.
It is noted in Equation (4) that the results for rij - d = 0 is meaningless when the contact between the two particles occur, and they will adhere to each other or rebound back. Thus, we treat the condition with as the situation when the two particles are separated, so that Equation (4) works. Otherwise, it results in coalescence between the two particles. Once the coalescence between colliding particles takes place, the clusters start growing up.
In this study, three cases with different diameters of particles (Case 1), different volume fractions (Case 2), and different zeta potentials (Case 3) are simulated, respectively (Table 1). For these cases, the parameters of the material, as well as other relevant parameters, are illustrated in Table 2.
Table 1. Three cases with different diameters of particles, volume fractions, and zeta potentials
Table 2. Parameters used in this simulation
In this simulation, the boundary conditions in the x and y directions (Table 2) in the horizontal plane are both periodic, whereas the top and bottom walls of the simulation domain in the z-direction are treated as adhesive walls to which the particles adhere immediately once they come into contact with them. It is reasonable to conclude thus, since the agglomerated particle clusters always adhere to the bottom walls or the top interfaces.
Initially, a random distribution is given to the particles. As time advances, the possible movements of particles are computed through solution of the governing equations.
Case 1: effect of particle sizes
This section deals with the effect of particle sizes on the coagulation and sedimentation. Figure 1 shows the simulation results at t = 0, 5, 10, and 50 μs for Case 1 due to the effects of different sizes of particles. Figure 1a, b, c, d shows the results of d = 10 nm at t = 0, 5, 10, and 50 μs, respectively. Similarly, Figure 1e, f, g, h shows the results for d = 25 nm, whereas Figure 1i, j, k, l shows them for d = 50 nm. It is seen that the coagulation takes place the most intensively and rapidly for the smallest size of particles (Figure 1a, b, c, d), moderately for the intermediate size of particles (Figure 1e, f, g, h), and weakly and slowly for the largest size of particles (Figure 1e, f, g, h, i, j, k, l). More importantly, the results of the intermediate sizes are due to coagulation but with weak sedimentation, whereas the results of the smallest sizes are due to both the effects of coagulation and sedimentation. For the largest size of the particles, it is neither due to coagulation nor sedimentation. It looks complicated. As is known, the larger particles bear the major effect of gravity, and they are the most prone to sedimentation. However, it is only true of the single particle without coagulation. With the superposition of the coagulation effect, it can amplify or augment the trend of sedimentation through coagulation. Owing to the increasing agglomeration of the particles, the gravity effect may play an important and even a dominant role, which causes possible sedimentation of the whole agglomeration (the upper part of the agglomeration in Figure 1d is due to the adhesive boundary on the upper wall). In other words, there exists a balance between the sedimentation effect of the large-sized individual particles and the sedimentation effect of small-sized aggregated particles. The former is caused solely by the gravity effect, whereas the latter is caused by the superposition of the coagulation and the gravity effects, i.e., the amplification and augmentation of the gravity effects of the aggregated particles due to coagulation.
Figure 1. Snapshots of simulation results. (Case 1 for d = 10 nm (a-d), d = 25 nm (e-h), d = 50 nm (i-l) at t = 0, 5, 10, and 50 μs respectively).
It is necessary to mention that Figure 1i, j, k, l does not indicate the stability of the nanofluids. Alternatively, it indicates a relatively stable characteristic compared to Figure 1a, b, c, d. After the evolution over a long time, possible coagulation or sedimentation may also occur.
In order quantify the degree of coagulation, we need to define some functions. Let us divide the simulation domain Lx × Ly × Lz by Nx × Ny × Nz cubic meshes by the cell volume (δx × δy × δz). The mean number concentration of particles is the mean number of particles within each mesh volume (δx × δy × δz). Then, the rms value of the particle concentration R1 and the flatness factor of the number concentration R4 are formulated as follows:
The rms of concentration means the fluctuation of the number concentration of particles, and it is closely related to the formation of particle clusters due to coagulation. The flatness factor means the intensity of fluctuation of the number concentration, thereby indicating the intensity of coagulation. Thus, these two functions are helpful in enabling the quantification of particle coagulation.
Figure 2a, b shows the R1 and R4 for Case 1. It is seen from Figure 2a, b that the coagulation of the small particles is the fastest. They become almost totally coagulated immediately even at the beginning. Comparatively, the coagulation of the larger particles takes place slowly and increasing steadily. However, the final level of coagulation of the larger particles is greater than that of smaller particles. In addition, the degree as well as the rapidity of the coagulation of the intermediate particles is intermediate between that of the smaller and the larger particles.
Figure 2. The flatness factor of the concentration distribution of nanoparticles. (The R1(a) and R4(b) for Case 1).
Case 2: effects of volume fractions
In this section, the effect of volume fraction, i.e., the concentration of particles, is studied. As aforementioned, the smaller particles are more prone to coagulate than the larger particles, under the same condition of the volume fractions. However, the process of coagulation is also closely related to the number of particles contained in it.
For example, Figure 3a, b, c, d shows one of the results of Case 2 where only np = 400 particles are simulated. Compared to Figure 1a, b, c, d, it is seen that the coagulation does not takes place at t = 50 μs. It says that the coagulation is to be regarded appropriately as a process with an excess of particle content. When the particle numbers go beyond the superior limit of the resolvent, then the coagulation will certainly take place. Thus, when the np = 4200 particles are simulated, more intensive coagulations are observed correspondingly (Figure 3e, f, g, h).
Figure 3. Snapshots of simulation results. (Case 2 for np = 400 (a-d) and np = 4200 (e-h) at t = 0, 5, 10, and 50 μs respectively).
In addition, Figure 4a, b shows the R1 and R4 of Case 2. It is seen from Figure 4 that the R1 and R4 for np = 400 are always relatively of small value, indicating a stable status almost without coagulation and sedimentation, although R4 is slightly fluctuated when t < 0.06. Moreover, compared to np = 1200, it is seen that the concentration fluctuation R1 and flatness factor of concentration R4 for np = 400 are relatively of lower values. It validates the conclusions derived from the observation of Figure 3a, b, c, d, in comparison with Figure 1e, f, g, h.
Figure 4. The flatness factor of the concentration distribution of nanoparticles. (The R1(a) and R4(b) for Case 2).
With the increased number of particles, it is seen that the R1 and R4 are increased too (np = 1200 and 2100, respectively, Figure 4). However, when the particle number is extremely large, all the spaces are almost stuffed with particles, leading to a homogeneous distribution and a low fluctuation in the number concentration (np = 4200, in Figure 4).
Case 3: effects of ζ potentials
The previous sections showed the results with ζ = 0 eV. As seen from Equation (5), no repulsive effect has been considered between the particles since fe = 0. Thus, this section will focus on the effect of the repulsive effect by varying the ζ potentials.
Comparing with Figure 1e, f, g, h, it is seen that the degree of coagulation with ζ = 0.01 eV is attenuated (Figure 5a, b, c, d), and it almost disappears with ζ = 0.05 eV (Figure 5e, f, g, h). It indicates that the repulsive effect induced by the ζ potentials is beneficial to the stability of nanofluids, since it acts against the coagulation process.
Figure 5. Snapshots of simulation results. (Case 3 for ζ = 0.01 eV (a-d) and ζ = 0.05 eV (e-h) at t = 0, 5, 10, and 50 μs, respectively).
This conclusion is also validated by the variations of R1 and R4 (Figure 6a, b, respectively). It is seen that the fluctuation of the number concentration of particles with ζ = 0.05 eV increases much more slowly than the cases with smaller ζ. Although the R1 and R4 for ζ = 0.05 eV are still increasing with time, which indicates that the coagulation process may still take place after a fairly long time, their effects against the formation of coagulation are very clear. In other words, it is positively beneficial for the stability of nanofluids.
Figure 6. The flatness factor of the concentration distribution of nanoparticles. (The R1(a) and R4(b) for Case 3).
The findings of this study are briefly summed up as follows:
1. A complicated superposition of the coagulation and sedimentation effects for small particle is observed. The mechanisms of sedimentation for the larger and the smaller particles are different. The former is caused mainly by the great gravity effect of any individual particle, whereas the latter is mainly due to the coagulation process, and the superposition of coagulation causes the sedimentation of the whole agglomeration of particles.
2. There exists a superior limit of the fluid for particle content. When the volume fraction is below the limit, it is hard for the coagulation to occur. In contrast, the coagulation will certainly take place when the concentration of nanoparticles is beyond the capacity of "resolution" of the fluids.
3. The effect of ζ potentials is beneficial for the stability of nanofluid, since it resists the formation of coagulation. In other words, increase in the value of ζ potentials is helpful to make the nanofluid more stable.
The authors declare that they have no competing interests.
All authors contributed equally.
This study is supported by the National High Technology Research and Development of China 863 Program (2007AA05Z254), for which the authors are grateful.
J Thermophys Heat Transfer 1999, 13:474-480. Publisher Full Text
Trans ASME 1999, 121:280-289. Publisher Full Text
Mater Today 2005, 8:36-44. Publisher Full Text
Annu Rev Mater Res 2004, 34:219-246. Publisher Full Text
Heat Transfer Eng 2008, 29:432-460. Publisher Full Text
Appl Phys Lett 2003, 83:3374-3376. Publisher Full Text
Appl Phys Lett 2007, 90:184107. Publisher Full Text
Int J Heat Mass Transfer 2010, 53:1542-1553. Publisher Full Text
Wear 1998, 216:92-96. Publisher Full Text
J Chem Phys 1986, 85:4079-4086. Publisher Full Text
J Chem Phys 1992, 96:6996-7004. Publisher Full Text
Mol Phys 1999, 96:893-903. Publisher Full Text
J Chem Phys 1978, 69:1352. Publisher Full Text
J Chem Phys 1969, 50:4831. Publisher Full Text