Entropy Optimization and Slant MHD Mixed Convection Hybrid Nanofluid Flow Within an Oblique Irregular Lid-Driven Enclosure Contains Baffles: Local Thermal Non-equilibrium Model

This paper explores the heat transfer by mixed convection within double lid-driven enclosures filled with Cu–Al2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}O3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document} hybrid nanofluids. The domain contains T-shaped baffles and is filled with local thermal non-equilibrium (LTNE) porous medium. The wavy boundaries are partially cooled, while the inner baffles have constant heat flux conditions. Mathematical formulations for the system irreversibility in cases of LTNE, constant heat flux conditions and double diffusive are presented and analyzed. The solution methodology is depending on the finite volume scheme in case of non-orthogonal grids. The major results revealed that the alteration of the undulation number λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda $$\end{document} from 1 to 5 gives an augmentation in values of θf,max\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _\text {f,max}$$\end{document} up to 4.57%. Also, the increase in the lengths of the baffles causes a reduction in the flow, while θf,max\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _\text {f,max}$$\end{document} is rising. Furthermore, the alteration in Ri from 0.5 to 10 gives an augmentation in θf,max\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _\text {f,max}$$\end{document} up to 12.97%.


Introduction
Recently, a dependable technique to promote the thermal conductivity of conventional liquids such as ethylene glycol, water, and engine oil has been to load nanometer-sized elements of metals, oxides, carbides, nitrides or nanotubes in the regular liquid. Miscellaneous applications of nanoliquid are observed in nanotechnology, particularly in thermal engineering, including the automobile, microfluidics, electronics applications, power stations, microelectronics and processing plants to obtain a highly effective heat transmission performance [1][2][3][4][5][6][7]. Further efforts of scientists and researchers on nanoliquids led to a new sort, namely hybrid nanofluids. This type of fluids have better thermal and heat transmission, rheological, and synergistic features than the traditional mono nanoliquids. Hybrid nanofluids represent an expansion of nanofluids that involve two or more variant sorts of solid nanometer-sized particles merged together in a pure fluid. Now, that hybrid nanoliquid has become the primary desired topic among the scientists to realize the mechanism behind such fluid well. Devi and Devi [8] reported a comparison between hybrid and traditional nanoliquids in terms of the heat transmission attitude. The authors detected that the hybrid nanofluid's heat transmission rate is greater than the regular nanofluid, even in the existence of a magnetic field. Mahdy et al. [9] simulated numerically free-forced convection flow of a micropolar hybrid nanofluid with an oblique magnetic field in a lid-driven of odd-shaped porous enclosure with radiation impact and discrete heating. A numerical treatment of flow and heat transfer of Ag-MgO-H 2 O hybrid nanoliquid along a slant permeable stretching shrinking plate was exhibited by Anuar et al. [10]. As observed if the angle of slant and suction factors rise, the gradient of temperature boosts. The hybrid nanoliquid flow in a differentially heated containers loaded by aluminum metal foam or glass spheres as porous media was inspected by Mehryan et al. [11]. The utilizing regular liquid was H 2 O, and the hybrid nanoelements were Cu-Al 2 O 3 . Their outcomes illustrated that the nanoparticles in porous media lessen the rate of heat transfer. MWCNT-Fe 3 O 4 -water hybrid nanometer-sized particles have synthesized by Sundar et al. [12], whereas a Ti 2 O-Cu hybrid nanometer elements have been formed by Madhesh et al. [13]. TiO 2 -Cu nanosize particles were suspended in water as a base fluid to be a hybrid nanofluid. Radiative attitude of nanoelements with various shapes of nanomaterials inside a permeable domain was simulated using a new technique. To see the characteristics of alumina, water has been developed by Chu et al. [14]. Ramzan et al. [15] introduced an innovative new mathematical pattern to discuss the flow of 3-D tangent hyperbolic nanofluid due to Hall current and Ion slip imposed by a strong magnetic field through a stretching plate. In addition, Chu with other co-authors investigated flow of nanofluid in different cases. [16][17][18][19][20][21].
Fluid flow and heat transmission are irreversible operations, which can be quantified in relationships of entropy optimization. The optimal design of heat transmission operation in variant industries is determined with precision computations of entropy optimization since it illustrates energy losses in a system distinctly. Therefore, entropy optimization is exhibited widely by many authors [22][23][24][25][26][27]. Induced entropy also has been investigated on MHD nanofluid, and Mahmoudi et al. [28] investigated the entropy analysis and improvement of heat transmission in free convection flow using copper CuH 2 O nanofluid in the existence of a constant magnetic field in a two-dimensional trapezoidal cavity. The outcomes illustrate that at Ra = 10 4 and 10 5 , the strength of the Nusselt number due to nanometer elements boosts with the Hartman number, while at higher values of Rayleigh number, a weakness was noticed. Additionally, it was found that the induced entropy reduces as the nanoparticles are present, but the magnetic field generally rises the magnitude of the induced entropy. Entropy generation has been investigated in different cases by many authors [29][30][31][32].
Besides, as the temperature of the fluid phase T f equals that of solid phase T s through the flow in porous media, i.e., (T f = T s ), hence, the model has one energy equation, this case is known as local thermal equilibrium model LTEM, while the local thermal non-equilibrium model LTNEM is observed if there occurs a contrast among the temperature of the fluid and solid phases, i.e., there is a net heat transfer from one phase to the other, and in this situation, the model has two energy equations. A number of scientists and researchers interested in the case of LTEM [33][34][35][36]. Some of investigations have dealt with LTNEM, Izadi et al. [37] imposed the LTNEM to examine the free convection of CuO-H 2 O micropolar nanofluids within an enclosure filled by an isotropic porous medium. The authors detected that an increment in the thermal conductivity ratio weakens the thermal resistance of the mixture, and therefore, the rate of heat transfer is boosted. Ahmed and Abd El-Aziz [38] scrutinized the boundary layer nanoliquid flow past a vertical wavy surface immersed in a porous medium utilizing LTNEM. They concluded that as Nield number increases, both of the fluid and porous phases local Nusselt numbers enhance.
Prediction of heat transmission through irregular shapes and baffles in cavities represents a topic of massive significance, and it is noticed in miscellaneous engineering processes. Adis [39] explored the heat transmission within an enclosure with heated base attached to fins to illustrate the impacts of configuration parameters of the fins on the heat transmission of laminar flow. Numerically, Ahmed et al. [40] examined the baffles in exhaust mufflers to improve their transmission loss. The calculations clarify that the temperature of outlet gases reduces by 15% with the change in baffle cut ratio from 75 to 25%. Eva et al. [41] scrutinized combined conduction, convection and radiation heat transmission in a closed horizontal airy enclosure, where its thickness and inner surface emissivity are tested with unchanged heat flux. The authors depicted that the rise in cavity thickness reduces the heat flux. The heat transmission is improved by around 120% utilizing short, thin fins compared with that of cavity with absence of fins. The heating block position within a cavity is filled with nanoliquids. The computations clarify that the solid volume fraction of nanofluid and Rayleigh number boosts the heat transmission that is presented by Sannad et al. [42]. Muthtamilselvan et al. [43] analyzed free-forced convection flow of Cu-H 2 O nanoliquid in a lid-driven container for variant aspect ratios employing the control volume approach. Finite difference technique for analyzing the natural convection in a partially heated wavy container loaded by nanofluid is introduced by Pop et al. [44]. Abundant published papers depend on simple or irregular enclosures [45][46][47][48].
In the presented contribution, our emphasis is to formulate the mathematical modeling of entropy optimization through double diffusive mixed convection hybrid nanofluid flow and heat transmission attitude in a slant enclosure of wavy walls containing baffles placed in its middle. A slant magnetic field with of inclination is applied. The LTNEM is imposed between the phases of hybrid nanofluid and porous medium, while the LTEM between the regular fluid (water) and the nanometer particles (alumina and copper) is considered. The impacts of dimensionless factors are interpreted graphically and numerically. The main object of this investigation is to answer the subsequent core questions:

