Effect of Straight, Inclined and Curved Fins on Natural Convection and Entropy Generation of a Nanofluid in a Square Cavity Influenced by a Magnetic Field

In this paper, the free convective heat transfer of nanofluids in a square cavity is simulated using a numerical method. The angle of the cavity could be changed in the horizontal axis from 0 to 90 degrees. The cavity is exposed under a constant magnetic field. Two opposite walls of the cavity are cold and warm, and the rest of the walls are insulated. On the hot wall, there are two fins with the same wall temperature. The equations were discretized by the finite volume method (FVM) and then solved using the SIMPLE algorithm. Three different fin configurations (straight, inclined and curved) were studied in terms of heat transfer rate and generation of entropy. According to the simulation results, the heat transfer rate was improved by tilting the fins toward the top or bottom of the cavity. At Ra = 105 and Ha = 20, the maximum heat transfer rate was achieved at a cavity inclination of 90° and 45°, respectively, for straight and curved fins. In the horizontal cavity, heat transfer rate could be improved up to 6.4% by tilting the fins and up to 4.9% by warping them. Increasing the Hartmann number from 0 to 40 reduced the Nusselt number and entropy generation by 37.9% and 33.8%, respectively.


Introduction
Natural convection heat transfer has extensive applications in several industries. This type of heat transfer occurs without any force applied by a fan or similar equipment, for example, in ovens and furnaces. Heat transfer by natural convection in closed enclosures has been used in various industries, such as aerospace and solar industries [1,2]. Accordingly, there has been a growing interest in this type of enclosure in recent years [3][4][5]. Additionally, numerous studies have been recently conducted on nanofluids in closed enclosures as an attractive research area.
Nanoparticles have many applications, and researchers have used nanoparticles in various fields, including solar, heat exchangers, medicine, etc. [6][7][8]. The enhancement of heat transfer rate (HTR) by adding nanofluids is one motivation behind these studies. Several types of research have focused on the HTR of nanofluids [9][10][11][12][13][14][15][16][17]. 2 of 23 Fins can be used to improve HTR in cavities by increasing heat transfer surface area. The effectiveness of fins has been confirmed in forced convection. However, due to reduced velocity in natural convection, the size and shape of the fin determine whether it improves or diminishes HTR, particularly in cavities [18][19][20][21][22][23][24]. Alnaqi et al. [25] placed a fin of constant thermal conductivity in an inclined square cavity to investigate its effect on natural convection heat transfer. Their results indicated an intensification in the Nusselt number (Nu) with the Rayleigh number (Ra). Alrashed [26] compared the effect of two fins with different inclinations (angles) on the HTR. The fins with adjustable inclinations were placed on the sides of the cavity. This study revealed the significant impact of fin inclination on HTR in the cavity.
Magnetic fields may also influence HTR. According to the literature on forced convection, HTR may be enhanced by applying a magnetic field [27]. In contrast, natural convection heat transfer may be diminished by applying a magnetic field due to a reduction in the fluid velocity [28][29][30][31][32][33][34][35][36][37]. For example, Pordanjani et al. [38] considered the effect of fins on HTR in an enclosure under a magnetic field of different intensities. According to their results, the magnetic field reduced HTR by decreasing fluid velocity. The HTR further decreased at higher field intensities.
Generation of entropy (GOE) can be calculated as a better way of studying the efficiency of thermal equipment. GOE is a measure of irreversibility and energy loss in the system. Accordingly, in addition to HTR, GOE has also been investigated in numerous studies [39][40][41][42][43][44][45]. Zhou et al. [46] studied HTR and GOE in a nanofluid-filled cavity at a Grashof number (Ga) of 10 4 -10 6 and a nanoparticle volume fraction of 0.01-0.15%. They found an increase in HTR and local GOE with an increasing Ga and nanoparticle volume fraction. Dutta et al. [47] investigated GOE and HTR of the Cu-water nanofluid in a rhombic enclosure exposed to a magnetic field. Their results showed a decrease in GOE with an increasing Hartmann number (Ha) at all Rayleigh number (Ra) values and cavity inclinations.
In another study, Li et al. [44] examined the free convection of nanofluids in a cavity with entropy generation despite a circular barrier in the center of the cavity. Their findings showed that increasing Ra could increase Nu by 4.5 times, but increasing the Hartmann number reduced GOE by 35%.
Fin-equipped cavities are widely used in the industry, including in furnaces and kitchen ovens. HTR may be improved by adding a nanofluid to the cavity. However, it should be noted that this can have different impacts on natural convection heat transfer due to increased viscosity. The flow in the cavity can be controlled by a magnetic field. Fins with different geometries and inclinations can significantly affect heat transfer in cavities. This study compares the effect of three different fins on heat transfer and GOE in a square cavity exposed to a magnetic field. The cavity was filled with a nanofluid and can assume different inclinations relative to the horizon. Two fins were installed on the hot wall of the cavity, and the effect of straight, inclined and curved fins was studied. As an innovative aspect of this study, the impact of the curved fin on the flow and temperature fields was investigated. GOE was also investigated by adjusting the fin inclination, curvature and cavity inclination.

