# Analysis of Heat Transfer in a Triangular Enclosure Filled with a Porous Medium Saturated with Magnetized Nanofluid Charged by an Exothermic Chemical Reaction

Authors: Raees-ul-Haq Muhammad
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

This paper intends to numerically study the steady-state free convection heat transfer in the presence of an exothermal chemical reaction governed by Arrhenius kinetics within a right-angled enclosure of triangular shape filled by porous media saturated with magnetized nanofluid. An approximation named as Darcy–Boussinesq approximation along with a nanofluid model mathematically propounded by Buongiorno has been implemented to model physical phenomenon representing fluid flow, heat transfer, and nanoparticle concentration. The mathematical equations in a dimensionless form describing the stream function for circulation of the fluid, the energy equation for heat, and nanoparticle volume fraction for concentration are solved using the finite difference method. The validity of the numerical procedure is established by comparing present results with the formerly available works in both statistical and graphical approaches. Streamlines, isotherms, and isoconcentrations are plotted and discussed for the various parametric regimes. The graphical description depicts that the average Nusselt and Sherwood numbers are the decreasing function of the Rayleigh number. The study revealed the accountable influence of model parameters such as thermophoresis and Brownian diffusion on the local Sherwood number, whereas a minimum impact on the local Nusselt number is observed.


---

## Body

## 1. Introduction

The flow of fluid and transfer of heat in a cavity filled by a nanofluid saturated with porous media are of vital importance in numerous practical situations. Some of these applications include geophysical applications, chemical reactors, building heating and cooling operations, electronic cooling equipment, and so on. Convective flows are studied extensively in enclosures filled with porous media as accounted in books by Pop and Ingham [1] and Bear [2]. Similarly, the enclosure of various geometric shapes is an essential aspect in engineering applications such as boilers or ovens with porous materials. Baytas and Pop [3] investigated the heat transfer characteristics for a steady-state flow in an inclined cavity filled with a porous medium. Sun and Pop [4] studied free convection in a triangular region filled with a porous medium. In their study, they placed a finite size flush mounted heater on the vertical wall of the cavity. They concluded that the average Nusselt number has a maximum value with a smaller aspect ratio of enclosure dimensions. It also has a lesser value by the lowest position of the heater with the highest magnitude of Rayleigh number and the maximum size of the heater. Izadi et al. [5] performed a numerical investigation of the natural convection heat transfer of a hybrid nanofluid into a porous cavity exposed to a variable magnetic field. Their findings indicate that at higher Rayleigh number, by increasing the Hartmann number, a significant decrease in Nusselt number occurs, which can be attributed to the decreased power of the flow. Arpino et al. [6] studied unsteady free convection in a partially porous annular enclosure and employed an adaptive local time-stepping procedure based on a properly stabilized full matrix-inversion scheme. They concluded that the position of the porous layer in the annulus strongly affected the stability of the free convection. Arpino et al. [7] have also studied porous medium-free fluid interface problems in the presence of high source terms by using a stable, accurate, and efficient artificial compressibility (AC) characteristic-based split (CBS) algorithm. In another paper, Arpindo et al. [8] presented an efficient three-dimensional algorithm based on the fully explicit matrix-inversion free finite element version of the characteristic-based split (CBS) scheme. The scheme was applied for the first time to the simulation of complex thermo-fluid-dynamic problems in domains containing a fluid and a porous layer. Massaroti et al. [9] adopted a dual time-stepping procedure based on the transient artificial compressibility version of the characteristic-based split algorithm to solve the transient equations of the generalized model for heat and fluid flow through porous media. Shermet and Pop [10] also studied steady-state free convection heat transfer in a right-angle triangular cavity filled with a porous medium using the mathematical nanofluid model proposed by Buongiorno [11].The subject of nanofluids has received considerable concentration from the researchers for its role in engineering and Biophysics fields. The presence of nanoscale particles in the fluid enhances heat transfer characteristics and has been used widely in the industry. For example, materials with sizes of nanometers possess unique physical and chemical properties (Oztop and Abu-Nada [12]). They tested Cu (copper), Al2O3 (alumina), and TiO2 (titania) as the base fluid to study the effect of nanoparticles on natural convection and thermal profiles. Choi [13] first introduced the term nanofluid which rooted from the addition of nanoscale particles to the fluids. The presence of nanoparticles in the fluid enhances the thermal conductivity within it. In a recent catalogue of articles by Kuznetsov and Nield [14] and Nield and Kuznetsov [15], a lot more have implemented Buongiorno’s mathematical model for nanofluids to study problems on Newtonian fluids and porous media filled by nanofluids. The authors reported that the use of surfactant or surface charge technology could cause suspension of nanoparticles in the nanofluid, hence preventing the particles from deposition on the porous matrix. Nanofluids have vital importance in the areas involving solar collectors to absorb solar thermal radiations. A wide set of nanofluid applications is recently reported in the books by Nield and Bejan [16] and YuS et al. [17]. Among several nanofluid mathematical models, Buongiorno’s model [11] received considerable attention since he noted that the effect of the velocity between nanoparticles and base fluids. Also, he has written down the conservation equations considering the two important effects which are the Brownian diffusion and the thermophoresis effects.The area of fluid mechanics that investigates the characteristics of fluids with electrical conductivity in the presence of an electromagnetic field is known as Magnetohydrodynamics (MHD). MHD convection flow has a wide range of applications in the area of science and engineering. These include geothermal energy extractions, petroleum reservoirs thermal insulation, plasma confinement, crystal growth, and nuclear reactor coolers. Chamkha [18] investigated the hydrodynamic convection flow in a vertical lid-driven cavity with absorption or internal heat generation. Malleswaram and Sivasankaran [19] presented a numerical analysis on magnetohydrodynamic-mixed convection in a lid-driven cavity when both the side-walls are provided with nonuniform heating. AlimM et al. [20] have presented a numerical study for the conjugate effect of Joule heating and MHD-mixed convection in a lid-driven square cavity. Rahman et al. [21] also studied the free convection within a tilted cavity filled with nanofluid in the presence of a porous medium and a sloping magnetic field energized by an exothermic chemical reaction. They concluded that the presence of the magnetic field causes the nanoparticles to separate from the fluid with small velocity. The process of free convection in the presence of an exothermic surface reaction has widely been studied by Merkin and Chaudhary [22], Mahmood and Merkin [23], Ikeda et al. [24], Gray and Merkin [25, 26], and Merkin and Pop [27]. Single first-order kinetic equation with Arrhenius [28] thermal energy dependence models the reaction on the surface and is defined as(1)A⟶B,rate=k0ae−E/RT,where A and B are the reactants, k0 denotes pre-exponential factor, a represents the reactant A′s concentration, and E and R represent the activation energy and fluid constant, respectively. Rahman [29] investigated a higher order chemical reaction along with internal heat generation on a non-Darcian forced convection for an incompressible viscous fluid flow with thermal conductivity across an elastic surface interlaced in a porous medium. He concluded that with the increase in Darcy number, the dynamic viscosity alongside the thermal conductivity reduces.It is apparent from the preceding literature study that most of the reviews are performed considering the nanofluids in cavities, but very less attention is paid to study the cavities filled with the porous medium saturated with nanofluids particularly in the presence of MHD and chemical reaction. So, in this paper, we study free convection within a triangular cavity filled with porous media saturated with nanofluid using Buongiorno’s model [11] incorporating Darcy’s law for the flow in the porous medium. We also employed the Boussinesq approximation to model the governing equations for the convective forces. The resulting partial differential equations of the advection-diffusion type are solved numerically using a finite difference method. The description of the scheme is presented in Appendix. To the best of authors’ knowledge, the problem of natural convection within a triangular cavity in the presence of a magnetic field driven by an exothermic chemical reaction using the mathematical nanofluid model has not been reported yet in the literature. As such, the focus of this paper is to study the detailed and attributed analysis on the flow and heat transfer due to the effects of pertinent parameters arising in the problem.


