Magnetohydrodynamic mixed convection of TiO2–Cu/water between the double lid-driven cavity and a central heat source surrounding by a wavy tilted domain of porous medium under local thermal non-equilibrium

The magnetohydrodynamic (MHD) mixed convection of heat and mass transfer is carried out using finite difference method applied inside a tilted porous cavity saturated with a hybrid nanofluid due to the presence of the double-moving lid and the heat sources. In contrast to the earlier research, various effects which are recognized by heat generation in the local thermal non-equilibrium case at the extended Brinkman Darcy model subjected to inclined magnetic field are thoroughly examined numerically. For instance, unusual observations of the cold mass surrounding the heat source emphasize that the maximum fluid temperature highly depends on the forced convection. Additionally, solid-phase temperature acts in accordance to the heat source location while fluid temperature is agitated by the moveable sides which points up the disparity at the thermal energy transportation. However, the transfer of heat and mass at the model requires a specific conduct due to the existence of damping factors. The magnetic field, for example, suppresses the fluid flow. Moreover, the thermal non-equilibrium condition deteriorates the global heat generation. The flow shear at the proposed model is more feasible for the heat and mass transfer than the buoyancy force. Local thermal non-equilibrium case creates inaccessible regions that disapprove the harmony in the flow behavior. The ratio of natural to forced convection is the most pertinent to the thermal nonequilibrium scenario. The flow shear at the proposed model is more feasible for the heat and mass transfer than the buoyancy force. Local thermal non-equilibrium case creates inaccessible regions that disapprove the harmony in the flow behavior. The ratio of natural to forced convection is the most pertinent to the thermal nonequilibrium scenario.