Problem Description
A 2D square cavity with unit length filled with a Al 2 O 3 -water nanofluid was studied. The cavity can tilt and make an angle of γ with the horizon. A constant, horizontal magnetic field was applied to the cavity. The right and left walls were respectively maintained at constant temperatures of T c and T H . Moreover, the top and bottom sides were insulated. Two fins were installed on the hot wall at the same temperature as the wall. As illustrated by Figure 1, the fins installed may be straight, inclined or curved. Natural convection was established in the cavity by Boussinesq approximation.
Processes 2021, 9, x FOR PEER REVIEW 3 of 24 wall. As illustrated by Figure 1, the fins installed may be straight, inclined or curved. Natural convection was established in the cavity by Boussinesq approximation.

Governing Equations and Boundary Condition
The equations governing the nanofluid flow are given below. The equations with the indices a and b are, respectively, dimensional and dimensionless. These equations include conservations of mass, momentum, and energy [38,44].
The equations with the indices a were nondimensionalized using the following relations and transformed to those with the index b [36,39].
According to the dimensionless Equation (5), the dimensionless boundary conditions were obtained, as shown in Figure 2.

Governing Equations and Boundary Condition
The equations governing the nanofluid flow are given below. The equations with the indices a and b are, respectively, dimensional and dimensionless. These equations include conservations of mass, momentum, and energy [38,44].
The equations with the indices a were nondimensionalized using the following relations and transformed to those with the index b [36,39].
According to the dimensionless Equation (5), the dimensionless boundary conditions were obtained, as shown in Figure 2.

Auxiliary Equations
GOE can be calculated from the following dimensionless equation [36,39].
The relations used for measuring the nanofluid properties are presented below [36].

Auxiliary Equations
GOE can be calculated from the following dimensionless equation [36,39].
Further, the total GOE is expressed as follows: The Nusselt number (Nu) shows the extent of heat transfer and is calculated as follows: Moreover, the average Nu is obtained from the following equation.
The relations used for measuring the nanofluid properties are presented below [36].
Processes 2021, 9, 1339 5 of 23 Various relations have been presented in the literature [48] for measuring the thermal conductivity and viscosity of nanofluids. This model considers the Brownian motion of nanoparticles for better simulation of the thermal conductivity and viscosity. The Brownian motion of nanoparticles in nanofluids refers to their random continuous motion in the fluid induced by collisions with liquid molecules, leading to dispersion of nanoparticles in the fluid [49].
The thermal model used in [48,50] is depicted in Figure 3. Further, this study uses the viscosity models presented by Vajjha et al. [48] and Brinkman [51]. The values of β and f (T‚ϕ) for the Al 2 O 3 -water nanofluid are given below [48].  and different grid sizes.

Validation
First, the results from the developed model were compared with those reported by Oztop and Abu-Nada [52], who investigated a simple nanofluid-filled cavity with vertical hot and cold walls. According to Figure 4, which plots the average Nu versus nanoparticle volume fraction, there is a slight discrepancy between the results.

Numerical Procedure
The equations were discretized before being solved by the finite volume method. The cavity geometry was meshed uniformly using identical quadrilateral elements. Next, a Fortran code was developed based on the SIMPLE algorithm. Each run took nearly 40 min on a computers with a 3.7 GHZ core i7 CPU and 16 GB DDR4 RAM. Further, the following convergence criterion was selected to terminate the solution.

Grid Analysis
To find the best grid size, the grid independence was investigated based on the average Nu over the cold wall and GOE at Ra = 10 5 , Ha = 20, ϕ = 0.03, γ = 0 • by changing the grid size from 3600 to 40,000. The variations of the average Nu and GOE for different grid sizes are shown in Figure 3. Eventually, a grid with 25,600 elements (160 × 160) was adopted.

Validation
First, the results from the developed model were compared with those reported by Oztop and Abu-Nada [52], who investigated a simple nanofluid-filled cavity with vertical hot and cold walls. According to Figure 4, which plots the average Nu versus nanoparticle volume fraction, there is a slight discrepancy between the results.

Validation
First, the results from the developed model were compared with those reported by Oztop and Abu-Nada [52], who investigated a simple nanofluid-filled cavity with vertical hot and cold walls. According to Figure 4, which plots the average Nu versus nanoparticle volume fraction, there is a slight discrepancy between the results.  Max. Error=5.34% Max. Error=3.01% Max. Error=1.91% Next, the model was evaluated for flow under a magnetic field. Accordingly, the model was compared with the results of Aminossadati et al. [53]. The average Nu for various Nu and Ha is indicative of the magnetic field are plotted in Figure 5. Obviously, the developed model exhibits excellent accuracy for natural convection in a nanofluid exposed to a magnetic field and, therefore, can be employed.
Processes 2021, 9, x FOR PEER REVIEW 7 of 24 the developed model exhibits excellent accuracy for natural convection in a nanofluid exposed to a magnetic field and, therefore, can be employed.

