- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Table of Contents

Advances in Mechanical Engineering

Volume 2013 (2013), Article ID 240438, 9 pages

http://dx.doi.org/10.1155/2013/240438

## Numerical Investigation on Heat and Flow Characteristics of Temperature-Sensitive Ferrofluid in a Square Cavity

^{1}Department of Mechanical Engineering, Dong-A University, Hadan 840, Saha-gu, Busan 604-714, Republic of Korea^{2}Graduate School of Mechanical Engineering, Dong-A University, Hadan 840, Saha-gu, Busan 604-714, Republic of Korea

Received 14 May 2013; Revised 1 October 2013; Accepted 2 October 2013

Academic Editor: Guy Lauriat

Copyright © 2013 Moo-Yeon Lee and Jae-Hyeong Seo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The objective of this paper is numerically to study the heat and flow characteristics of temperature-sensitive ferrofluid in the square cavity with and without the magnetic intensity. The numerical model was developed to predict the behavior of the ferrofluid using finite element method (FEM) and showed good agreement with the existing data within 5% at all Rayleigh number ranges from 10^{3} to 10^{6}. Natural convection and heat transfer characteristics of the ferrofluids within the tested cavity were found to depend on both magnetic intensity and magnetic volume fractions of magnetite. In addition, the mean Nusselt numbers and mean velocity of the ferrofluid in a square cavity were increased with the rise of the magnetic intensities and increased by 23.2% and 143.7%, respectively, at both magnetic intensity of and the elapsed time of seconds.

#### 1. Introduction

The ferrofluid, which means novel functional fluids using functional materials, has both characteristics of metals and fluids. The ferrofluid generally consisted of magnetite nanosized particle of around 1 to 100 nm and carrier liquid such as water, oils, and hydrocarbons with the aid of surfactants in a continuous carrier phase. It can be controlled by both magnitude and direction of an external magnetic field and temperature [1–12]. In addition, because of the nanosized magnetic particles consisting of the ferrofluids and the surfactant attached to magnetic particles, the ferrofluid could be prevented from particles sticking to each other or precipitate with Brownian motion [2]. Therefore, the ferrofluid could be applied for various industrial fields: enhancement and depression of the heat transfer of the thermal devices, magnetic sealing, damping and bearing of machines, energy conversion system, drug delivery for disease curing, optical devices, micro-/nanofluidic devices, and so forth [3–5].

Numerous numerical and experimental studies have been widely investigated scientifically for various scopes of application. Horng et al. [3] reported the structural pattern formation of the magnetic columns in the ferrofluid thin film under magnetic fields and framework of the application of the ferrofluid based on the remarkable optical properties caused by these magnetically induced structures. Wang et al. [6] studied the effects of magnetic force on the natural convection in porous enclosure under magnetic field. They explained the heat and flow characteristics on the effect of the inclination angle, Darcy number, and magnetic force parameter. Basak et al. [7] analyzed the heat flow due to natural convection within porous trapezoidal enclosures with hot bottom wall and cold side walls in presence of insulated top walls. They examined the thermal mixing within the cavity for various material processing applications and suggested the generalized non-Darcy model, neglecting the Forchheimer inertia term, to predict the flow in porous medium. Zhao et al. [8] introduced the novel optical devices based on the tunable refractive index of ferrofluid and their characteristics. They reported that the developed new method would play an important role in the fields of optical information communication and sensing technology. Wang et al. [9] investigated synthesis of Fe_{3}O_{4} ferrofluid used for magnetic resonance imaging and hyperthermia, which was to obtain a stable water based ferrofluid with fast magnetotemperature response for using it both in magnetic resonance imaging (MRI) and ferrofluid hyperthermia (MFH) with protecting magnetic precipitates from oxidation. Jin et al. [10] simulated a porous cavity filled with thermosensitive ferrofluids using the lattice Boltzmann method (LBM). They showed the magnetization, heat, and flow characteristics of ferrofluids on the effect of the Rayleigh number and the magnetic Rayleigh number. However, concise and precise studies on the thermal and its flow control of the ferrofluid in a cavity are few in open literature.

