Numerical Study of the Flow and Thermomagnetic Convection Heat Transfer of a Power Law Non-Newtonian Ferrofluid within a Circular Cavity with a Permanent Magnet

: The aim of this study is to analyze the thermo-magnetic-gravitational convection of a non-Newtonian power law ferroﬂuid within a circular cavity. The ferroﬂuid is exposed to the magnetic ﬁeld of a permanent magnet. The ﬁnite element method is employed to solve the non-dimensional controlling equations. A grid sensitivity analysis and the validation of the used method are conducted. The effect of alterable parameters, including the power law index, 0.7 ≤ n ≤ 1.3, gravitational Rayleigh number, 10 4 ≤ Ra T ≤ 10 6 , magnetic Rayleigh number, 10 5 ≤ Ra M ≤ 10 8 , the location of the hot and cold surfaces, 0 ≤ λ ≤ π /2, and the length of the magnet normalized with respect to the diameter of the cavity, 0.1 ≤ L ≤ 0.65, on the ﬂow and heat transfer characteristics are explored. The results show that the heat transfer rate increases at the end of both arcs compared to the central region because of buoyancy effects, and it is greater close to the hot arc. The location of the arcs does not affect the heat transfer rate considerably. An increase in the magnetic Rayleigh number contributes to stronger circulation of the ﬂow inside and higher heat transfer. When the Kelvin force is the only one imposed on the ﬂow, it enhances the heat transfer for magnets of length 0.2 ≤ L ≤ 0.3.


Introduction
Nanotechnology and nanoengineering have piqued the curiosity of many scientists and engineers in recent decades.The use of ferrofluids is one of the most effective and crucial ways to transfer heat in this new area of nanoscience, as such fluids have superior thermal conductivity compared to regular fluids [1].
Forming a ferrofluid, also known as magnetic fluid, by dispersing ferroparticles in a non-magnetic base fluid (typically water, kerosene, or oil) can significantly enhance the thermal and rheological properties when subjected to a magnetic field [2][3][4].The thermophysical properties are very sensitive to the magnetic field strength and direction [1].
The use of an external magnetic field to manage heat transfer fluxes in ferrofluids is known as thermomagnetic convection.This kind of convection can be significantly greater than natural convection and may be leveraged to improve the performance of cooling systems [1] The application of a magnetic field may either enhance or disrupt the fluid flow and heat transfer, resulting in a continuous convection of heat in a magnetized ferrofluid that has an uneven temperature distribution.The adequate use of magnetic convection in ferrofluids, where the magnetization changes in the ferrofluid due to temperature gradients, may increase the performance of small-scale cooling mechanisms [5].
With the development of better magnetic coils and permanent magnets, much effort has been made to approach effective cooling based on ferrofluid convection.Nakatsuka et al. [6] studied the heat transfer properties of a heat pipe containing a water-based magnetic fluid and found that a non-uniform magnetic field increased heat transfer efficiency by 13%.Kumar et al. [7] investigated the influence of the intensity of a self-induced magnetic field surrounding a current-carrying wire on thermomagnetic convection cooling in ferrofluid using both experimental and computational methods.This comparative study showed that the cooling was due to the self-induced magnetic field interacting with the magnetic fluid, rather than natural convection or other nanofluid-related mechanisms.
Combined natural and thermomagnetic convection in cavities, on the other hand, has gained a lot of attention from researchers because of the many potential applications, making it advantageous compared to natural convection alone [8][9][10].
The natural convection heat transfer properties of a diester-based ferrofluid under a permanent magnet effect in a channel with a constant cross-section, in which the bottom wall is heated and the top wall cooled, were studied by Cheng et al. [11].The measurement of the effect of cooling temperatures on the thermal behavior of ferrofluids was done in a simplified model of a magnetic fluid seal.The permanent magnetic field and temperature gradient influence the ferrofluid's convective heat transfer.As a result of the temperature gradient, both natural convection and thermomagnetic convection are generated.The heterogeneous distribution of the Kelvin force is caused by the non-uniform magnetic field, which, in turn, causes thermomagnetic convection.
Khedher et al. [12] numerically analyzed the behavior of a magnetic nanoliquid under the influence of Brownian diffusion inside a cavity with two local heaters and two magnetic sources.The findings demonstrated that an increase in the upper magnetic source power reflects an increase in the thermal diffusive zone in the upper half of the cavity as the ratio of the strengths of the magnetic sources rises.In a revolving enclosure with magnetic fins, the rate of heat conduction through a cooling fin is significantly influenced by the Kelvin body force produced on oxygen atoms by the fin [13].The Nusselt number increases by 5% in a magnetic field with a residual magnetic flux density of 5.0 T.However, gravitational convection dominates in a ferrofluid in the gap between two concentric cylinders, where the magnetic field formed by a current-carrying conductor is found to be inadequate to compete with gravitational convection [14].
Recently, the thermomagnetic cooling performance of many magnetic cooling devices has been examined.Pattanaik et al. [15] designed a novel magnetic cooling device for long-distance heat transfer.Experiments and numerical models were used to evaluate heat load cooling.
Wu et al. [16] used an external magnetic field to induce the sliding movement of ferric oxide nanoparticles in an aqueous lithium bromide solution in order to increase the vapor absorption rate.The experimental findings suggest that employing this strategy increases the vapor absorption rate significantly.The average absorption rates attained with the magnetic field effect were 1.58 times greater than those achieved without it.The findings suggest that using a nanofluid may improve the vapor absorption process.Particle image velocimetry data were used by Lian et al. [17] to quantify the velocity of a ferrofluid in a flow channel.
The development of novel cooling systems based on thermomagnetic convection demonstrates the relevance of prior research in this field in achieving this goal.However, further research is needed to better understand this issue and to maximize both the thermomagnetic and natural convections inside cooling systems of various shapes and configurations.This study aims to investigate the effects of the presence of a permanent magnet on the magnetic-gravitational heat transfer convection of a non-Newtonian ferrofluid in a cylindrical cavity.The goal of this research is to see how the non-Newtonian nature of a ferrofluid and the magnet dimension influence the flow and thermal fields, as well as the rate of heat transfer.The geometry used in this study is suggested as a way to examine and enhance the advantages of the magneto ferrofluid.Finally, in the studied geometry, the effects of the location and arc size of the hot and cold walls in the presence of the permanent magnet are investigated.