## 2. Problem Formulation

### 2.1. Physical Model

We consider an isosceles triangular-shaped cavity filled with a porous medium saturated with nanofluid consisting of water and nanoparticles. Thex¯-axis is along the bottom wall of width H, while the y¯-axis is along the vertical wall of the cavity with length L and maintained at the constant temperature To and the constant nanoparticle volume fraction Co. A similar assumption is also made for the inclined boundary, while the bottom wall along the x¯-axis is assumed to be adiabatic. A magnetic field Fb with strength B0 is imposed within the fluid along the positive x-axis. An exothermal surface reaction governed by Arrhenius kinetics [28] is also taken into account inside the cavity which is of order one and nonisothermal. Figure 1 represents the schematic diagram of problem geometry and the coordinate system. Water-based nanofluid inside the cavity is chosen for the analysis and is set to be laminar and incompressible. The base fluid (i.e., water) and the nanoparticles are assumed to be in thermal equilibrium. Hence, no slip occurs between them. The suspension of nanoparticles within the nanofluid is expected to be carried out using either surfactant or surface charge technology. Inside the triangular cavity, the fluid flow is considered to be weak, laminar, and incompressible with no heat transfer due to the viscous dissipation as well as the radiation. The nanofluid density depends on the volume fraction of nanoparticles and the temperature, while the other thermophysical properties of the nanofluid are expected to be constant. Compared to the initial temperature T0, the increase in temperature inside the flow domain is assumed to be smaller. As reported by CardosoT et al. [30], the rise in temperature before a chemical reaction is approximately below 10−50K, while general consideration of initial temperature T0 is below 632.2K. Hence, the Darcy–Boussinesq approximation is adapted for the present problem.Figure 1
Problem geometry.


### 2.2. Governing Equations

With the assumption that viscous dissipation effects are negligible, the Oberbeck–Boussinesq approximation is employed, and Darcy’s law is applicable. By following the nanofluid model proposed by Buongiorno [11], the governing equations for total mass, momentum, energy, and nanoparticles, respectively, are written as(2)∇⟶⋅V⟶=0,0=−∇⟶P∗−μKV⟶+F⟶b+ρgj^,ρcfV⟶⋅∇⟶T=km∇⟶2T+ερcpDB∇⟶T⋅∇⟶C+DTT0∇⟶T⋅∇⟶T+Qk0ae−E/RT,1εV⟶⋅∇⟶C=DB∇⟶2C+DTT0∇⟶2T,where Fb→ is given by Fb→=σu→×B→×B→, B→=B0,0,0 is the constant magnetic field strength along x¯−axis, and the buoyancy term ρg is approximated by the Boussinesq approximation:(3)ρg=−ρf0β1−C0∂T∂x¯g+ρp−ρf0∂C∂x¯g.Here,V⟶=u¯,v¯,0 is the velocity vector with u¯ and v¯ the nanofluid velocity components along x¯− and y¯−axis, respectively, T is the nanofluid temperature, C is the nanoparticle volume fraction, and other quantities of physical interest are mentioned in Nomenclature. To nondimensionalize the governing equations along with boundary conditions, we make use of the following similarity transformations:(4)x=x¯H,y=y¯H,u=Hαmu¯,v=Hαmv¯,θ=ET−T0RT02,ϕ=CC0,P=P∗Kμαm.Following the traditional way of defining the dimensionless stream functionψ from velocity components as follows:(5)u¯=∂ψ∂y¯,v¯=−∂ψ∂x¯.Using the dimensionless variables defined in equation (4) and the stream function in equation (5), the governing equations now take the form:(6)∇2ψ+Ha2∂2ψ∂x2+Ra∂θ∂x−RaNr∂ϕ∂x=0,(7)∇2θ−∂ψ∂y∂θ∂x−∂ψ∂x∂θ∂y+Nb∂ϕ∂x∂θ∂x+∂ϕ∂y∂θ∂y+Nt∂θ∂x2+∂θ∂y2+Fkeθ=0,(8)∇2ϕ−Le∂ψ∂y∂ϕ∂x−∂ψ∂x∂ϕ∂y+NtNb∇2θ=0.The boundary conditions are as follows:(9)ψ=0,θ=0,ϕ=1,atx=0,(10)ψ=0,θ=0,ϕ=1,atx+y=1,(11)ψ=0,∂θ∂y=0,∂ϕ∂y=0,aty=0,where Ha is the magnetic field parameter also known as the modified Hartmann number, Ra is the Rayleigh number, Nr is the buoyancy-ratio parameter, Nb is the Brownian motion parameter, Nt is the thermophoresis parameter, Fk is the Frank-Kamenetskii number, Le is the Lewis number, and the Laplacian operator ∇2 are all defined as(12)Ha=σB02Kμ,Ra=1−C0ρf0gβKRT02LμαmE,Nr=ρp−ρf0C0Eρf0β1−C0RT02,Nb=τC0DBαm,Nt=τDTRT0αmE,Fk=EQk0aexpE/RT0L2RT02αm,Le=αmεDB,∇2=∂2∂x2+∂2∂y2.For problems related to heat and mass transfer across the cavity, it is essential to determine rates of heat and mass transfer in terms of local and average Nusselt and Sherwood numbers, respectively, using the following formulas:(13)Nu=−∂θ∂xx=0,Sh=−∂ϕ∂xx=0,Nu¯=∫01Nudy,Sh¯=∫01Shdy.


