Next Article in Journal
Editorial for the Special Issue on MEMS/NEMS Sensors: Fabrication and Application
Next Article in Special Issue
Electrical Coupling and Simulation of Monolithic 3D Logic Circuits and Static Random Access Memory
Previous Article in Journal
A Review: Preparation, Performance, and Applications of Silicon Oxynitride Film
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermomagnetic Convection of Ferrofluid in an Enclosure Channel with an Internal Magnetic Field

1
Graduate School of Mechanical Engineering, Sungkyunkwan University, Suwon 16419, Korea
2
School of Mechanical Engineering, Sungkyunkwan University, Suwon 16419, Korea
*
Author to whom correspondence should be addressed.
Micromachines 2019, 10(9), 553; https://doi.org/10.3390/mi10090553
Submission received: 17 July 2019 / Revised: 6 August 2019 / Accepted: 16 August 2019 / Published: 21 August 2019
(This article belongs to the Special Issue NANO KOREA 2019)

Abstract

:
Ferrofluid is a colloidal liquid in which magnetic nanoparticles such as Fe3O4 are dispersed in a nonconductive solution, and the average diameter of the nanoparticles is 10 nm. When a magnetic field is applied, the ferrofluid generates magnetization, which changes the physical properties of the fluid itself. In this study, characteristics of the thermomagnetic convection of ferrofluid (Fe3O4) by the permanent magnet in the enclosure channel were studied. To effectively mix the ferrofluid (Fe3O4) and disturb the boundary layer, the heat dissipation of the heat source depending on the strength of the magnetic field and the shape of the enclosure channel was numerically studied. In particular, four different enclosure channels were considered: Square, separated square, circle, and separated circle. The hot temperature was set at the center of the enclosure channel. The ferrofluid was affected by the permanent magnet in the center of the channel. The magnetic field strength in the region close to the permanent magnet was enhanced. The magnetophoretic (MAP) force increased with increasing magnetic field strength. The MAP force generated a vortex in the enclosure channel, disturbing the thermal boundary. The vortex occurs differently, depending on the shape of the enclosure channel and affects the thermomagnetic convection. The temperature and velocity fields for thermomagnetic convection were described and the convective heat flux was calculated and compared. Results show that when the magnetic field strength was 4000 kA/m and the shape of the enclosure channel was a circle, the maximum convective heat flux of 4.86 × 105 W/m2 was obtained.

1. Introduction

