MHD Mixed Convection in a Lid-Driven Cavity with a Bottom Trapezoidal Body: Two-Phase Nanofluid Model

The current work examines the effects of a bottom trapezoidal solid body and a magnetic field on mixed convection in a lid-driven square cavity. The Al2O3-water nanofluid used is assumed to obey Buongiorno’s two-phase model. An isothermal heater is placed on the bottom base of the trapezoid solid body, while the cavity’s vertical walls are kept cold at temperature Tc. The top moving wall and the remaining portions of the cavity’s bottom wall are thermally insulated. The Galerkin weighted residual finite element method is employed to solve the dimensionless governing equations. The parameters of interest are the Richardson number (0.01 ≤ Ri ≤ 100), Hartmann number (0 ≤ Ha ≤ 50), nanoparticle volume fraction (0 ≤ φ ≤ 0.04), and the length of the bottom base of the trapezoidal solid body. The obtained results show that increasing the Richardson number or decreasing the Hartmann number tends to increase the heat transfer rate. In addition, both the thermophoresis and Brownian motion greatly improve the convection heat transfer. It is believed that the current work is a good contribution to many engineering applications such as building design, thermal management of solar energy systems, electronics and heat exchange.


Introduction
Mixed convection mechanism is important for controlling heat transfer in many applications, such as electronic systems and nuclear reactors.Mixed convection in cavities can be easily realised by moving surfaces [1].Lid-driven cavities appear in many engineering applications, such as oil extraction, design of heat exchangers, crystal growth and float glass production, and cooling of electronic equipment [2].Reference Torrance et al. [3] considered the effects of an upper moving wall on convection in a rectangular cavity.Moallemi and Jang [4] considered the effects of the Prandtl number on mixed convection in a lid-driven square cavity.They found that the temperature distribution and the flow field are clearly influenced by the Prandtl number.The steady mixed convection heat transfer problem in a double lid-driven square cavity was investigated numerically using the finite control volume by Oztop and Dagtekin [5].Reference Al-Amiri et al. [6] employed the Galerkin weighted residual method to study steady mixed convection in a lid-driven square cavity.Reference Sharif [7] studied laminar mixed convection in an inclined lid-driven cavity.The steady mixed convection in a water-filled square lid-driven cavity in the presence of partial slip condition was investigated by Ismael et al. [8] who showed that the heat transfer rate increases with the Richardson number.Reference Ismael [9] considered mixed convection in a cavity with an arc-shaped moving wall.
The use of nanofluids could help to enhance heat transfer and the process of heat removal.Nanofluids are defined as the addition of nanoparticles (<100 nm-diameter), like metal or ceramics, to base fluids that have a very low thermal conductivity, like water and oil.The problem of mixed convection in a double lid-driven square cavity filled with Cu-water nanofluid was reported by Tiwari and Das [10].The flow and heat transfer were shown to be affected by the Richardson number, the solid volume fraction, and the moving walls' directions.There are two models for nanofluids: the single-phase (homogenous) and two-phase models.In the single-phase model, both the nanoparticles and fluid phase are assumed to be in thermal equilibrium, i.e., the fluid-nanoparticles' slip velocity is negligible [1].Many studies have considered the effects of the single-phase nanofluid model on mixed convective heat transfer inside a lid-driven cavity.Reference Talebi et al. [11] numerically investigated the mixed convection of nanofluids in a lid-driven square cavity.They observed that increasing the Reynolds number tended to reduce the effects of the solid volume fraction.Reference Chamkha and Abu-Nada [12] studied the effects of two different viscosity models on the steady mixed convection of Al 2 O 3 -water nanofluid in a double lid-driven cavity.Reference Karimipour et al. [13] considered the mixed convectio of Cu-water nanofluid in an inclined lid-driven cavity using the lattice Boltzmann method.They found that increasing the nanoparticle volume fraction had the effect of increasing the average Nusselt number.The unsteady natural convection in a nanofluid-filled trapezoidal cavity with a non-uniform side-wall temperature was studied by Alsabery et al. [14].
The slip velocity between the base fluid and nanoparticles was shown experimentally to be non-negligible by Wen and Ding [15].Thermophoresis and Brownian diffusion are the two main nanofluid slip mechanisms of the two-phase nanofluid model of Buongiorno [16].A number of investigators have employed the two-phase nanofluid model to model and discuss nanofluid heat transfer.Considering the two-phase model, Alinia et al. [17] numerically studied the mixed convection of a nanofluid in an inclined double lid-driven square cavity.Adding the SiO 2 nanoparticles was shown to greatly enhance the convection heat transfer.Using the control volume based finite element method, Sheikholeslami et al. [18] studied convection heat transfer inside a semi-annulus enclosure based on the two-phase nanofluid model.They observed that the Nusselt number decreased with the Lewis number.Reference Sheremet and Pop [19] reported the problem of the steady laminar mixed convection and heat transfer of a nanofluid in a lid-driven square cavity.The mixed convection of alumina-water nanofluid in microchannels was investigated by Malvandi and Ganji [20] using a modified two-phase model of Buongiorno.Garoosi and Talebi [21] studied the effects of the two-phase model on the conjugate natural and mixed convection of a nanofluid in a square cavity.Recently, Alsabery et al. [22] examined the effects of a solid inner insert on the mixed convection of a two-phase nanofluid model in a double lid-driven square cavity.They observed various effects on the rate of heat transfer with the augmentation of nanoparticles at a low Reynolds number and high Richardson number.
Recently, more emphasis has been placed on studying the effect of a magnetic field on the convection heat transfer in cavities [23,24].Reference Chamkha [25] investigated the effects of magnetic field and internal heat generation on mixed convection in a lid-driven cavity.The results showed that the magnetic field reduced the average Nusselt number.MHD mixed convection and heat transfer in a lid-driven square cavity with a corner heater was reported by Oztop et al. [26].They found that increasing the Hartmann number reduced the heat transfer.Reference Chatterjee [27] investigated the effect of a heat source on MHD mixed convection in a lid-driven square cavity.Reference Selimefendigil et al. [28] studied the effect of an inclined magnetic field on MHD mixed convection and entropy generation in a nanofluid-filled lid-driven cavity.MHD mixed convection of a nanofluid in a partially-heated wavy-walled lid-driven cavity was considered by Öztop et al. [29].They found that the Hartmann and Richardson numbers affected the heat transfer rate with the augmentation of the nanoparticle volume fraction.Reference Chamkha et al. [30] considered the effects of a heat source/sink and a partial slip condition on the MHD mixed convection and entropy generation of Cu-water nanofluid in a lid-driven porous cavity.Recently, the MHD mixed convection of a two-phase nanofluid model in a lid-driven square cavity with an inner solid block and corner heater was studied by Alsabery et al. [31].They observed that increasing the concentration of the nanoparticles by more than 2% had negative effects on the rate of heat transfer.
To the best of the authors' knowledge, no study has been conducted on the MHD mixed convection of a nanofluid in a lid-driven cavity with a bottom trapezoidal body.We employ the two-phase nanofluid model of Buongiorno and the cavity is heated from below at the bottom base of the trapezoidal body.It is believed that the current work is a good contribution to many engineering applications such as building design, thermal management of solar energy systems, electronics, and heat exchange.