## 2.1. Physical Model

We consider an isosceles triangular-shaped cavity filled with a porous medium saturated with nanofluid consisting of water and nanoparticles. Thex¯-axis is along the bottom wall of width H, while the y¯-axis is along the vertical wall of the cavity with length L and maintained at the constant temperature To and the constant nanoparticle volume fraction Co. A similar assumption is also made for the inclined boundary, while the bottom wall along the x¯-axis is assumed to be adiabatic. A magnetic field Fb with strength B0 is imposed within the fluid along the positive x-axis. An exothermal surface reaction governed by Arrhenius kinetics [28] is also taken into account inside the cavity which is of order one and nonisothermal. Figure 1 represents the schematic diagram of problem geometry and the coordinate system. Water-based nanofluid inside the cavity is chosen for the analysis and is set to be laminar and incompressible. The base fluid (i.e., water) and the nanoparticles are assumed to be in thermal equilibrium. Hence, no slip occurs between them. The suspension of nanoparticles within the nanofluid is expected to be carried out using either surfactant or surface charge technology. Inside the triangular cavity, the fluid flow is considered to be weak, laminar, and incompressible with no heat transfer due to the viscous dissipation as well as the radiation. The nanofluid density depends on the volume fraction of nanoparticles and the temperature, while the other thermophysical properties of the nanofluid are expected to be constant. Compared to the initial temperature T0, the increase in temperature inside the flow domain is assumed to be smaller. As reported by CardosoT et al. [30], the rise in temperature before a chemical reaction is approximately below 10−50K, while general consideration of initial temperature T0 is below 632.2K. Hence, the Darcy–Boussinesq approximation is adapted for the present problem.Figure 1
Problem geometry.


## 2.2. Governing Equations

With the assumption that viscous dissipation effects are negligible, the Oberbeck–Boussinesq approximation is employed, and Darcy’s law is applicable. By following the nanofluid model proposed by Buongiorno [11], the governing equations for total mass, momentum, energy, and nanoparticles, respectively, are written as(2)∇⟶⋅V⟶=0,0=−∇⟶P∗−μKV⟶+F⟶b+ρgj^,ρcfV⟶⋅∇⟶T=km∇⟶2T+ερcpDB∇⟶T⋅∇⟶C+DTT0∇⟶T⋅∇⟶T+Qk0ae−E/RT,1εV⟶⋅∇⟶C=DB∇⟶2C+DTT0∇⟶2T,where Fb→ is given by Fb→=σu→×B→×B→, B→=B0,0,0 is the constant magnetic field strength along x¯−axis, and the buoyancy term ρg is approximated by the Boussinesq approximation:(3)ρg=−ρf0β1−C0∂T∂x¯g+ρp−ρf0∂C∂x¯g.Here,V⟶=u¯,v¯,0 is the velocity vector with u¯ and v¯ the nanofluid velocity components along x¯− and y¯−axis, respectively, T is the nanofluid temperature, C is the nanoparticle volume fraction, and other quantities of physical interest are mentioned in Nomenclature. To nondimensionalize the governing equations along with boundary conditions, we make use of the following similarity transformations:(4)x=x¯H,y=y¯H,u=Hαmu¯,v=Hαmv¯,θ=ET−T0RT02,ϕ=CC0,P=P∗Kμαm.Following the traditional way of defining the dimensionless stream functionψ from velocity components as follows:(5)u¯=∂ψ∂y¯,v¯=−∂ψ∂x¯.Using the dimensionless variables defined in equation (4) and the stream function in equation (5), the governing equations now take the form:(6)∇2ψ+Ha2∂2ψ∂x2+Ra∂θ∂x−RaNr∂ϕ∂x=0,(7)∇2θ−∂ψ∂y∂θ∂x−∂ψ∂x∂θ∂y+Nb∂ϕ∂x∂θ∂x+∂ϕ∂y∂θ∂y+Nt∂θ∂x2+∂θ∂y2+Fkeθ=0,(8)∇2ϕ−Le∂ψ∂y∂ϕ∂x−∂ψ∂x∂ϕ∂y+NtNb∇2θ=0.The boundary conditions are as follows:(9)ψ=0,θ=0,ϕ=1,atx=0,(10)ψ=0,θ=0,ϕ=1,atx+y=1,(11)ψ=0,∂θ∂y=0,∂ϕ∂y=0,aty=0,where Ha is the magnetic field parameter also known as the modified Hartmann number, Ra is the Rayleigh number, Nr is the buoyancy-ratio parameter, Nb is the Brownian motion parameter, Nt is the thermophoresis parameter, Fk is the Frank-Kamenetskii number, Le is the Lewis number, and the Laplacian operator ∇2 are all defined as(12)Ha=σB02Kμ,Ra=1−C0ρf0gβKRT02LμαmE,Nr=ρp−ρf0C0Eρf0β1−C0RT02,Nb=τC0DBαm,Nt=τDTRT0αmE,Fk=EQk0aexpE/RT0L2RT02αm,Le=αmεDB,∇2=∂2∂x2+∂2∂y2.For problems related to heat and mass transfer across the cavity, it is essential to determine rates of heat and mass transfer in terms of local and average Nusselt and Sherwood numbers, respectively, using the following formulas:(13)Nu=−∂θ∂xx=0,Sh=−∂ϕ∂xx=0,Nu¯=∫01Nudy,Sh¯=∫01Shdy.


