Comparison of particles.

## Abstract

Nanobiotechnology has huge number of applications in medical science thereby improving health care practices. Keeping in mind the applications of nanoparticles and the convection patterns in biological fields, behaviour of nanofluids is explored for small temperature difference in the layer. The flow of nanofluids is usually described by system of differential equations. A mathematical model for the system based on conservation laws of mass, momentum and energy is formed. To get the insight of the problem, complex equations are simplified wherever needed to get interesting results without violating the necessary physics. The influence of physical properties such as density and conductivity of metallic/non-metallic nanoparticles is examined on the onset of convection currents in the fluid layer.

### Keywords

- nanofluids
- natural convection
- conservation equations
- metallic and non-metallic nanoparticles

## 1. Introduction

In 1959, the celebrated physicist and Nobel laureate Richard Feynman presented an idea of nanotechnology in his talk “There is a plenty of room at the bottom-An invitation to enter a new field of physics” by emphasizing on the fact that the laws of physics allow us to arrange the atoms the way we want. Almost a century ago, Maxwell [1] initiated working on this issue theoretically and unveiled that the particles of size of micrometer and millimeter, if used in traditional fluids can resolve the motive in a more efficient manner. Yet they had few drawbacks like clogging, erosion in micro channel and settling down which were curbed with the evolution of better substitute; nanosized particles (called as nanoparticles). The suspension of nanoparticles in the regular fluids comprised the nanofluids [2]. Nanofluids have also shown many interesting properties, and the distinctive features (refer Table 1) resulting in unprecedented potential for many applications particularly in biological, medical and biomedical applications.

Properties | Microparticles | Nanoparticles |
---|---|---|

Stability | Not stable | Stable |

Surface/Volume ratio | One | One thousand times more than that of microparticles |

Conductivity | Less | High |

Clogging | More | Negligible |

Erosion | Yes | No |

Nanoscale phenomenon | No | Yes |

The catalytic role of nanoparticles in intensifying the thermal conductivity of nanofluids is analyzed by many researchers: Masuda et al. [3], Eastman et al. [4], Das et al. [5] and others. In 2006, Buongiorno [6] pioneered the formulation of conservation equations of nanofluids by incorporating the impacts of diffusion due to Brownian motion and thermophoresis of nanoparticles. He made an observation that the velocity of nanoparticles can be perceived as a sum of base fluid and relative (slip) velocities. To prosecute his research, he considered seven slip mechanisms; inactivity, magnus effect, Brownian motion, diffusiophoresis, thermophoresis, gravitational settling and fluid drainage. Throughout his investigation, he agreed that from all these seven techniques, Brownian diffusion and thermophoresis have a significant part in the absence of turbulent effects. Choi et al. [7] found that carbon nanotubes provide highest thermal conductivity enhancement of nanofluids. There are ample number of evaluations on thermal conductivity of nanofluids [8, 9, 10, 11] in which they discussed and analyzed the theoretical as well as experimental results. Heat transfer in nanofluids because of convection has been examined and contemplated by Das and Choi [12], Ding et al. [13] and Das et al. [14].

The ballistic character of heat transfer within nanoparticles has been studied by Chen [15]. Abnormal increase in viscosity is generally observed in relation to the base fluid. The presence of nanoparticles has found to enhance thermal conductivity [4, 7, 16, 17, 18, 19]. At very low nanoparticle volume fractions (<0.1%), a heat transfer enhancement up to 40% has been reported [8] and this percentage is found to enhance with temperature [5] and concentration of nanoparticles [16]. The results of Choi et al. [7] established the unexpected non-linear character of measured thermal conductivity with nanotube loadings at low concentration while all theoretical studies concluded a linear relationship. Also, it was discovered that thermal conductivity strongly depends on temperature [5] and particle size [20]. Pak and Cho [21] in their study also reported the heat transfer data for turbulent flow of nanofluids having nanoparticles as aluminum and titanium in circular tubes. They found that Nusselt number is up to 30% more than that of base fluid. Nowadays nanofluids are also used in drug delivery systems [22] and advanced nuclear systems [23] due to enriched thermal properties. The nanofluid technology is still in its early stage and various researchers are using nanofluids as a tool to solve technological riddles of the modern society. Figure 1 establishes big impact of small particles in view of the diverse applications of nanofluids in fields of industrial, residential, biomedical and transportation.