Introduction
Mixed convection with magnetic fields is successfully used in polymers, metallurgy, geothermal energy extraction, and fusion reactors. Due to its properties, it is utilized in a wide range of industries, including mechanical engineering, the petroleum industry, chemical engineering, the medical profession, and other physiological frameworks [1][2][3][4][5]. This idea is especially helpful for applications like magnetic resonance imaging, blood pumping, radiofrequency ablation, cancer tumor treatment, nuclear reactor cooling, magnetic electro-catalysis, metallurgical processes, thermal protection, receiving and transmitting antennas, magnetic electro-catalysis [6][7][8][9][10], etc.
Additionally, obstructions existence at the convection have been compared by Scott et al. [11] to the opening function, where the blockages are like flaps. They split motion into two main categories: bulk density variations and motion pressure gradient, which may both be utilized to explain natural convection flow in multizone enclosures. Partitions should span the ground and roof of a square cavity, according to Chang et al. [12], such that heat transfer reduces as the thickness of the partition climbs. Moreover, Eshaghi et al. [13] have studied the doublediffusive natural convection in an H-shaped cavity with a baffle. They reported that the best Nusselt number and Sherwood number have been grasped. According to Nansteel and Greif [14] and Lin and Bejan [15], impediments between zones tend to keep the boundary layer's thermal equilibrium condition stable, which suggests a general drop in convective heat transmission. As per Sahoo et al. [16], the rectangular blockage at the vertical channel causes a "fluid flow step" that is "backward-facing" and can alter the maximum average temperature. In addition, heat transmission in vertical channels with semi-circular barriers has been studied by Said and Krane [17]. Depending on the Rayleigh number, the obstruction resulted in a drop in the average Nusselt number at a uniform wall temperature (it increases as Ra decreases). In the same context, Seyyedi et al. [18] have investigated heat transfer characteristics as well as entropy generation in a hexagonal enclosure. While, An economic analysis was carried out [19] for heat transfer between two non-similar circular cylinders. However, Dogonchi et al. [20] have proved that the suction process has a dual effect on heat transfer and fluid flow using the Duan-Rach approach.
Shahriari et al. [21] reported their findings in a research on Magnetohydrodynamic natural convection within the angled undulating cavity packed with CuO-water nanofluid. They observed that the magnetic field's inhibitory impact causes the mean Nusselt number and Bejan number to become decreasing functions as the Hartmann number increases, while entropy formation is enhanced. Another illustration is Sheremet et al.'s [22] study of the Magnetohydrodynamic non-driven convection of a nano-fluid in an open porous "high cavity" with a "corner heater". Moreover, Dogonchi et al. [23] reported that "Rayleigh number could overshoot the velocity gradient" for the natural convection at a porous medium saturated between two square cylinders filled with magnetic nanofluids. In addition, Mandal et al. [24] studied the magnetohydrodynamic mixed bioconvection.
On the other hand, the larger amplitude and higher frequency of the moving wall can effectively mix the viscous fluid, according to a significant study by Li et al. [25].
Moreover, Zou et al. [26] concluded that the flow direction should be taken into account to promote convective mixing. To examine the effects of the Lorentz force on the free convection inside the wavy triangular cavity, Shekaramiz et al. [27] carried out a numerical research. They discovered that entropy generation first rises with rising Hartmann numbers before remaining constant.
For instance, in a square cavity containing Cu-water nanofluids, Rashad et al. team 's [28] numerically examined the MHD mixed convection flow. In a porous cavity with a lid with a "Cattaneo-Christov" heat flux pattern by mixed convection, Jakeer et al. [29] numerically investigated the effect of the heated barrier location on the magneto-hybrid nanofluid flow. For MHD natural convection, Sheikholeslami and Sadoughi [30] investigated the "lattice Boltzmann" technique inside a porous cavity that was filled with a nanofluid enclosing that contained 4 square-heaters. In addition, CuO-MWCNT-oil hybrid nanofluid's turbulent flow field, heat transfer, and entropy formation in a trapezoidal container under the influence of a magnetic field were all quantitatively investigated by Aghaei et al. [31] using the finite volume approach. According to their research, the average Nusselt numbers for different Rayleigh numbers and Hartmann numbers decrease with the increased "volume fraction" of nanoparticles. Additionally, the efficiency of heat exchange is dramatically altered by the nanoparticles dispersed, especially those made of metals (Cu, Al 2 O 3 , ZnO, Ag, and SiO 2 ) inside the fluid matrix as compared to regular liquid flows (Khanafer et al., Chamkha et al., and Aly et al., [32][33][34][35]). Tiwari and Das [36] looked into the supplementary effect for a lid-driven cavity, while Billah et al. [37] looked into the buoyancy-driven force at a slanted triangular compartment. Moreover, The Nusselt number greatly change based on the shape of the used nanoparticles [38].
In accordance with the size, form, and composition of the nanoparticles that are dispersed inside them, nanofluids have different thermal conductivities (Suresh et al. [39]). For instance, metallic particles have a "higher thermal conductivity than non-metallic particles". Nevertheless, as a result of their size, cylindrical nanoparticles function better thermally than spherical nanoparticles ones. Consequently, when the atomic structure of the fluid matrix is precisely treated, significant changes take place. For businesses looking for a product that exhibits excellent stability and superior chemical inertness, Al 2 O 3 is the solution. However, this will oftentimes be inadequate for nanofluid applications because they will lose their benefit of heat conductivity. This problem can be handled by adding a few metal particles to an aluminum framework to improve its thermal characteristics. By chemically preparing CuO-Al 2 O 3 combinations and hydrogen-reducing them, Jeena et al. [40] created the Cu-Al 2 O 3 nanocomposites.
Biswas et al. [41] were able to rise the thermal performance at a tilted porous cavity using Cu-Al 2 O 3 -water nanofluid subjected to the partially active magnetic field. Sung-Tag Oh et al. [42] demonstrated the fabrication of an Al 2 O 3 -Cu nanocomposite using powder mixtures of fine Al 2 O 3 and copper oxide sizes. As an Al 2 O 3 -Cu/water hybrid nanofluid flowed through a circular tube with a constant heat flux, Suresh et al. [39] assessed the "Nusselt number and friction factor". Davood Toghraie et al. [43] concluded from their experimental studies on ZnO-TiO 2 /EG hybrid nanofluids that the significant increase in heat transmission occurs either at a warmer temperature proportionate to solid volume fraction or vice versa.
However, it is reported that when carbon nanotubes (CNTs) were attached to the surface of spherical oxide nanoparticles (Huang et al. [44]) then the thermal conductivity was significantly enforced (Han et al. [45]) due to the strong thermal action among the CNTs attached to the alumina/iron oxide particle. On the other hand, nanoencapsulated phase change material has been utilized [46] for free convection in a porous, wavy enclosure. For a better comprehension, Suresh et al. [39] discovered that for a hybrid Al 2 O 3 -Cu/deionized water nanofluid with a "0.1% volume concentration", the Nusselt number "was 13.56% greater than that of water at Reynold number (Re) = 1730". Furthermore, Labib et al. [47] reported that improved heat transmission in mixes of Al 2 O 3 nanoparticles in CNT/ water nanofluids is enabled by the shear shining features of CNT nanofluids. We contend that because the increment is independent of any discernible chemical change in the research area, qualitative advancement can happen at the heat transfer. The distribution of Nusselt numbers at cavities, on the other hand, is said to follow the geometric pattern of the wall by Arefmanesh et al. [48], which is one of the key purposes of wavy walls at cavities. Engineers can create thermal access points depending on flexible usage at various locations using these models. Additionally, it reduced the development of entropy (Alsabery et al. [49]). Alongside, Das and Mahmud [50] discovered that the heat exchange rate is significantly influenced by the amplitude of the wavy wall and the amount of undulations, particularly in the convective mode (Kumar [51]). A further important aspect is that when the hollow wall is wave shaped, the fluid domain is vulnerable to some disruption, as is the case in storm courses, earthquake systems, and maritime transportation.
Nevertheless, because various mediums have varying temperatures, heat transportation to trace energy dynamics follows a path along a temperature gradient. The viewpoint of thermal equilibrium is used to analyze the above phenomenon. In contrast, the term "local thermal non-equilibrium" (LTNE) is used to describe situations in which there is a "significant seepage velocity, a significant temperature difference between the fluid and solid matrix of a porous medium, or when a heat source exists in either a solid domain or a pure fluid domain within the porous medium" as was stated by Alsabery et al. [52]. Accordingly, when the heat transmission coefficient is low, the temperature of the solid phase is lower than the temperature of the fluid phase as a result of Baytaş's [53] application of a heat source to the solid phase. Aside from that, Sheremet et al. [54] show that a lower interface heat transfer coefficient implies to decrease in the heat exchange rate, which inhibits the LTNE. In addition, Zargartalebi et al. [55] found that the buoyancy ratio has a negative impact on the heat exchange ratio while positively affecting the Rayleigh number. This finding indicates that the Rayleigh number enhances the thermal equilibrium because the thermal transport via diffusion is greater than that via convection. The system's increasing entropy state, however, causes the buoyancy ratio to induce a non-equilibrium state to rise. In light of the influence of Buongiorno's concept, Tahmasebi et al. [56] discovered that speeding up the exchange of heat at the interface entails slowing down the interchange of the micro liquid phase. Regarding Brownian motion, this model. On the basis of our knowledge, Brownian motion prevents heat exchange at the contact. However, Nield et al. [57] found that the local Nusselt number is poorly dependent on the Darcy number or the fluid-solid heat exchange parameter and relies entirely on the Peclet number and the solid/fluid conductivity. This could be due the advection/diffusion regimes are so effective since the velocity profile overtakes heat transfer in boundary layer flows. Nonetheless, we may conclude that there is a connection between the Darcy number and the heat exchange parameter in that scenario to understand why the heat exchange is ineffective. The reliability of the local thermal equilibrium for all porous media when the conduction mode is predominate has been shown by Kim and Jang [58].
Furthermore, when the interstitial heat transfer coefficient is proportional to Re exponent>1 (such as sintered metals and cellular ceramics) and the convection mode is dominant, the local thermal equilibrium is applicable. For the purpose of choosing an appropriate analytical interpretation, Amiri et al. [59] made the supposition that the constant heat flux at boundaries is satisfied when either heat propagation is based on temperature gradients (effective conductivities) or each phase receives the same amount of the flux at the interface. In addition, Mahmoudi and Karimi [60] have shown that the Nusselt number for the first method is higher than the other one for the second. The Darcy number was negatively associated with the thickness of the local thermal non-equilibrium porous zone, which they were able to measure. Aside from that, Izadi et al. [61] have observed that the thermal nonequilibrium case deteriorates as the thermal conductivity ratio and heat transfer parameter increase. Moreover, the reduction in the thermal resistance [76] of fluid phase can enhance the heat transmission rate due to the increase in the thermal conductivity ratio. In addition, permeability influences the Nu number (the dimensionless quantity) of both solid/fluid phases due to nanoparticles rotating around the center of gravity. Sivasankaran et al. [62] have found that exchange speed of thermal energy is enhanced by increasing values of the modified conductivity ratio and the porosity of the media. Entirely, Tahmasebi et al. [63] have found that the average Nu is supported by Rayleigh number, Darcy number, and thermal conductivity ratios at solid and fluid phases. Furthermore, the convection interaction regime at pores increases the Nu for the solid phase whereas decreases it for the nano-fluid.
Regarding the running study, the control parameters of the model encompasses various heat and mass dynamics have been analyzed numerically in line with the MHD mixed convection around a central circular heat source and the other cold wavy boundaries. The manuscript is organized as following; first section is the introduction where we shed light on the different styles of heat transfer mechanisms due to fluid flow within miscellaneous media under different conditions. In the second section, the mathematical formulation for the problem is designated. Then the fundamentals of the numerical investigation are demonstrated in the third section. Therefore, at section four, the results are displayed in the form of the contours of the temperature and fluid flow which disclose the active regions within the cavity while the Nusselt number profiles monitor the heat flux behavior at the boundaries. Ultimately, the fifth section sum up the findings of the proposed work.  Figure 1 shows the preliminary geometry of an undulating inclined cavity. The involved wavy cavity is filled by porous media and hybrid nanofluids. The flow of porous medium is subjected to the extension of Darcy's law using Dupuit-Forchheimer relationship. The right and left side wavy-walls are cold ( T c ), the bottom and top sides are adiabatic. Moreover, the fluid domain is inclined with angle α and the magnetic field has an inclination angle Φ . The hybrid nanofluid convection is not in a local thermodynamic equilibrium condition. The normal direction and constant value are considered for the gravity acceleration. Dirichlet type applied on all boundaries (no-slip condition). According to the previously specified hypotheses, continuity, momentum, concentration, and energy equations are formulated for hybrid nanofluids, incompressible, laminar, and single-phase flows in steady-state as follows [64,65]: The boundary conditions imposed on the flow field are taken as: On the left wall where u and v are the velocity components, T is temperature, s-index refers to solid phase, f-index refers to fluid phase, C stands for the concentration of nanoparticles, hnf is the density of the hybrid nanofluid, hnf is kinematic viscosity. g is a gravity, p is a pressure, hnf is a dynamic viscosity, β* = (− 1/ρ)(dρ/dC) T,p , the buoyancy rate due to nanoparticle concentration currents, and Q 0 is the heat generation.