The objective of this study is to numerically investigate the heat and flow characteristics including the control of the heat transfer of the ferrofluid in the square cavity with variations of the magnetic intensity and temperature difference between heated and cooled walls. In addition, the mean Nusselt numbers on heated and cooled walls were calculated with boundary conditions.

#### 2. Numerical Method

##### 2.1. Physical Model

Figure 1 shows the physical model of the square cavity filled with temperature-sensitive ferrofluid. Because the magnetic force was caused by the interactions between the external magnetic field and magnetization [11], in order to observe the effects on the temperature of walls and magnetic intensity of the ferrofluid at the square cavity, temperatures of the left and right walls were maintained to be 290 K and temperatures of the bottom and top walls were 320 K and 290 K, respectively. Magnetic intensities ranged from 0 to 15000 A/m with the interval of 1000 A/m and were applied uniformly along *y-*direction at bottom wall, respectively. Velocities of *x-* and *y-*directions at initial conditions were assumed to be 0 and with no-slip boundary conditions at all walls. Table 1 showed the thermophysical properties of the ferrofluids used in the numerical model. The thermophysical properties of the ferrofluids supplied by ferrofluid manufacturing company were used. The nominal diameter of magnetite particles was 10 nm and the base solution of the used ferrofluids was water. The volume fractions of magnetite of the ferrofluid were varied from 1.2% (FF1) to 4.5% (FF4).

##### 2.2. Governing Equations

The conservation equations for mass, momentum, and energy are expressed using Equations from (1) to (3) as transport equations of a natural convection problem and the governing equations for magnetization and Maxwell are expressed using (4) and (5) suggested by Sawada et al. [16]:

is the magnetization vector of the ferrofluid and expressed as a function of magnetic intensity (), differential magnetic susceptibility (), inertia moment of particle per unit mass (), and fluid temperature [17], where is the density of ferrofluid, is the velocity vector, is the rotational angular velocity vector of the particle, is the couple stress tensor, is the Levi-Civita symbol, is the total stress tensor, is the volume couple vector per unit mass, is the absolute temperature, is internal angular momentum per unit mass, is the heat flux vector, is the viscous dissipation function, is the intensity of heat source per unit mass, is the volume magnetic susceptibility, is the idemfactor, and is the relaxation time of a rotational motion by fluid friction. Consider where magnetic induction vector is expressed as the combination of magnetic intensity and magnetization vector .

Several physical assumptions as well as the well-known Boussinesq approximation are used to solve the above equations [18]. The thermophysical properties of ferrofluid are only function of the temperature and there is no internal heat generation by chemical reaction (). Also, density difference is the only function of temperature difference. Then, the density difference in the mass conservation could be negligible. The particles in the momentum conservation do not have the couple of volumes () and have a rigid body rotation (). Magnetic intensity is perpendicular to the rotational angular velocity vector () of the particle (). And the viscous dissipation function in the energy conservation could be negligible due to the low flow velocity [19]. The governing Equations (1) to (5) could be rearranged into (7) with the above-mentioned assumptions: That is, the governing equations describing a natural convection of a ferrofluid are rearranged as adding the magnetic term to the momentum equation, where is the representative density of ferrofluids, is the pressure, is the dynamic viscosity, is the eddy viscosity, is the equilibrium magnetization strength, is the cubical expansion coefficient, is the fluid temperature, is the representative temperature, is the gravitational acceleration vector, is the specific heat at constant pressure, is the thermal conductivity, and is the rotational angular velocity vector of fluid.

The above governing equations could be rearranged as dimensionless forms like (8) using the dimensionless analysis technique [20]: where is the Reynolds number, is the polarity effect parameter, is the magnetic effect parameter, is the rotational Peclet number, is the Rayleigh number, is the Prandtl number, and is the unit vector of direction of the gravitational acceleration [21, 22]. Dimensionless parameters are expressed as follows: where is the intensity of applied magnetic field, is the relaxation time of a Brown motion, and subscript is the representative value of each parameter in the dimensionless process. Also, the representative velocity, , was calculated using by the approximation of because the flow velocity of the ferrofluid for natural convection was defined as the ratio of the buoyancy force and the viscosity force [23].

