Entropy production and mixed convection within trapezoidal cavity having nanofluids and localised solid cylinder

The entropy production and mixed convection within a trapezoidal nanofluid-filled cavity having a localised solid cylinder is numerically examined using the finite element technique. The top horizontal surface moving at a uniform velocity is kept at a cold temperature, while the bottom horizontal surface is thermally activated. The remaining surfaces are maintained adiabatic. Water-based nanofluids (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {Al}_2\text {O}_3$$\end{document}Al2O3 nanoparticles) are used in this study, and the Boussinesq approximation applies. The influence of the Reynolds number, Richardson number, nanoparticles volume fraction, dimensionless radius and location of the solid cylinder on the streamlines, isotherms and isentropic are examined. The results show that the solid cylinder’s size and location are significant control parameters for optimising the heat transfer and the Bejan number inside the trapezoidal cavity. Furthermore, the maximum average Nusselt numbers are obtained for high R values, where the average Nusselt number is increased by 30% when R is raised from 0 to 0.25.

Mixed convective heat transfer inside cavities acts as an extensive rule in various engineering applications such as solar panels, material processors, solar ponds, heat exchangers and many more. Several researchers have investigated the combination of the shear effect and buoyancy force. For instance, mixed convection in square or rectangle cavities have been investigated in several studies by Roy et al. 1 , Abbasian Arani et al. 2 , Khorasanizadeh et al. 3 , Sebdani et al. 4 and Nayak et al. 5 . Additionally, some researchers have performed numerical investigations on mixed convection inside trapezoidal enclosures. Arefmanesh et al. 6 concluded that the nanofluid's effective thermal conductivity for a variable-property form had a significant performance toward the entropy production and heat transfer within the employed enclosure. Bhattacharya et al. 7 stated that the non-isothermal bottom wall plays an important role in multiple steady elements into either natural convection controlled regime or mixed convection regime for convection dominated heat transport regime. Kareem et al. 8 investigated mixed convection in a trapezoidal cavity using different types of nanofluid without an inlet opening site. The results indicated that a wide range of nanofluids had a higher Nusselt number than pure water. Furthermore, the Nusselt number increased as the volume fraction increased, and the rotational angle declined as the measurement of the nanoparticles increased. Selimefendigil and Öztop 9 investigated mixed convection heat transfer into a lid-driven trapezoidal cavity filled with nanofluid and having an inclined magnetic field. The results revealed that the Richardson number, the magnetic field force, and nanoparticle volume fraction were enhanced, while a disparity was observed within the average Nusselt number rises concerning systems by various electrical conductivity criteria. Aghaei et al. 10 examined the impacts of the magnetic field on entropy production and heat transfer of Cu-water nanofluid mixed convection in a trapezoidal cavity. The authors showed that the convection of nanofluids and the intensity of the flow decreased, while the flow-oriented across natural convection and lastly toward pure conduction by introducing and enhancing the magnetic field. In contrast, the effects of entropy production due to friction were insignificant, and the entire entropy production was mainly because of the irreversibility associated with heat transfer. The average Nusselt number grows with the developing Reynolds Scientific Reports | (2021) 11:14700 | https://doi.org/10.1038/s41598-021-94238-w www.nature.com/scientificreports/ number for all Darcy numbers, aspect ratios, and nanoparticle volume fractions 11 . The influence of the cavity's tilt angle on mixed convection and heat transfer within two separate lid-driven trapezoidal enclosures was previously investigated by Hasib et al. 12 . It did note that the features of heat transfer and fluid flow inside the cavities were entirely conditional on Reynolds and Grashof numbers. Chamkha and Ismael 13 explored the mixed convection of the right heat-inclined sidewall within a lid-driven nanofluid-filled trapezoidal cavity with the impact of the uniform magnetic field. The authors showed that the Nusselt number's action varies from the Richardson number associated with the lid's direction. The numerical outcomes of the heat transfer into mixed convection in a liddriven flow inside a trapezoidal cavity with a steady magnetic field were examined by Javed et al. 14 . The findings revealed that the effects of the moving lid became negligible for Ra = 10 5 , although the increasing Rayleigh number led to more excellent streamline circulation and dominant convection effects within the enclosure. While in another investigation, the results indicated that the average Nusselt number decreased as the Lewis number increased, while the total entropy production increased 15 . Al-Rashed et al. 16 examined mixed convection and entropy generation in a cubical heated up by a central isothermal block. Natural convection heat transfer and total entropy generation of nanofluid in a square cavity in presence of a refrigerant rigid body was considered by Rahimi et al. 17 . Recently, Rashidi et al. 18 investigated a three-dimensional free convection and entropy generation of Al 2 O3 water-based nanofluid in a cylindrical cavity. Natural convection and entropy generation in an air-filled cubical cavity with active lateral walls was studied by Alnaqi et al. 19 . Several studies have been performed using either the base fluid or nanofluid or the combination of both fluids. For instance, Alsabery et al. 20 studied natural convective heat transfer in a nanofluid filled-non-Darcian porous and wavy-walled domain under the local thermal non-equilibrium condition. Islam et al. 21 investigated the MHD convection flow within a prism shape cavity filled with Cu-H 2 O nanofluid. Numerical simulations of a hybrid nanofuid flow through a permeable curved enclosure was given by Shah et al. 22 . Gibanov et al. 23 investigated the conjugate mixed convection and entropy production of nanofluid (alumina-water) in a liddriven cavity with a solid bottom surface. It was shown that the bottom wall thickness and volume fraction of nanoparticles perform an important function in enhancing heat transfer. Astanina et al. 24 examined the mixed convection of nanofluids (alumina-water) into a lid-driven cavity containing two porous layers with different thermal properties, permeability, and porosity on the bottom wall. The authors discovered that the extension of nanoparticles contributed to the enhancement of heat transfer in natural convection, while the addition in the nanoparticles volume fraction led toward an apparent decline in heat transfer for mixed and forced convection. Mehmood et al. 25 investigated the mixed convection within nanofluid-filled (alumina-water) lid-driven square cavity containing square blockage and magnetic field. It was observed that as the nanoparticle volume fraction improved, the average Nusselt number and internal energy were increased. Their findings also showed that the rise in the magnetic field directed to reduce the average Nusselt number, and average entropy generation was due to heat transfer. Goodarzi et al. 26 investigated the two-phase mixture rule laminar and turbulent mixed convective heat transfer of water-Cu nanofluids in a shallow rectangular cavity. By increasing the volume fraction of nanoparticles, the coefficient of convective heat transfer increased together with the Nusselt number for specific Grashof and Richardson numbers. The laminar mixed convection flow inside lid (single and double) square cavities filled with nanofluids (alumina-water) was also investigated numerically by Chamkha and Abu-Nada 27 . Due to nanoparticles' existence, it was observed that notable heat transfer enhancement did accomplish, which was strengthened by growing the nanoparticle volume fraction on moderate and high Richardson numbers. Alinia et al. 28 examined the mixed convection of a nanofluid having water and SiO 2 into a tilted enclosure. The authors concluded that nanoparticles' inclusion significantly improved heat transfer in the cavity and induced significant alterations within the flow distribution. Shariat et al. 29 investigated the two-phase mixture model on laminar mixed convection alumina-water nanofluid flow inside elliptic ducts with a fixed heat flux boundary. The results indicated that increasing the nanoparticle volume fraction improved the Nusselt number for a given Reynolds and Richardson number while decreasing the skin friction.
In addition to studies using nanofluids and various types of cavities, studies on heat transfer and fluid flow inside cavities with an inner body have experienced significant awareness during recent years due to their practical engineering purposes. Free convection in cavities having a square inner solid was previously discussed in studies by Alsabery et al. 30 , Alsabery et al. 31 , Mahmoodi and Sebdani 32 , and Sheremet et al. 33 . The influence of nanoparticles' addition on the heat transfer rate was necessary for the low Rayleigh number and number of undulations 34 . Alsabery et al. 35 investigated the time-dependent mixed convection of alumina-water nanoliquid inside the differentially-heated chamber having a wavy upper surface and central solid cylinder. The findings discovered that the average Nusselt number of the hot surface improved on the moving parameter and inner solid cylinder diameter. Liao and Lin 36 investigated mixed convection within the square cavity which contains an isothermally-rotating cylinder. It was observed that the rate of heat transfer depends on the Rayleigh number, aspect ratio, and Prandtl number. The study on the MHD-mixed convection of nanofluids (Cu-water) inside a triangular cavity including a rotating cylinder was conducted by Shariat et al. 37 . It was observed that heat transfer and total entropy production increased as the nanoparticles' solid volume fraction improved. Billah et al. 38 investigated the mixed convection into a lid-driven cavity containing a heated cylinder placed at the centre of the cavity. The results indicated that the flow field and temperature distribution were remarkably dependent on the used cylinder's diameter. Recently, Alhashash 39 employed the nonhomogeneous two-phase Buongiorno's model to study the effect of an inner cylinder on natural convection in a square porous cavity with all walls kept cold. Chatterjee et al. 40 were probably the first to examine the simultaneous effects of a moving lid and an inner rotating cylinder on the mixed convective transport of nanofluids ( Cu-H 2 O ) in a bottom-heated square enclosure. Alsabery et al. 41 considered the effect of an inner rotating cylinder on entropy generation and convective heat transfer in a wavy porous cavity partially-heated from below. The entropy production and mixed convection inside a wavy-walled, nanofluid-filled cavity which is partially-heated from below and contains an inner rotating conductive cylinder was investigated by Alsabery et al. 42  Based on the studies mentioned above and to the best of our knowledge, no investigation has been performed on entropy generation and mixed convection of nanofluids within a lid-driven trapezoidal cavity containing an inner solid cylinder. The heat and fluid flow can be affected significantly by the inclusion of an inner body (like a square or a cylinder) without consuming extra energy. Therefore, this investigation intends to examine the influence of solid cylinder's location on the entropy production and mixed convection inside a lid-driven trapezoidal cavity having water-Al 2 O 3 nanofluids. The outcomes shall be reported for various Reynolds number, Richardson number, nanoparticles volume fraction, the solid cylinder's dimensionless radius and location toward streamlines, isotherms, isentropic, local and average Nusselt numbers, Bejan number, and the global entropy generation.

