# Cooling a Hot Semiannulus with Constant Heat Flux by UsingFe3O4-Water Nanofluid and a Magnetic Field: Natural Convection Mechanism

Authors: Zhengqiang Yang; S. M. Bouzgarrou; Riadh Marzouki; Fatma Aouaini; M. A. El-Shorbagy; Mahidzal Dahari; Said Anwar Shah; D. L. Suthar
Journal: Journal of Nanomaterials
 (2022)

---

## Abstract

In this simulation, the nanoparticle distribution and entropy generation were studied using the Buongiorno’s developed two-phase model and magnetic field inside a porous semiannulus cavity. The influence of three terms was considered in the Buongiorno’s developed two-phase model such as Brownian motion, thermophoresis, and magnetophoresis effects. In addition, the entropy generation was assessed due to temperature and velocity gradient. The evidence showed that the effects of the magnetic field in high porosities and volume fraction of nanoparticles were remarkable on the Nusselt number and entropy generation. Also, irreversibility due to heat transfer is much greater than fluid friction.


---

## Body

## 1. Introduction

Nowadays, heat transfer processes have widely exerted for many applications such as cooling battery, CPU, forging, radiator in vehicles, heating the home, and powerplant [1–5]. One method that is popular between thermal engineering is injecting various nanoparticles such as metal, nonmetal, and oxide nanoparticles and nanoencapsulated phase change material (NEPCM) to host fluid [6–8], because mixing nanoparticles into host fluid can changed the thermophysical properties, specifically, thermal conductivity of host fluid [9–11].In recent years, the heat transfer rates by computational fluid dynamics (CFD) and entropy generation are significantly investigated by researchers and companies due to having accurate with experimental studies [12–15]. Moreover, owing to its large specific area and higher solid thermal conductivity, there has been growing interest in heat transfer inside the porous media. According to previous studies, porous structure improves thermal performance of nanofluid [16–19]. A numerical investigation of heat transfer and fluid flow in a parabolic trough solar receiver with internal annular porous structure and synthetic oil-Al2O3 nanofluid was carried out by Bozorg et al. [20]. According to their results, utilization of porous structure and nanofluids enhances heat transfer coefficient 7% and 20%, respectively. Jamal-Abad et al. [21] experimentally studied the thermal efficiency of a solar parabolic trough collector filled with porous media; they illustrated an enhancement in efficiency of the collector by increasing the mass flow rate.It is confirmed that in nanofluid simulation, two-phase approach provides better accuracy compared to single-phase approach [22–24]. A two-phase model presented by Buongiorno [25] has received significant attention; in this model, it is suggested that among seven slip mechanisms, Brownian motion and thermophoresis diffusion play vital role in nanofluid distribution. Natural convection of nanofluid in an inclined cavity and inside porous medium considering two-phase approach is studied by studies [6] and [15]. Their results demonstrate good agreement with experimental studies. Thermophoresis and Brownian motion effect on boundary layer flow of nanofluid in presence of thermal stratification due to solar energy is analyzed by Anbuchezhian et al. [26]. Their results highlighted that Brownian motion and thermophoresis distribution can affect the heat transfer properties; they reported a substantial impact on the boundary layer flow field by Brownian motion in the presence of thermal stratification. Kaloudis et al. [27] numerically investigated on parabolic trough solar collector with nanofluid using a two-phase model. They reported that two-phase simulation of nanofluids in solar studies shows better agreement with experimental studies. Also, the presence of nanoparticles improves collector’s efficiency.To obtain the optimal configuration of solar collectors, it is essential to analyze the entropy generation. Farshad and Sheikholeslami [28] scrutinize exergy loss and heat transfer of mixture of aluminum oxide and H2O through a solar collector. Thermal performance and entropy generation analysis of a high concentration ratio parabolic trough solar collector was studied by Mwesigye et al. [29]. They reported a decrement in entropy generation with augmentation of nanofluid volume fraction for some ranges of Re number. Verma et al. [30] experimentally analyzed exergy efficiency and entropy generation in flat plate solar collectors for different types of nanofluids. They highlighted that rise of Bejan number towards unity illustrates the improvement of system performance due to efficient conversion of the available energy into useful functions. Also, Sheikholeslami et al. [31] scrutinized impact of Lorentz forces on magnetic nano fluid of Fe3O4 with entropy and exergy analyzing inside a semiannulus. Afrand et al. [32] studied free convective heat transfer and entropy generation of Al2O3-water nanofluid in a triangular enclosure. They illustrated that the Bijan number increases by decreasing the Ra and increasing the Ha. The maximum heat transfer rate takes place at the enclosure angle of 60°.To the best knowledge of the authors, there has been no detailed investigation of impacts of nonuniform magnetic fields on PTC thermal performance using nanofluid, considering Brownian motion and thermophoresis distribution as well as entropy generation. In the present work, a coil is wrapped around semiannulus to produce a variable magnetic field. Also, it focuses on the local distribution of nanoparticles, entropy generation due to fluid friction and heat transfer, and Nusselt number variation.