The local Nusselt numbers on the heated wall were defined as (10) and the mean Nusselt numbers on the heated wall were calculated using local Nusselt numbers like the following:

##### 2.3. Algorithm

The governing equations were calculated using generalized simplified maker and cell method (GSMAC) suggested by Kawai et al. [24] and the finite element method (FEM) suggested by Vinogradova and Zarubin [25]. And they were finally expressed as forms of the discretization equations from (12) to (13) where Based on the research of Seo [26], the numerical analysis was iteratively calculated until the tolerance is converged. In order to decide the proper grid size and number, mesh sensitivity analysis with different numbers of mesh size (, , and ) was performed under , , aspect ratio = 1, cold wall temperature of 300 K (left, top, and bottom walls), and hot wall temperature of 320 K (Right wall), as shown in Figure 2. Reasonably good agreement of 3.32% deviation was found between the and grids at the value of maximum velocity and therefore the grid used was for considering a proper iteration time and the accuracy of the results. In other words, the case with grids was used for deciding both the accuracy and convergence rate taken into account. Figure 2 shows the comparison of the velocity distributions of the ferrofluid in the square cavity at different number of grids. In addition, the steady state of the ferrofluid in the square cavity was decided at temperature difference of between the initial state and final state during 0.5 second. The interval of isotherms was 0.05.

#### 3. Results and Discussion

##### 3.1. Validation

Validation of numerical solution was carried out and compared with that of Barakos at al. [14] and that of Markatos and Pericleous [15]. Temperatures of both top and bottom walls of the cavity maintain adiabatic state. Temperature of the left wall maintains hot state and temperature of the right wall maintains cold state. All conditions except the above-mentioned conditions were similarly maintained with conditions considered in the previous researches. Table 2 shows the comparison of the mean Nusselt numbers with the existing data. The predictions of the present model showed good agreement within 5% of the data of other researchers. The mean and average deviations between the present data and the existing data were 1.0% and −0.77%, respectively, for Barakos at al. [14], 2.0% and −2.0%, respectively, for Markatos and Pericleous [15], at all Rayleigh numbers ranges from 10^{3} to 10^{6}, as shown in Figure 3.

Comparisons were performed keeping a temperature difference of 20 K between the vertical walls of the square cavity and establishing the reference temperature at 293 K. Therefore, from Figure 3, we could be confident that the present analysis and model used are correct.

##### 3.2. Nature of Isotherms

We present the numerical results for isotherms. The isotherm characteristics of the ferrofluids were analyzed under both nonmagnetic and magnetic intensities and also compared with the data of the Newtonian fluid. Newtonian fluid used in this study was water. And numerical analysis was conducted at and took 7 minutes at each condition.

Figure 4 shows the isotherm characteristics of the Newtonian fluid and ferrofluid with the magnetic volume fractions of 2.0% (FF2) and 4.5% (FF4), respectively, of the square cavity under nonmagnetic intensity and the elapsed time of 10000 seconds. Heat transfer of the Newtonian fluid in the tested cavity was slightly activated more than that of the ferrofluid because the magnetite of the ferrofluid could not increase the heat transfer under nonmagnetic intensity. In addition, the heat transfer decreased with the rise of the magnetic volume fractions from 2.0% (FF2) to 4.5% (FF4). This is because the number of particles of magnetite per unit volume increased by 118.2%. It could not be helpful to activate the heat transfer of the ferrofluid in the tested cavity under nonmagnetic intensity.