Mathematical formulation
The two-dimensional geometric design regarding mixed convection flow and heat transfer inside a trapezoidal cavity with a bottom wall of length L and top wall with range L/2, containing an internal solid cylinder with dimensional radius r that is in various locations (D1-D5) is described schematically in Figure 1. The bottom surface is assumed to be heated with a constant temperature of T h , and the top surface slide has uniform velocity from left to right at +U and is preserved at a constant cold temperature, T c . Both sloping surfaces with length 0.65L and inclination angle ϕ are maintained under adiabatic processes. The fluid flow in the trapezoidal cavity was examined to hold steady, laminar, and loaded with alumina-water nanofluids. The governing equations of Navier-Stokes and energy equations toward the viscous incompressible flow are indicated in the dimensional model as the following 43 : The heat equation of the solid inner body remains as: where x and y are the Cartesian coordinates measured in the horizontal and vertical directions respectively, g is the acceleration due to gravity, ρ nf is the density of the nanofluid and ν nf is the kinematic viscosity of the nanofluid.
The following relations can describe the thermophysical features of the applied nanofluids 44 : The following non-dimensional variables are introduced as follows: This then yields the dimensionless governing equations as follows: where Ri = Gr Re 2 is the Richardson number and Gr = is the Grashof number. The dimensionless boundary conditions as indicated in Eqs. (15) and (19) are shown below: On the bottom heated surface(AB) : On the top moving cold surface(DC) : The entropy generation relation as indicated by 41,45 is shown in dimensionless form below: is the irreversibility distribution ratio and S GEN = S gen The terms of Eq. (27) can be separated based on the following formula: where S θ and S are the entropy generation due to heat transfer irreversibility (HTI) and fluid friction irreversibility (FFI), respectively, By integrating Eq. (28) over the domain, the global entropy generation (GEG) for the present two-dimensional study is obtained as follows: It is necessary to mention the Bejan number in determining which is dominant, either heat transfer or fluid friction irreversibility. Bejan number is defined as: The HTI is dominant when Be > 0.5 , while the FFI is dominant when Be < 0.5.