Figure 5.
A comparison of the developed model with the results of previous works for a nanofluid exposed to a magnetic field.
According to the two comparisons made above, the highest error of the results of this article with previous articles is 5.34% and 4.82% for the results of Oztop and Abu-Nada [52] and Aminossadati et al. [53], respectively. From the above comparison, it was found that there is a slight difference between the results of the present work and the results of the references [52,53]. Therefore, the code used for this research has acceptable accuracy.  Figure 6a shows that as the fins are angled, the number of stream functions increases at the center of the vortex. This can be caused by the tip of the fins approaching the walls and thus increasing the space in the middle of the cavity. Increasing the value in the center of the contours means increasing the maximum value of the stream function, which also means increasing the spin speed of the vortex. In all the above cases, the vortex rotates clockwise. Vortex penetration in the middle of the fins can lead to an increase in fluid collision with the wall in that area. In these areas, the rate of heat transfer also increases. In Figure 6b, in temperature counters, the presence of more constant temperature lines in an area means an increase in the temperature difference in that area. The highest density of constant temperature lines is in the two areas of the upper part of the right wall and the lower part of the left wall. In these areas, due to the direction of rotation of the vortex, the temperature difference is high, and the density of the lines is higher. The heat transfer rate is also high in these areas. When the Max. Error=2.34%

Effect of Fin Inclination
Max. Error=2.12% Max. Error=4.82% Figure 5. A comparison of the developed model with the results of previous works for a nanofluid exposed to a magnetic field.
According to the two comparisons made above, the highest error of the results of this article with previous articles is 5.34% and 4.82% for the results of Oztop and Abu-Nada [52] and Aminossadati et al. [53], respectively. From the above comparison, it was found that there is a slight difference between the results of the present work and the results of the references [52,53]. Therefore, the code used for this research has acceptable accuracy.  Figure 6a shows that as the fins are angled, the number of stream functions increases at the center of the vortex. This can be caused by the tip of the fins approaching the walls and thus increasing the space in the middle of the cavity. Increasing the value in the center of the contours means increasing the maximum value of the stream function, which also means increasing the spin speed of the vortex. In all the above cases, the vortex rotates clockwise. Vortex penetration in the middle of the fins can lead to an increase in fluid collision with the wall in that area. In these areas, the rate of heat transfer also increases. In Figure 6b, in temperature counters, the presence of more constant temperature lines in an area means an increase in the temperature difference in that area. The highest density of constant temperature lines is in the two areas of the upper part of the right wall and the lower part of the left wall. In these areas, due to the direction of rotation of the vortex, the temperature difference is high, and the density of the lines is higher. The heat transfer rate is also high in these areas. When the fin is lower, line densities may occur on the fin tip, so heat transfer from the tip of the fin increases in this area. In Figure 6c, the appearance of entropy production lines in locally produced entropy contours means the formation of entropy in that area. By comparing the temperature and entropy contours produced, it can be seen that the density areas of the lines are similar in both. The dependence of the generated entropy on the temperature gradient has caused the density of the entropy field lines to be very similar to the temperature field. It can also be seen that with different variations of the fin angle, the amount of local entropy is always high in the lower corner of the lower fin. Increasing the angle of the fins causes the place of maximum entropy generation to be closer to the lower wall.