## 2. Physical Model

In the present study, 2-dimensional and steady-state natural convection is simulated in a semiannulus enclosure. Schematic presentation of the problem is presented in Figure1. In Figure 1, two horizontal walls are thermally insulated, the inner semicircle wall at constant heat flux q", and the outer semicircle wall is fixed at constant temperature Tc. The working fluid is Fe3O4-water nanofluid. Natural convection fluid flow is simulated based on Boussinesq’s approximation. Also, heat transfer, nanoparticles distribution, and entropy generation in the presence of a nonuniform magnetic field are investigated. Therefore, it is assumed that nanoparticle distribution is based on Buongiorno two-phase model. According to the aforementioned assumption, the governing equations are presented in the next section.Figure 1
Schematic presentation of the problem.


## 3. Relations and Hypothesis

The nondimensional parameters are as follows:X=x/L,Y=y/L,T∗=kfT−Tc/q′′L,V∗=VL/vf,P∗=PL2/ρfvf2,H∗=H/H0,M∗=M/M0,φ∗=φ/φAve,u∗=u/αfL,v∗=v/αfL,∇∗=L∇,δ=kfTc/q′′L,DB∗=DB/DB0,DT∗=DT/DT0,DT0=γμf/ρfφAve,H0=I/2πλ,DB0=KBTc/3πμfdnp,andM0=χφAve,TfH0forφAve=0.02,andTm=Tc+q′′L/2kf.The dimensionless forms of continuity equations, momentum, energy, and volume fraction are as follows:Nondimensional continuity equation:(1)∇∗.V∗=0.Nondimensional momentum equation:(2)1ε2ρnfρfV∗.∇∗V∗=−∇∗P∗−μnfμf1DaV∗+∇∗.μnfμf∇∗V∗+ρβnfρfβfRafPrT∗.e^+MnM∗.∇∗H∗.Nondimensional heat transfer equation:(3)ρCPnfρCPfV∗.∇∗T∗=1Pr∇∗keqkf∇∗T∗+ε1Pr1LeDB∗∇∗φ∗.∇∗T∗+DT∗NBT∇∗T∗.∇∗T∗1+T∗/δ−DB∗ξLξφ∗∇∗H∗.∇∗T∗H∗.Nondimensional mass transfer equation:(4)1εV∗.∇∗φ∗=1Sc∇∗.DB∗∇∗φ∗+DT∗NBT∇∗T∗1+T∗/δ−DB∗ξLξφ∗∇∗H∗H∗.The dimensionless numbers in the above relationships are defined as:(5)Da=KL2,Pr=vfαf,αf=kfρcpf,Raf=gBfq′′L4kfαfvf,Rap=RafDa,e^=gg,Mn=μ0H0M0L2ρfvf2,Le=kfρCpnpDB0φAve,NBT=φAveDB0δDT0,Sc=vfDB0.Mn is the magnetic number, which is defined as the ratio of the Kelvin force to the kinematic viscosity. Moreover, Da,Ra,Pr,Le,andSc denote the Darcy, Rayleigh, Prandtl, Lewis, and Schmitt numbers, respectively.The average Nusselt on the constant heat flux wall (inner cylinder wall) is calculated as follows:(6)Nuloc=keqkeff1T∗,Nuave=1π∫0πNulocζdζ.In this study, the entropy generation is considered due to the irreversibility of the velocity gradients and temperature gradients. According to Shavik et al. [33], the entropy generation is obtained as
(7)ss=μnfTm2∂u∂x2+2∂v∂y2+∂u∂y+∂v∂x2+knfTm2∂T∂x2+∂T∂y2.In the mentioned equation, the first and second terms are local entropy generation due to the fluid friction (SL,FF) and the heat transfer (SL,HT) irreversibility, respectively. Also, the dimensionless equation of entropy generation is as follows:
(8)SL,HT=knfkf∂T∗∂X2+∂T∗∂Y2,SL,FF=Χμnfμf2∂u∗∂X2+∂v∗∂Y2+∂u∗∂Y+∂v∗∂X2,Χ=μfkfTmαfL2q"2,SL,s=ss×kfTm2q"2=SL,FF+SL,HT.Bejan number is the ratio of the entropy generation due to the heat transfer to the total entropy generation. This is as follows:(9)BeL=SL,HTSL,s.Total entropy generation and the total Bejan number are obtained by integrating the relations of local entropy generation and Bejan number.(10)ST,HT=∫ASL,HTdA,ST,FF=∫ASL,FFdA,ST,s=∫ASL,sdA,Beave=∫ABeLdA∫AdA.