Numerical method and validation
The Galerkin weighted residual with finite element methods was employed to investigate the control equations (15)- (19) subject to the boundary conditions shown in Eqs. (20)- (24). The finite element analysis of the momentum equations (16) and (17) is represented by the following procedure: Firstly, the penalty finite element method was applied by excluding the pressure (P) and including a penalty parameter ( ) as follows: (22) On the left and right sloping surfaces(AD and BC) : θ = θ s , at the outer solid cylinder surface, The weak (or weighted-integral) formulation regarding the momentum equations was obtained by multiplying the equation by an internal domain ( ) and integrating it over the computational domain which is discretised towards small triangular elements as revealed in Fig. 2. The following weak formulations are obtained: The selection of interpolation functions as an approximation towards velocity distribution and temperature distribution is as follows: The non-linear residual equations for the momentum equations obtained from the Galerkin weighted residual finite-element method are as follows: where the superscript k is the relative index, subscripts i, j, and m are the residual number, node number, and iteration number, respectively. For clarifying the non-linear terms into the momentum equations, a Newton-Raphson iteration algorithm was employed. The convergence of the solution was allowed using the relative error to any of the variables that satisfy the resulting convergence criteria as follows: To assure the confidence of the existing numerical solution at the grid dimension of the adopted domain, different grid sizes have been applied to calculate the average Nusselt number and Bejan number concerning the case of Ri = 1 , Re = 100 , φ = 0.02 , R = 0.15 and D1. The outcomes are exhibited in Table 1 which show insignificant variations for the G5 grids and above. Hence, for all calculations in the current problem, the G5 uniform grid remains used.
For data validation, the results obtained in this study were compared with previous numerical findings presented by Khanafer and Aithal 43 for free convective flow and heat transfer in a square cavity filled with pure fluid and partial heating from below as shown, in Figs. 3 and 4. In addition, a comparison was made between the resulting figures and the one provided by Ilis et al. 45 for entropy generation and natural convection in a square cavity fully heated from the side walls as shown in Fig. 5. Figure 6 shows alternative comparisons for the enhancement in the thermal conductivity due to the addition of Al 2 O 3 nanoparticles with two different experimental results and the numerical results obtained by Corcione et al. 46 . These results provide support for the accuracy obtained using the numerical method illustrated in this study.