On the top wall
On the bottom wall: On the inner circular: The local Nusselt number of the fluid can be formulated as: The local Nusselt number of the solid can be written as: The average Nusselt number of the fluid and solid are

Thermophysical properties of nanofluid and hybrid nanofluid
While some research has looked at determining the thermophysical properties of nanoparticles (Table 1), classical models do not necessarily apply to nanofluids. Experimental results facilitate the selection of an appropriate model for a particular property. Defining the effective properties of TiO 2 /water nanofluids and TiO 2 -Cu/water hybrid nanofluids (Table 1), respectively, are as follows: which determine density of nanofluid, the density of hybrid nanofluid is specified as [66]: where ϕ is the overall volume concentration of two different types of nanoparticles dispersed in hybrid nanofluid and is calculated as; = TiO 2 + Cu , and the heat capacitance of the nanofluid given is as, where bf and p are the coefficients of thermal expansion of the fluid and of thesolid fractions respectively.
Hence, for hybrid nanofluid, thermal expansion can be defined as follows: Thermal diffusivity, nf of the nanofluid isdefined as: k nf is the thermal conductivity of the nanofluid which for spherical nanoparticles, according to the Maxwell-Garnetts model [67], is: Thus, thermal diffusivity, hnf of the hybrid nanofluid can be defined as; If the thermal conductivity of hybrid nanofluid is defined according to Maxwell model, The effective dynamic viscosity of the nanofluid based on the Brinkman model [68] is given by where μ bf is the viscosity of the fluid fraction, then the effective dynamic viscosity of the hybrid nanofluid is; According to (6), heat capacity of hybrid nanofluid can be determined as follows: The thermal expansion coefficient of the nanofluid can be determined by:  geometrical complexity of wavy boundaries, stagnation points, and critical points which often cause disturbances in the approximation solution. This issue could be more realistic from the physical domain view but it sacrifices the stability of the solution environment. Therefore, such an approved computational domain is used with the finite volume method presented in [72,73] which is extended here to the non-orthogonal grids. This technique starts with definitions of the new coordinates (ξ, η) as following: Using Eq. (38), the irregular physical model is mapping to a rectangular computational model as it is shown in Fig. 2. The resulting transformed system is treated using the finite volume method. Here, the advection terms are evaluated using "the second upwind scheme", whilst the Laplace operators are handled using the central difference methods. Iteratively utilizing ADI (alternating direction implicit), the resultant algebraic structure is solved. This study's convergence criteria were selected to be systematic. Additionally, the size of the grid is discovered to be appropriate for all computations after doing several grid tests. Additionally, Fig. 3 compares the outcomes of the current simulation to those attained by Cheong et al. [74] in certain instances (ϕ = φ = Ha = Q = 0). It is discovered that (38) the outcomes are in good accordance with one another. Furthermore, another comparison test is performed with the results of Khanafer and Chamkha [75] where the mixed convection due to the moving of the upper wall of a square enclosure has been examined. Table 2 shows an excellent agreement between the results are obtained.