These days, nanoparticles are used in almost every biomedical application. Recent usage of nanotechnology in medicine and cancer therapy has attracted a lot of interest in thermal properties of nanofluid such as blood with nanoparticles suspension. Researchers have made the efforts to construct a mathematical model that shows the physical system or phenomenon nearly exact behaviour [24, 25]. Motivated by their work, we also intended to form an analytical model for the analysis of the convection currents in a horizontal nanofluid layer which is in accordance with the physical laws. Consequently, the onset of convection currents in the nanofluid layer is investigated mathematically with the help of partial differential equations. To begin with, equations are non-dimensionalized to get Rayleigh number in the system. Then small disturbances are added to the initial flow and new set of equations are obtained. Further PDE’s are converted into ordinary differential equations using normal modes and expression for Rayleigh number is obtained. It is found that density and conductivity of nanoparticles are important parameters in deciding the stability of the system.

## 2. Instability of fluids under small temperature gradient: Rayleigh Bênard convection

The convective motions occur in a fluid layer heated underside in which a small temperature gradient is maintained across its boundaries. The maintained temperature across the boundaries must surpass a certain value before the instability can manifest itself. This Phenomenon was discovered by Bénard [26] in 1900. In most of his experiments, he found that if a fluid layer is heated underside, the layer at the bottom expands due to higher temperature. This makes the fluid density lighter at the bottom than that on the top making the system top heavy. Here viscosity and thermal diffusivity tend to oppose the convective motions but with the application of higher temperature gradient across the fluid layer, the thermal convection process gets initiated showing the pattern of cellular motions (called Bénard convection). Bénard [27] performed an experiment with metallic plate with a thin non-volatile liquid layer of 1 mm depth maintained under constant temperature.

Keeping the upper layer of fluid exposed to free air, he observed that the fluid layer was decomposed into number of cells (showing cellular motion) called Bénard cells. Thus in the standard Bénard problem, density difference due to variation in temperature across the upper and lower boundaries of the fluid becomes the main reason for the occurrence of instability. Figure 2 shows the schematic representation of Rayleigh-Bénard convection. Rayleigh [28] was the first person who gave an analytical treatment of the problem related to identifying the conditions responsible for breakdown of basic state. As a subsequent work carried out by Rayleigh and Bénard, thermal instability of fluids is known as Rayleigh Bénard convection. The condition for convective motions (depends on temperature gradient) can be represented in dimensionless form by the critical Rayleigh number. He figured out the condition for the instability of free surfaces by showing that the instability would occur on a large temperature gradient

## 3. Conservation equations for a nanofluid layer

We start this section with the description of Boussinesq approximation which is used to write the conservation equations of nanofluids in simplified form.

As is the case of regular fluid [32], equations of nanofluids are difficult to solve because of their non-linear character. Therefore some mathematical approximations are to be used to simplify the basic equations without violating the physical laws. The contribution of Boussinesq [40] in the solution of thermal instability problems is in the form of approximations which is after his name. This approximation has been used by a many researchers for solving different problems of fluids. Boussinesq suggested that inertial effects of density variations can be neglected as compared to its gravitational effects as such situations exist in the domain of meteorology and oceanography. So, density is assumed to be constant everywhere in the equations of motion except in the term with external force. Therefore, we change

Anoop et al. [41] explained various experimental techniques using which nanoparticles can be suspended in the base fluid and that suspension remain stable for several weeks. Buongiorno [6] adopted the formalism of Bird et al. [42] and Chandrasekhar [32] to write conservation equations for nanofluids by considering nanoscale effects; Brownian diffusion and thermophoresis. A model for convective transport in regular fluids was reformulated for nanofluids to accommodate these nanoscale effects as follows.

The random motion of nanoparticles is called Brownian motion and results into the continuous collisions with the base fluid molecules. The Brownian diffusion coefficient due to Brownian motion is given by

where

where

Thermophoresis is the phenomenon in which particles diffuse due to temperature gradient and the effect is similar to one of well-known effects of solute; Soret effect. The thermophoretic velocity is defined as.

Here,

where

The nanoparticles mass flux due to Brownian diffusion (Eq. (1)) and thermophoresis (Eq. (4)) are used to develop a two-component model for convective transport in nanofluids with the following assumptions:

The nanofluid flow is incompressible.

There are no chemical reactions in the fluid layer.

The external forces are negligible.

The mixture is dilute with nanoparticle volume fraction less than 1%.

The viscous dissipation is negligible in the fluid.

The radiative heat transfer is negligible.

The nanoparticles and base fluid are locally in thermal equilibrium.

The seven equations based on basic conservation laws with the above mentioned assumptions are given as follows.

Equation of state (one).

Equation of continuity (one).

Equation of nanoparticles (one).