Results and discussion
The numerical results for the streamlines, isotherms, and isentropic lines (the local dimensionless entropy generation) along with various values for Richardson number (0.01 ≤ Ri ≤ 10) , Reynolds number ( Table 2. Figure 7 shows the distributions of streamlines, isotherms, and isentropic lines for different values of Richardson number for φ = 0.02 , R = 0.15 , Re = 100 , and the inner solid located at the centre of the solid. Therefore, the increase in Ri is due to the increase in the Grashof number. Natural convection is negligible if the Richardson number is less than 1.0. Due to the motion direction of the upper wall, one small cell appeared at the top of the solid wall in Fig. 7(a,b). As Ri increase, the streamlines extend downwards as shown in Fig. 7(c) since both natural and forced convection exist. When Ri increases to 10, multi-cellular patterns appear at the top of the wall as well as the bottom of the wall. Since the bottom wall is heated and the upper wall is cooled, the isotherm patterns appear as curved lines. It can be observed that the curved lines appear at the top left and the bottom right of the solid as in Fig. 7(a). No change in the isotherm patterns was observed as Ri increased. However it was observed that the isotherm curved lines elongated to the right of the solid cylinder at Ri = 10 . The entropy regions are produced by the heat transfer irreversibility (HTI) and the nanofluid flow irreversibility (NFI) that occur at high temperatures and the velocity gradients. The isentropic lines depict a concentrated entropy generation close to three regions, namely at the top and right bottom of the enclosure and the surface of the solid cylinder. No change was observed for the isentropic lines at a smaller Ri number. However, the isentropic lines changed slightly at Ri = 10 as the lines appeared on the bottom (centre) of the solid and the sloping surfaces of the wall. Figure 8 shows the effect of various Reynold numbers on streamlines, isotherm, and isentropic with Ri = 1 , φ = 0.02 , and R = 0.15 . In Figure 8(a), the streamlines appear with two rotating cells at the top of the wall and one rotating cell at the bottom of the wall. This is due to the shear force and buoyant force that dominate at a smaller Reynolds number.As Re increases, the cell at the bottom of the wall and the top right shrink, while the cell at the top left expands as shown in Figure 8(b). Furthermore, the cells at the top right and bottom do not exist, while the cell at the top left expands at higher Re as shown in Figure 8(c) and (d). This is because the velocity of the top lid increases as the Reynolds number increases. Since the velocity of the top lid has increased, forced convection is more dominant as compared to the buoyant force. Since the cooling activity occurs at the top wall and the heating activity occurs at the bottom wall, the thick thermal boundary layers are clustered at the top left and bottom right of the cavity wall. The isotherms are approximately horizontal. The increase in Re disrupts the isotherm inside the cavity. At Re = 50 , it can be observed that the intensity of the isotherm patterns increased at the top left and bottom right walls and decreased at the right side of the cavity. Furthermore, the intensity of the isotherm pattern decreased in the cavity as Re increases. The isentropic lines indicate a concentrated entropy generation inside the cavity and around the solid cylinder. These entropy regions result from heat transfer irreversibility (HTI) and nanofluid irreversibility (NFI) that arise from shear force and buoyant force. As Re increases, the intensity of isentropic lines decreases, especially around the solid cylinder. At Re = 500 , the isentropic lines only appear at the top and bottom of the wall as the intensity of isentropic lines was high at the top of the wall due to the dominant shear force as Re increased. Figure 9 depicts the local Nusselt number interface with X and the local entropy generation at the outer cylinder surface for different Richardson numbers. The local Nusselt number was maximum for all Richardson number at X = 0.4 . The highest local Nusselt number was obtained when Ri = 10 as the system became dominant with natural convection. At the lowest Richardson number, the local entropy generation was highest at both the vertical walls. Figure 10 depicts the distribution of the local Nusselt number, in which the maximum heat transfer was exhibited at X = 0.4 for all the Reynolds number. A higher Reynolds number provides maximum heat transfer as compared to a smaller Reynolds number. Figure 11 illustrates the effect of different nanoparticle volume fractions on the average Nusselt number, Bejan number, and global entropy generation with the Richardson number at Re = 100 , R = 0.15 , and the location of the inner solid at the centre of the enclosure. As Ri increases, the average Nusselt number also increases for all www.nature.com/scientificreports/ the nanoparticle volume fractions. It was observed that the best performance of the average Nusselt number was at higher nanoparticle volume fractions. Since the value of Be was close to 1, the heat transfer irreversibility was dominant as compared to nanofluid irreversibility in the system. The value of Be was shown to decreases for all nanoparticle volume fractions. The system is considered to have dominant forced convection if the Richardson number is small. As the system changes to natural convection, the GEG was observed to increase. A higher nanoparticle volume fraction provided the best results for global entropy generation as compared to a smaller nanoparticle volume fraction. Figure 12 depicts the effect of different nanoparticle volume fractions on the average Nusselt number, Bejan number, and the global entropy generation with the Reynolds number at Ri = 1 , R = 0.15 , and the location of the inner solid at the centre of the enclosure. The increase in Re is indicative of the increase in the velocity of the top lid. As Re increases, the average Nusselt number and GEG increases. The change in the nanoparticle volume fraction also has an effect on the average Nusselt number, Bejan number, and global entropy generation. The best performance for the average Nusselt number and GEG was observed when the nanoparticle volume fraction was at 0.04. From the graph in Fig. 12(a), it was observed that the system was HTI dominant for all nanoparticle volume fractions. As Re increased to 50, a corresponding increase was observed in the graph. However, a decrease was observed in the graph when Re was greater than 50, possibly due to the combined effects of natural and forced convection.
The effects of various radius sizes of the cylinder on the streamlines, isotherms, and isentropic are depicted in Figure 13 with Ri = 1 , Re = 100 , and φ = 0.02 . As the radius of the cylinder increases, the passage width between the cylinder and the top and bottom wall decreases. Therefore, the size of the cell decreases as the radius of the cylinder increases. It was observed that two cells appeared when the radius of the cylinder was greater than 0.2. As a result, the increase in the solid cylinder showed a high intensity of isotherm lines at the bottom and top of the cavity. For small radius sizes, the intensity of the isentropic was similar to the isentropic lines appearing at the top and bottom of the wall. At a higher radius, the isentropic lines appeared at the solid cylinder as well. Figure 14 depicts the local Nusselt number interfaces with X and the local entropy generation at the outer solid cylinder surface with different radius sizes of the solid cylinder. Figure 14(a) shows that the maximum heat transfer occurs mostly at the centre of the cavity. Figure 14(b) illustrates that the local entropy generation is not altered at a smaller radius of the solid cylinder. The local entropy generation depicts a sinusoidal shape for a higher radius of the solid cylinder.
It was also shown that the radius size changes of the solid cylinder affected the average Nusselt number, Bejan number, and global entropy generation as shown in Figure 15. The average Nusselt number (Figure 15(a)) and global entropy generation (Figure 15(b)) were found to increase when the Re number increased. This observation was probably due to the increase in the velocity of the top lid as Re increases, thus affecting thermal performance. Furthermore, the Bejan numbers were high for all the radius sizes, thus indicating that the system is Figure 6. Comparison of (a) thermal conductivity ratio with results obtained by Chon et al. 47 and Corcione et al. 46 and (b) dynamic viscosity ratio with results obtained by Ho et al. 48 and Corcione et al. 46 .  Figure 15(b). The graph illustrates that the Bejan number showed an increasing trend for all radius sizes except for R = 0.25 up to Re = 50 and decreased from that point onwards. This is due to the increase in fluid friction as Re increases, thus resulting in the decreased Bejan number. Figure 16 illustrates the effect of various positions of the solid cylinder on the streamlines, isotherm, and isentropic evolution based on the Richardson number, Reynolds number, nanoparticle volume fraction and radius of the solid cylinder with values of 1, 100, 0.02, and 0.15, respectively. In Fig. 16(a) and (b), the position of the solid cylinder was kept constant at y = 0.5 , while the value of x was indicated by x = 0.35 and x = 0.65 respectively. In Fig. 16(c) and (d) the x position was kept constant at 0.5, while the y was indicated by y = 0.35 and y = 0.65 , respectively. It can be observed that the streamlines appeared according to the position of the solid cylinder. However, the streamlines also appeared at the top of the enclosure for all positions of the solid cylinder due to the movement of the top lid wall. In addition, there was no streamline protruding in the solid cylinder. Likewise it can be observed that the isotherm lines appeared in the enclosure in a wavy shape. However, the It can also be observed that the isentropic lines appeared around the solid cylinder, in high density when the solid cylinder was at the bottom of the enclosure as shown in Figure 16(d).
Various positions of the solid cylinder based on the local Nusselt number interfaces with X and the local entropy generation at the outer solid cylinder surface is depicted in Figure 17. The local Nusselt number was minimum at the vertical position of the cavity wall and maximum between X = 0.2 and X = 0.6 based on the position of the solid cylinder. The best thermal performance was obtained when the position of the solid cylinder was at the centre of the enclosure. In addition, the local entropy generation at the outer solid cylinder displayed a sinusoidal shape. It was observed that the local entropy generation was high when the position of the solid cylinder was at the bottom of the enclosure. www.nature.com/scientificreports/ Figure 18 illustrates the average Nusselt number, Bejan number, and global entropy generation with nanoparticle volume fractions for different positions of the solid cylinder at Ri = 1 , Re = 1 , and R = 0.15 . Different solid cylinder positions are thought to affect the average Nusselt number, average entropy generation, and global entropy generation. It can be observed that the average Nusselt number and global entropy generation were high when the position of the solid cylinder was at the centre of the enclosure. However, both the average Nusselt number and global entropy generation increased slightly as the nanoparticle volume fraction increased for all positions except when the solid was located at the left side of the enclosure. The graph in Figure 18(b) shows that the system was HTI dominant for all positions of the solid cylinder since the value of Be was greater than 0.5. As

Conclusions
In this study, the finite difference method was used to analyse entropy generation and mixed convection in a trapezoidal cavity with various positions of the solid cylinder filled with Al 2 O 3 -water nanofluid. The detailed computational results for the streamlines, isotherms, and the local entropy generation for different values of Ri, Re, nanoparticle volume fractions, radius, and location of the solid cylinder were investigated. Several important conclusions from the study are shown below: 1. The distinct location of the solid cylinder appears to influence flow behaviour, whereby temperature distribution and local entropy generation are both affected by resistance due to the conductive heat transfer in solid cylinder. 2. Increasing the Richardson number and Reynolds number increases the rate of heat transfer as well as nanoparticles volume fraction due to the increment of buoyancy and viscous forces. 3. The maximum average Nusselt numbers are gained for high R values, where the average Nusselt number is increased by 30% when R is raised from 0 to 0.25. 4. The GEG increases with the rise of Ri and Re, while a counteractive behaviour is observed for Be. The Be number also appears to decrease as the Re and Ri increase. 5. In the case of the intense convection regime (Re), the Bejan number is a decreasing function of the solid cylinder radius, whereas the GEG is an increasing function of the solid cylinder radius of the same regime. 6. The heat transfer rate and GEG are high when the solid cylinder is located at the centre of the cavity, while the Bejan number is high when the solid cylinder is located at the top of the cavity.