Results and discussion
In this section, the obtained results are presented stepby-step in terms of key parameters such as the heat generation parameter (Q) is varied from 0.1 to 1, the effects of Hartmann number (Ha) is varied from 0 to 150, the It is important to glance on the prominent terms at the governing equations. For example for momentum conservation equations, the surface forces are represented by the pressure components at the porous medium, while the body forces come from gravitational term which the linear combination of both heat and nanoparticles concentration and the magnetic field. The term H * (θ s − θ f ) in energy equation points to the internal heat exchange between the solid and fluid phases which characterizes the non-thermal equilibrium at porous medium.
For streamlines as in Fig. 4, there are two vortices of the utmost stream function that are concentrated adjacent to the mobile surfaces (top and bottom). The unnormal issue is the stagnant fluid around the heat source which means that the forced convection dominates more regions than the natural convection. It could be back to the value of Richardson number where "Ri (→ 1)" in this study which is relatively small. The effect of heat generation parameter "Q" is weak due to the non-thermal equilibrium. For the fluid isotherms, two vortices are established near to those of streamlines. Still the amazing event is the cold fluid around the heat source region meaning that the supremum fluid temperature is created basically from forced convection. The solid isotherms are different than fluid isotherms at the location. The isotherms are two adjoining vortices and their centers (maximum solid temperature values) are away from those of streamlines because the solid isotherms approaches to the heat source. The two vortices are united at the left half may be due to the numerical calculations or the effect of the cavity inclination angle. The two solid vortices appeal to make a circle around the heat source. For the interpretation; the nonthermal equilibrium enhances the inconsistency between the fluid and solid isotherms distributions. We can see that fluid temperature is greatly affected by the movable sides while the solid temperature is influenced in the heat source position. the mutual effect between solid and fluid phase could not be neglected but it doesn't represents the basic factor of heat transmission here. Moreover, the isoconcentrations contours are concentrated at the bottomleft corner due to the effect of inclination angle "α".
The fluid Nusselt number along the bottom side (Fig. 5) scores nonzero values at the origin point of the cavity which is the conjunction point of movable bottom side and wavy left side. Besides, Nusselt number completely vanishes when arrives the right side so the left side is so active region while the right corner is semi isolated region at the cavity may be due to the effect of the inclination angle. Normally, As "Q" increases the Nusselt number increases (Fig. 6). On the other hand, the solid Nusselt number along the bottom side is represented by a curve of the second order which is symmetric around the heat source location. Consequently, this represents an evidence of the robust dependency between the heat source position and the heat transfer at the solid phase which has been mentioned above at the contours discussion. In addition, the average Nusselt number (Fig. 7) for the fluid phase is larger (~ 100%) than that of the solid phase, which reflex the strong fluid response to the Q parameter more than the solid phase. The effect of Hartmann number appears as a feedback of the inclination angle of the magnetic field used in the model as it is shown at Fig. 8 for the streamlines. Consequently, the magnetic field enhances the contours deformation along the second diagonal of the square cavity. Accordingly, as "Ha" increases, the supreme isotherms extend through the left side of the cavity (because of the cavity inclination angle effect too) and join to each other. Therefore, the magnetic field supports the thermal equilibrium case where the fluid and solid isotherms are congruent. Iso-concentrations contours show a slightly reasonable response to the isotherm's dynamics may be because both of them is added to momentum conservation equations as a linear component. Fluid Nusselt number along the bottom boundary (Fig. 9) shows variations at the left side of the cavity as "Ha" changes. By the way, the difference between left and right halves is discussed above (Fig. 6). The noticeable point here is that Hartmann number suppresses the heat transfer (more than ~ 18%); meaning that the magnetic field has a negative effect on the forced convection which clearly appears at Fig. 9 for the average fluid Nusselt number. Therefore, the magnetic field works for the equilibrium between the forced and natural convection regimes.
Solid Nusselt number (Fig. 10) along the bottom shows a somewhat change near to the midpoint as the radiation parameter increases while it keeps as a latent at the limits. The fluid average Nusselt number (Fig. 11) discovers the inhibitory role of the radiation which could be due to the Titanium nanoparticles characters while for the solid mean   Figure 12 reflects again the opposite behavior of heat transfer at the solid and fluid phases. It's found that the solid average Nusselt number is enhanced by "H * " while it decays for the fluid phase. The reason could be referred to some factors such as the thermophysical properties of the used nanoparticles, the dominant convection regime, or the boundary condition effects.
Streamlines increase as the porosity increases (Fig. 13) which agrees with the natural behavior. On the other hand, Fig. 14 shows that the average Nusselt number decreases as the porosity increases, because of the forced convection decay. Exhaustively, the average fluid Nusselt falls down ~ 60% while the average solid Nusselt descends only ~ 6%. In the same context, Fig. 15 shows the Darcy number effect, where "Da" tracks the permeability at the unit area intersection. Decreasing "Da" means weak forced convection which clearly appears from streamlines. Consequently, the supreme isotherms vortices transfer from the top/bottom boundaries to the circular source at the center. The amazing phenomenon is that the concentration descends all over the boundaries while grows up around the circle heat source.
It's known that the augmentation for the Richardson number supports the natural convection [76], therefore the streamlines and isotherms are concentrated around the heat source (Fig. 16). As a consequence the fluid Nusselt number is depressed ~ 23% while the solid Nusselt number arises ~ 7% (Fig. 17). Th interpretation is for the large "Ri", the natural convection dominates the fluid flow regime, so the fluid temperature difference between the heat source and the other sides decreases in general, but for solid temperature difference is still great due to the non-thermal equilibrium condition. Ultimately, for the cavity geometry, Fig. 18 shows the effect of wave amplitude. As shown, for small A, the contours present uniform vortices as the fluid mass increases. As "A" increases, the  (Fig. 19) decreases ~ 15% because as "A" increases, the fluid flow is hampered. Therefore, the convection regime is deficient and fluid Nusselt number descends.

