Numerical Analysis of Nanofluids in Differentially Heated Enclosure Undergoing Orthogonal Rotation

Natural convection heat transfer in a rotating, differentially heated enclosure is studied numerically in this paper. The rotating enclosure is filled with water-Ag, water-Cu, water-Al 2 O 3 , or water-TiO 2 nanofluids. The governing equations are in velocity, pressure, and temperature formulation and solved using the staggered grid arrangement togetherwithMACmethod.The governing parameters considered are the solid volume fraction, 0.0 ≤ φ ≤ 0.05, and the rotational speeds, 3.5 ≤ Ω ≤ 17.5 rpm, and the centrifugal force is smaller than the Coriolis force and both forces were kept below the buoyancy force. It is found that the angular locations of the local maximums heat transfer were sensitive to rotational speeds and nanoparticles concentration. The global quantity of heat transfer rate increases about 1.5%, 1.1%, 0.8%, and 0.6% by increasing 1% φ of the nanoparticles Ag, Cu, Al 2 O 3 , and TiO 2 , respectively, for the considered rotational speeds.


Introduction
The fluid flow and heat transfer characteristics of a rotating enclosure were treated in detail according to its functional and practical importance in the design of guided missiles and a rotating electrical machine.Consistently too high temperature occurring in the rotating machine will reduce the lifetime of the machine and may lead to serious failure, so that the efficient cooling system is needed.A new smart generation cooling system can be achieved by utilizing nanofluids.Nanofluids are dilute liquid suspension of nanoparticles with at least one of their principal dimensions smaller than 100 nm.Nanofluids are proposed by Choi [1].Nanoparticles are basically metal (Cu, Ni, Al, etc.), oxides (Al 2 O 3 , TiO 2 , CuO, SiO 2 , Fe 2 O 3 , Fe 3 O 4 , BaTiO 3 , etc.), and some other compounds (AlN, SiC, CaCO 3 , graphene, etc.) and base fluids usually include water, ethylene glycol, propylene glycol, engine oil, and so forth.Due to very small sizes and large specific surface areas of the nanoparticles, nanofluids have superior properties like high thermal conductivity, minimal clogging in flow passages, long-term stability, and homogeneity.Application of nanofluids had been reviewed intensively by Yu and Xie [2].Some theoretical and experimental studies can be found in the literature related to convective heat transfer enhancement by utilizing nanofluids, but for a configuration not closely similar to this work.The review of the rotating cubic and rectangular enclosures problem will be given first.
Buehler and Oertel [3] studied thermal cellular convection in a rotating rectangular box and found that the roll cells changed orientation with increasing Taylor number.The fluid flow and heat transfer characteristics of a rotating square enclosure were studied experimentally and numerically by Hamady et al. [4].They concluded that the Coriolis force arising from rotation may have a remarkable influence on heat transfer when compared with nonrotating results and correlations of Nusselt number as functions of Taylor and Rayleigh number were built.Lee and Lin [5] and Ker and Lin [6,7] studied a differentially heated rotating cubic enclosure.A significant flow modification was obtained when the rotational Rayleigh number is greater than the Rayleigh 2 Advances in Mathematical Physics number or when the Taylor number is greater than the Rayleigh number.They also examined the effect of rotation on flow stabilization.A significant increase or decrease in heat transfer in a rotating and differentially heated square enclosure could be achieved due to rotational effects as reported by Baig and Masood [8] and Baig and Zunaid [9].Jin et al. [10] studied numerically the rectangular enclosure with discrete heat sources and found that rotation results in imbalance of clockwise and counterclockwise circulations, increases heat transfer in the worst stage, reduces the oscillation of Nusselt number, and improves or reduces mean performance in each cycle.The effects of Coriolis force, centrifugal force, and thermal buoyancy force were segregated numerically by Tso et al. [11] on a differentially heated square enclosure.The effects of the Coriolis and centrifugal forces were found small and differentiated from those of other forces.
Natural convection studies of nanofluids in a square enclosure had been investigated by Khanafer et al. [12].They showed that heat transfer coefficient increases with increasing Cu nanoparticles concentration.Oztop et al. [13] observed the heat transfer coefficient increases by adding CuO nanoparticles concentration in an inclined nonuniformly heated enclosure.On the other hand, Hwang et al. [14], Santra et al. [15] and Abu-Nada [16] reported a decrease in the heat transfer by increasing the nanoparticles concentration for a particular Rayleigh number.Ho et al. [17] investigated the effects of uncertainties due to adopting different formulas for the ratio thermal conductivity and dynamic viscosity for water Al 2 O 3 nanofluids.They found that the heat transfer could be enhanced or mitigated depending on the formulas used for the estimated dynamic viscosity of the nanofluids.Recently, Qi et al. [18] and He et al. [19] found that the flow and heat transfer characteristics nanofluid were more sensitive to viscosity than to thermal conductivity.
In the present study, the problem of natural convection heat transfer in a rotating, differentially heated enclosure filled with nanofluids water-Ag, water-Cu, water-Cu, water-Al 2 O 3 , and water-TiO 2 is studied numerically.The effects of the nanoparticles types and concentration as well as the rotational speeds on characteristics of fluid and temperature field as well as heat transfer rate are considered.The numerical results of this study are thought to be useful in the design of an effective cooling system for rotating electrical machine to help ensure effective and safe operational conditions.