## 3. Numerical Method

We employ the finite difference method in combination with an iterative procedure [31] to solve the governing equations (6)–(8) along with boundary conditions (9). A five-point cell-centered stencil with a central difference formula is applied to obtain the system of linear equations which are then addressed in MATLAB. The step size is chosen to be uniform in both x and y directions. For the convergence criteria, 10−6 is chosen as a tolerance error for all dependent variables. The iterative procedure is terminated once the following convergence criteria are satisfied:(14)∑i=1M∑j=1Nχi,jn+1−χi,jn∑i=1M∑j=1Nχi,jn+1≤10−6,n≥1,where χ represents the dependent variables ψ, θ or ϕ, while M and N are the number of grid points along x and y directions, respectively, and n is the iteration index.


### 3.1. Test for Grid Independence

For correct usage of the finite difference method, it is essential to perform grid sensitivity analysis on the converged solutions. Thus, eight cases of a uniform grid varying from8×8 to 128×128 with an increment of 32×32 are tested. The numerical simulations for the grid independence test are shown for selected values of the model parameters Ra=10, Ha=2, Nr=0.1, Nb=Nt=0.4, Le=1, and Fk=0.1. The test is conducted by plotting the average Nusselt number against the number of grid points. As can be seen in Figure 2, the value of the average Nusselt number almost becomes independent of grid size when the grid points are higher than 128×128. Hence, a regular grid of 128×128 points has been selected for the following analysis.Figure 2
Average Nusselt number for different number of grid points.


### 3.2. Corroboration of the Numerical Code

We performed graphical as well as statistical tests for the corroboration of a numerical code. One of the two tests compares the average Nusselt numberNu¯ for a regular fluid with different values of Rayleigh number Ra=10−104. The results are reported in Table 1 for the parameters Nt=Nr=Nb=0,Ha=0,Le=0, and Fk=0, which reduce the problem to one mentioned in references [3, 10, 32–37].Table 1
Average Nusselt numberNu¯ for various values of Ra.
Previous results
R
a
10
100
1000
10000
Sheremet and Pop [10]
1.071
3.104
13.839
49.253
Walker and Homsy [32]
—
3.097
12.96
51.0
Bejan [33]
—
4.2
15.8
50.8
Beckermann et al. [34]
—
3.113
—
48.9
Gross et al. [35]
—
3.141
13.448
42.583
RamosS et al. [36]
1.065
2.801
—
—
Baytas and Pop [3]
1.079
3.16
14.06
48.33
Manole and Lage [37]
—
3.104
13.637
48.117
Present results
1.022
3.4
13.366
43.74The second test is performed to compare the obtained results for streamlines, isotherms, and isoconcentrations with the numerical data of [10]. They have studied the steady free convection flow and heat transfer in a triangular enclosure filled with a porous medium which is of the same geometry as our problem. The numerical simulations are performed by setting Ha=0 and Fk=0, which reduce the present study to the reference problem. Figures 3 and 4 compare the present results with previous results and are approximately the same compared to those of Sheremet and Pop [10]. These tests for the comparison of results obtained for the present study and their agreement boost the confidence of the current numerical scheme.Figure 3
Comparison of results forNr=0.1, Nb=0.1, and Nt=0.1 when Ra=100 and Le=1.0. (a) Streamlines ψ of Sheremet and Pop [10]. (b) Streamlines ψ of present results. (c) Isotherms θ of Sheremet and Pop [10]. (d) Isotherms θ of present results. (e) Isoconcentrations ϕ of Sheremet and Pop [10]. (f) Isoconcentrations ϕ of present results.
(a)
(b)
(c)
(d)
(e)
(f)Figure 4
Comparison of results forNr=0.4, Nb=0.4, and Nt=0.4 when Ra=100 and Le=1.0. (a) Streamlines ψ of Sheremet and Pop [10]. (b) Streamlines ψ of present results. (c) Isotherms θ of Sheremet and Pop [10]. (d) Isotherms θ of present results. (e) Isoconcentrations ϕ of Sheremet and Pop [10]. (f) Isoconcentrations ϕ of present results.
(a)
(b)
(c)
(d)
(e)
(f)


## 3.1. Test for Grid Independence

For correct usage of the finite difference method, it is essential to perform grid sensitivity analysis on the converged solutions. Thus, eight cases of a uniform grid varying from8×8 to 128×128 with an increment of 32×32 are tested. The numerical simulations for the grid independence test are shown for selected values of the model parameters Ra=10, Ha=2, Nr=0.1, Nb=Nt=0.4, Le=1, and Fk=0.1. The test is conducted by plotting the average Nusselt number against the number of grid points. As can be seen in Figure 2, the value of the average Nusselt number almost becomes independent of grid size when the grid points are higher than 128×128. Hence, a regular grid of 128×128 points has been selected for the following analysis.Figure 2
Average Nusselt number for different number of grid points.