Physical Configuration
A two-dimensional wavy square enclosure of length H is taken for this investigation with inclination α which filled by an anisotropic porous medium. A steady mixed convection heat transfer of hybrid nanofluid with heat flux is examined with entropy analysis. The bottom wall of the cavity is assumed to be the x-axis with T c temperature, while the left vertical wall is considered y-axis. A part of the wavy wall is kept at temperature T c and the other parts of the wavy are adiabatic. A heated T-shaped baffles are placed in the middle of the cavity as observed in Fig. 1. The top and bottom walls move with constant velocity λ d U 0 . MHD is imposed through our contribution to control flow changes in a free-forced convection system. A slant magnetic field of intensity B 0 is imposed on left wall of the container with an inclination angle along the positive horizontal trend, whereas the gravity force direction is downward. Hybrid nanoparticles of copper (Cu) and aluminum dioxide (Al 2 O 3 ) are immersed in the cavity, and the hybrid nanofluid is simulated as a single-phase homogeneous fluid. The fluid's Joule heating and viscous dissipation impacts are however ignored. Two energy equations are formulated, namely hybrid nanofluid and porous phase energy equations (LTEM). Table 1 presents the thermophysical properties of the base fluid and the nanoparticles. Under all these mentioned assumptions, the coupled boundary layer continuity, momentum and energy equations, including the boundary conditions, are formulated as [11,27,36]: in where, the proper imposed conditions for the porous wavy enclosure are On the left side On the right side where notations appear in governing flow Eqs. (1)-(5) are defined as x and y stand for Cartesian coordinates scaled along the vertical and horizontal walls of the container, respectively, u and v symbolize the hybrid nanofluid velocities along the x-and y-axes, respectively, T denotes the hybrid nanofluid temperature,P indicates the hybrid nanofluid pressure, g indicates the gravity acceleration, K means the permeability and Q 0 means the volumetric heat generation/absorption rate. The dependent factors of the Cu, Al 2 O 3 and water hybrid nanoliquid, i.e., viscosity, density, electrical conductivity, thermal conductivity, heat capacity and effective electrical conductivity are stated as (Mahdy et al. [9,48]); where ϕ 1 = ϕ Al 2 O 3 and ϕ 2 = ϕ Cu . To non-dimensionalize the aforementioned governing Eqs. (1)-(6), the following dimensionless variables are inserted which mutate the flow governing equations into the nondimensional subsequent equations in the form Re Pr In the previous non-dimensional equations, the emerging physical factors are expressed mathematically as In addition, the mutated boundary conditions are On the left side On the right side On the top side On the bottom side Due to a science and engineering perspective, the physical quantities of curiosity, such as rate of heat transfer in terms of Nusselt number, have plentiful applications. The mathematical expressions for the physical significant amounts of the local Nusselt number of the fluid phase on the vertical and horizontal obesticals are stated as Nu ssv = 1 θ s v Nu ssh = 1 θ s h (24) and the average Nusselt number of the fluid phase is given by