Effect of Fin Inclination
Processes 2021, 9, x FOR PEER REVIEW 8 of 24 fin is lower, line densities may occur on the fin tip, so heat transfer from the tip of the fin increases in this area. In Figure 6c, the appearance of entropy production lines in locally produced entropy contours means the formation of entropy in that area. By comparing the temperature and entropy contours produced, it can be seen that the density areas of the lines are similar in both. The dependence of the generated entropy on the temperature gradient has caused the density of the entropy field lines to be very similar to the temperature field. It can also be seen that with different variations of the fin angle, the amount of local entropy is always high in the lower corner of the lower fin. Increasing the angle of the fins causes the place of maximum entropy generation to be closer to the lower wall. The variations of the dimensionless velocity (a) and dimensionless temperature (b) over the line X = 0.5 are plotted in Figure 7 at different fin inclinations for the specified constants. In the dimensionless velocity diagram, tilting the fin toward the top and bottom of the cavity causes an increase in the maximum and minimum dimensionless velocities. This can be attributed to the increased distance of the fins from the center of the cavity, as well as the increased vortex velocity. Nonetheless, it must be noted that the orientation of the fins also plays an important role in the loci of maximum and minimum velocities n the cavity. When the fin is tilted up, the maximum velocity occurs at the bottom of the cavity, while the highest velocity is found at the top when the fins are tilted down. According to the dimensionless temperature diagram, fluid temperature increases toward the top of the cavity due to the upward convection of hot fluid by buoyancy. Temperature changes with cavity inclination as a result of changes in the vortex velocity and fin location in the cavity. The variations of the dimensionless velocity (a) and dimensionless temperature (b) over the line X = 0.5 are plotted in Figure 7 at different fin inclinations for the specified constants. In the dimensionless velocity diagram, tilting the fin toward the top and bottom of the cavity causes an increase in the maximum and minimum dimensionless velocities. This can be attributed to the increased distance of the fins from the center of the cavity, as well as the increased vortex velocity. Nonetheless, it must be noted that the orientation of the fins also plays an important role in the loci of maximum and minimum velocities n the cavity. When the fin is tilted up, the maximum velocity occurs at the bottom of the cavity, while the highest velocity is found at the top when the fins are tilted down. According to the dimensionless temperature diagram, fluid temperature increases toward the top of the cavity due to the upward convection of hot fluid by buoyancy. Temperature changes with cavity inclination as a result of changes in the vortex velocity and fin location in the cavity.
The variations of (a) the local Nu over the cold wall of the cavity for negative and (b) the local Nu over the cold wall of the cavity for positive trigonometric fin inclination angles for Ra = 10 5 , Ha = 20, ϕ = 0.03, γ = 0 • are plotted in Figure 8 for the specified conditions. Figure 8a shows that the most significant amount of Nu occurred at the top of the wall. In that area, the fluid inside the vortex, which is heated by the fins, strikes the cold wall for the first time, resulting in the highest HTR. Increasing the fin angle has slightly increased the Nu. Figure 8b also shows that the largest amount of Nu occurs in the upper area of the wall. Furthermore, increasing the angle in the negative direction of the fin increases the amount of heat transfer rate similar to that seen above. In the lower part of the wall, where the fluid has lost its temperature by moving along the side of the wall, the amount of heat transfer is less, and as a result, the lowest Nu is seen in that area.
The average Nu across the cold wall and GOE at different fin inclinations are plotted in Figure 9 for Ra = 10 5 , Ha = 20, ϕ = 0.03, γ = 0 • . Increasing the fin inclination either in upward or downward directions enhances HTR and GOE. The vortex rotates faster with increases the distance of the tip of the fin from the center of the cavity. Furthermore, the faster vortex rotation causes an increase in the temperature difference, which, in turn, enhances the HTR. Considering the fact that GOE depends on the temperature and velocity gradients, GOE is also increased. The variations of (a) the local Nu over the cold wall of the cavity for negative and (b) the local Nu over the cold wall of the cavity for positive trigonometric fin inclination angles for = 10 , = 20, = 0.03, = 0 are plotted in Figure 8 for the specified conditions. Figure 8a shows that the most significant amount of Nu occurred at the top of the wall. In that area, the fluid inside the vortex, which is heated by the fins, strikes the cold wall for the first time, resulting in the highest HTR. Increasing the fin angle has slightly increased the Nu. Figure 8b also shows that the largest amount of Nu occurs in the upper area of the wall. Furthermore, increasing the angle in the negative direction of the fin increases the amount of heat transfer rate similar to that seen above. In the lower part of the wall, where the fluid has lost its temperature by moving along the side of the wall, the amount of heat transfer is less, and as a result, the lowest Nu is seen in that area.   The variations of (a) the local Nu over the cold wall of the cavity for negative and (b) the local Nu over the cold wall of the cavity for positive trigonometric fin inclination angles for = 10 , = 20, = 0.03, = 0 are plotted in Figure 8 for the specified conditions. Figure 8a shows that the most significant amount of Nu occurred at the top of the wall. In that area, the fluid inside the vortex, which is heated by the fins, strikes the cold wall for the first time, resulting in the highest HTR. Increasing the fin angle has slightly increased the Nu. Figure 8b also shows that the largest amount of Nu occurs in the upper area of the wall. Furthermore, increasing the angle in the negative direction of the fin increases the amount of heat transfer rate similar to that seen above. In the lower part of the wall, where the fluid has lost its temperature by moving along the side of the wall, the amount of heat transfer is less, and as a result, the lowest Nu is seen in that area.

Effect of Fin Curvature
The (a) flow, (b) temperature and (c) entropy fields are depicted in Figure 10 for the specified constants. Figure 10a shows that increasing the curvature of the fin increases the strength of the flow field. As mentioned, increasing the curvature of the fin increases the value of the maximum stream function in the middle of the vortex due to the increase in the distance between the tip of the fins and the middle of the cavity. Increasing the space for fluid movement has increased the speed of fluid movement. Thus, as can be seen in the streamlines contour, the speed of the vortex has increased with increasing fin curvature. Figure 10b also shows slight changes in the isotherms with increasing fin curvature. These changes can be seen as a slight increase in the density of the lines in the areas under the fins. In Figure 10c, which is related to the entropy field, the compactness of the entropy production lines has increased with increasing curvature, which could indicate more entropy production, especially in those areas. Therefore, increasing the fin curvature can increase local entropy production.
in Figure 9 for = 10 , = 20, = 0.03, = 0 . Increasing the fin inclination either in upward or downward directions enhances HTR and GOE. The vortex rotates faster with increases the distance of the tip of the fin from the center of the cavity. Furthermore, the faster vortex rotation causes an increase in the temperature difference, which, in turn, enhances the HTR. Considering the fact that GOE depends on the temperature and velocity gradients, GOE is also increased.