Mathematical Formulation
Consider a lid-driven square cavity of length L filled with a water-based nanofluid with Al 2 O 3 nanoparticles.The cavity contains a bottom trapezoidal body with a bottom base of length h, a top base of length h/2, and height d, as shown in Figure 1.The bottom base of the trapezoidal body is heated at a higher temperature, T h .Both the left and right walls are kept at a cold temperature, T c , while the top moving wall and the remaining parts of the bottom wall are thermally insulated.The acceleration due to gravity is assumed to be vertically downward.The applied magnetic field is uniform and inclined at angle γ.All four boundaries are taken to be impermeable.The fluid physical properties are considered to be constant except for the density (Boussinesq approximation).The flow is steady and laminar.The governing equations are as follows: where v is the velocity vector, p is the pressure, T is the temperature, B is the magnetic field, g is the gravitational acceleration vector, and ϕ and J p are, respectively, the local volume fraction of nanoparticles and the nanoparticles' mass flux.The trapezoidal body's energy equation is The nanoparticles' mass flux J p can be written as (Buongiorno's model): where The nanofluid's heat capacitance, effective thermal diffusivity, effective density, and thermal expansion coefficient are given as, respectively, (ρβ The electrical conductivity, dynamic viscosity and thermal conductivity ratios of Al 2 O 3 -water nanofluid (33 nm particle-size in the ambient condition) are, respectively (cf.[32,33]), where ) is the Boltzmann constant, and l f = 0.17 nm is the mean path of the fluid particles.The water molecular diameter [33] is where M, N and ρ f are, respectively, the molecular weight of water, the Avogadro number, and the density of the base fluid at standard temperature (310 K).
The governing equations can be made dimensionless if we take the following non-dimensional variables Hence, the dimensionless governing equations are subject to the dimensionless boundary conditions.On the bottom wall: On the cold left and right vertical walls: On the horizontal adiabatic moving top wall: On the trapezoidal solid body's walls: Here, V = (U 0 , V 0 ) and B* = (Ha 2 sin γ, Ha 2 cos γ).The reference Brownian diffusion and thermophoretic diffusion coefficients are In (27), the thermal conductivity ratio K r is k w /k n f .Both the local and average Nusselt numbers at the heated part of the bottom wall can be evaluated, respectively, using Nu n f dX.(28)