## 3.2. Corroboration of the Numerical Code

We performed graphical as well as statistical tests for the corroboration of a numerical code. One of the two tests compares the average Nusselt numberNu¯ for a regular fluid with different values of Rayleigh number Ra=10−104. The results are reported in Table 1 for the parameters Nt=Nr=Nb=0,Ha=0,Le=0, and Fk=0, which reduce the problem to one mentioned in references [3, 10, 32–37].Table 1
Average Nusselt numberNu¯ for various values of Ra.
Previous results
R
a
10
100
1000
10000
Sheremet and Pop [10]
1.071
3.104
13.839
49.253
Walker and Homsy [32]
—
3.097
12.96
51.0
Bejan [33]
—
4.2
15.8
50.8
Beckermann et al. [34]
—
3.113
—
48.9
Gross et al. [35]
—
3.141
13.448
42.583
RamosS et al. [36]
1.065
2.801
—
—
Baytas and Pop [3]
1.079
3.16
14.06
48.33
Manole and Lage [37]
—
3.104
13.637
48.117
Present results
1.022
3.4
13.366
43.74The second test is performed to compare the obtained results for streamlines, isotherms, and isoconcentrations with the numerical data of [10]. They have studied the steady free convection flow and heat transfer in a triangular enclosure filled with a porous medium which is of the same geometry as our problem. The numerical simulations are performed by setting Ha=0 and Fk=0, which reduce the present study to the reference problem. Figures 3 and 4 compare the present results with previous results and are approximately the same compared to those of Sheremet and Pop [10]. These tests for the comparison of results obtained for the present study and their agreement boost the confidence of the current numerical scheme.Figure 3
Comparison of results forNr=0.1, Nb=0.1, and Nt=0.1 when Ra=100 and Le=1.0. (a) Streamlines ψ of Sheremet and Pop [10]. (b) Streamlines ψ of present results. (c) Isotherms θ of Sheremet and Pop [10]. (d) Isotherms θ of present results. (e) Isoconcentrations ϕ of Sheremet and Pop [10]. (f) Isoconcentrations ϕ of present results.
(a)
(b)
(c)
(d)
(e)
(f)Figure 4
Comparison of results forNr=0.4, Nb=0.4, and Nt=0.4 when Ra=100 and Le=1.0. (a) Streamlines ψ of Sheremet and Pop [10]. (b) Streamlines ψ of present results. (c) Isotherms θ of Sheremet and Pop [10]. (d) Isotherms θ of present results. (e) Isoconcentrations ϕ of Sheremet and Pop [10]. (f) Isoconcentrations ϕ of present results.
(a)
(b)
(c)
(d)
(e)
(f)


## 4. Results and Discussion

This section discusses the numerical simulation results and presents the graphical results of streamlines, isotherms, isoconcentrations, and average Nusselt number and average Sherwood number for various model parameters in a triangular enclosure. The effects of different values of Rayleigh numberRa=10−104, modified Hartmann number Ha=0−10, Brownian motion parameter Nb=0.1−0.4, thermophoresis parameter Nt=0.1−0.4, buoyancy-ratio parameter Nr=0.1−0.4, Lewis number Le=1−10, and Frank-Kamenetskii number Fk=0.1−0.3 on the flow and heat transfer characteristics are analyzed.


### 4.1. Explanation of Physical Parameters

It is useful to present the physical reasons for variation in each vital parameter. In fluid mechanics, Rayleigh numberRa is the property of a fluid that determines how heat is transferred throughout the fluid. The Rayleigh number is closely related to Grashof number Gr, and both numbers are used to describe natural convection and heat transfer by natural convection. The Rayleigh number is simply defined as the product of the Grashof number, which represents the relationship between buoyancy and viscosity within a fluid, and the Prandtl number, which describes the relationship between momentum diffusivity and thermal diffusivity. A smaller value of Rayleigh number represents the laminar flow, while a higher range indicates a turbulent flow. The buoyancy-ratio parameter describes an effect of the buoyancy force due to the concentration difference in comparison with the buoyancy force due to temperature difference. Nb, the Brownian motion parameter describes the movement of nanoparticles and enlightens the thermal effects of the nanofluid and nanoparticle diameter, whereas the thermophoresis parameter Nt explains the temperature gradient effect on nanoparticle diffusion depending on the thermal conductivity of the fluid and particle materials. The Lewis number is a dimensionless number, named after Warren K. Lewis, defined as the ratio of thermal diffusivity and mass diffusivity. It is used to characterize fluid flows where there is simultaneous heat and mass transfer. The Lewis number is, therefore, a measure of the relative thermal and concentration boundary layer thicknesses. Ha is a dimensionless number known as Hartmann number, which gives a measure of the relative importance of drag forces resulting from magnetic induction and viscous forces in Hartmann flow and determines the velocity profile for such flow.


### 4.2. Effects of Frank-Kamenetskii and Rayleigh Numbers