Problem Definition
This work deals with the convection flow of a power law non-Newtonian fluid which is influenced by the buoyancy forces resulting from a permanent magnet and gravity.Figure 1 depicts the computational domain along with the defined Cartesian coordinate system.The domain is a circular enclosed medium having a diameter, d.The higher and lower temperatures T h and T c are imposed on the left and right arcs of the circular enclosure with a length of G = δd/2.The diameter passing through the middle of the cold sector forms an angle λ with respect to the defined x-axis.A permanent magnet of length l applies the magnetic field of strength H d to the ferrofluid.The surfaces of the magnet, located in the center of the enclosure, are adiabatic.The ferrofluid flow is considered to be steady, incompressible and laminar.The electronic current due to the movement of the magnetic flow may be neglected.Density changes are ignored except in the buoyancy term, which is estimated by the Boussinesq linear approximation.The development of novel cooling systems based on thermomagnetic convection demonstrates the relevance of prior research in this field in achieving this goal.However, further research is needed to better understand this issue and to maximize both the thermomagnetic and natural convections inside cooling systems of various shapes and configurations.This study aims to investigate the effects of the presence of a permanent magnet on the magnetic-gravitational heat transfer convection of a non-Newtonian ferrofluid in a cylindrical cavity.The goal of this research is to see how the non-Newtonian nature of a ferrofluid and the magnet dimension influence the flow and thermal fields, as well as the rate of heat transfer.The geometry used in this study is suggested as a way to examine and enhance the advantages of the magneto ferrofluid.Finally, in the studied geometry, the effects of the location and arc size of the hot and cold walls in the presence of the permanent magnet are investigated.

Problem Definition
This work deals with the convection flow of a power law non-Newtonian fluid which is influenced by the buoyancy forces resulting from a permanent magnet and gravity.Figure 1 depicts the computational domain along with the defined Cartesian coordinate system.The domain is a circular enclosed medium having a diameter, .The higher and lower temperatures  h and  c are imposed on the left and right arcs of the circular enclosure with a length of  =  2 ⁄ .The diameter passing through the middle of the cold sector forms an angle  with respect to the defined  -axis.A permanent magnet of length  applies the magnetic field of strength  d to the ferrofluid.The surfaces of the magnet, located in the center of the enclosure, are adiabatic.The ferrofluid flow is considered to be steady, incompressible and laminar.The electronic current due to the movement of the magnetic flow may be neglected.Density changes are ignored except in the buoyancy term, which is estimated by the Boussinesq linear approximation.The variables listed in Equation ( 1) are used to normalize the controlling equations, where  is the component of velocity in the -direction,  is the component of velocity in the -direction,  is thermal diffusivity,  is pressure and  is the density of the fluid. is magnetic field strength, B is magnetic flux density and  is magnetization (magnetic moment per unit volume) of the fluid.The subscripts n, d, r and s mean dimensionless, with dimensions, remanent and saturation, respectively. is the permeability of free space.The variables listed in Equation ( 1) are used to normalize the controlling equations, where u is the component of velocity in the x-direction, v is the component of velocity in the y-direction, α is thermal diffusivity, p is pressure and ρ is the density of the fluid.H is magnetic field strength, B is magnetic flux density and M is magnetization (magnetic moment per unit volume) of the fluid.The subscripts n, d, r and s mean dimensionless, with dimensions, remanent and saturation, respectively.µ 0 is the permeability of free space.
Applying the dimensionless variables shown in Equation (1), i.e., P, H n , H X , H Y and M n , the controlling equations in dimensionless coordinates may be expressed as in Equations ( 2)-(4).Equations ( 5)-( 8) define the variables .γ and the Prandtl (Pr), Rayleigh (Ra), Mach (Ma) and temperature (ε) numbers for the magnetic ferrofluid.K is the thermal conductivity of the fluid, which has a specific heat capacity at constant pressure of C p .Other variables are defined in subsequent equations.
The components of the magnetic field in dimensionless coordinates are given in the X-direction by Equation ( 9) and in the Y-direction by Equation (10), where L = l/d in both cases.The dimensionless form of the magnetic field properties is given by the expressions in Equation (11), where ε is the temperature number of the fluid defined in Equation (8) and N is the number of magnetic dipoles per unit volume.