Effect of Fin Curvature
The (a) flow, (b) temperature and (c) entropy fields are depicted in Figure 10 for the specified constants. Figure 10a shows that increasing the curvature of the fin increases the strength of the flow field. As mentioned, increasing the curvature of the fin increases the value of the maximum stream function in the middle of the vortex due to the increase in the distance between the tip of the fins and the middle of the cavity. Increasing the space for fluid movement has increased the speed of fluid movement. Thus, as can be seen in the streamlines contour, the speed of the vortex has increased with increasing fin curvature. Figure 10b also shows slight changes in the isotherms with increasing fin curvature. These changes can be seen as a slight increase in the density of the lines in the The variations of the dimensionless velocity (a) and dimensionless temperature (b) over the line X = 0.5 are plotted in Figure 11 for different fin curvatures and the specified constants. As seen, more significant maxima and minima dimensionless horizontal velocities are obtained with increasing fin curvatures. According to the flow field, this can be related to the higher vortex velocity due to an increase in the velocity in the cavity by increasing the curvature of the fins. Furthermore, with increasing the fin curvature, as their tips approach the upper parts of the cavity, velocity contours draw away from the left-side wall. The locus of maximum velocity is displaced toward the middle of the cavity. As seen in the dimensionless temperature diagram, the temperature decreases in the middle of the cavity but increases at the top by increasing the fin curvature. The tips of fins pulling away from the middle of the cavity and the displacement of the end parts of fins toward the top of the cavity are the main reasons behind these temperature variations.
The local Nu is plotted in Figure 12 for the specified constants. Apparently, the local heat transfer first increases at low curvatures, and the Nusselt number (Nu) changes slightly at higher curvatures. The increased Nu can be attributed to the higher fluid velocity, particularly at lower fin curvatures. The final reduction of the local Nu at the top of the wall is the result of the vortex failing to reach the sharp corners of the cavity, which reduces fluid-wall collisions and, consequently, HTR. Nevertheless, the high-temperature gradient in this area still allows for a considerable HTR at the upper corner of the cavity. areas under the fins. In Figure 10c, which is related to the entropy field, the compactness of the entropy production lines has increased with increasing curvature, which could indicate more entropy production, especially in those areas. Therefore, increasing the fin curvature can increase local entropy production. The variations of the dimensionless velocity (a) and dimensionless temperature (b) over the line X = 0.5 are plotted in Figure 11 for different fin curvatures and the specified constants. As seen, more significant maxima and minima dimensionless horizontal velocities are obtained with increasing fin curvatures. According to the flow field, this can be related to the higher vortex velocity due to an increase in the velocity in the cavity by increasing the curvature of the fins. Furthermore, with increasing the fin curvature, as their tips approach the upper parts of the cavity, velocity contours draw away from the left-side wall. The locus of maximum velocity is displaced toward the middle of the cavity. As seen in the dimensionless temperature diagram, the temperature decreases in the middle of the cavity but increases at the top by increasing the fin curvature. The tips of fins pulling away from the middle of the cavity and the displacement of the end parts of fins toward the top of the cavity are the main reasons behind these temperature variations. The average Nu and GOE are presented in Table 2. Increasing the curvature, at first, enhances HTR and GOE. However, at β = 0.48, increasing the curvature of the two fins causes a decrease in both parameters. A further increase in the fin curvature causes an increase in GOE and HTR. HTR and GOE strongly depend on the vortex velocity. Moving the ends of the fins away from the middle of the cavity accelerates the vortex. On the other hand, the limited reach of the flow toward the walls, particularly their middle parts, may reduce fluid-wall collisions and, consequently, the HTR. However, the changes in Nu and GOE with the fin curvature were found to be insignificant in the studied range (4.4 and 5.3, respectively).  The local Nu is plotted in Figure 12 for the specified constants. Apparently, the local heat transfer first increases at low curvatures, and the Nusselt number (Nu) changes slightly at higher curvatures. The increased Nu can be attributed to the higher fluid velocity, particularly at lower fin curvatures. The final reduction of the local Nu at the top of the wall is the result of the vortex failing to reach the sharp corners of the cavity, which reduces fluid-wall collisions and, consequently, HTR. Nevertheless, the high-temperature gradient in this area still allows for a considerable HTR at the upper corner of the cavity. The average Nu and GOE are presented in Table 2. Increasing the curvature, at first, enhances HTR and GOE. However, at β = 0.48, increasing the curvature of the two fins The local Nu is plotted in Figure 12 for the specified constants. Apparently, the local heat transfer first increases at low curvatures, and the Nusselt number (Nu) changes slightly at higher curvatures. The increased Nu can be attributed to the higher fluid velocity, particularly at lower fin curvatures. The final reduction of the local Nu at the top of the wall is the result of the vortex failing to reach the sharp corners of the cavity, which reduces fluid-wall collisions and, consequently, HTR. Nevertheless, the high-temperature gradient in this area still allows for a considerable HTR at the upper corner of the cavity. The average Nu and GOE are presented in Table 2. Increasing the curvature, at first, enhances HTR and GOE. However, at β = 0.48, increasing the curvature of the two fins