With advances in technology, electronic devices are getting smaller and the amount of power used is increasing. The increased power generates heat in the electronics, which reduces the performance of the devices. For this reason, the importance of thermal management in the industrial field is emphasized. Recently, ferrofluid has been highlighted with cooling technology. It has higher thermal conductivity than conventional heat transfer fluid and has the advantage of controlling the fluid by using the magnetic field. Ferrofluid is a material that contains nanoparticles with nonconductive liquid (water, oil, etc.) as a solvent. The nanoparticles in the ferrofluid are Fe3O4, Al2O3, FePt, CoPt, etc. and the average diameter of the particles is 10 nm. The ferrofluid is entirely superparamagnetic, so that the flow of ferrofluid can be controlled by the magnetic field. A magneto-hydrodynamic (MHD) method using a magnetic field is used primarily for flow control of ferrofluid, in which a magnetic force moves magnetic particles when a magnetic field is externally applied to a ferrofluid [1,2]. The flow of ferrofluid can be controlled by magnetophoresis and the temperature gradient changes, depending on the relationship between the intensity of the applied magnetic field and the flow of the ferrofluid. For these reasons, many researchers of engineering fields proposed new design parameters to various applications using fluid element [3,4,5]. Blennerhassett et al. [6] found that the ratio of convective and conduction heat transfer can be up to 10% greater than in the absence of a magnetic field through nonlinear analysis. Strek et al. [7] numerically studied the effect of a magnetic dipole located under a 2D channel on the ferrofluid in the channel. They observed that ferrofluid flows in the direction of the magnetic field at cold temperature. Aminfar et al. [8] numerically investigated the hydro-thermal characteristics of ferrofluid (water and Fe3O4) exposed to a non-uniform magnetic field, caused by the current passing through a wire placed in parallel under the duct. It was shown that heat transfer is improved by the increase of Nusselt number and friction factor when the magnetic field is applied. Song et al. [9] numerically studied the effects of permanent magnet position and intensity on thermomagnetic convection of oxygen in a rectangular enclosure. The thermomagnetic convection was largely influenced by the relative position of the permanent magnet. Tangthieng et al. [10] studied finite element simulations of heat transfer to a ferrofluid in a box in the presence of a magnetic field. Simulations showed improved heat transfer due to magnetic field gradient. Yamaguchi et al. [11] conducted experiments and numerical studies on the natural convection of ferrofluids in a 2D cell with an aspect ratio of 1. Vertically imposed magnetic fields have unstable effects and have shown that the flow is different from a state without magnetic fields. Furthermore, it has been shown that the magnetic field causes transition flow to improve heat transfer. Ghasemi et al. [12] studied the natural convection of enclosures filled with water–Al2O3 nanoparticles affected by magnetic fields. It was shown that the heat transfer rate was increased with increasing Rayleigh number but decreased with increasing Hartmann number. Jiang et al. [13] investigated numerically the thermomagnetic convection of air in a 2D square enclosure under the magnetic field on four sides. They presented the flow and temperature distributions for the thermomagnetic convection of air. The local and average Nusselt numbers on the walls were calculated and compared. The results showed that the magnetic field strength, Darcy number, and Rayleigh number affect the heat transfer in the square enclosure. Szabo et al. [14] experimentally analyzed thermodynamic convection using infrared thermography and numerically analyzed it with the finite element model. The temperature dependence of the Kelvin body force was explained by relating the heat flux and temperature profile.
Many researchers studied the thermomagnetic convection characteristics of ferrofluid in a square enclosure when applied to an external magnetic field. In this study, the thermomagnetic convection characteristic of ferrofluid (Fe3O4), depending on the enclosure channel shape (square, separated square, circle, separated circle) with the internal magnetic field, was studied using COMSOL Multiphysics 5.3a (COMSOL Inc, Burlington, MA, USA). Temperature and velocity distributions were graphically depicted with various geometrical configurations and operating conditions.

2. Governing Equations

Since there is no electromagnetic free current, ferrofluid can be assumed to be a nonconductor. It was assumed that there is no magnetic field effect due to the temperature gradient. Therefore, the Maxwell–Amperes law in the magnetostatics condition can be expressed as follows:
· B = 0
× H = 0
where B   is the magnetic flux density vector and H is the magnetic field vector.
H = 1 μ 0 1 + χ m B
where χ m   is the magnetic susceptibility and μ 0   is the magnetic permeability in vacuum.
The viscous torque of nanoparticles can be ignored. Ferrofluids need not consider magneto- dissipation [15]. To describe the flow of ferrofluid, the following governing equations are used:
Continuity:
ρ t + · ρ v = 0
Momentum:
ρ v t + v · v = p + · τ i j + β ρ g T T 0 k ^ + M · B
Energy:
ρ C p T t + v · T = k 2 T + η Φ μ 0 T M T v · H
where ρ is the density of ferrofluid, v is the velocity vector, p is the pressure, T is the temperature of ferrofluid, τ i j is the viscous stress tensor, k ^ is the unit vector of gravity force, β is the thermal expansion coefficient of the ferrofluid, k is the thermal conductivity, η is the viscosity of ferrofluid, and η Φ is the viscous dissipation term in the energy equation.
The magnetic induction vector can be expressed as follows:
B = μ 0 1 + χ m H
The variation of the total magnetic susceptibility χ m is treated as being dependent on temperature [16].
χ m = χ m T = χ m , 0 1 + β T T 0

3. Numerical Analysis