Entropy Optimization
Physically, entropy is defined as the measure of irreversibility and signifies a state of unrest in the system and its surroundings. If the heat is not utilized for work completely, the amount of entropy it creates is determined. Entropy optimization of the fluid and solid is formulated as [9,24] Due to the dimensionless factors mentioned in Eq. (20), the formulation of the dimensionless entropy optimization may be stated as Herein, the irreversibility rations can be stated by Also, the following Nusselt ration, total entropy generation ration and thermal performance criteria reported by Ilis et al. [26] are defined as In order to provide the impact of nanometer particles, magnetic field and difference of temperature on the average Nusselt number, total induced entropy and Began number, the following Nusselt ration, total induced entropy ration and Began number ration are expressed as Ha=0 , Ha=0 , (34)

Results and Discussion
Finite volume procedure in case of non-orthogonal grids has been imposed to get the numerical computations through this paper.  computed against with the aforementioned values of the controlling parameters. Fig. 2 illustrates the flow features, thermal and concentrations distributions for the alteration of the undulation number λ. The results disclosed that ψ is represented by two forced eddies near the moving walls with non-active middle area due to the presence of the baffles. The alteration of λ causes a weakness of the flow and a reduction in the absolute values of ψ. The mixture temperature θ f shows that the convective mode is dominance on their distributions and the growing in λ causes an enhancement in its maximum temperature θ f,max . However, the solid-phase temperature θ s shows that the conduction mode is dominance and the higher values of λ leads to a clear diminishing in values of θ s,max . Additionally, the features of the mass concentration revealed that it is concentrated near the top-left and bottom-right corners due to the moving of the upper and lower walls; those pushed the mass in opposite direction of the movements. Here, the growing in λ enhances the complexity of the domain, and hence, the suspension flow faces a higher obstruction. Fig. 3 displays the features of ψ, θ f , θ s and φ for the variations in the baffles lengths B. It is remarkable that the increase in B enhances the obstruction of the flow, and hence, a clear reduction in values of ψ is given. Also, the upper eddy is shifted toward the outer boundaries as B is increased. However, a clear augmentation in either values of θ f,max , θ s,max or the distributions of θ f and θ s is seen as B is growing. Physically, the inner baffles have the constant heat flux conditions, and hence, the increase in their lengths enhances the heat flux in the flow domain, and as results, the rising in the isotherms is obtained. On the contrary, features of φ disclosed that the mass concentration is diminishing as B is altered. In the fact, the increase in lengths of the baffles causes a weakness of the buoyancy-induced flow, and hence, the mass concentration is reduced. It is seen that for all values of B, Nu mf is rising as ϕ is altered due to the increase in the mixture thermal conductivity. Also, the heat transfer rate is augmented as the lengths of the baffles are increased due to the increase in the heat flux. On the contrary, a rising in NP concentration causes a reduction in values of Be f due to the decrease in the velocity gradients, while the growing in B makes the heat transfer irreversibility more dominance. Figs. 6, 7 and 8 illustrate the influences of the Richardson number (Ri = Gr/Re 2 ) on the features of the streamlines, isotherms for both the fluid and solid phases and iso-concentration together with the Nusselt number for the fluid and solid phases. The flow features show a reduction in values of ψ min as Ri grows due to the decrease in the forced flow, while the values of θ f,max are rising. Also, the average heat transfer rate for the fluid phase Nu mf is diminishing as either Ri or Ha is altered due to the decrease in the convective flow. Furthermore, the Nusselt number along the vertical baffle is reduced, gradually, as Ri is enhanced due to the decrease in the Reynolds number. In the same context, Figs. 9, 10 display the system entropy for the fluid S f and for the solid S s for the considered range of Ri and ϕ. Two inverse behaviors are noted for S s and S f when Ri is altered. Here, S f is diminishing as Ri is increased due to the invers relation between S f and θ f , while as Ri increases, S s get higher features. Fig. 11 illustrates the contours of ψ, θ f , θ s and φ for the variations in the Darcy number Da. The figure discloses that the flow has a great reduction as Da decreases. In this context, ψ min is varied from −0.0591891 at Da = 0.1 to −0.0016694 at Da = 0.00001. This reduction is due to the decrease in the medium's permeability which makes the flow faces an obstruction due to the medium solid elements. In fact, there is a well-known inverse relation between the fluid speed and θ f,max . Therefore, it can be observed that θ f,max is enhanced, clearly, as Da is altered. Further, the iso-concentration lineaments revealed that the alteration of Da enhances the mass concentration within the flow area.
Figs. 12, 13 and 14 depict the ratio of the fluid Bejan number values in the presence/absence of magnetic field Be + f , ratio of the total entropy to the Nusselt number e + f and values of the Nusselt number at the vertical baffle Nu fsv for the variations in , Ha, H and Q. It is noted that Be + f is reduced as the medium porosity or the Hartmann number is increased due to the reduction in the convective case.
Also, e + f is diminishing as H or Ha is altered due to the Lorentz force. Furthermore, the growing in the heat generation parameter decreases the temperature differences within the flow dominant, and hence, Nu fsv gets lower features.

Conclusions
Conclusion parametric investigations have been carried out for the mixed convection and mass transfer due to the movements of the two-sided edges of two-sided wavy enclosures filled with hybrid nanofluids. Also, T-shaped baffles with constant heat flux conditions were placed in the middle part of the domain. The flow area was filled by LTNE porous elements and an inclined magnetic field together with a heat source/sink has been considered. The finite volume method (FVM) in case of the non-orthogonal grids was applied to solve the mathematical forms, and the irreversibility of the system for the fluid and solid phases was computed. The current analyses led to the following major results: when λ is altered from 1 to 5. 2. The growing in the lengths of the baffles enhances the flow obstruction, and hence, values of ψ min are diminishing, while the average heat transfer rate for the suspension is augmented. 3. The alteration in Ri from 0.5 to 10 gives an augmentation in θ f,max up to 12.97% 4. As Da decreases, the convective transport is reduced while θ f,max is enhanced. 5. The increase in the undulation number λ from 1 to 5 gives a reduction in the flow up to 5.66%. 6. Be + f and e + f are diminishing, clearly, as Ha is growing. 7. The rising in the heat generation parameter causes a reduction in the Nusselt number along the vertical baffles 8. The LTNE model gives more realistic results comparing to LTE cases. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.