Numerical Method and Validation
We used the Galerkin weighted residual finite element method to solve problems ( 17)-( 27).The reader is referred to Alsabery et al. [34] for a brief description of the method.For convergence of the solution, we used the convergence criterion (Γ i+1 − Γ i )/Γ i+1 ≤ 10 −6 , where Γ represents the variable and i is the iteration number.The ranges of values of the Reynolds number and Richardson number considered in the current work held the nanofluid flow to be incompressible and laminar.
First we validated the present code against the experimental and numerical work of Calcagni et al. [35].Figure 2 presents the results of streamlines and isotherms for the case where there was no solid body in the cavity.Very good agreement was demonstrated between our results and the experimental and numerical results of Calcagni et al. [35].Next, we considered the case of natural convection in a Buongiorno's nanofluid-filled cavity.Figure 3 shows a comparison of our calculated average Nusselt number to that of the experimental work of Ho et al. [36] and the numerical results of Sheikhzadeh et al. [37] and Motlagh and Soltanipour [38].We can see that a reasonably good agreement was reached.Figure 4 further shows a good agreement between the present method and the numerical results of Corcione et al. [39] for the streamlines, isotherms, and nanoparticle distribution of natural convection.Finally, we can see from Figure 5 that our calculated thermal conductivity and dynamic viscosity compare reasonably well with the results of Chon et al. [40] and Corcione et al. [39] and Ho et al. [36] and Corcione et al. [39], respectively.