## 4. CFD Setting

The presented nonlinear governing PDE equations are solved based on finite volume method (FVM). The continuity and momentum equations are coupled and solved in an algorithm termed SIMPLE. The energy and concentration equations simultaneously solved. The advection terms in the governing equations are discretized based on first-order upwind schemes, and diffusion terms are solved based on second-order central schemes. The numerical convergence criterion was residual values. The residual values at convergence for velocity and pressure fields were10−5, and for temperature and φ were 10−6. Under relaxation factors for velocity, pressure, temperature, and φ were 0.4, 0.6, 0.2, and 0.01, respectively. Square uniform mesh is selected for presented study, based on nanoparticle distribution on the hot wall 40,000 (200×200) numbers of mesh selected. Moreover, for validation of the entropy generation, the results of the present study is compared with the work of Shavik et al. [33] (Table 1). In all cases, the results of this study are in good agreement.Table 1
Comparison of total entropy generation (Sg,tot) and Bejan number (Beave), for Pr=0.71 and irreversibility (Χ=10−4).
Ra=103Ra=105Sg,totBeaveSg,totBeavePresent work1.1540.9723.20.193Shavik et al. [33]1.150.9723.270.194


## 5. Results and Discussion

Extended Buongiorno’s two-phase model is considered for the distribution of nanoparticles. The numerical simulation has investigated for porous Rayleigh numberRap=10and1000, the volume fraction of nanoparticles φAve=0,0.01and0.03, porosity number ε=0.4and0.7, and magnetic number 0≤Mn≤8×107. Constant values included Pr=4.623, Sc=3.55×104, Tc=310K, q′′=48.01w/m2, 1.71×105<Le<6.84×105, Da=10−3K=0.625×10−6, δ=161, NBT=0.245, Χ=3×10−11, and Tm=Tc+Lq′′/2kf=310.96K.The effects of the mentioned parameters have studied the distribution of nanoparticles and the entropy generation contours due to fluid friction and heat transfer, and Bejan.Figure2 shows the nondimensional distribution of nanoparticles on the lower and upper cylinder walls in φAve=0.03 and Rap=10and1000 for the different magnetic numbers. In Figure 2(a) and porous Rayleigh number 10, the density of nanoparticles is higher on the top wall (cold) than the down wall (hot); that is due to the thermophoresis term in the volume fractional equation. With the rising of porous Rayleigh number, the dimensional distribution of nanoparticles is almost identical due to the increasing flow velocity on the inner and outer walls. By increasing the magnetic number to Mn=2×107 for all of the cases, the density of nanoparticles for both walls is almost equal to the unit value, and no change occurs. But with the increasing magnetic number to Mn=8×107, the peak of nanoparticle density appears near the wires, which is due to the absorption of nanoparticles by the magnetic field. The above arguments are also true for Figure 2(b), which the porosity has increased to 0.7. Besides, in all graphs, the peak of nanoparticle density in the inner wall is wider than the outer wall.Figure 2
Nondimensional distribution of nanoparticles on the outer and inner cylinder walls forφAve=0.03, Rap=10and1000, and different magnetic numbers at (a) ε=0.4 and (b) ε=0.7.
(a)(b)Figure3 shows a better comparison of the effect of different parameters on the dimensional distribution of nanoparticles. Figure 3(a) is shown for magnetic numbers and different volume fractions in Rap=1000 and ε=0.7. The dimensionless nanoparticle density is the same for the volume fraction 0.01 and 0.03 in a constant magnetic number. Therefore, the volume fraction does not affect the dimensionless distribution of nanoparticles. Also, in a constant volume fraction, the density of nanoparticles increases near the wires due to increasing the magnetic number. Figure 3(b) is plotted for different porosity and porous Rayleigh number in φAve=0.03 and Mn=8×107. Porosity and porous Rayleigh numbers do not affect on the dimensionless distribution of nanoparticles at the high magnetic field.Figure 3
Nondimensional distribution of nanoparticles on inner and outer walls (with angle ofζ) for (a) different magnetic numbers and volume fraction at Rap=1000 and ε=0.7, (b) different porous Rayleigh number and porosity at φAve=0.03 and Mn=8×107.
(a)(b)Figure4 illustrates the dimensionless local entropy generation due to heat transfer (SL,HT) and fluid friction (SL,FF) and summation (SL,s) inside a semiannulus for φAve=0.03, ε=0.7, and Mn=0and8×107 for different porous Rayleigh number. In Figure 4(a) in the absence of a magnetic field and Rap=10, the SL,HT contours are similar to temperature contours. The highest value of SL,HT is on the inner cylinder wall (hot) due to extreme temperature gradients. In Figure 4(b), as increasing porous Rayleigh number to 1000, the shape of the contours SL,HT changes completely, in such a way that the densities of the contours increase near the inner cylinder wall and a core is created near the outer cylinder wall. By adding a magnetic field, the density of SL,HT contours increase near the wires, and, with the rise of porous Rayleigh number, there is no change in the shape of the contours. In all cases, due to gradients of high velocity on the walls and near the wires, the maximum value of SL,FF are in these regions. By comparing the values of SL,HT and SL,FF, it can be seen that the effects of entropy generation due to heat transfer are much greater than the entropy generation due to fluid friction. As a result, the contours SL,s are very similar to SL,HT contours.Figure 4
The dimensionless local entropy generation due to heat transfer (SL,HT), fluid friction (SL,FF) and summation (SL,s) for φAve=0.03, ε=0.7, and Mn=0and8×107at (a) Rap=10 and (b) Rap=1000.
(a)(b)Figure5 presents total entropy generation due to the heat transfer (ST,HT) and fluid friction (ST,FF) for the ε=0.7 and different volume fractions of nanoparticle in porous Rayleigh numbers 10 and 1000. According to Figure 5(a), by intensifying the magnetic field for both porous Rayleigh numbers, ST,HT decreases linearly for the volume fraction of nanoparticles 0.01 and 0.03. In Figure 5(b), by intensifying the magnetic field at both porous Rayleigh numbers 10 and 1000, ST,FF decreases linearly and nonlinearly for volume fractions 0.01 and 0.03, respectively. Also, in all cases and the absence of a magnetic field, the value of ST,FF is almost zero. Figure 6(a) shows the dimensionless summation total entropy generation (ST,s) at ε=0.7 and different volume fraction of nanoparticles. According to the charts, the charts of ST,s are similar to ST,HT charts, because of the values of ST,HT dominates the values of ST,FF. Also, according to Figure 6(b), with the increasing magnetic number for both porous Rayleigh numbers, Bejan number decreases linearly and nonlinearly for volume fractions 0.01 and 0.03, respectively.Figure 5
Dimensionless total entropy generation due to heat transfer (ST,HT) and fluid friction (ST,FF) for ε=0.7 and different volume fraction of nanoparticle.
(a)(b)Figure 6
Dimensionless total entropy generation due to summation (ST,s) and Bijan number for ε=0.7 and different volume fraction of nanoparticle.
(a)(b)


## 6. Conclusion

In present work, effects of porous Rayleigh number, the volume fraction of nanoparticles, porosity, and magnetic number are investigated on nondimensional distribution of nanoparticles and entropy generation. The main findings can be condensed following point:(i)
Adding a magnetic field increases the distribution of nanoparticles near the wires and causes the formation of vortices and increasing the flow velocity(ii)
In the presence of the nonuniform magnetic field, with the increasing porosity and porous Rayleigh number, the distribution of nanoparticles becomes uniform. But, in the absence of a magnetic field, porosity and porous Rayleigh number do not affect the dimensionless distribution of nanoparticles(iii)
By increasing magnetic number and volume fraction of nanoparticles, Bejan number, entropy generation due to heat transfer, and summation decrease but entropy generation due to fluid friction increases(iv)
Entropy generation due to heat transfer is much greater than the entropy generation due to fluid friction(v)
With increasing the magnetic number, entropy generation due to fluid friction increases near the wires



---

*Source: 4761989-2022-05-26.xml*