The numerical models of this study are shown in Figure 1. The channels consisted of squares and circles with the same area. The size of the enclosure channel was set based on the conventional Central Processing Unit (CPU) chip. The polydimethylsiloxane (PDMS) was used as the material of the body and ferrofluid (Fe3O4) was used as the working fluid. The hot temperature (Th) in the center of the channel was 313.15 K. To apply a magnetic field, a permanent magnet was placed at the center of the channel. The magnetic field strength was changed from H = 1000 to 4000 kA/m to analyze the effect of the magnetic field on the ferrofluid in the enclosure channel. The grid applied to the channel is shown in Figure 2. The grid system was set up with dense free triangular for the whole model and an inflation layer was set in the vicinity of wall to check the change according to the channel. Grid dependency was tested by the magnetic flux density in the center of the channel. All cases set the mesh to 3 × 104 (see Figure 3). Detailed boundary conditions are given in Table 1 and Table 2.

4. Results & Discussion

The applied magnetic field has a different magnetic energy distribution, depending on the shape of the channel. This affects the heat transfer characteristics of the ferrofluid in the channel. To verify the heat transfer characteristics of the ferrofluid in the channel, the influence of the magnetophoretic force due to the magnetic field was analyzed and the temperature distribution was analyzed. When the magnetic field strength is 2000 kA/m, the magnetophoretic force distributions according to the geometry are shown in Figure 4. The MAP force can be expressed as follows [17]:
F MAP = μ m M e f f · H
F MAP = 2 π r p 3 μ m ReK μ p , μ m H 2
where M e f f is the magnetic dipole moment by external force, r p is the radius of magnetic nanoparticle, and R e K μ p , μ m is the Clausius–Mossoti factor using magnetic permeability.
The MAP force represented the behavior of the nanoparticle, due to its interaction with the magnetic field gradient. The MAP force decreased as the distance from the permanent magnet (center of the channel) increased and depended on the shape of the enclosure channel. When the magnetic field strength is H = 2000 kA/m, the streamlines and temperature distributions, according to the geometry, are shown in Figure 5. The temperature distribution was consistent with the streamline. The vortex was responsible for spreading the heat from the center outside. The non-separated shape showed a high temperature in a particular area and the vortex influenced the edge of the enclosure channel when the shape of the channel was circular. As the magnetic field intensity increased, a large vortex was generated in a circular channel, and the central heat shifted to the outside (see Figure 6). Figure 7 shows the convective heat flux for all cases. The heat in the center of the enclosure channel was spread by the MAP force. Therefore, the convective heat flux increased as the magnetic field strength increased. The highest convective heat flux was obtained when the channel shape was a circle. The average temperature of the channel according to the convective heat flux is shown in Figure 8. As the convective heat flux increased, the average temperature increased. However, when the convective heat flux exceeded 4 W/m2, the increasing slope of the average temperature value tended to decrease.

5. Conclusions

In this study, the thermomagnetic convection of the ferrofluid depending on the internal magnetic field was studied numerically by changing the shape of the enclosure channel. When the magnetic field was applied through the permanent magnet, the MAP force was found to occur strongly near the permanent magnet. The flow distribution of the ferrofluid formed a vortex field around the position where the MAP force was strongly generated. The MAP force plays an important role in promoting heat transfer by forming the recirculation region. The temperature distribution corresponds with the formed vortex. As the MAP force increased, the size of the vortex increased. According to the shape of the enclosure channel, a different vortex occurred in the same MAP force. Convective heat flux was larger in the circle- than in the square-shaped channel in the same area and external condition. When the magnetic field strength was 4000 kA/m and the shape of the enclosure channel was a circle, the maximum convective heat flux of 4.86 × 105 W/m2 was obtained. Through study, it was confirmed that the enclosure channel shape and magnetic field affect the thermodynamic convection characteristics of the ferrofluid. It is possible to embody the effective heat transfer performance through the magnetic field and channel design.

Author Contributions