Conclusion
The mixed convection phenomenon of hybrid nanofluids within an undulating porous cavity has been investigated within this paper. The contours of the streamlines, isotherms of both fluid and solid phases and isoconcentrations have been inspected as well as the profiles of local and average Nusselt number of the fluid/solid phase under local thermal non-equilibrium case at various key parameters such as the coefficient of heat generation/absorption Q, Hartmann number Ha, porosity parameter ε, an inter-phase heat transfer coefficient H * , undulation parameter λ, Darcy parameter Da, magnetic field inclination angle Φ, and hybrid nanofluid parameter ϕ. The remarkable points at the thermodynamic model under the study could be concluded as: • The hybrid nanofluid around the circular heat source is somewhat clotted where there is low temperature and very slow motion, • This infers that the isolated wavy sides put down the buoyancy force which is the main pillar of natural convection, • The local thermal non-equilibrium case (LTNE) plays a prominent role as it could create some isolated regions through the entire cavity domain. Therefore, the fluid isotherms are pulled toward the movable surfaces while the solid isotherms concentrated close to the central heat source, • Consequently, the heat transfer is not accompanied with the mass transfer. Moreover, there is somewhat an isolation between the nanofluid mass and the solid matrix, • The driven lids mounting at top and bottom surfaces enforces the forced convection dominancy, • In this arrangement, the circular heat source serves as a thermal barrier and flow resistance. • The fluid response to heat generation "Q" is much stronger (~ 100%) than the solid response, • The magnetic field strengthens the thermal equilibrium state which satisfies the balance between the natural convection and the forced convection, • The thermal effect of the parameters {Ha, R d , H*, Ri} is reverse when the phase changes (solid/fluid). While the effect of {Q, ε} doesn't change. This because the last two factors don't change the convection regime. However, the former parameters are either directly related to the LTNE case or change the dominant regime obviously. • The average fluid Nusselt number descent due to the porosity is ~ 60% which represents ten times the of that at the solid descent. • Iso-concentrations contours show a slightly reasonable response to the isotherm's dynamics when the magnetic field varies. • The magnetic field has a negative effect on the forced convection. • The average Nusselt number decreases as the porosity increases.
Acknowledgements Authors appreciate the great efforts provided by Prof. Sameh E. Ahmed. His help has greatly improved the article status.

Author contributions
Prof. Mansour has established the theoretical basis of the model and has conducted the numerical investigation, Bakier has introduced for the study, discussed the results, and organized the sections through the paper.
Funding Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB). During the preparation of this manuscript, no funds, grants, or other support was received by the authors.

Data availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Conflict of interest
The authors have no relevant financial or nonfinancial interests to disclose.
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.