Numerical Method, Grid Test, and Verification
The program COMSOL Multiphysics 5.6 is used to solve the governing equations and the boundary conditions described in Section 2. To examine the dependence or otherwise of the results on the grid (mesh) size in the numerical simulation, a mesh sensitivity test is performed, in which several cases of different grid element densities are compared.In this study, the average Nusselt number and the maximum strength of the fluid inside the cavity (for two cases of Newtonian fluid, where n = 1 by definition, and a non-Newtonian fluid in which n = 0.7) are considered due to the importance of the heat transfer rate and power of fluid circulation.The parameters used for this test are Ra T = 10 5 , Ra M = 10 7 , Pr = 100, Ma = 10 −4 , β c = 60, L = 0.3, λ = 0 and ε = 0.5.The simulation grids tested contained 3851, 9515, 23723 and 41242 elements.The grid with 23723 elements is the coarsest mesh that gives results that do not change with further refinement (Figure 2).
the cavity (for two cases of Newtonian fluid, where  = 1 by de tonian fluid in which  = 0.7) are considered due to the import rate and power of fluid circulation.The parameters used for this t 10 , Pr = 100, Ma = 10 ,  c = 60,  = 0.3,  = 0 and  = 0. tested contained 3851, 9515, 23723 and 41242 elements.The grid the coarsest mesh that gives results that do not change with furth In order to ensure that the solution is accurate, the results o compared with the those reported in published studies [18,19].the streamlines and isotherms of a Newtonian ferrofluid under t magnet are compared using the conditions Ra T = 0, Ra M = 1  = 1.0,  = 0.2,  c = 60 and  = 5.0.The current results agree w previously (Figure 3a,b) [18].There is also good agreement betw sults and the equivalent values calculated in this study (  In order to ensure that the solution is accurate, the results of the present study were compared with the those reported in published studies [18,19].In the first comparison, the streamlines and isotherms of a Newtonian ferrofluid under the effect of a permanent magnet are compared using the conditions Ra T = 0, Ra M = 10 9 , Ma = 0, Pr = 1000, n = 1.0, l = 0.2, β c = 60 and ε = 5.0.The current results agree well with those published previously (Figure 3a,b) [18].There is also good agreement between other published results and the equivalent values calculated in this study (Figure 3c) [19].Furthermore, to verify the outcomes for a non-Newtonian fluid, the average Nusselt number under the buoyancy-driven flow of a non-Newtonian fluid in a square enclosure is compared to similar published results [20].There is a good match between the results for different values of the power law index (Table 1).buoyancy-driven flow of a non-Newtonian fluid in a square enclosure is compared to similar published results [20].There is a good match between the results for different values of the power law index (Table 1). (a)

Results and Discussion
The following parameters were studied in the ranges specified: 10 4 ≤ Ra T ≤ 10 6 , 10 5 ≤ Ra M ≤ 10 8 , 0.7 ≤ n ≤ 1.3, 0 ≤ λ ≤ π/2, 0.1 ≤ L ≤ 0.65.The other parameters, which were kept fixed throughout the study, were Pr = 100, Ma = 10 −4 , β c = 60, and ε = 0.2. Figure 4 shows the variation of the streamlines for power law indices, n = 0.7, 1.0 and 1.3 and of the angle λ when Ra T = 10 4 , Ra M = 10 7 , Pr = 100, Ma = 0, β c = 60, δ = 30 • L = 0.3, and ε = 0.2.All subfigures show two vorticities formed inside the circular cavity.For λ = 0, it is seen that the fluid velocity is greater for fluids with a lower power law index, n (Figure 4).This is attributed to the fact that the viscosity of a fluid with a lower power law index is lower for a given state of flow, which allows the ferrofluid to flow faster under given conditions.Also, the center of the vorticities is closer to the cold wall, showing the tendency of the movement of the ferrofluid away from the hot region.As the angle λ increases from 0 to λ = π/2 the center of the vorticities moves counterclockwise with the hot and cold regions of the wall.For the lower power-low index fluids, the centers of the vortices are closer together and a strong flow circulation is seen.For a high power-low index, such as n = 1.3, the movement of the vorticities reduces as the friction force of the fluid increases.
Figure 5 shows the variations of the isotherm lines inside the circular cavity for different values of the angle, λ, and power law index, n, when Ra T = 10 4 , Ra M = 10 7 , Pr = 100, Ma = 10 −4 , β c = 60, δ = 30 • and L = 0.3.For each case, the isotherms are near vertical close to the hot sector of the cavity due to the dominance of the heat conduction in this area.Similar vertical isotherms are seen for all values of λ when n = 1.3.This is because the viscosity of the ferrofluid is higher in this case, and so more energy is dissipated due to the flow of the fluid.In turn, this leads to less flow overall and, therefore, less convection inside the circular cavity.Consequently, heat transfer is primarily by conduction for n = 1.3.In the other cases, the circulation of the fluid is more apparent and thermo-magnetic convective heat transfer is significant.

Results and discussion
The following parameters were studied in the ranges specified: 10 ≤ Ra T ≤ 10 , 10 ≤  ≤ 10 , 0. For  = 0, it is seen that the fluid velocity is greater for fluids with a lower power law index,  (Figure 4).This is attributed to the fact that the viscosity of a fluid with a lower power law index is lower for a given state of flow, which allows the ferrofluid to flow faster under given conditions.Also, the center of the vorticities is closer to the cold wall, showing the tendency of the movement of the ferrofluid away from the hot region.
As the angle  increases from 0 to  = π 2 ⁄ the center of the vorticities moves counterclockwise with the hot and cold regions of the wall.For the lower power-low index fluids, the centers of the vortices are closer together and a strong flow circulation is seen.For a high power-low index, such as  = 1.3, the movement of the vorticities reduces as the friction force of the fluid increases.Figure 5 shows the variations of the isotherm lines inside the circular cavity for different values of the angle, , and power law index, , when Ra T = 10 , Ra M = 10 , Pr = 100, Ma = 10 ,  c = 60,  = 30° and  = 0.3.For each case, the isotherms are near vertical close to the hot sector of the cavity due to the dominance of the heat conduction in this area.Similar vertical isotherms are seen for all values of  when  = 1.3.This is because the viscosity of the ferrofluid is higher in this case, and so more energy is dissipated due to the flow of the fluid.In turn, this leads to less flow overall and, therefore, less convection inside the circular cavity.Consequently, heat transfer is primarily by conduction for  = 1.3.In the other cases, the circulation of the fluid is more apparent and thermo-magnetic convective heat transfer is significant.The local Nusselt number is higher in the hot region of the cavity than the cold region (Figure 6, where the "hot arc" and "cold arc" are those regions of the circumference of the cavity which have fixed temperatures T h and T c , respectively).The local Nusselt number represents the heat transfer at these zones.The rate of heat transfer close to the hot arc is higher than the cold one.For both arcs, the heat transfer rate increases considerably toward the ends of arcs, as in the middle region of each arc there is not much fluid displacement and so heat transfer is dominated by conduction.Conversely, in the end regions of the arcs, the movement of the fluid base and buoyancy effects causes the heat transfer to increase considerably.The greater the power law index, the lower the rate of heat transfer.This is the result of greater energy dissipation and correspondingly lower fluid flow and, so, convective heat transfer.Figure 5 shows the variations of the isotherm lines inside the circular cavity for different values of the angle, , and power law index, , when Ra T = 10 , Ra M = 10 , Pr = 100, Ma = 10 ,  c = 60,  = 30° and  = 0.3.For each case, the isotherms are near vertical close to the hot sector of the cavity due to the dominance of the heat conduction in this area.Similar vertical isotherms are seen for all values of  when  = 1.3.This is because the viscosity of the ferrofluid is higher in this case, and so more energy is dissipated due to the flow of the fluid.In turn, this leads to less flow overall and, therefore, less convection inside the circular cavity.Consequently, heat transfer is primarily by conduction for  = 1.3.In the other cases, the circulation of the fluid is more apparent and thermo-magnetic convective heat transfer is significant.The local Nusselt number is higher in the hot region of the cavity than the cold region (Figure 6, where the "hot arc" and "cold arc" are those regions of the circumference of the cavity which have fixed temperatures  h and  c , respectively).The local Nusselt number represents the heat transfer at these zones.The rate of heat transfer close to the hot arc is higher than the cold one.For both arcs, the heat transfer rate increases considerably toward the ends of arcs, as in the middle region of each arc there is not much fluid displacement and so heat transfer is dominated by conduction.Conversely, in the end regions of the arcs, the movement of the fluid base and buoyancy effects causes the heat transfer to increase considerably.The greater the power law index, the lower the rate of heat transfer.This is the result of greater energy dissipation and correspondingly lower fluid flow and, so, convective heat transfer.

Hot arc
Cold Arc  The local Nusselt number is higher in the hot region of the cavity than the cold region (Figure 6, where the "hot arc" and "cold arc" are those regions of the circumference of the cavity which have fixed temperatures  h and  c , respectively).The local Nusselt number represents the heat transfer at these zones.The rate of heat transfer close to the hot arc is higher than the cold one.For both arcs, the heat transfer rate increases considerably toward the ends of arcs, as in the middle region of each arc there is not much fluid displacement and so heat transfer is dominated by conduction.Conversely, in the end regions of the arcs, the movement of the fluid base and buoyancy effects causes the heat transfer to increase considerably.The greater the power law index, the lower the rate of heat transfer.This is the result of greater energy dissipation and correspondingly lower fluid flow and, so, convective heat transfer.

Hot arc
Cold Arc The average Nusselt number, Nu, inside the circular cavity does not change significantly as a function of the location of cold and hot arcs of the cavity wall, as given by the variable λ (Figure 7).However, for power law indices, n ≥ 1, the average Nusselt number shows a slight decrease as an increasing function of λ, and for the power-low indices n = 0.7 and 0.8, it exhibits a slight increase.The average Nusselt number, Nu avg , has its maximum value at λ = 45 • for all power law index values except n = 0.7, which has a maximum at 90 • .As mentioned previously, the lower power law indices result in less viscous energy dissipation, and so faster fluid flow.Therefore, the fluid with n = 0.7 (the lowest in this study) experiences the highest average Nusselt number, while n = 1.3 (greatest value) has the minimum value.For clarity, Figure 7b shows the relative change in the average Nusselt number that are shown in Figure 7a, compared to the base condition (λ = 0, at which the arcs are located on opposite sides of the cylinder to each other with their centers on the x-axis).This change is denoted as the variable P and is expressed as a percentage in Figure 7b.As λ increases up to λ = 30 • , P increases up to around +5%, indicating an increase in the rate of heat transfer.After this angle, fluids with a low power law indexes, i.e., n = 0.7 and 0.8, keep increasing in heat transfer rate, albeit more slowly than at lower angles; however, those fluids with power index values of n ≥ 0.9 show a decrease of average Nusselt number with increasing angle, λ, up to P = −15%.It worth mentioning that although P decreases above λ = 30 • , it is never negative in the case of n = 0.9, indicating that the lowest average Nusselt number for that case occurs when λ = 0. maximum at 90°.As mentioned previously, the lower power law indices result in less viscous energy dissipation, and so faster fluid flow.Therefore, the fluid with  = 0.7 (the lowest in this study) experiences the highest average Nusselt number, while  = 1.3 (greatest value) has the minimum value.For clarity, Figure 7b shows the relative change in the average Nusselt number that are shown in Figure 7a, compared to the base condition ( = 0, at which the arcs are located on opposite sides of the cylinder to each other with their centers on the -axis).This change is denoted as the variable  and is expressed as a percentage in Figure 7b.As  increases up to  = 30°,  increases up to around +5%, indicating an increase in the rate of heat transfer.After this angle, fluids with a low power law indexes, i.e.,  = 0.7 and 0.8, keep increasing in heat transfer rate, albeit more slowly than at lower angles; however, those fluids with power index values of  ≥ 0.9 show a decrease of average Nusselt number with increasing angle, , up to  = −15%.It worth mentioning that although  decreases above  = 30°, it is never negative in the case of  = 0.9, indicating that the lowest average Nusselt number for that case occurs when  = 0. Figures 8 and 9 demonstrate the effects of the thermal and magnetic Rayleigh numbers on the streamlines and temperature field inside the circular cavity when n = 0.8, λ = 0, Pr = 100, Ma = 10 −4 , β c = 60, δ = 30 • and L = 0.3.For a specific magnetic Rayleigh number, as the buoyancy Rayleigh number increases, the buoyancy force imposed on the fluid flow increases.This leads to stronger fluid flow and enhances the strength of the vortices inside the cavity.As the buoyancy Rayleigh number increases, the conductive area close to the hot and cold arcs diminishes, and the effect of the convective term increases due to the greater buoyancy effects.At a buoyancy Rayleigh number of Ra T = 10 4 , the increase in the magnetic Rayleigh number leads to the increase of fluid flow circulation by Kelvin forces.The increase in the magnetic Rayleigh number enhances flow rotation and heat transfer.In terms of isotherms, an increase of the buoyancy Rayleigh number contributes to the smaller vertical area of the isotherms close to the hot and cold arcs, which implies that conduction is less important and free convection dominates.Also, the effect of the magnetic Rayleigh number on the isotherms in a specific buoyancy Rayleigh number is not significant.
Table 2 presents the effects of the magnetic and buoyancy Rayleigh numbers on the average Nusselt number inside the cavity.It is evident that both the magnetic and buoyancy Rayleigh numbers cause the average Nusselt number within the circular cavity to increase.The mechanism by which the magnetic Rayleigh number manifest itself is the circulation of the fluid flow, while the buoyancy Rayleigh number is concerned with the upward movement of the fluid due to buoyancy, causing intense fluid flow.When buoyancy effects are neglected (i.e., Ra T = 0), only the Kelvin force is imposed on the fluid flow, and the average Nusselt number increases by approximately 100% as the magnetic Rayleigh number increases from 10 7 to 10 8 .There is a similar increase when the buoyancy Rayleigh number increases from 10 5 to 10 6 at constant magnetic Rayleigh number.However, for higher buoyancy Rayleigh numbers, such as Ra T = 10 6 , the average Nusselt number decreases slightly when magnetic Rayleigh number increases.
vortices inside the cavity.As the buoyancy Rayleigh number increases, the conductive area close to the hot and cold arcs diminishes, and the effect of the convective term increases due to the greater buoyancy effects.At a buoyancy Rayleigh number of Ra T = 10 , the increase in the magnetic Rayleigh number leads to the increase of fluid flow circulation by Kelvin forces.The increase in the magnetic Rayleigh number enhances flow rotation and heat transfer.In terms of isotherms, an increase of the buoyancy Rayleigh number contributes to the smaller vertical area of the isotherms close to the hot and cold arcs, which implies that conduction is less important and free convection dominates.Also, the effect of the magnetic Rayleigh number on the isotherms in a specific buoyancy Rayleigh number is not significant.Figure 10 shows the dependence of the flow streamlines and isotherms on the magnet size for two different Rayleigh numbers when Ra M = 10 7 , n = 0.8, λ = π/3, Pr = 100, Ma = 10 −4 and β c = 60.The streamline at Ra T = 0 (i.e., when buoyancy is neglected and only the Kelvin force is imposed on the fluid flow) shows that a stronger magnet helps the rotation of the fluid flow inside the circular cavity.There are two vorticities for L = 0.1 and L = 0.3, while four vorticities are formed at L = 0.5.However, as the space available for the fluid flow diminishes, since the (larger) magnet takes up more of the space in the cavity, the strength of the vorticities reduces at L = 0.5.In terms of isotherms at this buoyancy Rayleigh number, the bigger magnet causes the dominant cooling mechanism to change from convection to conduction.At Ra T = 10 5 , as in the previous case, a larger magnet limits the space for the flow rotation and movement, leading to the weaker vorticities.However, the buoyancy effect due to the higher Rayleigh number, causes the strength of the fluid circulation to increase considerably, compared to the case where there is no buoyancy force.This enhances the buoyancy-driven convective heat transfer inside the circular cavity.Table 2 presents the effects of the magnetic and buoyancy Rayleigh numbers on the average Nusselt number inside the cavity.It is evident that both the magnetic and buoyancy Rayleigh numbers cause the average Nusselt number within the circular cavity to increase.The mechanism by which the magnetic Rayleigh number manifest itself is the circulation of the fluid flow, while the buoyancy Rayleigh number is concerned with the upward movement of the fluid due to buoyancy, causing intense fluid flow.When buoyancy effects are neglected (i.e., Ra T = 0), only the Kelvin force is imposed on the fluid flow, and the average Nusselt number increases by approximately 100% as the magnetic Rayleigh number increases from 10 to 10 .There is a similar increase when the buoyancy Rayleigh number increases from 10 to 10 at constant magnetic Rayleigh number.However, for higher buoyancy Rayleigh numbers, such as Ra T = 10 , the average Nusselt number decreases slightly when magnetic Rayleigh number increases.Figure 11 shows the variations of the average Nusselt number due to different magnet sizes for the entire range of power law indexes at Ra T = 0 and Ra T = 10 6 .For the case of Ra T = 0, where buoyancy is not possible and so only the Kelvin force is imposed on the flow, the average Nusselt number varies from 0 to 6.It increases initially, as the magnet is made bigger, up to some size between L = 0.2 and L = 0.3, depending on the power law index.Above this peak value, the average Nusselt number decreases.This is due to the limited space available for fluid flow when the magnet becomes large, which restricts heat transfer.Also, as is expected, low power law index values reduce the energy dissipation in the fluid, leading to more fluid flow and a high Nusselt number.The average Nusselt number when Ra T = 10 6 is constant for magnet sizes up to l = 0.55.Above this, fluids with a power law index of n ≤ 1 exhibit an upward trend of average Nusselt number, while other fluids show a fall.The reason of this constant trend is that the buoyancy force is dominant compared to the Kelvin force for smaller magnets, while the influence of the Kelvin force dominates the buoyancy force at magnet sizes larger than l = 0.55.The increase in average Nusselt number for power law index ≤ 1 is due to the low energy dissipation, even when the magnet is so large as to restrict the space for fluid flow.For fluid with the power-low exponent n > 1, the high viscosity leads to significant energy dissipation, and therefore, low heat transfer and a corresponding low average Nusselt number.Figure 11 shows the variations of the average Nusselt number due to different magnet sizes for the entire range of power law indexes at Ra T = 0 and Ra T = 10 .For the case of Ra T = 0, where buoyancy is not possible and so only the Kelvin force is imposed on the flow, the average Nusselt number varies from 0 to 6.It increases initially, as the magnet is made bigger, up to some size between  = 0.2 and  = 0.3, depending on the power law index.Above this peak value, the average Nusselt number decreases.This is due to the limited space available for fluid flow when the magnet becomes large, which restricts heat transfer.Also, as is expected, low power law index values reduce the energy dissipation in the fluid, leading to more fluid flow and a high Nusselt number.The average Nusselt number when Ra T = 10 is constant for magnet sizes up to  = 0.55.Above ence of the Kelvin force dominates the buoyancy force at magnet sizes larger than  = 0.55.The increase in average Nusselt number for power law index ≤ 1 is due to the low energy dissipation, even when the magnet is so large as to restrict the space for fluid flow.
For fluid with the power-low exponent  > 1, the high viscosity leads to significant energy dissipation, and therefore, low heat transfer and a corresponding low average Nusselt number.
Ra T = 0.0 Ra T = 10 Figure 12 shows streamlines and isotherms within the cavity for different angular sizes of the heated and cooled arcs, .Two connected vorticities are formed at the left and right sides of the magnet, and as  increases, the strength of the flow increases within the cavity.In fact, when the arc size increases, more heat reaches the ferrofluid, and consequently, the flow strength grows.This flow strength is higher at the right side close to the cold arc of the cavity.Also, streamlines closer to the magnet show greater strength due to the positive effects of the magnet on the flow movement.The flow velocity is higher close to the colder side as the flow starts from the hot section (on the left of the domain) toward the cold one (on the right side).The numbers represented by the isotherms indicate the greater importance of convective heat transfer as δ increases.For example, when δ = 45°, the horizontal isotherms representing convective term occupy more area of the cavity, while for δ = 15°, this area is smaller, and the conductive term plays a major role in the heat transfer rate.The average Nusselt number for δ = 15°, 30°, 45°, and 60° is 13.18, 11.32, 10.22, and 9.37, respectively.When the heated surface expands (i.e., a larger value of ), the temperature gradient between the hot wall and the ferrofluid close to it decreases, leading to a lower rate of the heat transfer and a lower average Nusselt number.Figure 12 shows streamlines and isotherms within the cavity for different angular sizes of the heated and cooled arcs, δ.Two connected vorticities are formed at the left and right sides of the magnet, and as δ increases, the strength of the flow increases within the cavity.In fact, when the arc size increases, more heat reaches the ferrofluid, and consequently, the flow strength grows.This flow strength is higher at the right side close to the cold arc of the cavity.Also, streamlines closer to the magnet show greater strength due to the positive effects of the magnet on the flow movement.The flow velocity is higher close to the colder side as the flow starts from the hot section (on the left of the domain) toward the cold one (on the right side).The numbers represented by the isotherms indicate the greater importance of convective heat transfer as δ increases.For example, when δ = 45 • , the horizontal isotherms representing convective term occupy more area of the cavity, while for δ = 15 • , this area is smaller, and the conductive term plays a major role in the heat transfer rate.The average Nusselt number for δ = 15 • , 30 • , 45 • , and 60 • is 13.18, 11.32, 10.22, and 9.37, respectively.When the heated surface expands (i.e., a larger value of δ), the temperature gradient between the hot wall and the ferrofluid close to it decreases, leading to a lower rate of the heat transfer and a lower average Nusselt number.

Conclusions
In this study, the thermofluidic characteristics of a non-Newtonian ferrofluid were investigated numerically.A circular cavity was considered, and two arcs located on opposite sides if its perimeter were considered hot and cold, while the other parts were fully insulated.COMSOL Multiphysics commercial code was utilized to perform numerical calculations.The effects of different parameters on the thermofluidic behavior of the non-Newtonian ferrofluid including thermal and magnetic Rayleigh numbers, power law index of the (non-Newtonian) fluid, the position and size of the hot and cold arcs at the cavity walls and the dimension of the magnet were surveyed in this study.It was seen that an increase in the power law index inhibited the fluid movement due to stronger

Conclusions
In this study, the thermofluidic characteristics of a non-Newtonian ferrofluid were investigated numerically.A circular cavity was considered, and two arcs located on opposite sides if its perimeter were considered hot and cold, while the other parts were fully insulated.COMSOL Multiphysics commercial code was utilized to perform numerical calculations.The effects of different parameters on the thermofluidic behavior of the non-Newtonian ferrofluid including thermal and magnetic Rayleigh numbers, power law index of the (non-Newtonian) fluid, the position and size of the hot and cold arcs at the cavity walls and the dimension of the magnet were surveyed in this study.It was seen that an increase in the power law index inhibited the fluid movement due to stronger friction.The vorticities formed inside the cavity tended to move toward the cold arc location.Also, as the power law index increased, vorticities form further from each other.The isotherm pattern in the high power law index regime demonstrated the relative dominance of conductive heat transfer over convection.The heat transfer rate increased at the end of both arcs compared to the middle part of each arc because of the greater buoyancy effects at the ends of the arcs and was greater close to the hot arc.The location of the arcs did not affect the heat transfer rate considerably.The heat transfer rate was greater for lower power law indices due to the increased contributions of flow movements and free convection.As the buoyancy Rayleigh number increases, stronger vorticities form inside the cavity and enhance the effects of convection.An increase in the magnetic Rayleigh number contributed to stronger circulation of the flow inside the cavity and also led to higher heat transfer.When the Kelvin force was the only one imposed on the flow, a larger magnet increased the heat transfer for magnets of size 0.2 ≤ L ≤ 0.3; however the heat transfer decreased for larger magnets, as the magnet occupied more space inside the cavity and hindered fluid

Figure 1 .
Figure 1.Computational domain along with the Cartesian coordinate system used in this study.

Figure 1 .
Figure 1.Computational domain along with the Cartesian coordinate system used in this study.

Figure 2 .
Figure 2. Dependence of the average Nusselt number, Nu avg , and max inside the cavity,  max , on the number of elements in the mesh used for Figure verify the outcomes for a non-Newtonian fluid, the average Nu

Figure 2 .
Figure 2. Dependence of the average Nusselt number, Nu avg , and maximum strength of the fluid inside the cavity, ψ max , on the number of elements in the mesh used for numerical simulation.

Figure 3 .Figure 3 .
Figure 3. (a,b) Streamlines patterns (left column) and temperature (temperature) of Newtonian ferro-fluid; (a) Ashoori et al. [18] and (b) present investigation (Ra T = 0, Ra M = 10 , Ma = 0, Pr = 1000,  = 1.0,  = 0.2,  c = 60 and  = 5.0); (c) comparison of local average Nusselt number calculated in the current study and values available in literature [19].Permission to use the figures contained in (a) is granted via the terms of the International Association of Scientific, Technical & Medical Publishers, to which both the copyright holder of the original work and MDPI are signatories.

Figure 5 .
Figure 5. Dependence of the isotherms on the power law index, , and angle of the hot and cold arcs relative to the horizontal, , when Ra T = 10 , Ra M = 10 , Pr = 100, Ma = 10 ,  c = 60,  = 30° and  = 0.3.

Figure 7 . 3 λ 3 Figure 7 .
Figure 7. Dependence for different values of power law index, , on the angle of the location of the centers of the hot and cold arcs of the cavity wall,  of (a) the average Nusselt number in the entire cavity, Nu avg , and (b) changes in the average Nusselt number, compared to the base case of  = 0, when Ra T = 10 , Ra M = 10 , Pr = 100, Ma = 10 ,  = 60,  = 30° and  = 0.3

Table 1 .
[20]arison between the average Nusselt number, Nu avg in a magnetic, in the work conducted Khezar et al.[20]and this study when Ra T = 10 5 and Pr = 100.

Table 1 .
[20]arison between the average Nusselt number, Nu avg in a magnetic, in the work conducted Khezar et al.[20]and this study when Ra = 10 and Pr = 100.