The impact of the Frank-Kamenetskii number and the Rayleigh number on the evolution of streamlines, isotherms, and isoconcentrations is described in Figure5 when Ra=10 and Fk=0.1 and in Figure 6 when Ra=100 and Fk=0.1,1,and3, respectively. From these figures, it is observed that in the case of streamlines, two opposite rotating vortices appear which are independent of the values of Ra and Fk. The values for other model parameters are Le=10, Nb=0.4, Nr=0.1, Nt=0.4, and Ha=1. A decrease in thermal and concentration boundary layer thicknesses is observed with the increase of Ra and Fk. The value of θmax decreases with the rise in Ra, while it increases with the rise of Fk. The exothermic reaction produces more heat with growing value for Fk, hence accelerating the convection phenomenon within the cavity. For the smaller value of Ra, the isoconcentrations become concentrated towards the vertical wall where the heat is low. This concentration is due to the thermophoresis phenomenon since the nanoparticles flee towards cool regions from the hot regions. However, a vortex appears near a bottom wall for a higher value of Ra, which shifts towards the center of the cavity as the value Fk increases. Increase in Ra and Fk also intensifies the convective flow in the cavity accounting by the maximum ψ.Figure 5
Stream functionψ, isotherms θ, and isoconcentrations ϕ for Ra=10, Ha=1, Le=10, Nr=0.1, Nb=Nt=0.4, and Fk=0.1 for Figures (a)–(c), Fk=1 for Figures (d)–(f), and Fk=3 for Figures (g)–(i). (a) ψ. (b) θ. (c) ϕ. (d) ψ. (e) θ. (f) ϕ. (g) ψ. (h) θ. (i) ϕ.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)Figure 6
Stream functionψ, isotherms θ, and isoconcentrations ϕ for Ra=100, Ha=1, Le=10, Nr=0.1, Nb=Nt=0.4, and Fk=0.1 for Figures (a)–(c), Fk=1 for Figures (d)–(f), and Fk=3 for Figures (g)–(i). (a) ψ. (b) θ. (c) ϕ. (d) ψ. (e) θ. (f) ϕ. (g) ψ. (h) θ. (i) ϕ.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)The monotonic effects of the Rayleigh number and the Frank-Kamenetskii number on heat transfer in terms of local Nusselt are depicted in Figure7. An accretion of Ra with Fk leads to an increase in Nu due to an interaction of the thermal boundary layers close to the vertical walls. Similar results are obtained for local Sherwood number and plotted in Figure 8. For the higher value of Ra, the local Sherwood number decreases with y¯, and particles become concentrated towards the vertical wall with a less significant effect of increasing Frank number Fk, resulting in a decrease in the nanoparticle volume fraction (Figure 9(b)). Figure 9 depicts that average Sherwood and Nusselt numbers decrease with the increasing value of Ra. This shows that the rate of heat and mass transfer decreases along the vertical wall, and the heat transfer is dominated by the conduction. Figure 10 is the mesh diagram for stream function ψ depicting the formation of two vortices. This pattern can also be viewed in Figures 5 and 6, where we noticed that the nanofluid flows inside the triangular cavity which occurs in two major vortices. One vortex has a clockwise rotation in the left, while the second vortex has an anticlockwise rotation in the right half of the cavity. So, the Lorentz force increases the motion in the left half, and it decelerates the motion in the other half because the flow and the field act in the opposite direction. Since the applied magnetic field prevents the convection current, as a result, the average Nusselt number and Sherwood number decreases.Figure 7
Variation of local Nusselt number whenLe=10, Ha=1, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)Figure 8
Variation of local Sherwood number whenLe=10, Ha=1, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)Figure 9
Variation of average Nusselt number (a) and average Sherwood number (b) at the vertical wall with varyingRa when Le=10, Ha=1, Nb=0.4, Nr=0.1, and Nt=0.4.
(a)
(b)Figure 10
Mesh diagram forψ.


### 4.3. Effects of Hartmann and Rayleigh Numbers

In Figures11 and 12, the effects of the Hartmann number on local Nusselt and Sherwood numbers with the variation of the Rayleigh number are plotted against the vertical wall of the cavity. For a weak flow (Ra = 10), the effect of Hartmann number on local Nusselt number is negligible, as shown in Figure 11(a). This is because when the Rayleigh number is small, the flow convection is insignificant and the heat transfer in the cavity is dominated by conduction. However, in the case of higher Rayleigh number, the Nu increases but depicts a small variation with Ha, as observed in Figure 11(b). Due to high Rayleigh number, the flow convection becomes stronger which in turn increases the heat transfer in the cavity. Similar to Figure 12(a), the effect of the Hartmann number on Sherwood number is minimal, and it decreases for both the lower and higher values of Rayleigh number, as shown in Figure 12(b).Figure 11
Variation of local Nusselt number whenLe=10, Fk=1, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)Figure 12
Variation of local Sherwood number whenLe=10, Fk=1, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)It is also noted that a magnetic field with high intensity prevents the convection and may cause separation among nanoparticles. They are hence breaking the linearity of the gradient of nanoparticle concentration. To further check this aspect, we have plotted the concentration gradients∂ϕ/∂y along the center line (x=0.5) of the enclosure for both cases of Rayleigh number, i.e., Ra=10 and Ra=100, for different values of Hartmann number, as displayed in Figure 13. For higher Rayleigh number, Ra=100, the nanoparticles gather and start concentrating along the boundaries of the cavity away from the center, as can also be seen from Figure 6. Also, the deviations in gradients for nanoparticle concentration are nonlinear for a strong flow (Ra=100), while for a weak flow (Ra=10), ∂ϕ/∂y is completely linear. Figure 14 shows the effects of the Hartmann number and Rayleigh number on local Nusselt and local Sherwood numbers. As the Hartmann number increases, the Nusselt number decreases. Similar is the trend observed for the Rayleigh number. The existence of the maximum and minimum points in the local Nusselt number profile corresponds to the presence of thermal plumes.Figure 13
Gradients for nanoparticle concentration alongx=0.5 for different values of Le=10, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)Figure 14
Variation of local Nusselt number (a) and local Sherwood number (b) whenLe=10, Fk=1, Nb=0.4, Nr=0.01, Nt=0.1, and Ra=1000.
(a)
(b)


### 4.4. Effects of Brownian Diffusion and Thermophoresis Parameter

