A numerical study on heat transfer of a ferrofluid flow in a square cavity under simultaneous gravitational and magnetic convection

Abstract

Thermomagnetic convection is based on the use of external magnetic fields to better control heat transfer fluxes in ferrofluids, finding important applications in engineering and many related areas. The improvement of such methods relies on fundamentally understanding the flow of ferrofluids under temperature gradients and external magnetic fields. However, the underlying physics of this phenomenon is very complex and not yet well characterized. The problem we analyze in this paper consists of a ferrofluid confined in a square cavity heated from the top and subjected to an external magnetic field applied in the horizontal direction. Differently from earlier investigations, this scenario leads to a clear competition between stabilizing gravitational forces and destabilizing magnetic forces; the unbalance between them drives the flow and dictates the mechanisms of heat transfer in the system. The problem is described by the equations of conservation of mass, momentum, and energy; both gravitational and magnetic effects are accounted for through the definition of the body force following the well-known Boussinesq approximation to express the ferrofluid’s density and magnetization as functions of temperature. The resulting set of fully coupled, nonlinear equations of the model is solved with a fully implicit finite element method. The results show that thermomagnetic convection increases convective heat transfer fluxes in the flow, but its net effects essentially depend on a complex balance among viscous, magnetic, and gravitational forces which determines the pattern of recirculation regions inside the cavity. There are two critical values associated with the external field’s intensity: the first marks the onset of thermomagnetic convection when the destabilizing magnetic effects become comparable to the stabilizing gravitational ones and the second corresponds to the external field’s strength that suffices to make the effects of gravity nearly negligible in the flow dynamics. The latter regime is dictated by a balance between viscous and magnetic forces only, and the corresponding numerical predictions agree notoriously well with a scaling analysis which suggests that $${\overline{\mathrm{Nu}}} \sim \mathrm{Ra}_\mathrm{m}^{1/4}$$, where $${\overline{\mathrm{Nu}}}$$ is the average Nusselt number at the isothermal walls and $$\mathrm{Ra}_\mathrm{m}$$ is the magnetic Rayleigh number, a dimensionless parameter associated with the intensity of the external field.