Supervision, Y.-J.K.; Writing original draft, M.L.; Writing review & editing, M.L. and Y.-J.K.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Ministry of Education, Science and Technology (2017R1D1A1B03036263).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rosensweig, R.E. Ferrohydrodynamics; Cambridge University Press: Cambridge, UK, 1985. [Google Scholar]
  2. Yu, P.X.; Qiu, J.X.; Qin, Q.; Tian, Z.F. Numerical investigation of natural convection in a rectangular cavity under different directions of uniform magnetic field. Int. J. Heat and Mass Transf. 2013, 67, 1131–1144. [Google Scholar] [CrossRef]
  3. Bibik, E.E.; Matgullin, B.Y.; Raikher, Y.L.; Shliomis, M.I. The magnetostatic properties of colloidal magnetite. Magnetohydrodynamics 1973, 9, 58–62. [Google Scholar]
  4. Janocha, H. Application potential of magnetic field driven new actuators. Sens. Actuators A 2001, 91, 126–132. [Google Scholar] [CrossRef]
  5. Afifah, A.N.; Syahrullail, S.; Sidik, N.A.C. Magnetoviscous effect and thermomagnetic convection of magnetic fluid: A review. Renew. Sustain. Energy Rev. 2016, 55, 1030–1040. [Google Scholar] [CrossRef]
  6. Blennerhassett, P.J.; Lin, F.; Stiles, P.J. Heat Transfer through strongly magnetized ferrofluids. Proc. R. Soc. Lond. A 1991, 433, 165–177. [Google Scholar] [CrossRef]
  7. Strek, T.; Jopek, H. Computer simulation of heat transfer through a ferrofluid. Phys. Status Solidi B 2007, 244, 1027–1037. [Google Scholar] [CrossRef]
  8. Aminfar, H.; Mohammadpourfard, M.; Zonouzi, S.A. Numerical study of the ferrofluid flow and heat transfer through a rectangular duct in the presence of a non-uniform transverse magnetic field. J. Magn. Magn. Mat. 2013, 327, 31–42. [Google Scholar] [CrossRef]
  9. Song, K.W.; Tagawa, T. Thermomagnetic convection of oxygen in a square enclosure under non-uniform magnetic field. Int. J. Therm. Sci. 2018, 125, 52–65. [Google Scholar] [CrossRef]
  10. Tangthieng, C.; Finlayson, B.A.; Maulbetsch, J.; Cader, T. Heat transfer enhancement in ferrofluids subjected to steady magnetic fields. J. Magn. Magn. Mat. 1999, 201, 252–255. [Google Scholar] [CrossRef]
  11. Yamaguchi, H.; Zhang, Z.; Shuchi, S.; Shimada, K. Heat transfer characteristics of magnetic fluid in a partitioned rectangular box. J. Magn. Magn. Mat. 2002, 252, 203–205. [Google Scholar] [CrossRef]
  12. Ghasemi, B.; Aminossadati, S.M.; Raisi, A. Magnetic field effect on natural convection in a nanofluid-filled square enclosure. Int. J. Therm. Sci. 2011, 50, 1748–1756. [Google Scholar] [CrossRef]
  13. Jiang, C.; Feng, W.; Zhong, H.; Zeng, J.Y.; Zhu, Q.M. Effects of a magnetic quadrupole field on thermomagnetic convection of air in a porous square enclosure. J. Magn. Magn. Mat. 2014, 357, 53–60. [Google Scholar] [CrossRef]
  14. Szabo, P.S.B.; Bekovic, M.; Früh, W.G. Using infrared thermography to investigate thermomagnetic convection under spatial non-uniform magnetic field Peter. Int. J. Therm. Sci. 2017, 116, 118–128. [Google Scholar] [CrossRef]
  15. Muller, H.W.; Engel, A. Dissipation in ferrofluid: Mecroscopic versus hydrodynamic theory. Phys. Rev. E. 1999, 60, 7001–7009. [Google Scholar] [CrossRef] [PubMed]
  16. Ganguly, R.; Sen, S.; Puri, I.K. Heat transfer augmentation using a magnetic fluid under the influence of a line dipole. J. Magn. Magn. Mat. 2004, 271, 63–73. [Google Scholar] [CrossRef]
  17. James, C.D.; McClain, J.; Pohl, K.R.; Reuel, N.; Achyuthan, K.E.; Bourdon, C.J.; Rahimian, K.; Galambos, P.C.; Ludwig, G.; Derzon, M.S. High-efficiency magnetic particle focusing using dielectrophoresis and magnetophoresis in a microfluidic device. J. Micromech. Microeng. 2010, 20, 04515. [Google Scholar] [CrossRef]