This section elaborates the effects of thermophoresisNt and Brownian diffusion Nb parameters on the rate of heat and mass transfer which are observed along the vertical wall (x¯=0). Figure 15 illustrates the variation of local Nusselt and local Sherwood numbers for Brownian diffusion parameter. It is noticed from Figure 15(a) that the impact of Nb is quite less on the Nusselt number which distinctively agrees with the conclusions of Saghir and Mohamed [38]. However, the effect of these parameters on the mass transfer rate (Sherwood number) is quite noticeable. From Figure 15(b), it can be observed that as the particle size is smaller, the mass transfer is quite increased, while an increase in particle size diminishes the Brownian diffusion effect resulting in no change in the rate of mass transfer.Figure 15
Deviations in local Nusselt number (a) local Sherwood number (b) whenLe=10, Ha=1, Fk=1, Nr=0.1, Nt=0.4, and Ra=100.
(a)
(b)Figure16 depicts the effects of thermophoresis parameter (Nt) on heat and mass transfer. Physically, the nanoparticles move from the hotter region towards the colder region due to the thermophoresis effect. From Figure 16(a), we found that the effect of Nt is almost negligible on the rate of heat transfer for the nanofluid flow inside the cavity filled with the porous medium. Once again the results reconfirm the conclusions of Saghir and Mohamed [38]. Figure 16(b) shows that the rate of mass transfer decreases along the vertical wall. But, the rate of the mass transfer becomes prominent with the growing Nt effect which leads to the homogeneous distribution of particles within a cavity.Figure 16
Deviations in local Nusselt number (a) local Sherwood number (b) whenLe=10, Nb=0.4, Nr=0.1, and Ra=100.
(a)
(b)


## 4.1. Explanation of Physical Parameters

It is useful to present the physical reasons for variation in each vital parameter. In fluid mechanics, Rayleigh numberRa is the property of a fluid that determines how heat is transferred throughout the fluid. The Rayleigh number is closely related to Grashof number Gr, and both numbers are used to describe natural convection and heat transfer by natural convection. The Rayleigh number is simply defined as the product of the Grashof number, which represents the relationship between buoyancy and viscosity within a fluid, and the Prandtl number, which describes the relationship between momentum diffusivity and thermal diffusivity. A smaller value of Rayleigh number represents the laminar flow, while a higher range indicates a turbulent flow. The buoyancy-ratio parameter describes an effect of the buoyancy force due to the concentration difference in comparison with the buoyancy force due to temperature difference. Nb, the Brownian motion parameter describes the movement of nanoparticles and enlightens the thermal effects of the nanofluid and nanoparticle diameter, whereas the thermophoresis parameter Nt explains the temperature gradient effect on nanoparticle diffusion depending on the thermal conductivity of the fluid and particle materials. The Lewis number is a dimensionless number, named after Warren K. Lewis, defined as the ratio of thermal diffusivity and mass diffusivity. It is used to characterize fluid flows where there is simultaneous heat and mass transfer. The Lewis number is, therefore, a measure of the relative thermal and concentration boundary layer thicknesses. Ha is a dimensionless number known as Hartmann number, which gives a measure of the relative importance of drag forces resulting from magnetic induction and viscous forces in Hartmann flow and determines the velocity profile for such flow.


## 4.2. Effects of Frank-Kamenetskii and Rayleigh Numbers

The impact of the Frank-Kamenetskii number and the Rayleigh number on the evolution of streamlines, isotherms, and isoconcentrations is described in Figure5 when Ra=10 and Fk=0.1 and in Figure 6 when Ra=100 and Fk=0.1,1,and3, respectively. From these figures, it is observed that in the case of streamlines, two opposite rotating vortices appear which are independent of the values of Ra and Fk. The values for other model parameters are Le=10, Nb=0.4, Nr=0.1, Nt=0.4, and Ha=1. A decrease in thermal and concentration boundary layer thicknesses is observed with the increase of Ra and Fk. The value of θmax decreases with the rise in Ra, while it increases with the rise of Fk. The exothermic reaction produces more heat with growing value for Fk, hence accelerating the convection phenomenon within the cavity. For the smaller value of Ra, the isoconcentrations become concentrated towards the vertical wall where the heat is low. This concentration is due to the thermophoresis phenomenon since the nanoparticles flee towards cool regions from the hot regions. However, a vortex appears near a bottom wall for a higher value of Ra, which shifts towards the center of the cavity as the value Fk increases. Increase in Ra and Fk also intensifies the convective flow in the cavity accounting by the maximum ψ.Figure 5
Stream functionψ, isotherms θ, and isoconcentrations ϕ for Ra=10, Ha=1, Le=10, Nr=0.1, Nb=Nt=0.4, and Fk=0.1 for Figures (a)–(c), Fk=1 for Figures (d)–(f), and Fk=3 for Figures (g)–(i). (a) ψ. (b) θ. (c) ϕ. (d) ψ. (e) θ. (f) ϕ. (g) ψ. (h) θ. (i) ϕ.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)Figure 6
Stream functionψ, isotherms θ, and isoconcentrations ϕ for Ra=100, Ha=1, Le=10, Nr=0.1, Nb=Nt=0.4, and Fk=0.1 for Figures (a)–(c), Fk=1 for Figures (d)–(f), and Fk=3 for Figures (g)–(i). (a) ψ. (b) θ. (c) ϕ. (d) ψ. (e) θ. (f) ϕ. (g) ψ. (h) θ. (i) ϕ.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)The monotonic effects of the Rayleigh number and the Frank-Kamenetskii number on heat transfer in terms of local Nusselt are depicted in Figure7. An accretion of Ra with Fk leads to an increase in Nu due to an interaction of the thermal boundary layers close to the vertical walls. Similar results are obtained for local Sherwood number and plotted in Figure 8. For the higher value of Ra, the local Sherwood number decreases with y¯, and particles become concentrated towards the vertical wall with a less significant effect of increasing Frank number Fk, resulting in a decrease in the nanoparticle volume fraction (Figure 9(b)). Figure 9 depicts that average Sherwood and Nusselt numbers decrease with the increasing value of Ra. This shows that the rate of heat and mass transfer decreases along the vertical wall, and the heat transfer is dominated by the conduction. Figure 10 is the mesh diagram for stream function ψ depicting the formation of two vortices. This pattern can also be viewed in Figures 5 and 6, where we noticed that the nanofluid flows inside the triangular cavity which occurs in two major vortices. One vortex has a clockwise rotation in the left, while the second vortex has an anticlockwise rotation in the right half of the cavity. So, the Lorentz force increases the motion in the left half, and it decelerates the motion in the other half because the flow and the field act in the opposite direction. Since the applied magnetic field prevents the convection current, as a result, the average Nusselt number and Sherwood number decreases.Figure 7
Variation of local Nusselt number whenLe=10, Ha=1, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)Figure 8
Variation of local Sherwood number whenLe=10, Ha=1, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)Figure 9
Variation of average Nusselt number (a) and average Sherwood number (b) at the vertical wall with varyingRa when Le=10, Ha=1, Nb=0.4, Nr=0.1, and Nt=0.4.
(a)
(b)Figure 10
Mesh diagram forψ.


## 4.3. Effects of Hartmann and Rayleigh Numbers

In Figures11 and 12, the effects of the Hartmann number on local Nusselt and Sherwood numbers with the variation of the Rayleigh number are plotted against the vertical wall of the cavity. For a weak flow (Ra = 10), the effect of Hartmann number on local Nusselt number is negligible, as shown in Figure 11(a). This is because when the Rayleigh number is small, the flow convection is insignificant and the heat transfer in the cavity is dominated by conduction. However, in the case of higher Rayleigh number, the Nu increases but depicts a small variation with Ha, as observed in Figure 11(b). Due to high Rayleigh number, the flow convection becomes stronger which in turn increases the heat transfer in the cavity. Similar to Figure 12(a), the effect of the Hartmann number on Sherwood number is minimal, and it decreases for both the lower and higher values of Rayleigh number, as shown in Figure 12(b).Figure 11
Variation of local Nusselt number whenLe=10, Fk=1, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)Figure 12
Variation of local Sherwood number whenLe=10, Fk=1, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)It is also noted that a magnetic field with high intensity prevents the convection and may cause separation among nanoparticles. They are hence breaking the linearity of the gradient of nanoparticle concentration. To further check this aspect, we have plotted the concentration gradients∂ϕ/∂y along the center line (x=0.5) of the enclosure for both cases of Rayleigh number, i.e., Ra=10 and Ra=100, for different values of Hartmann number, as displayed in Figure 13. For higher Rayleigh number, Ra=100, the nanoparticles gather and start concentrating along the boundaries of the cavity away from the center, as can also be seen from Figure 6. Also, the deviations in gradients for nanoparticle concentration are nonlinear for a strong flow (Ra=100), while for a weak flow (Ra=10), ∂ϕ/∂y is completely linear. Figure 14 shows the effects of the Hartmann number and Rayleigh number on local Nusselt and local Sherwood numbers. As the Hartmann number increases, the Nusselt number decreases. Similar is the trend observed for the Rayleigh number. The existence of the maximum and minimum points in the local Nusselt number profile corresponds to the presence of thermal plumes.Figure 13
Gradients for nanoparticle concentration alongx=0.5 for different values of Le=10, Nb=0.4, Nr=0.1, Nt=0.4, and Ra=10 (a) and Ra=100 (b).
(a)
(b)Figure 14
Variation of local Nusselt number (a) and local Sherwood number (b) whenLe=10, Fk=1, Nb=0.4, Nr=0.01, Nt=0.1, and Ra=1000.
(a)
(b)


## 4.4. Effects of Brownian Diffusion and Thermophoresis Parameter

This section elaborates the effects of thermophoresisNt and Brownian diffusion Nb parameters on the rate of heat and mass transfer which are observed along the vertical wall (x¯=0). Figure 15 illustrates the variation of local Nusselt and local Sherwood numbers for Brownian diffusion parameter. It is noticed from Figure 15(a) that the impact of Nb is quite less on the Nusselt number which distinctively agrees with the conclusions of Saghir and Mohamed [38]. However, the effect of these parameters on the mass transfer rate (Sherwood number) is quite noticeable. From Figure 15(b), it can be observed that as the particle size is smaller, the mass transfer is quite increased, while an increase in particle size diminishes the Brownian diffusion effect resulting in no change in the rate of mass transfer.Figure 15
Deviations in local Nusselt number (a) local Sherwood number (b) whenLe=10, Ha=1, Fk=1, Nr=0.1, Nt=0.4, and Ra=100.
(a)
(b)Figure16 depicts the effects of thermophoresis parameter (Nt) on heat and mass transfer. Physically, the nanoparticles move from the hotter region towards the colder region due to the thermophoresis effect. From Figure 16(a), we found that the effect of Nt is almost negligible on the rate of heat transfer for the nanofluid flow inside the cavity filled with the porous medium. Once again the results reconfirm the conclusions of Saghir and Mohamed [38]. Figure 16(b) shows that the rate of mass transfer decreases along the vertical wall. But, the rate of the mass transfer becomes prominent with the growing Nt effect which leads to the homogeneous distribution of particles within a cavity.Figure 16
Deviations in local Nusselt number (a) local Sherwood number (b) whenLe=10, Nb=0.4, Nr=0.1, and Ra=100.
(a)
(b)


## 5. Conclusions

A numerical analysis has been carried out to study the free convection within the triangular cavity filled with a porous medium saturated with nanofluids in the presence of magnetic field initiated by a chemical reaction influenced by Arrhenius kinetics. The finite difference method has been applied to solve the governing equations. The numerical results have been obtained and discussed for a wide range of physical parameters appearing in the problem. It is found that increasing the value of Rayleigh number decreases both the average Nusselt number and average Sherwood number. Moreover, the local Nusselt number increases with the variation of Frank-Kamenetskii number, Hartmann number, Brownian diffusion, and thermophoresis parameter for the increasing value of Rayleigh number. But, local Sherwood number decreases with the variation of Frank-Kamenetskii number, Hartmann number, Brownian diffusion, and thermophoresis parameter for the increasing value of Rayleigh number.



---

*Source: 7451967-2019-12-06.xml*