Figure 5 shows the isotherm characteristics of ferrofluids under magnetic intensity with the variation of the magnetic volume fractions from 1.2% (FF1) to 4.5% (FF4) of the square cavity after the elapsed time of 10000 seconds. As shown in Figures 4 and 5, isotherms of the ferrofluid under the magnetic intensity expanded more rapidly than that of ferrofluid under the nonmagnetic intensity at the same time due to the increased magnetic volume force at the given magnetic intensity. In addition, the heat transfer of the ferrofluid in the tested cavity under magnetic intensity increased greatly with the rise of the magnetic volume fractions from 1.2% (FF1) to 4.5% (FF4) due to the increased heat transfer. Namely, this is because the diffusive concentration gradients lead to the coupling effect between heat and mass transport as mentioned in Rinaldi et al. [2] Also, the presence of the magnetic intensity tends to accelerate the fluid motion inside the tested cavity as mentioned in Grosan et al. [27].

Figure 6 shows the effect of the isotherms and mean Nusselt numbers on the temperature difference between top and bottom walls of the tested cavity under the elapsed time of seconds, magnetic intensity of , and magnetic volume fractions of 2.0% (FF2). Temperature difference in the tested cavity was changed with 30°C to 50°C. The expansion speed of the isotherms counter of the ferrofluid in the tested cavity increased with the rise of the temperature difference between top and bottom walls because it caused an increased heat flow by natural convection with the rise of the density difference, as shown in Figure 6(a). Namely, increased temperature difference between top and bottom walls of the square cavity means the increased superficial magnetic volume force of the ferrofluid at the given magnetic intensity. In addition, the mean Nusselt number was increased from 12.70 to 12.84 with the rise of the temperature difference from 30°C to 50°C due to the increased superficial magnetic volume force, as shown in Figure 6(b). Also, increased temperature gradient causes a concentration gradient of magnetic particles and leads to diffusive concentration gradients known as Soret effect by Lange [28].

##### 3.3. Nusselt Numbers

Figure 7 shows the mean Nusselt numbers and velocity vectors in the square cavity with the volume fractions of the magnetite and magnetic intensities at the elapsed time of seconds at . And it took 21 minutes at each condition. Figure 7(a) shows the mean Nusselt numbers of the present model. The mean Nusselt numbers of the ferrofluid under magnetic intensities of and , respectively, increased by 9.6% and 23.2%, respectively, with the rise of the volume fraction of magnetite from 1.2% (FF1) to 4.5% (FF4) due to the increased heat transfer with the rise of the magnetic intensities. The increased heat transfer means the increased isotherms expansion speed of ferrofluids in the tested cavity with the rise of the mean velocities as shown in Figure 7(b). The mean velocity increased with the rise of the volume fraction of magnetite and magnetic intensities. Figure 7(b) shows the mean velocity of the present model. The mean velocity increased by 49.4% and 143.7%, respectively, at the magnetic intensities of and , respectively.

#### 4. Conclusions

This study investigated numerically the heat and flow characteristics including the control of both suppression and enhancement of the heat transfer of the ferrofluid in the square cavity with variations of the magnetic intensity and temperature difference between heated and cooled walls. The mean and average deviations between the predicted data by developed model and the existing data were 1.0% and −0.77%, respectively, for Barakos at al. [14], 2.0% and −2.0%, respectively, for Markatos and Pericleous [15], at all Rayleigh number ranges from 10^{3} to 10^{6}. Heat and flow characteristics of the ferrofluids within the tested cavity were found to depend on both magnetic intensity and magnetic volume fractions of magnetite. Isotherms of the ferrofluid under the magnetic intensity of expanded more rapidly than those of the ferrofluid under the nonmagnetic intensity at the same time due to the increased magnetic volume force. And the mean Nusselt number and the expansion speed of the isotherm contours of the ferrofluid increased with the rise of the temperature difference between top and bottom walls in the tested cavity. In addition, the mean Nusselt numbers and mean velocity of the ferrofluid were increased with the rise of the magnetic intensities and increased by 23.2% and 143.7%, respectively, at both magnetic intensity of and the elapsed time of seconds.

#### Acknowledgment

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (no. 2013R1A1A1062152), and Business for Cooperative R&D between Industry, Academy, and Research Institute funded Korea Small and Medium Business Administration in 2013 (Grant no. C0102915).

#### References