Mathematical Formulation
A schematic diagram of a square enclosure of side ℓ executes a steady uniform counterclockwise angular velocity about horizontal axis as shown in Figure 1, with the geometric layout and the Cartesian coordinates (, ) rotating with the enclosure.The fluid in the enclosure is a water-based nanofluids containing Ag, Cu, Al 2 O 3 , or TiO 2 nanoparticles.The surface at  = ℓ/2 has constant hot temperature ( ℎ ) and the surface at  = −ℓ/2 has a constant cold temperature (  ).The temperatures along the lateral wall are assumed to be linearly distributed between  ℎ and   to consider conjugate heat transfer in the lateral wall of the experiments.
The  shown in Figure 1 is defined as an angular position.The Boussinesq approximation is used for natural convection and its validation in the rotation situation can be seen in [3].The terms representing the thermal and rotational buoyancies and Coriolis force are, respectively, equal to  nf g nf ( −   ), − nf g nf ( −   )Ω ×(Ω × r), and −2 nf [1 −  nf ( −   )] ⋅ Ω × V.The continuity, momentum, and energy equations for the laminar flow development can be described as follows: where   is the motion pressure defined as +  nf g sin (Ω) +  nf g cos (Ω) . ( Hence, The effective density of the nanofluids,  nf , is given as and  is the solid volume fraction of nanoparticles.Thermal diffusivity of the nanofluids is where the heat capacitance of the nanofluids given is The thermal expansion coefficient of the nanofluids can be determined by The governing equation ( 1) can be converted to nondimensional forms using the following nondimensional parameters: The Coriolis buoyancy force is neglected, because | nf (−   )| ≪ 1 in the present studies.The nondimensional continuity, momentum, and energy equations are written as follows: where  =  = 0 on the walls and the boundary conditions for the nondimensional temperatures are Θ =  + 0.5 at  = −0.5,Θ =  + 0.5 at  = 0.5, The nondimensional governing equations contain three independent dimensionless parameters, namely, Prandtl number, Rayleigh number, and Taylor number.In addition, rotational Rayleigh number reflects the effect of centrifugal buoyancy force, but it depends on the other dimensionless parameters too [11].
The ratio dynamic viscosity of the nanofluids given by Brinkman [21] is The ratio thermal conductivity of nanofluids restricted to spherical nanoparticles is approximated by the Maxwell-Garnetss (MG) [20] model The viscosity and conductivity of the nanofluids are integrated in Figure 2 and the physical properties of the water and nanoparticles are given in Table 1.
The fluid motion is displayed using the stream function Ψ obtained from velocity components  and .The relationships between the stream function and the velocity components are  = Ψ/ and  = −Ψ/, which yield a single equation, where Ψ = 0 at all walls of the enclosure.The physical quantities of interest in this problem are the average Nusselt number, representation of the average heat transfer rate on the heated or cooled plates that is defined by The time-integrated of the average heat transfer rate in one cycle is

Numerical Method and Validation
Staggered grid arrangement, together with Marker and Cell (MAC) method [22], is adopted to solve the governing equation ( 9) subject to the boundary conditions (10).Due to lack of boundary conditions for pressure, the use of the staggered grid and MAC formulation provided an advantage.That is, one may locate the secondary grid along the boundaries of the domain where only specification of velocity boundary conditions is required but not of the pressure.The fictitious values of velocity outside the domain are obtained by extrapolation of the interior points as given by Hoffmann and Chiang [23].A second-order central difference approximation is used for the space discretization and a first-order approximation is used for temporal derivative.The solution of the Poisson pressure equation is obtained by applying an iterative Gaussian-SOR method.The velocities are then computed by the projection method.
In this study, the convergence criterion for the Poisson equation is set at  = 10 −5 and the time stepping is chosen, Δ = (1/4)(Δ) 2 Pr bf to meet stability criteria where  is the safety factor with a value between 0 and 1. Uniform grid distribution is used for the whole enclosure.The effect of grid resolution was examined in order to select the appropriate grid density as demonstrated in Figure 3(a) for water-Cu nanofluids with volume fraction of 5% and Ta bf = 4.7 × 10 3 at  = /2.The results indicate that an 120 × 120 mesh can be used in the final computations.A typical computation process for the present periodic oscillation problem is shown in Figure 3(b) where, after two rotations, the final periodic oscillation of the Nu is obtained.Therefore, all the computations reported in this work were carried out beyond two rotations.As a validation, our results for the isotherms compare well with those obtained by Hamady et al. [4] and Tso et al. [11] for a pure fluid, Pr = 0.7, Ra = 1.2 × 10 5 , Ta = 8.9 × 10 4 , and Ω = 8.5 rpm as shown in Figure 4.

Results and Discussion
The analyses in the undergoing numerical investigation are performed in the following range of the associated dimensionless groups: the solid volume fraction, 0.0 ≤  ≤ 0.05, and the Taylor number for the base fluid, 8 × 10 2 ≤ Ta bf ≤ 2 × 10 4 as tabulated in Table 2.The Prandtl number and the Rayleigh number are fixed at Pr bf = 6.2 and Ra bf = 4.62 × 10 5 , respectively.The rotational Rayleigh number, (Ra  ) bf , reflecting centrifugal buoyancy force was not discussed except to specify it explicitly.This is due to the fact that the (Ra  ) bf is not being an independent parameter but depending on Pr bf , Ra bf , and Ta bf .It worth mentioning that in the present research, the buoyancy force is stronger than Coriolis forces and the Coriolis force is greater the centrifugal force.Present work will also compare utilization   1.The flow and temperature fields of nanofluids and base fluid during rotation will be shown in Figures 5 and 6.The average Nusselt number over one period will be presented in Figures 7 and 8. Finally, plots of the time-integrated average Nusselt number will be displayed in Figures 10 and 9.
Figures 5(a)-5(h) show the effects of rotation on the flow field for water-Cu nanofluids with volume fraction 5% and the pure water with Ta bf = 4.7 × 10 3 .There are two basic flow structures, clockwise (negative sign) and counterclockwise (positive sign) circulations, depending on the heated position as well as the body forces.The fully clockwise circulations were obtained in Figures 5(a)-5(h).These circulations are induced by the buoyancy force.The weak counterclockwise flows due to the Coriolis force then appear in the upperleft and lower-right corners of the enclosure at Figure 5(e), heating from below.When  > /2, the gravitational buoyancy force now induces a counterclockwise circulation.However the flow inertia from the clockwise circulation still persists and puts off the transition angle for flow reversal and this flow is aided by the Coriolis force on the other hand.Therefore, the clockwise and the counterclockwise flow compete with each other as displayed at Figure 5

(f).
The large vortex driven is by Coriolis force and a smaller vortex is driven by thermal buoyancy force.In Figure 5(g), the clockwise flow totally disappears and replaces the counterclockwise, unicellular motion flow as the angle of the gravitational force increases.This is due to the fact that buoyancy force balances the effect of the Coriolis force in the clockwise circulation.Finally, in Figure 5(h), the transition from clockwise to counterclockwise motion is completed, where the Coriolis force is directed outwards from the core that suppressed the flow.It should have made clear that this sequence repeats itself, the fully clockwise flow then transition flow, and totally counterclockwise flow.It observed that, during rotation, the strength of the flow circulation of nanofluids is always stronger than the base fluid.It observed the flow patterns between the nanofluids and pure water were different enough in some angular locations.It also observed that the streamlines of the nanofluids and pure water in all of Figures 5(a  circulation and then higher temperature gradient near the walls.Later, the temperature gradient near the walls decreases by increasing the rotational angle.This due to the fact that the Coriolis force is directed outwards from the core and tends to stabilize, thus reducing the heat transfer on the wall.The heat transfer rate of the nanofluids is always higher than the base fluid at any angular position.This is due to the fact that, at any locations, nanofluid has higher thermal conductivity leading to higher energy, which accelerates the flow.Actually, this phenomenon is also accompanied by an undesirable effect promoted by higher viscosity that suppresses the flow but is small compared to the favorable effect driven by the presence of the high thermal conductivity. Figure 7 depicts the variation of the average Nusselt number, Nu, over one cycle for several rotational speeds (Taylor numbers) at  = 0.05 and water-Cu.The average Nusselt number profiles display different characteristics for the various Taylor numbers.It observed that increasing the Taylor number reduces the number of local maximum Nu, that is, four, three, two, and one local maximum for Ta bf = 8.0 × 10 2 , Ta bf = 2.0 × 10 3 , Ta bf = 4.7 × 10 3 , and Ta bf = 2.0 × 10 4 , respectively.The local maximums of the heat  (Ψ min ) bf = −29.4(Ψ min ) nf = −30.2 (Ψ min ) bf = −38.9(Ψ min ) nf = −39.9The stagnant Nu by increasing the fluid conductivity is also found in some angular locations.This due to at these locations and the flow is stabilized by the Coriolis force.The Coriolis force is directed outwards from the core that would decelerate the higher velocities promoted by the higher fluid conductivity, thus stagnating the heat transfer on the wall.Figure 9 shows the time-integrated average Nusselt number of the water-Cu against the volume fraction for different Taylor numbers.The global heat transfer performance increases linearly by increasing the nanoparticles concentration.Higher rotational speeds of the nanofluids or pure water lead to better heat transfer rate.This finding agrees with the previous experimental result by Lin and Yan [24].It should be better if we simulate for Ω > 30 rpm, but higher rotational speeds at the higher Rayleigh number bring about difficulties in numerical simulation because of the instability of fluid motion.The Nu in Figure 9 displays monotone curves, where the Nu increases about 1.1% by increasing 1%  for each of rotational speeds.
Figure 10 displays the time-integrated average Nusselt number against the volume fraction for different nanoparticles types at Ta bf = 4.7 × 10 3 .The results demonstrated in this figure help compare the heat transfer appearance of the argentum (Ag), copper (Cu), alumina (Al 2 O 3 ), and titania (TiO 2 ) suspended in the rotating water.The lowest heat transfer rate was obtained for the TiO 2 since the TiO 2 has lower value of thermal conductivity ratio (see Figure 2 and Table 1) that leads to lower thermal energy and smaller enhancements in the heat transfer performance.As it can be seen in Figure 2, increasing  will increase the ratio thermal conductivity monotonically that eventually leads to monotone increasing in the heat transfer rate as shown in Figure 10.The Nu increases about 1.5%, 1.1%, 0.8%, and 0.6% by increasing 1%  of the Ag, Cu, Al 2 O 3 , and TiO 2 , respectively.

Conclusions
The present work has studied natural convection heat transfer characteristics and flow fields in a nanofluids-filled, differentially heated square enclosure.The enclosure was rotated above its longitudinal horizontal axis.The dimensionless forms of the governing equations were solved using the staggered grid arrangement together with MAC method.Detailed computational results for flow and temperature fields and the heat transfer performance have been presented in graphical forms.The periodic oscillation of the flow and temperature fields as well as the heat transfer was obtained.The main conclusions of the present analysis are as follows.
(a) A single primary flow circulation is obtained for the heating from the top and the sides.The flow circulation patterns between the nanofluids and pure water were different enough in some angular locations and it presented centrosymmetry shape during rotation.
(b) The angular locations of the local maximums heat transfer are sensitive to rotational speeds and nanoparticles concentration.
(c) The higher nanoparticles concentration leads to better heat transfer rate for any rotational angles considered here.The global quantity of heat transfer rate increases about 1.5%, 1.1%, 0.8%, and 0.6% by increasing 1%  of the Ag, Cu, Al 2 O 3 , and TiO 2 , respectively, for any rotational speeds.
The theoretical prediction in this paper is expected to be a useful guide for the experimentalists to study the effectiveness of utilizing nanofluids in the cooling system of the rotating equipment.

Nomenclature
Cp: Specific heat capacity (J/kg K) : Gravitational acceleration (m/s 2 ) : Th e r m a lc o n d u c t i v i t y( Wm −1 K −1 ) ℓ: Width and height of enclosure (m) Nu: Nusselt number Ra: Rayleigh number Ra  : Rotational Rayleigh number Ta: Taylor number : Temperature (K) , V:

Figure 1 :
Figure 1: Schematic representation of the model.
)-5(h) present centrosymmetry shape.Figures6(a)-6(h) exhibit the effects of rotation on the temperature field for water-Cu nanofluids with volume fraction of 5% and the pure water with Ta bf = 4.7 × 10 3 .The heated wall positions move according to the angular position.Starting from  = /4 where the heated wall is above the cold wall, the fluid temperatures adjoining the hot wall rise.The fluid moves due to buoyancy force from the left region the right region and creates a clockwise circulation cell as shown in the previous figure.Denser boundary layers near the hot and cold walls were formed as the rotational angle takes higher up to  = .This due to the fact that the Coriolis is directed into the core region that leads to stronger flow
V e l oc i t yc o m po n e n t si nth e-and -direction (m/s) , &, : Space coordinates (m) & dimensionless space coordinates.Greek Symbols : Thermal expansion coefficient (1/K) : Dimensionless time   : Dimensionless time for one rotation Θ: Dimensionless temperature Ω: Magnitude of angular rotation rate (rpm, rad/s) : Density (kg/m 3 ) : Solid volume fraction Ψ: Dimensionless stream function : Angular position of enclosure (rad) : Dynamic viscosity (N s/m 2 ).
bf Figure 2: The ratio thermal conductivity of water-Ag, water-Cu, water-Al 2 O 3 , and water-TiO 2 using Maxwell-Garnett model and the ratio dynamical viscosity using Brinkman model versus volume fraction.

Table 2 :
Fixed variables and constants.
of different nanoparticles that is, Ag, Cu, Al 2 O 3 , or TiO 2 as specified in Table