Equations of motion (three).

Equation of energy (one).

### 3.1 Equation of state

Variables of state depend only upon the state of a system. The physical quantities: *T*; the temperature and *ρ*; the density are the variables of state. We have three thermodynamic variables and a relation between them is given as.

For substances with which we shall be principally concerned, the equation of state can be written as

where *T0* is the temperature at which *ρ = ρ0*.

### 3.2 Equation of continuity-conservation of mass

The equation of continuity for nanofluids is

where

For an incompressible flow (using equation of state)

so that the Eq. (7) reduces to

and in vector form continuity equation for nanofluid is expressed as

### 3.3 Equation of nanoparticles-conservation of mass

The conservation equation for nanoparticles in absence of chemical reactions is

where

Combining Eqs. (11) and (12), nanoparticles conservation equation becomes

Eq. (13) reveals that the nanoparticles move consistently with fluid (second term of left-hand side) and possess velocity relative to fluid (right-hand side) due to Brownian diffusion and thermophoresis.

### 3.4 Equations of motion-conservation of momentum

The equation of motion is derived from Newton’s second law of motion which states that

The momentum equation for nanofluid with negligible external forces is

where

Thus Eq. (14) becomes

Note that in the absence of nanoparticles, Eq. (16) reduces to momentum equation for regular fluid.

### 3.5 Equation of energy-conservation of energy

The thermal energy equation for nanofluid with the assumptions (i)–(v) is

where

Substituting Eq. (18) in Eq. (17), we get

with assumption of negligible external forces

Note that if

Thus, Eqs. (10), (13), (16), (20) constitute the convective transport model for nanofluids which further can be solved for different parameters once the initial and boundary conditions are known. It is interesting to note that all the equations are strongly coupled meaning thereby that the one parameter depends on various other parameters.

Let us introduce non-dimensional variables to get the expression for thermal Rayleigh number as.

where

Using Eqs. (21), (10), (13), (16), (20) after dropping the dashes are

where thermal Rayleigh number

## 4. Initial and perturbed flow

At the initial state, it is assumed that nanoparticle volume fraction is constant and fluid layer is still while temperature and pressure vary in horizontal direction. We get initial solution of Eqs. (22)–(25) using the fact that thermal diffusivity is very large as compared to Brownian diffusion coefficient (refer Buongiorno [1]) as

Let us add perturbations to initial solution and write

The Eq. (27) in Eqs. (22)–(25) give

Making use of the identity

where

## 5. Method of normal modes

To change PDE’s to ODE’s, Eqs. (30)–(32) are solved using normal mode analysis and perturbed variables are written as

Thus above mentioned equations reduce to

where

We write

using of the orthogonality to the functions; gives eigenvalue equation as

where

## 6. Results and discussion

## 7. Discussions on analytical results using various metallic/non-metallic nanoparticles

Table 2 shows the ratios of density to conductivity of various metallic/non-metallic nanoparticles and density 997.1 and conductivity 0.613 of water is used.

Physical properties | Al | Cu | Ag | Fe | Al_{2}O_{3} | SiO_{2} | CuO | TiO_{2} |
---|---|---|---|---|---|---|---|---|

2700 | 9000 | 10,500 | 7900 | 3970 | 2600 | 6510 | 4250 | |

237 | 401 | 429 | 80 | 40 | 10.4 | 18 | 8.9 | |

11.3 | 22.4 | 24.47 | 98.7 | 99.25 | 250 | 361.6 | 477.5 |

It is observed that ratio of density to conductivity is accountable for hastening the onset of convection in the system. The ratio is more for non-metals than metals establishes the lesser stability of non-metallic nanoparticles than metals. Alumina is most stable and titanium oxide is least stable among the nanoparticles under consideration. Density is found to be more influential than conductivity towards deciding the onset of convection in the layer.

## 8. Conclusions

Tremendous applications of nanofluids in pharmaceutical industry with respect to drug invention and cancer imaging motivated the scientists to study convection currents in fluid layer mathematically as well as experimentally. In the present work, the onset of instability in layer is studied under small temperature difference with the help of equations based on conservation laws. The expression of non-dimensional number Rayleigh number is found analytically which decides the instability of the system. Approximations are made whenever needed without violating the necessary physics to get the useful results. Analysis reveals that lesser the ratio of density to conductivity, higher is the stability of the layer. It is found that convection currents majorly depends on density and conductivity and precisely concluding density is more pronounced property than conductivity of nanoparticles. Metallic oxides make the system more stable than metallic nanoparticles in the fluid.