Figure 1. Configurations for numerical analysis.
Figure 1. Configurations for numerical analysis.
Micromachines 10 00553 g001aMicromachines 10 00553 g001b
Figure 2. Schematic of the grid systems.
Figure 2. Schematic of the grid systems.
Micromachines 10 00553 g002
Figure 3. Grid dependency test.
Figure 3. Grid dependency test.
Micromachines 10 00553 g003
Figure 4. Magnetophoretic force distribution with the same magnetic field strength, H = 3000 (kA/m).
Figure 4. Magnetophoretic force distribution with the same magnetic field strength, H = 3000 (kA/m).
Micromachines 10 00553 g004
Figure 5. Streamlines and temperature distributions with the same magnetic field strength, H = 2000 (kA/m).
Figure 5. Streamlines and temperature distributions with the same magnetic field strength, H = 2000 (kA/m).
Micromachines 10 00553 g005
Figure 6. Streamlines and temperature distribution in a circle channel with the different magnetic field strength.
Figure 6. Streamlines and temperature distribution in a circle channel with the different magnetic field strength.
Micromachines 10 00553 g006
Figure 7. Convective heat flux in all cases.
Figure 7. Convective heat flux in all cases.
Micromachines 10 00553 g007
Figure 8. Average temperature due to convective heat flux in all cases.
Figure 8. Average temperature due to convective heat flux in all cases.
Micromachines 10 00553 g008
Table 1. Physical properties of the ferrofluid (EFH-1).
Table 1. Physical properties of the ferrofluid (EFH-1).
Density1221 kg/m³
Relative permeability of ferrofluid1.01
Relative permeability of magnetite (Fe3O4)16
Magnetic susceptibility1.552
Dynamic viscosity0.00727 Pa∙s
Thermal conductivity0.19W/m∙K
Heat capacity at static pressure1840 J/kg∙K
Thermal expansion coefficient8.6 × 10−4 (1/K)
Table 2. Boundary conditions applied in this study.
Table 2. Boundary conditions applied in this study.
Working fluidFerrofluid
Geometric configurationsSquare, Separated square, Circle, Separated circle
Magnetic field strength, H1000–4000 (kA/m)
Directions of magnetic energyHorizontal
Initial temperature293.15 K
Hot temperature (Heat source), Th313.15 K
Wall conditions of fluid domainNo-slip condition

Share and Cite

MDPI and ACS Style

Lee, M.; Kim, Y.-J. Thermomagnetic Convection of Ferrofluid in an Enclosure Channel with an Internal Magnetic Field. Micromachines 2019, 10, 553. https://doi.org/10.3390/mi10090553

AMA Style

Lee M, Kim Y-J. Thermomagnetic Convection of Ferrofluid in an Enclosure Channel with an Internal Magnetic Field. Micromachines. 2019; 10(9):553. https://doi.org/10.3390/mi10090553

Chicago/Turabian Style

Lee, Myoungwoo, and Youn-Jea Kim. 2019. "Thermomagnetic Convection of Ferrofluid in an Enclosure Channel with an Internal Magnetic Field" Micromachines 10, no. 9: 553. https://doi.org/10.3390/mi10090553

APA Style

Lee, M., & Kim, Y. -J. (2019). Thermomagnetic Convection of Ferrofluid in an Enclosure Channel with an Internal Magnetic Field. Micromachines, 10(9), 553. https://doi.org/10.3390/mi10090553

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop