Entropy Generation and MHD Convection within an Inclined Trapezoidal Heated by Triangular Fin and Filled by a Variable Porous Media

Analyses of the entropy of a thermal system that consists of an inclined trapezoidal geometry heated by a triangular fin are performed. The domain is filled by variable porosity and permeability porous materials and the working mixture is Al2O3-Cu hybrid nanofluids. The porosity is varied exponentially with the smallest distance to the nearest wall and the permeability is depending on the particle diameter. Because of using the two energy equations model (LTNEM), sources of the entropy are entropy due to the transfer of heat of the fluid phase, entropy due to the fluid friction and entropy due to the porous phase transfer of heat. A computational domain with new coordinates (ξ,η) is created and Finite Volume Method (FVM) in case of the non-orthogonal grids is used to solve the resulting system. Various simulations for different values of the inclination angle, Hartmann number and alumina-copper concentration are carried out and the outcomes are presented in terms of streamlines, temperature, fluid friction entropy and Bejan number. It is remarkable that the increase in the inclination angle causes a diminishing of the heat transfer rate. Additionally, the irreversibility due to the temperature gradients is dominant near the heated fins, regardless of the values of the Hartmann number.


Introduction
Various industrial applications depend on the properties and nature of porous media. Thermal insulation, grain stocking and drying processes, the petroleum reservoir, compacted beds for the chemical industry and geophysical systems are good examples of these applications. This importance of these materials makes researchers focus on understanding their nature and properties [1][2][3][4][5]. In fact, the first category of porous medium investigations is concerned with situations where the properties (porosity, permeability and thermal conductivity) of the medium are constants [6][7][8][9][10]. Moreover, the presence of nanoparticles inside a thermal application solvent is receiving increasing attention as it allows one to control the performance in terms of heat transfer and entropy [11][12][13]. Therefore, researchers are targeting the analysis of situations where the two (porous materials and nanofluids) interact. Recent trends of such studies can be found in Ahmed [14] who discussed nanofluids motion in an inclined geometry filled with anisotropic porous elements. His outcomes mentioned that an increase in nanoparticles concentration by 4% results in an enhancement of the maximum temperature of the liquid by 12.27%. A porous medium with constant properties was considered by Rashed et al. [15] as a flow medium within a rectangular domain filled by nanofluids. The thermal dispersion impacts are of interest in this study. Their findings refer to a support in the values of the Nusselt (Nu) coefficient given as the Brownian parameter is growing. The two-energy model (LTNEM) in case of the constants porosity and permeability is used by Ahmed [16] to survey non-Darcy flow within trapezoidal geometries. The flow is due to the presence of a triangular fin that is located at the bottom. The domain is divided into nanofluid layer and porous layer. The outcomes uncovered that variations of the Neild number from 1 to 1000 results in an enhancement in the Nu coefficient by 95.6% at the value of the thermal conductivity ratio Kr = 0.1. The second category of porous media studies is focused on the variable properties of the medium, i.e., variable porosity, variable permeability and variable thermal conductivity. Examples of these studies are those of Al-Weheibi et al. [17], Abelman et al. [18], Saif et al. [19], Nithiarasu et al. [20], EL-Kabeir et al. [21] and Amiri and Vafai [22]. In all the mentioned studies, the porosity is varied exponentially with the smallest distance to the nearest wall and the permeability depends on the particle diameter.
On the other side, the dynamical behaviors of nanofluids within irregular geometries (trapezoidal, triangle, wavy, etc.) have been examined by various investigators in recent years. Treatments of such domains may be made using the body-fitted method, i.e., converting the complex real physical model into a rectangular computational model. While cavities with regular shapes and fins have been studied extensively in many research works, such as [23][24][25], non-regular geometries have received less focus and the aforementioned method is effective in simulations of the transfer of heat problems concerning them. The real coordinates (x, y) are expressed as functions in new coordinates (ξ, η), and based on these functions, all the partial derivatives of the dependent variables are transformed into the new coordinates [26,27]. Alsabery et al. [28] used a trapezoidal geometry as a flow domain to discuss the transfer of heat by Darcian nanofluid flow. They considered that the domain is inclined and divided into two layers (porous layer and nanofluid layer). They noted a remarkable enhancement in the rate of heat transfer in the case of silverwater mixture. Alsabery et al. [29] discussed a trapezoidal domain that is divided into nanofluid porous layer and non-Newtonian fluid layer. Visualization of the heat lines within the system is of interest. Their results reveal that for fixed values of the power-law index (n = 0.7), activity of the mixture flow is enhanced as the Prandtl number grows. Nanofluids flow within triangular domains heated from inside/outside (two heating modes) is examined by Ahmed et al. [30]. The geometry is filled by copper as nanoparticles and includes heated/cooled fins with a variable length. The outcomes disclose that the boost in the height of the fins augments the transfer of heat rate. Additionally, these irregular geometries with various thermal conditions have been handled by different scientific researchers [31][32][33][34][35]: particularly, magneto-convective flow of nanofluids confined open trapezoidal enclosures by Miroshnichenko et al. [36] and Astanina et al. [37]. Finally, it should be noted that many research articles that deal with different shapes and geometries along with many different nanofluids and hybrid nanofluids for conventional fluids or rarefied gases under the influence of many parameters can be found in [38][39][40][41][42].
As stated previously, using triangular fins within the containers as controlling factors for the nanofluids flows and the variable properties of the porous medium have not been considered widely. Thus, the main objective of this simulation is to study the magnetohydrodynamic hybrid nanofluid flow in trapezoidal containers heated by triangular fins. Unusually, the magnetic field is in a horizontal direction while the flow domain is inclined. The properties of the porous medium, namely, the porosity and permeability, as well as values of the Da, are considered as variables. Entropy analyses in the case of the two-energy model (LTNEM) are examined. The real physical domain is transformed to a rectangular computational domain, then the finite volume method based on non-orthogonal grids is applied to solve the resulting system. The novelty significance of the current study appears in the following points:

•
Using an irregular flow domain heated by irregular fin that was not presented before.

•
Properties of the porous medium are considered variables and this assumption makes the work more attractive to readers.
• Physically, the local thermal non-equilibrium state (two-energy equations model) is more realistic than the local thermal equilibrium case.

•
Most of the published works in this field consider the magnetic field as inclined and the geometry as non-inclined, and hence, the formulation of the governing system when the magnetic force is horizontal and the geometry is inclined is unusual and novel.

•
Analyses of the second law of thermodynamics for such kinds of complex geometries (triangular fin within a trapezoidal enclosure) have not been presented before.

•
The numerical methodology which depends on mapping between the real-irregular domain and regular-rectangular computational domain has received less focus. • Finally, various practical applications for the current simulations can be found, e.g., air conditioning systems in buildings, furnace and home heating, electronic equipment cooling, drying foods and double pane windows.

Formulation of the Problem
A horizontal magnetic field with strength β 0 is influencing on an inclined trapezoidal domain heated by a triangular fin, as illustrated in Figure 1. The inclined walls of the containers are referred to with a low temperature T c while the upper wall is insulated. The inclination angle is γ and the gravity and magnetic vectors are, respectively, (−g sin γ, −g cos γ), (β 0 cos γ, −β 0 sin γ). The trapezoidal angle is Φ and the geometry is filled by a variable properties porous medium. The non-equilibrium thermal model between the hybrid nanofluid and porous phases is valid. The mixture consists of alumina and copper nanoparticles suspended in a base fluid (water). Properties of the nanoparticles and water are given in Table 1. In addition, the porous elements are considered as a compacted spheres bed. Therefore, the model proposed by Ergun [43] for permeability is introduced. One of the features of this model is that the permeability depends on the particle size (particle diameter dp) and the porosity of the medium ε. It is expressed as: Appl. Sci. 2020, 10, x FOR PEER REVIEW 3 of 20 • Physically, the local thermal non-equilibrium state (two-energy equations model) is more realistic than the local thermal equilibrium case. • Most of the published works in this field consider the magnetic field as inclined and the geometry as non-inclined, and hence, the formulation of the governing system when the magnetic force is horizontal and the geometry is inclined is unusual and novel. • Analyses of the second law of thermodynamics for such kinds of complex geometries (triangular fin within a trapezoidal enclosure) have not been presented before. • The numerical methodology which depends on mapping between the real-irregular domain and regular-rectangular computational domain has received less focus. • Finally, various practical applications for the current simulations can be found, e.g., air conditioning systems in buildings, furnace and home heating, electronic equipment cooling, drying foods and double pane windows.

Formulation of the Problem
A horizontal magnetic field with strength is influencing on an inclined trapezoidal domain heated by a triangular fin, as illustrated in Figure 1. The inclined walls of the containers are referred to with a low temperature while the upper wall is insulated. The inclination angle is and the gravity and magnetic vectors are, respectively, (− sin , − cos ), ( cos , − sin ). The trapezoidal angle is Φ and the geometry is filled by a variable properties porous medium. The non-equilibrium thermal model between the hybrid nanofluid and porous phases is valid. The mixture consists of alumina and copper nanoparticles suspended in a base fluid (water). Properties of the nanoparticles and water are given in Table 1. In addition, the porous elements are considered as a compacted spheres bed. Therefore, the model proposed by Ergun [43] for permeability is introduced. One of the features of this model is that the permeability depends on the particle size (particle diameter ) and the porosity of the medium . It is expressed as:    In addition, the variable porosity of the medium is taken as a function of the particle size dp and the smallest distance to the nearest wall |y|. This relation is given as: Values of the empirical constants a 1 and a 2 are introduced by Cheng and Hsu [44] as 1.4 and 5, respectively, and the value of ε ∞ is considered as 0.36. Using this correlation, the porosity near the bottom wall is equal to 0.864 and this value decreases exponentially as the distance from the bottom wall is increased. The Da is also considered a variable in this study.
When all these assumptions are taken into account, the two-energy model governing this physical case is given as: ε ρc p hn f The symbols in Equations (4)-(8) are the velocity fields u p , v p , the temperature for the hybrid nanofluid/solid phase T f p , T p , the porosity ε, the permeability K, the electrical conductivity σ hn f , the thermal conductivity k hn f , the independent variables (Cartesian coordinates and time) x, y, t and the gravity acceleration g. In addition, the following dimensionless quantities are submitted: Substituting Equation (9) into Equations (4)-(8), the following combined dimensionless system is given: The subjected conditions to the boundaries of the geometry are: On the inner fin: In Equations (10)- (15), is the Nield number and is the thermal conductivity ratio.

Correlations of the Hybrid Nanofluids
Correlations of the hybrid nanofluids were presented recently by Ahmed [45]. Following this investigation, the density, heat capacitance, thermal expansion, thermal diffusivity, thermal conductivity, dynamic viscosity and electrical conductivity are, respectively:

Heat Transfer Coefficient
The heat transfer rate is calculated at both the left and right walls. Since these walls are inclined, these quantities are given as: For the fluid phase: For the solid phase:

Entropy Generation Analysis
The entropy equations can be written in the following form: where T 0 = T h +T c 2 , using the dimensionless variables and the characteristics of entropy In the above equation, Θ = is the ratio of the irreversibility distribution. In addition, the local and average Bejan number (Be) is expressed as: Using the dimensionless quantities, the entropy generation is given as: The total dimensionless entropy generation is the sum of these components (Buonomo et al. [46]), and it is written as: The entropy generation due to the entire domain can be obtained by integrating the previous component over the entire domain as: The local Be is the ratio between the local entropy generation due to the transfer of heat (sum of the entropy generation due to the heat transfer for the fluid phase and for the solid phase) and the entropy generation due to the fluid friction; thus, it can be expressed as: Here, it should be mentioned that the ratio of the irreversibility distribution Θ is fixed at 10 −4 for all the computations and this value is selected according to Ilis et al. [47].

Numerical Treatments
The numerical treatments used to solve the governing equations are the non-orthogonal grids FVM. The real physical domain (X, Y) is converted to a rectangular domain (ξ, η) using the following functions: Using (32), the advection terms are given as: where Similar, the diffusive terms are expressed as: In Equation (38), α 11 , α 22 are expressed as: The upwind scheme is used to discretize the transformed convective terms while the diffusive terms are treated using the central differences scheme. The alternating direction implicit method (ADI) is introduced to solve the obtained algebraic system. The iterative technique is employed until the convergence criterion (10 −6 ) is fullfiled. Additionally, after many grid tests being made (as it is presented in Table 2), the grid size 101 × 101 was found suitable for all computations. Further, Figure 2 shows the mesh used in all computations, which consists of 101 × 101 nodes.  Similar, the diffusive terms are expressed as: In Equation (38), , are expressed as: The upwind scheme is used to discretize the transformed convective terms while the diffusive terms are treated using the central differences scheme. The alternating direction implicit method (ADI) is introduced to solve the obtained algebraic system. The iterative technique is employed until the convergence criterion (10 ) is fullfiled. Additionally, after many grid tests being made (as it is presented in Table 2), the grid size 101 × 101 was found suitable for all computations. Further, Figure 2 shows the mesh used in all computations, which consists of 101 × 101 nodes.   Figure 3 illustrates a comparison between the current code results and the results obtained by Ilis et al. [47] for the case of a square cavity in which the left wall is kept at a hot temperature, the right wall is kept at a cold temperature while the lower and upper walls of the cavity are considered adiabatic. The comparison is implemented for the case   [47] for the case of a square cavity in which the left wall is kept at a hot temperature, the right wall is kept at a cold temperature while the lower and upper walls of the cavity are considered adiabatic. The comparison is implemented for the case where Ra = 10 5 , Pr = 0.71, Θ = 10 −4 . Contours for the streamlines, isotherms, local Bejan number and local entropy generation due to heat transfer are constructed. Comparison shows a very good agreement between contours resulted from the current code to those presented in Ilis et al. [47].

Ilis et al. [47]
Streamlines Isotherms Local Bejan number Local entropy generation due to heat transfer   Figure 4 elucidates the streamlines and isotherms of the liquid phase and the isotherms of the porous phase for the cases where Ra = 10 6 , Ha = 15, H = 0.1, φ Al = φ Cu = 1%, Φ = π 3 and different inclination angles. As far as the streamlines are concerned, plots show a formation of two symmetric cells for an inclination angle of zero. The cell on the left hand side starts to dominate the region as the inclination angle increases. This is mainly due to the buoyancy forces. For the isotherms related to the fluid phase, plots show a formation of a hot thermal region near the fin. Moreover, plots show that for an inclination angle of zero, isotherms are symmetric throughout the centerline of the cavity with a maximum area covered by higher values of isotherms. As the angle increases, distortion of the isotherms happens as they shift counter clockwise. Additionally, it is obvious that as the inclination angle increases, less area of the domain is covered by higher values of the isotherms. Finally, the graph shows that for the porous phase, symmetric isotherms are observed for all inclination angles. Furthermore, as the inclination angle increases, the impact on the temperature distribution is very minimal as the porous temperature is not coupled with the velocity fields. show a formation of two symmetric cells for an inclination angle of zero. The cell on the left hand side starts to dominate the region as the inclination angle increases. This is mainly due to the buoyancy forces. For the isotherms related to the fluid phase, plots show a formation of a hot thermal region near the fin. Moreover, plots show that for an inclination angle of zero, isotherms are symmetric throughout the centerline of the cavity with a maximum area covered by higher values of isotherms. As the angle increases, distortion of the isotherms happens as they shift counter clockwise. Additionally, it is obvious that as the inclination angle increases, less area of the domain is covered by higher values of the isotherms. Finally, the graph shows that for the porous phase, symmetric isotherms are observed for all inclination angles. Furthermore, as the inclination angle increases, the impact on the temperature distribution is very minimal as the porous temperature is not coupled with the velocity fields.

Streamlines
Fluid Phase Temperature Solid Phase Temperature    Figure 5 illustrates features of the local entropy generation due to the fluid friction and local Be for the variation of the inclination angle γ (γ = 0, π 6 , π 3 , π 2 ) at Ra = 10 6 , Ha = 15, The graph shows that the S f decreases as the γ increases. The local entropy generation for γ equals zero appears to have the maximum values near the side walls. As the γ increases, this maximum value takes place at the upper wall of the cavity. Moreover, the graph shows that for γ, Be less than 0.5 contours occupies most of the domain of the cavity; this is mainly due to the friction. Additionally, Be contours greater than 0.5 are observed at or near the sidewalls of the cavity as well as the fin walls. For these cases where Be is greater than 0.5, the entropy generation is mainly due to the transfer of heat.
. 2020, 10, x FOR PEER REVIEW 11 of 20 15, = 0.1, = = 1%, Φ = . The graph shows that the Sf decreases as the increases. The local entropy generation for equals zero appears to have the maximum values near the side walls. As the increases, this maximum value takes place at the upper wall of the cavity. Moreover, the graph shows that for , less than 0.5 contours occupies most of the domain of the cavity; this is mainly due to the friction. Additionally, Be contours greater than 0.5 are observed at or near the sidewalls of the cavity as well as the fin walls. For these cases where is greater than 0.5, the entropy generation is mainly due to the transfer of heat.
Local Entropy Generation Due to the Fluid Friction Local Bejan Number        Figure 8. The graph shows that as the decreases, Nu decreases. This is because as the angle increases, Variations of the local Nu for the cases where nanoparticle volume fraction φ Al = φ Cu = 1% at Ha = 15, H = 0.1, γ = π 6 , Φ = π 3 and different Ra are plotted in Figure 8. The graph shows that as the γ decreases, Nu decreases. This is because as the angle increases, the thermal boundary layers near the inclined walls decrease. This will result in a reduction in the transfer of heat rate, and consequently, lower Nu is achieved.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 14 of 20 the thermal boundary layers near the inclined walls decrease. This will result in a reduction in the transfer of heat rate, and consequently, lower Nu is achieved.  The graph shows that as increases, the Nu increases. This is mainly due to the fact that as increases, the dominant mode of transfer of heat becomes convection, and consequently, higher Nu is given. Moreover, the graph shows that as the s increases for the case of = 10 and = 10 , better Nu is achieved. As increases, better effective conductivity is achieved and consequently better transfer of heat. At = 10 , the activity of the hybrid nanofluids flow is much affected by the increase in , . In this case, the overall viscosity is augmented, and thus, the convective process and likewise the Nu are decreased. In Figure 10,   Figure 8. Profiles of the Nu f for the variation of the inclination angle γ (γ = 0, π 6 , π 3 , π 2 ) at Ra = 10 6 , Ha = 15, H = 0.1, φ Al = φ Cu = 1%, Φ = π 3 . Figure 9 presents profiles of the local Nu for the variation of the Ra and nanoparticle volume fraction φ Al = φ Cu at Ha = 15, H = 0.1, γ = π/6, Φ = π 3 . The graph shows that as Ra increases, the Nu increases. This is mainly due to the fact that as Ra increases, the dominant mode of transfer of heat becomes convection, and consequently, higher Nu is given. Moreover, the graph shows that as the s φ increases for the case of Ra = 10 4 and Ra = 10 5 , better Nu is achieved. As φ increases, better effective conductivity is achieved and consequently better transfer of heat. At Ra = 10 6 , the activity of the hybrid nanofluids flow is much affected by the increase in φ Al , φ Cu . In this case, the overall viscosity is augmented, and thus, the convective process and likewise the Nu are decreased.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 14 of 20 the thermal boundary layers near the inclined walls decrease. This will result in a reduction in the transfer of heat rate, and consequently, lower Nu is achieved.  The graph shows that as increases, the Nu increases. This is mainly due to the fact that as increases, the dominant mode of transfer of heat becomes convection, and consequently, higher Nu is given. Moreover, the graph shows that as the s increases for the case of = 10 and = 10 , better Nu is achieved. As increases, better effective conductivity is achieved and consequently better transfer of heat. At = 10 , the activity of the hybrid nanofluids flow is much affected by the increase in , . In this case, the overall viscosity is augmented, and thus, the convective process and likewise the Nu are decreased. In Figure 10,   In Figure 10, profiles of the entropy generation due to the transfer of heat of the hybrid nanofluid phase in the entire domain for the variation of the Ra and nanoparticle volume fraction φ Al = φ Cu at Ha = 15, H = 0.1, γ = π 6 , Φ = π 3 are constructed. The graph shows that as Ra increases, the S T increases. Moreover, the graph shows that as the φ increases, S T increases as well. Higher volume fraction is associated with better conductivity and hence a higher entropy generation due to heat is achieved.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 15 of 20 shows that as increases, the ST increases. Moreover, the graph shows that as the increases, ST increases as well. Higher volume fraction is associated with better conductivity and hence a higher entropy generation due to heat is achieved.  Figure 11 shows profiles of the total entropy generation in the entire domain for the variation of the Ra and nanoparticle volume fraction = at = 15, = 0.1, = /6, Φ = . The graph shows that the Stotal increases as Ra increases. Furthermore, the graph shows for the cases of = 10 and = 10 , as increases, the Stotal increases. Finally, the graph shows that for the case where = 10 , as the increases, the Stotal increases. This is mainly because the entropy generation due to the heat increases with the volume fraction while the entropy generation due to friction decreases.   Figure 11 shows profiles of the total entropy generation in the entire domain for the variation of the Ra and nanoparticle volume fraction φ Al = φ Cu at Ha = 15, H = 0.1, γ = π/6, Φ = π 3 . The graph shows that the S total increases as Ra increases. Furthermore, the graph shows for the cases of Ra = 10 4 and Ra = 10 5 , as φ increases, the S total increases. Finally, the graph shows that for the case where Ra = 10 6 , as the φ increases, the S total increases. This is mainly because the entropy generation due to the heat increases with the volume fraction while the entropy generation due to friction decreases.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 15 of 20 shows that as increases, the ST increases. Moreover, the graph shows that as the increases, ST increases as well. Higher volume fraction is associated with better conductivity and hence a higher entropy generation due to heat is achieved.  Figure 11 shows profiles of the total entropy generation in the entire domain for the variation of the Ra and nanoparticle volume fraction = at = 15, = 0.1, = /6, Φ = . The graph shows that the Stotal increases as Ra increases. Furthermore, the graph shows for the cases of = 10 and = 10 , as increases, the Stotal increases. Finally, the graph shows that for the case where = 10 , as the increases, the Stotal increases. This is mainly because the entropy generation due to the heat increases with the volume fraction while the entropy generation due to friction decreases.   Figure 12 illustrates profiles of the entropy generation due to the transfer of heat of the solid phase in the entire domain for the variation of the Ra and nanoparticle volume fraction φ Al = φ Cu at Ha = 15, H = 0.1, γ = π/6, Φ = π 3 . The graph shows that the S TS increases as Ra increases. Additionally, the graph shows that as the φ increases, the S TS decreases. Physically, the increase in the concentration of the nanoparticles (φ Al , φ Cu ) enhances the temperature gradients for the fluid phase which results in a diminishing in the temperature gradients of the solid phase, and hence S Ts is reduced.  The graph shows that as increases, t Sf increases. Moreover, for low Ra values, the entropy generation is almost constant with the ; while this value decreases as the increases at Ra= 10 6 . In fact, high values of ( = 10 ) make the velocity fields more affected by variations of . In addition, an increase in augments the overall viscosity which in turn decreases the gradients of the velocity, and hence, the total entropy is reduced.

Conclusions
Magneto convection within the trapezoidal enclosures heated by triangular fins and filled by variable properties porous media has been carried out. The Darcy number is varied from a point to another in the flow domain based on a variable porosity of the medium. The magnetic force is taken in the horizontal direction while the geometry is inclined. The LTNEM occurs between the porous medium and the worked fluid. The hybrid nanofluid that consists of water, copper and alumina is assumed to be a worked fluid. The FVM in the case of the non-orthogonal grids is used to solve the governing equations. The main findings of this simulation are: 1 The irreversibility due to the transfer of heat is dominant along the fin boundary in the case of a horizontal domain, while in the case of a vertical geometry, it is dominant near the boundaries of the trapezoidal.   Figure 13 demonstrates profiles of the S f in the entire domain for the variation of the Ra and nanoparticle volume fraction φ Al = φ Cu at Ha = 15, H = 0.1, γ = π/6, Φ = π 3 . The graph shows that as Ra increases, t S f increases. Moreover, for low Ra values, the entropy generation is almost constant with the φ; while this value decreases as the φ increases at Ra= 10 6 . In fact, high values of Ra (Ra = 10 6 ) make the velocity fields more affected by variations of φ. In addition, an increase in φ augments the overall viscosity which in turn decreases the gradients of the velocity, and hence, the total entropy is reduced.  The graph shows that as increases, t Sf increases. Moreover, for low Ra values, the entropy generation is almost constant with the ; while this value decreases as the increases at Ra= 10 6 . In fact, high values of ( = 10 ) make the velocity fields more affected by variations of . In addition, an increase in augments the overall viscosity which in turn decreases the gradients of the velocity, and hence, the total entropy is reduced.

Conclusions
Magneto convection within the trapezoidal enclosures heated by triangular fins and filled by variable properties porous media has been carried out. The Darcy number is varied from a point to another in the flow domain based on a variable porosity of the medium. The magnetic force is taken in the horizontal direction while the geometry is inclined. The LTNEM occurs between the porous medium and the worked fluid. The hybrid nanofluid that consists of water, copper and alumina is assumed to be a worked fluid. The FVM in the case of the non-orthogonal grids is used to solve the governing equations. The main findings of this simulation are: 1 The irreversibility due to the transfer of heat is dominant along the fin boundary in the case of a horizontal domain, while in the case of a vertical geometry, it is dominant near the boundaries of the trapezoidal.

Conclusions
Magneto convection within the trapezoidal enclosures heated by triangular fins and filled by variable properties porous media has been carried out. The Darcy number is varied from a point to another in the flow domain based on a variable porosity of the medium. The magnetic force is taken in the horizontal direction while the geometry is inclined. The LTNEM occurs between the porous medium and the worked fluid. The hybrid nanofluid that consists of water, copper and alumina is assumed to be a worked fluid. The FVM in the case of the non-orthogonal grids is used to solve the governing equations. The main findings of this simulation are: 1 The irreversibility due to the transfer of heat is dominant along the fin boundary in the case of a horizontal domain, while in the case of a vertical geometry, it is dominant near the boundaries of the trapezoidal. 2 Activity of the hybrid nanofluid as well as the fluid friction entropy are diminished as the Ha grows. 3 The transfer of heat entropy is augmented as the concentration of the nanoparticles is boosted. The increase in the γ reduces the Nusselt number while it increases as the Rayleigh number or the volume fraction parameter are growing. 4 From the obtained results, it is recommended to use a non-inclined irregular domain to enhance the heat transfer rate. Additionally, the variable-properties porous medium is more realistic than the porous medium with constant properties. 5 It is recommended to use hybrid nanofluids for enhancement of the heat transfer instead of mono nanofluids.