- H. Yamaguchi, I. Kobori, Y. Uehata, and K. Shimada, “Natural convection of magnetic fluid in a rectangular box,”
*Journal of Magnetism and Magnetic Materials*, vol. 201, no. 1–3, pp. 264–267, 1999. View at Scopus - C. Rinaldi, A. Chaves, S. Elborai, X. He, and M. Zahn, “Magnetic fluid rheology and flows,”
*Current Opinion in Colloid and Interface Science*, vol. 10, no. 3-4, pp. 141–157, 2005. View at Publisher · View at Google Scholar · View at Scopus - H. E. Horng, C.-Y. Hong, S. Y. Yang, and H. C. Yang, “Novel properties and applications in magnetic fluids,”
*Journal of Physics and Chemistry of Solids*, vol. 62, no. 9-10, pp. 1749–1764, 2001. View at Scopus - S. Laohalertdecha, A. S. Dalkilic, and S. Wongwises, “A review on the heat-transfer performance and pressure-drop characteristics of various enhanced tubes,”
*International Journal of Air-Conditioning and Refrigeration*, vol. 20, no. 4, Article ID 1230003, 20 pages, 2012. View at Publisher · View at Google Scholar - A. Kaneko, M. Katsuta, and Y. Hamano, “Research Of heat transfer and pressure drop characteristic of concavo-convex plate,”
*International Journal of Air-Conditioning and Refrigeration*, vol. 21, no. 2, Article ID 1350012, 9 pages, 2013. View at Publisher · View at Google Scholar - Q. W. Wang, M. Zeng, Z. P. Huang, G. Wang, and H. Ozoe, “Numerical investigation of natural convection in an inclined enclosure filled with porous medium under magnetic field,”
*International Journal of Heat and Mass Transfer*, vol. 50, no. 17-18, pp. 3684–3689, 2007. View at Publisher · View at Google Scholar · View at Scopus - T. Basak, S. Roy, A. Matta, and I. Pop, “Analysis of heatlines for natural convection within porous trapezoidal enclosures: effect of uniform and non-uniform heating of bottom wall,”
*International Journal of Heat and Mass Transfer*, vol. 53, no. 25-26, pp. 5947–5961, 2010. View at Publisher · View at Google Scholar · View at Scopus - Y. Zhao, Y. Zhang, R. Lv, and Q. Wang, “Novel optical devices based on the tunable refractive index of magnetic fluid and their characteristics,”
*Journal of Magnetism and Magnetic Materials*, vol. 323, no. 23, pp. 2987–2996, 2011. View at Publisher · View at Google Scholar · View at Scopus - Y. M. Wang, X. Cao, G. H. Liu et al., “Synthesis of Fe3O4 magnetic fluid used for magnetic resonance imaging and hyperthermia,”
*Journal of Magnetism and Magnetic Materials*, vol. 323, no. 23, pp. 2953–2959, 2011. View at Publisher · View at Google Scholar · View at Scopus - L. Jin, X. Zhang, and X. Niu, “Lattice oltzmann simulation for temperature-sensitive magnetic fluids in a porous square cavity,”
*Journal of Magnetism and Magnetic Materials*, vol. 324, no. 1, pp. 44–51, 2012. View at Publisher · View at Google Scholar · View at Scopus - S. Kamiyama, K. Koike, and T. Oyama, “Pipe flow resistance of magnetic fluids in a nonuniform transverse magnetic field,”
*Journal of Magnetism and Magnetic Materials*, vol. 39, no. 1-2, pp. 23–26, 1983. View at Scopus - L. Liang, J. Zhu, and X. Xuan, “Three-dimensional diamagnetic particle deflection in ferrofluid microchannel flows,”
*Biomicrofluidics*, vol. 5, no. 3, Article ID 034110, 2011. View at Publisher · View at Google Scholar · View at Scopus - J. Liu, S.-H. Tan, Y. F. Yap, M. Y. Ng, and N.-T. Nguyen, “Numerical and experimental investigations of the formation process of ferrofluid droplets,”
*Microfluidics and Nanofluidics*, vol. 11, no. 2, pp. 177–187, 2011. View at Publisher · View at Google Scholar · View at Scopus - G. Barakos, E. Mitsoulis, and D. Assimacopoulos, “Natural convection flow in a square cavity revisited: laminar and turbulent models with wall functions,”
*International Journal for Numerical Methods in Fluids*, vol. 18, no. 7, pp. 695–719, 1994. View at Scopus - N. C. Markatos and K. A. Pericleous, “Laminar and turbulent natural convection in an enclosed cavity,”
*International Journal of Heat and Mass Transfer*, vol. 27, no. 5, pp. 755–772, 1984. View at Scopus - T. Sawada, H. Kikura, and T. Tanahashi, “Visualization of wall temperature distribution caused by natural convection of magnetic fluids in a cubic enclosure,”
*International journal of Applied Electromagnetics in Materials*, vol. 4, no. 4, pp. 329–335, 1994. View at Scopus - S. M. Snyder, T. Cader, and B. A. Finlayson, “Finite element model of magnetoconvection of a ferrofluid,”
*Journal of Magnetism and Magnetic Materials*, vol. 262, no. 2, pp. 269–279, 2003. View at Publisher · View at Google Scholar · View at Scopus - S. H. Jin, J. W. Park, and L. S. Seo, “Numerical analysis on the characteristic benard flow and control of magnetic fluid,” in
*Proceedings of the SAREK Winter Conference*, pp. 280–285, 2001. - J. W. Park, S. H. Jin, C. H. Jun, L. S. Seo, and Y. S. Kim, “Numerical analysis on the heat transfer characteristics of Benard flow in a magnetic fluids,”
*Journal of the Korean Magnetics Society*, vol. 13, no. 1, pp. 41–46, 2003. - J. W. Park, C. H. Jun, S. O. Ryu, and L. S. Seo, “Characteristics of heat transfer of natural convection for magnetic fluids in annular pipes,”
*Journal of the Korean Magnetics Society*, vol. 12, no. 2, pp. 73–79, 2002. - H. Okanaga, K. Shizawa, N. Yashima, and T. Tanahashi, “Numerical analysis of magnetic fluid flow in a square cavity (GSMAC finit element method of magnetic fluid),”
*Transactions of the Japan Society of Mechanical Engineers B*, vol. 53, no. 492, pp. 2414–2422, 1987. View at Scopus - Y. Oki and T. Tanahashi, “Numerical analysis of natural convection of thermo-electrically conducting fluids in a square cavity under a constant magnetic field (1st Report, estimation of induced heating term),”
*Transactions of the Japan Society of Mechanical Engineers B*, vol. 59, no. 562, pp. 1835–1842, 1993. View at Scopus - J. W. Park,
*A study on the heat transfer characteristics and fluid flow of magnetic fluids in concentric double pipe annuli [Ph.D. thesis]*, Dong-A University, 2002. - H. Kawai, Y. Katoh, T. Sawada, and T. Tanahashi, “GSMAC-FEM for incompressible viscous flow analysis (a modified GSMAC method),”
*The Japan Society of Mechanical Engineers*, vol. 33, no. 1, pp. 17–25, 1990. View at Scopus - P. Vinogradova and A. Zarubin, “A study of Galerkin method for the heat convection equations,”
*Applied Mathematics and Computation*, vol. 218, no. 2, pp. 520–531, 2011. View at Publisher · View at Google Scholar · View at Scopus - J. H. Seo,
*A study on the natural convection of magnetic fluid in a cubic cavity [M.S. thesis]*, Dong-A University, 2011. - T. Grosan, C. Revnic, I. Pop, and D. B. Ingham, “Magnetic field and internal heat generation effects on the free convection in a rectangular cavity filled with a porous medium,”
*International Journal of Heat and Mass Transfer*, vol. 52, no. 5-6, pp. 1525–1533, 2009. View at Publisher · View at Google Scholar · View at Scopus - A. Lange, “Magnetic Soret effect: application of the ferrofluid dynamics theory,”
*Physical Review E*, vol. 70, no. 4, Article ID 046308, pp. 1–46308, 2004. View at Publisher · View at Google Scholar · View at Scopus