Effect of Cavity Inclination
The flow field is depicted in Figure 13 at different cavity inclinations for Ra = 10 5 , Ha = 20, ϕ = 0.03 with (a) straight fins, (b) angular fins and (c) curved fins. Figure 13a shows that the angle of the cavity causes a large change in the strength of the vortex as well as the number of vortices. Comparing Figure 13a-c, it can be seen that the maximum stream function for angular and curved fins has greater values than the straight fin for any different angles of the cavity. In particular, the angular fin has a higher stream function value than the straight fin. It can also be seen that in the angular fin, the vortex has more penetration in the lower part of the fins and gets closer to the side wall. Additionally, for cases of curved and angular fins, less flow penetrates above the fin, so in this area, less fluid hits the hot wall with the upper wall of the fins compared to the straight fin. Changing the cavity angle can have a significant effect on the strength of the vortex because the position of the hot and cold walls changes relative to the gravity vector. In the horizontal cavity, the vortex strength is low, which increases with increasing angle. As when the cavity is horizontal (zero angle), the buoyancy force is less, so the amount of fluid displacement is less.
13a shows that the angle of the cavity causes a large change in the strength of the vortex as well as the number of vortices. Comparing Figure13a, b and c, it can be seen that the maximum stream function for angular and curved fins has greater values than the straight fin for any different angles of the cavity. In particular, the angular fin has a higher stream function value than the straight fin. It can also be seen that in the angular fin, the vortex has more penetration in the lower part of the fins and gets closer to the side wall. Additionally, for cases of curved and angular fins, less flow penetrates above the fin, so in this area, less fluid hits the hot wall with the upper wall of the fins compared to the straight fin. Changing the cavity angle can have a significant effect on the strength of the vortex because the position of the hot and cold walls changes relative to the gravity vector. In the horizontal cavity, the vortex strength is low, which increases with increasing angle. As when the cavity is horizontal (zero angle), the buoyancy force is less, so the amount of fluid displacement is less. The temperature field is depicted in Figure 14 at different cavity inclinations for = 10 , = 20, = 0.03 with (a) straight fins, (b) angular fins and (c) curved fins.
Comparing Figure 14a,c, it can be seen that for the horizontal cavity in cases of angular and curved fins, the density of isotherms is higher than that of the straight fin at the bottom of the fins. In this area, due to higher penetration, the amount of temperature change in the angular and curved fins is more than the straight fin. However, in these two fins, the flow penetration in the middle of the fins is less than the straight fin, and fewer isotherms are seen in this area. In the vertical cavity, the penetration of isotherms between the angular and curved fins is greater than in the straight case. In fact, for cases of angular and curved fins, the fluid moves from the top of the first fin to the cold wall, while for the case of straight fins, the fluid moves through the two fins to the cold wall. By changing the angle of the cavity, the place of density of the isotherms also changes and their curvature is in the direction of fluid motion. In the horizontal cavity, the fluid moves along the walls, resulting in greater curvature of the isotherms near the walls. In the vertical cavity, in all three cases, the fluid rises from the middle of the cavity and descends from its sides. Hence, the curvature of isotherms in the middle of the cavity is greater than that The temperature field is depicted in Figure 14 at different cavity inclinations for Ra = 10 5 , Ha = 20, ϕ = 0.03 with (a) straight fins, (b) angular fins and (c) curved fins. Comparing Figure 14a,c, it can be seen that for the horizontal cavity in cases of angular and curved fins, the density of isotherms is higher than that of the straight fin at the bottom of the fins. In this area, due to higher penetration, the amount of temperature change in the angular and curved fins is more than the straight fin. However, in these two fins, the flow penetration in the middle of the fins is less than the straight fin, and fewer isotherms are seen in this area. In the vertical cavity, the penetration of isotherms between the angular and curved fins is greater than in the straight case. In fact, for cases of angular and curved fins, the fluid moves from the top of the first fin to the cold wall, while for the case of straight fins, the fluid moves through the two fins to the cold wall. By changing the angle of the cavity, the place of density of the isotherms also changes and their curvature is in the direction of fluid motion. In the horizontal cavity, the fluid moves along the walls, resulting in greater curvature of the isotherms near the walls. In the vertical cavity, in all three cases, the fluid rises from the middle of the cavity and descends from its sides. Hence, the curvature of isotherms in the middle of the cavity is greater than that in other areas. Additionally, in the vertical cavity, the highest density of isotherms is in the middle of the cold wall, while in other cases, this density was seen on the upper part of the cold wall. As in the vertical cavity (at an angle of 90), the temperature gradient is high in the middle of the wall, but in other cases it is higher in the upper areas of the wall. In the vertical cavity, isotherms are also seen in the direction of straight fins, but in the other two cases, isotherms are seen asymmetrically due to the inclination of the fins to one side. The entropy field is depicted in Figure 15 at different cavity inclinations for = 10 , = 20, = 0.03 with (a) straight fins, (b) angular fins and (c) curved fins.
Comparing Figure 15a to 15c it is observed that the highest local entropy production occurs at the corners of the straight fin, while in the angular fin, the highest local entropy production is observed on the lower wall of the lower fin. In Figure 15c, which is related to the curved fin, the highest entropy production is observed on the lower wall and the vertical wall of the lower fin. In these places, due to the high temperature gradient, the amount of entropy production is also high. In the vertical cavity, the maximum local entropy production occurs in the middle of the cold wall, while for other angles, the maximum occurs in the upper region of the cold wall. For straight fins in the vertical cavity, entropy production lines can also be seen to be symmetrical with respect to the middle line of the cavity, while in curved and angled fins, the lines are different and asymmetric in different regions. The entropy field is depicted in Figure 15 at different cavity inclinations for Ra = 10 5 , Ha = 20, ϕ = 0.03 with (a) straight fins, (b) angular fins and (c) curved fins. Comparing Figure 15a to 15c it is observed that the highest local entropy production occurs at the corners of the straight fin, while in the angular fin, the highest local entropy production is observed on the lower wall of the lower fin. In Figure 15c, which is related to the curved fin, the highest entropy production is observed on the lower wall and the vertical wall of the lower fin. In these places, due to the high temperature gradient, the amount of entropy production is also high. In the vertical cavity, the maximum local entropy production occurs in the middle of the cold wall, while for other angles, the maximum occurs in the upper region of the cold wall. For straight fins in the vertical cavity, entropy production lines can also be seen to be symmetrical with respect to the middle line of the cavity, while in curved and angled fins, the lines are different and asymmetric in different regions. The dimensionless horizontal velocity is plotted in Figure 16 over the line X = 0.5 for different cavity inclinations and = 10 , = 20, = 0.03 with (a) straight fins (b) inclined fins, and (c) curved fins. The velocity plots corresponding to the curved and inclined fins are similar but differ from those for the straight fins. In all three cases, higher velocity maxima and minima were achieved at a cavity angle of 90°. The velocity increases at higher inclinations-particularly in the middle of the cavity. For the straight fin, the velocity minima near the walls are also higher. Meanwhile, for the other two fin configurations, the flow direction changes near the top of the cavity. As a result, near the top of the vertical cavity, the velocity is minimum instead of maximum, suggesting the reversal of the direction of velocity. The dimensionless horizontal velocity is plotted in Figure 16 over the line X = 0.5 for different cavity inclinations and Ra = 10 5 , Ha = 20, ϕ = 0.03 with (a) straight fins (b) inclined fins, and (c) curved fins. The velocity plots corresponding to the curved and inclined fins are similar but differ from those for the straight fins. In all three cases, higher velocity maxima and minima were achieved at a cavity angle of 90 • . The velocity increases at higher inclinations-particularly in the middle of the cavity. For the straight fin, the velocity minima near the walls are also higher. Meanwhile, for the other two fin configurations, the flow direction changes near the top of the cavity. As a result, near the top of the vertical cavity, the velocity is minimum instead of maximum, suggesting the reversal of the direction of velocity. The dimensionless temperature is plotted in Figure 17 for the specified conditions. The dimensionless temperature plots also show similar trends for the curved and inclined fins that are different from those for the straight fin. Similarly, the temperature rises considerably in the middle of the vertical cavity. The hot fluid flow through the middle of the cavity instead of the sides causes an increase in the temperature in the middle and a significant reduction in the temperature near the upper parts since the cold fluid moves toward the hot wall near the top of the vertical cavity. Increasing the inclination of the cavity from 0° increases the temperature in its lower parts, indicating a higher fluid temperature near the bottom wall. The local Nu is plotted in Figure 18 for the specified conditions. All local Nu plots are similar, showing smaller values at the bottom of the wall and higher ones on top for all cavity inclinations, except for 90°, due to the considerable or small temperature difference in that area. Contrary to the other cases, the Nu is higher in the middle of the wall and lower at its beginning or end of the vertical cavity. The temperature difference and, consequently, HTR are enhanced in this case due to the hot fluid coming in contact with the middle parts of the wall. Increasing the inclination, overall, improved HTR at the bottom of the wall and reduced it on top. These changes can be attributed to temperature difference variations in different cases that affect the local Nu. The dimensionless temperature is plotted in Figure 17 for the specified conditions. The dimensionless temperature plots also show similar trends for the curved and inclined fins that are different from those for the straight fin. Similarly, the temperature rises considerably in the middle of the vertical cavity. The hot fluid flow through the middle of the cavity instead of the sides causes an increase in the temperature in the middle and a significant reduction in the temperature near the upper parts since the cold fluid moves toward the hot wall near the top of the vertical cavity. Increasing the inclination of the cavity from 0 • increases the temperature in its lower parts, indicating a higher fluid temperature near the bottom wall. The dimensionless temperature is plotted in Figure 17 for the specified conditions. The dimensionless temperature plots also show similar trends for the curved and inclined fins that are different from those for the straight fin. Similarly, the temperature rises considerably in the middle of the vertical cavity. The hot fluid flow through the middle of the cavity instead of the sides causes an increase in the temperature in the middle and a significant reduction in the temperature near the upper parts since the cold fluid moves toward the hot wall near the top of the vertical cavity. Increasing the inclination of the cavity from 0° increases the temperature in its lower parts, indicating a higher fluid temperature near the bottom wall. The local Nu is plotted in Figure 18 for the specified conditions. All local Nu plots are similar, showing smaller values at the bottom of the wall and higher ones on top for all cavity inclinations, except for 90°, due to the considerable or small temperature difference in that area. Contrary to the other cases, the Nu is higher in the middle of the wall and lower at its beginning or end of the vertical cavity. The temperature difference and, consequently, HTR are enhanced in this case due to the hot fluid coming in contact with the middle parts of the wall. Increasing the inclination, overall, improved HTR at the bottom of the wall and reduced it on top. These changes can be attributed to temperature difference variations in different cases that affect the local Nu. The local Nu is plotted in Figure 18 for the specified conditions. All local Nu plots are similar, showing smaller values at the bottom of the wall and higher ones on top for all cavity inclinations, except for 90 • , due to the considerable or small temperature difference in that area. Contrary to the other cases, the Nu is higher in the middle of the wall and lower at its beginning or end of the vertical cavity. The temperature difference and, consequently, HTR are enhanced in this case due to the hot fluid coming in contact with the middle parts of the wall. Increasing the inclination, overall, improved HTR at the bottom of the wall and reduced it on top. These changes can be attributed to temperature difference variations in different cases that affect the local Nu.
GOE and HTR are plotted in Figure 19 for the specified conditions. As clearly seen, the maximum HTR is achieved in the inclined cavity using curved and inclined fins, while HTR is maximized with straight fins in the vertical cavity. Apparently, the best configuration for heat transfer through curved fins is one with a curved cavity. According to the entropy diagram, the inclined cavity offers the highest GOE in all cases, and the maximum GOE is achieved with curved fins. GOE is drastically reduced at a cavity angle of 90 • as a result of splitting of the vortex into two vortices. This, in turn, affects the velocity gradient, limiting GOE. GOE and HTR are plotted in Figure 19 for the specified conditions. As clearly seen, the maximum HTR is achieved in the inclined cavity using curved and inclined fins, while HTR is maximized with straight fins in the vertical cavity. Apparently, the best configuration for heat transfer through curved fins is one with a curved cavity. According to the entropy diagram, the inclined cavity offers the highest GOE in all cases, and the maximum GOE is achieved with curved fins. GOE is drastically reduced at a cavity angle of 90° as a result of splitting of the vortex into two vortices. This, in turn, affects the velocity gradient, limiting GOE.     GOE and HTR are plotted in Figure 19 for the specified conditions. As clearly seen, the maximum HTR is achieved in the inclined cavity using curved and inclined fins, while HTR is maximized with straight fins in the vertical cavity. Apparently, the best configuration for heat transfer through curved fins is one with a curved cavity. According to the entropy diagram, the inclined cavity offers the highest GOE in all cases, and the maximum GOE is achieved with curved fins. GOE is drastically reduced at a cavity angle of 90° as a result of splitting of the vortex into two vortices. This, in turn, affects the velocity gradient, limiting GOE.      Figure 20 shows the (a) dimensional vertical velocity and (b) dimensional temperature at Y = 0.5 for Ra = 10 5 , γ = 0, ϕ = 0.03 and straight fins for different Ha. Figure 20a shows that increasing the Hartmann number reduces the velocity values in the cavity. Increasing the Hartmann number increases the Lorentz force and thus decreases the strength of the vortices. Hence, the amount of velocity is reduced. As the velocity decreases, it can be seen that the amount of the velocity is tending to zero. Figure 20b also shows that increasing the Hartmann number reduces the curvature of the temperature graph due to the smaller vortex, which causes the temperature changes in the cavity to be less than before. In this case, the temperature graph is becoming linear. While the temperature lines are linear, the conduction heat transfer in the cavity is predominant.