Results and Discussion
We investigated the effects of the Richardson number (Ri), Hartmann number (Ha), nanoparticle volume fraction (φ), and the height of the trapezoidal body (D) on the streamlines, isotherms, nanoparticle distribution, and local and average Nusselt numbers.The other parameters were fixed as follows: Re = 100, Pr = 4.623, H = 0.5, Le = 3.5 × 10 5 , Sc = 3.55 × 10 4 , γ = π/4, N BT = 1.1, and δ = 155.The trapezoidal body was assumed to be made of fiberglass with k w = 0.045 W/(m K).Table 1 lists the thermophysical properties of the Al 2 O 3 -water nanofluid.Figure 6 describes the fluid flow (streamlines), temperature fields (isothermal lines), and nanoparticle distribution (volume fraction distribution) for several Ri, while the other parameters were fixed at Ha = 20, φ = 0.02 and D = 0.3.The flow within the square cavity takes a singular streamline cell that rotates clockwise next to the top moving wall.Most of the isotherms follow vertical lines close to the right cold wall; such a behaviour is clearly related to the dominant forced convection.This phenomenon can be obviously observed from the contours of the nanoparticle distribution; most of the nanoparticle concentration appears at the lower segment of the square cavity affected by the movement of the driven top wall.Increasing the Richardson number clearly affects the flow behaviour.A secondary streamline cell in the anticlockwise direction takes place at the lower part of the cavity near to the right wall.The strength of the streamlines is enhanced by the augmentation of Ri.This indicates that natural convection is dominant.The isotherm patterns transfer from the vertical shape into a plume-like pattern with high destruction within the middle of the cavity.The concentration of nanoparticles moves around the walls of the cavity which clearly shows the enhancement of the thermophoresis effects, as depicted in Figure 6c,d.The effect of the nanoparticle loading on the interface of the local Nusselt numbers with X for Ri = 10, Ha = 20, and D = 0.3 is described in Figure 11a.The local Nusselt number is an increasing function of the solid volume fraction due to the increment of the thermal gradient due to an addition of the nanoparticles.Such a behaviour is related to the enhancement of the Brownian diffusion and the thermophoresis effects.Figure 11b depicts the influence of the local heat transfer for different lengths of the trapezoidal body at Ri = 10, Ha = 20, and φ = 0.02.At the left segment of the cavity, the local Nusselt number shows a clear reduction with an increment in D, while this increment in D tends to enhance the local Nusselt number at the right segment of the cavity, which results in an increase of the thermal gradient.Figure 12 shows the variation in the average Nusselt number with various solid volume fractions (φ) and Richardson numbers (Ri) for Ha = 20 and D = 0.3.At a low value of the Richardson number (Ri = 0.01), when the forced convection is dominant, the convection heat transfer shows very clear enhancement with the addition of the solid volume fraction.This happens due to the increment of the thermal conductivity and thermophoresis.The augmentation of Ri leads to the enhancement of the thermophoresis effects, and as a result, the convection heat transfer increases.More importantly, when Ri is relatively high (Ri = 10 and 100), the addition of the solid volume fraction that exceeds 0.035 leads to a reduction in the average Nusselt number.This happens due to the fact that the buoyancy effects dominate the thermal conductivity and thermophoresis effects at this stage.The variation in the average Nusselt number with the Hartmann number (Ha) and the solid volume fraction (φ) for Ri = 10 and D = 0.3 are shown in Figure 13.Clearly, we can see the enhancement of the convection heat transfer with an increment in the nanoparticle loading.This is due to the increase in the thermal conductivity which leads to a rise in the thermal gradient.However, the increasing of the solid volume fraction that exceed 0.03 tends to decrease the average Nusselt number.Imposing an external magnetic field Ha > 0 leads to an obvious reduction in the rate of the heat transfer, which is due to the thermal gradient diminishing.
Figure 14a describes the influence of the solid volume fraction (φ on the average Nusselt number for various values of the Richardson number and Ha = 20 and D = 0.3.The convection heat transfer enhances with an augmentation of Ri, which is due to the increase in the buoyancy and the shear forces.The nanoparticle loading shows an improvement in the average Nusselt number due to the thermal conductivity and thermophoresis increment.However, the convection heat transfer shows a reduction behaviour when the augmentation of the solid volume fraction exceeds 0.03.The effects of φ on the average Nusselt number for various values of the Hartmann number (Ha) are described in Figure 14b for Ri = 10 and D = 0.3.A decreasing behaviour is observed for the average Nusselt number with an increment in the Hartmann number due the the reduction in the thermal gradient with as Ha increases.The effects of the length of the bottom base of the trapezoidal solid body (D) on the average Nusselt number for various values of the Hartmann number are shown in Figure 15a.The convection heat transfer reduces with an increase in the Hartmann number.Inserting a slight solid body within the square cavity leads to an enhancement in the rate of heat transfer.This behaviour is seen for a low magnetic field (Ha < 10).However, a different behaviour is obtained for the convection heat transfer when the magnetic field exceeds 25.The average Nusselt number enhances with as D increases.Figure 15b shows that the average Nusselt number appears to have an enhancement behaviour with as D increases for relatively low φ ≤ 0.025.

Conclusions
The current paper investigated the effects of a magnetic field on mixed convective heat transfer in an Al 2 O 3 -water nanofluid filled lid-driven square cavity with a bottom trapezoidal solid body.Buongiorno's two-phase model was adopted, i.e., the slip between the base fluid and the nanoparticles was taken into consideration.The governing equations were made dimensionless and then solved using the Galerkin weighted residual finite element method.The concluding remarks are as follows: 1.At a high Richardson number when the natural convection is dominant, thermophoresis can produce a homogeneous dispersion of nanoparticles.2. The convection heat transfer enhances with the nanoparticle volume fraction and is affected by increments in the thermal conductivity and thermophoresis.In contrast, at a high Richardson number, a boost in the nanoparticle volume fraction that exceeds 0.035 leads to a reduction in the average Nusselt number.3. A cavity with small trapezoidal solid body and low magnetic field tends to show more enhancement on the rate of heat transfer.However, imposing a high magnetic field in a cavity with a big trapezoidal solid body leads to the best improvement in convection heat transfer.4. A strong enhancement of the average Nusselt number occurs with a high nanoparticle volume fraction and a relatively small trapezoidal solid body. 5. Thermophoresis and Brownian motion effects play significant roles in the augmentation of convection heat transfer.

Figure 1 .
Figure 1.Sketch of the problem geometry.

Figure 7
Figure 7 examines the effect of the Hartmann number (Ha) on the streamlines, isotherms, and nanoparticle distribution for fixed values of Ri = 10, φ = 0.02, and D = 0.3.The flow within the square cavity appears with two streamline circular cells in different directions, one at the upper part and the other at the lower part next to the right wall for the case where Ha = 0.The isotherm patterns are destroyed at the left segment of the considered cavity.The nanoparticle distributions seem to be clearly affected by the dragging force that is supported by the convection heat transfer; the distribution of nanoparticles takes place around the trapezoidal body next to the walls of the cavity.The magnetic field reduces the buoyancy and inertial forces; as a result, the strength of the flow circulation decreases.This application of Ha leads to an inhibition of the thermal gradient, the isotherms pattern show less distortion, and more lines transfer to the right segment of the cavity.The increase in Ha shows an enhancement in the distribution of the nanoparticles which leads to a significant nanoparticle migration.Increasing the Hartmann number up to the higher value (Ha = 50) clearly affects the flow circulation; a third streamline cell appears at the lower right part of the cavity, as observed in Figure 7d.

Figure 7 .
Figure 7. Effects of the Hartmann number (Ha) on the streamlines (left), isotherms (middle), and nanoparticle distribution (right) for Ri = 10, φ = 0.02, and D = 0.3.The impacts of the solid volume fraction (φ) on the streamlines, isotherms, and nanoparticle distribution are described in Figure8for Ri = 10, Ha = 15 and D = 0.3.The flow circulation shows less enhancement as the solid volume fraction increases; however, the strength of the flow circulation increases, which is clearly related to the nanoparticle Brownian diffusion.The isotherm lines show more deterioration and more lines move to the left wall, which is evidence of the temperature gradient increment.The temperature gradient leads to an increment in the thermophoresis; the nanoparticles show high distribution next to the right cold wall.The Brownian movements increase with such a loading of nanoparticles; nanoparticle migration mostly occurs close to the trapezoidal solid body.

Figure 8 .
Figure 8. Effects of the solid volume fraction (φ) on the streamlines (left), isotherms (middle), and nanoparticle distribution (right) for Ri = 10, Ha = 15, and D = 0.3.The effects of the length of the bottom base of the trapezoidal body (D) on the streamlines, isotherms, and nanoparticle distribution are depicted in Figure9for Ri = 10, Ha = 20 and φ = 0.02.The lower length of the trapezoidal body (D = 0.15) shows a big area for the nanofluid circulation; as a result, a high buoyancy and inertial forces occur.The isotherm patterns appear mostly distorted close to the left vertical wall due to the high thermal gradient.Strong distribution occurs on the contours of the nanoparticles; most of the nanoparticles are distributed at the lower segment of the cavity.Increasing D tends to reduce the flow circulation area and the buoyancy effect which clearly tends to decrease the strength of the flow circulation, as demonstrated in Figure9c,d.The deterioration of the isotherms reduces and more heat moves to the right segment of the cavity.This is related to the reduction of the thermal gradient affected by the low thermal conductivity of the solid trapezoidal body.Imposing a big solid trapezoidal body tends to reduce the nanoparticle distribution which clearly leads to more nanoparticles migration within the square cavity.

Figure 10
Figure 10 shows the impact of the Richardson and Hartmann numbers on the interface of the local Nusselt numbers with X at the heated part of the bottom wall for φ = 0.02 and D = 0.3.The local heat transfer enhances with an augmentation of the Richardson number due to the high buoyancy effect which results in an increase in the thermal gradient, as shown in Figure 10a.Figure 10b presents the effect of the Hartmann number on the local heat transfer.As observed from this figure, an increase in Ha tends to reduce the local Nusselt number due to the reduced thermal gradient and Brownian diffusion.
Figure 10 shows the impact of the Richardson and Hartmann numbers on the interface of the local Nusselt numbers with X at the heated part of the bottom wall for φ = 0.02 and D = 0.3.The local heat transfer enhances with an augmentation of the Richardson number due to the high buoyancy effect which results in an increase in the thermal gradient, as shown in Figure 10a.Figure 10b presents the effect of the Hartmann number on the local heat transfer.As observed from this figure, an increase in Ha tends to reduce the local Nusselt number due to the reduced thermal gradient and Brownian diffusion.

Figure 10 .
Figure 10.Plots of the local Nusselt number vs. X for different (a) Ri and (b) Ha at φ = 0.02 and D = 0.3.

Figure 11 .
Figure 11.Plots of the local Nusselt number vs. X for different (a) φ and (b) D at Ri = 10 and Ha = 20.

Figure 12 .
Figure 12.Plots of the average Nusselt number vs. φ for different Ri at Ha = 20 and D = 0.3.

Figure 13 .Figure 14 .
Figure 13.Plots of the average Nusselt number vs. φ for different Ha at Ri = 10 and D = 0.3.

Figure 15 .
Figure 15.Variation in the average Nusselt number with (a) Ha and (b) φ for different D at Ri = 10.