Max NuM
Min NuM Max SG,T Min SG,T Straight Fins 2.42 ( = 90) 2.01 ( = 00) 7.28 ( = 90) 5.80 ( = 00) Angular Fins 2.43 ( = 45) 2.14 ( = 00) 8.12 ( = 45) 6.24 ( = 45) Curve Fins 2.46 ( = 45) 2.11 ( = 00) 8.04 ( = 60) 6.14 ( = 60) Figure 20 shows the (a) dimensional vertical velocity and (b) dimensional temperature at Y = 0.5 for = 10 , = 0, = 0.03 and straight fins for different Ha. Figure 20a shows that increasing the Hartmann number reduces the velocity values in the cavity. Increasing the Hartmann number increases the Lorentz force and thus decreases the strength of the vortices. Hence, the amount of velocity is reduced. As the velocity decreases, it can be seen that the amount of the velocity is tending to zero. Figure 20b also shows that increasing the Hartmann number reduces the curvature of the temperature graph due to the smaller vortex, which causes the temperature changes in the cavity to be less than before. In this case, the temperature graph is becoming linear. While the temperature lines are linear, the conduction heat transfer in the cavity is predominant.  Table 4 reports the effect of the magnetic field on the GOE and the Nu for the straight fin case at different Hartmann numbers and the volume fraction of nanoparticles of 3% for the horizontal cavity. It can be seen that increasing the Hartmann number has reduced both the values of the Nu and GOE, so that reaching the Hartmann number to 40 has reduced the value of the Nu by 37.9% and the GOE by 33.8%. The reduction in vortex velocity has been the main reason for this decrease because it has reduced both the amounts of HTR and GOE by velocity changes.   Table 4 reports the effect of the magnetic field on the GOE and the Nu for the straight fin case at different Hartmann numbers and the volume fraction of nanoparticles of 3% for the horizontal cavity. It can be seen that increasing the Hartmann number has reduced both the values of the Nu and GOE, so that reaching the Hartmann number to 40 has reduced the value of the Nu by 37.9% and the GOE by 33.8%. The reduction in vortex velocity has been the main reason for this decrease because it has reduced both the amounts of HTR and GOE by velocity changes. Finally, by comparing the present results in this work with the same work as in the past [25], it can be seen that by changing the conditions, the values of heat transfer in the cavity can be improved compared to the reference [25]. Alnaghi et al. [25] have studied a cavity equipped with two straight fins on a hot wall and have also considered the effects of a magnetic field. Compared to this work, it can be seen that by increasing the angle of the fin, the amount of heat transfer can be increased by 6.4% of the horizontal fin ratio. It is also possible to increase the amount of heat transfer by 4.9% compared to this article by bending the fins.

Conclusions
In this paper, the effect of installing three different types of fins on the wall of a square cavity on free convective heat transfer has been studied. These three types of fins include straight fins, inclined fins and curved fins. The cavity is affected by the magnetic field. By changing parameters such as fin angle, fin curvature, cavity angle and Hartmann number, the following results are obtained: