Next Article in Journal
Analytical Model of Interaction of an Oil Pipeline with a Support of an Overpass Built in a Mountainous Area
Previous Article in Journal
Determination of Geometrical Deviations of Large-Size Crankshafts with Limited Detection Possibilities Resulting from the Assumed Measuring Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Analysis of Magneto-Natural Convection Flows within a Partially Thermally Active Rectangular Enclosure

1
School of Mechanical and Mechatronic Engineering, University of Technology Sydney, Sydney, NSW 2007, Australia
2
Department of Mathematics, COMSATS University Islamabad, Islamabad 45550, Pakistan
3
Department of Science and Humanities, Sir Syed CASE Institute of Technology, Islamabad 44000, Pakistan
4
School of Energy and Power Engineering, Jiangsu University, Zhenjiang 212013, China
*
Author to whom correspondence should be addressed.
Energies 2023, 16(11), 4462; https://doi.org/10.3390/en16114462
Submission received: 16 March 2023 / Revised: 11 May 2023 / Accepted: 24 May 2023 / Published: 31 May 2023
(This article belongs to the Section J1: Heat and Mass Transfer)

Abstract

:
In this study, we numerically investigate heat transfer enhancement in a partially thermally active rectangular enclosure. The enclosure is filled with a ternary hybrid nanofluid (water, Carbon Nanotube, A l 2 O 3 , and Graphene). It is subjected to a magnetic field and uniform internal heat generation. The study also investigates the effect of magnetic field strength and direction on the natural convection flow, which arises from density fluctuations caused by partial heating of the left vertical wall. To solve the dimensionless governing equations, the finite element approach is employed. The parameters studied in detail include Rayleigh number ( R a ) , Hartmann number ( H a ) , nanoparticle volume fraction ( ϕ ) , and heat generation coefficient ( λ ) . The findings are presented graphically for the range of the parameters as follows: 10 3 R a   10 6 ,   0 H a 20 , 0.01 ϕ 0.05 , and 0 λ 15 . It is noted that these parameters have an impact on heat transfer enhancement, flow patterns, and temperature fields. The results show that the average Nusselt number ( N u ¯ ) increases with an increasing value of ϕ . Moreover, it has been noted that N u ¯ decreases as the value of H a increases, and the impact becomes more obvious at higher R a values. Finally, the influence of the heat generation coefficient on the heat transfer rate inside an enclosure is examined.

1. Introduction

Due to its diverse range of applications, including electronic appliances, microelectromechanical system appliances (MEMS), nuclear power plants, reactor cooling, buildings and thermally insulated systems, solar cells, the food storage industry, and geophysical fluid mechanics, the research of magnetohydrodynamic natural convection flows in an enclosure has lately acquired great attention [1]. Natural convection is influenced by magnetic fields in various practical instances, such as the formation of crystals in fluids, the casting of metals, the extraction of geothermal energy, and fusion reactors. Understanding the flow behavior and heat transport mechanism under the influence of magnetic field in an enclosure filled with electrically conducting fluid has become increasingly important [2,3,4,5]. When working with such fluids, two opposing forces, Lorentz and buoyancy forces (body forces), begin to interact with one another, thus altering the transport dynamics because of their mutual impacts.
External magnetic fields are increasingly used in the material manufacturing industry to control convection currents and improve product quality. Therefore, analysis of momentum and heat transport in such processes is essential. In the past, various researchers have investigated natural convection for different geometries and fluids subjected to magnetic fields [6]. For instance, Aydin and Yesiloz [7] investigated natural convection in a quadrantal enclosure filled with nanofluid with thermally active adjacent walls. Sathiyamoorthy and Chamkha [8] studied natural convection in a square cavity filled with liquid gallium and subjected to a magnetic field at an inclination angle ϕ to the horizontal plane. A numerical study of magneto-natural convection flow was conducted by Rudraiah et al. [9] in a rectangular enclosure with vertically isothermal walls and insulated horizontal walls.
Nanofluids have gained significance in heat transfer as the demand for heat transfer enhancement has grown. Choi and Eastman [10] proposed nanofluid, which several researchers are now considering. The use of nanofluids has been investigated in various engineering applications to enhance heat transfer performance. Researchers have extensively studied nanofluids with enhanced thermal properties to increase thermal efficiency in numerous engineering applications. Most researchers have found that incorporating high-thermal-conductivity nanoparticles into the base fluid enhances the thermal efficiency of the resulting nanofluid [11,12,13]. Some researchers believe that nanoparticles added to the base fluid can reduce heat transmission significantly. When it comes to enhanced enclosure design that is subjected to these factors, researchers should consider these factors. Thus, it is crucial to gain a better knowledge of the phenomenon.
A relatively new class of nanofluids (NFs), known as “hybrid” nanofluids, has recently appeared in parallel to the ongoing developments of the common nanofluids. This new class is produced by suspending many types of nanoparticles in a base fluid. A key objective of hybridization is compromised qualities between the advantages and disadvantages of the properties of specific nanoparticles. Additionally, the pricing of various nanoparticle varieties is substantially varied from one to another among providers of nanoparticles. For instance, copper nanoparticles cost around ten times more than alumina nanoparticles. Therefore, it is economical if one obtains the qualities of costly nanoparticles using a minimal amount. In light of this, scientists have developed modern-designed equipment by combining a mixture of three nanoparticles in the natural convection techniques. This has led to a significant improvement in thermal efficiency (heat transfer augmentation). Ternary hybrid nanofluid (THNF) is the term given to this new fluid. The creation of such a composite nanofluid is innovative since it significantly enhances thermophysical and heat transfer properties. Experimental studies show that these fluids’ thermal conductivities are higher than those of hybrid and mono nanofluids [14,15].
A numerical analysis of magnetohydrodynamic natural convection was performed by Raisi et al. [16] in a square enclosure containing nanofluid (water- A l 2 O 3 ). The study examined multiple parameters and demonstrated that heat transfer improved with increasing Rayleigh number ( R a ) but decreased with increasing Hartmann number ( H a ) . Kahveci and Oztuna [17] found that the heat transfer performance of an enclosure can be affected by the magnetic field’s direction. Their investigation indicated that the orientation of the magnetic field had a substantial impact on the heat transfer rate within the enclosure. Specifically, their findings showed that a magnetic field aligned in the x-direction was more effective at inhibiting heat transfer than a vertically oriented magnetic field. The numerical analysis carried out by P. Kandaswamy et al. [18] focused on magnetoconvection in a partially thermally active square cavity. The study found that low H a did not affect heat transfer rates at low R a . A numerical investigation by A. Mahmoudi et al. [19] examined the impact of a magnetic field and heat generation/absorption on heat transfer rates and entropy generation in a water- A l 2 O 3 nanofluid-filled square cavity using numerical methods. The study concluded that heat transfer rate was influenced by the heat production/absorption coefficient in the range of 10 < q < 5, but the entropy generation rate remained unaffected.
The present study aims to investigate flow behavior and heat transfer enhancement in a partially thermally active rectangular enclosure filled with a ternary hybrid nanofluid consisting of water (base fluid) and Carbon Nanotubes, A l 2 O 3 , and Graphene nanoparticles in natural convection flows. The study aims to analyze the influence of Rayleigh number ( R a ) , heat generation coefficient ( λ ) , and Hartmann number ( H a ) on flow patterns and temperature fields, as well as on the heat transfer rate within an enclosure. The magnetic field can alter the natural convection phenomena, which can be used in microgravity conditions as well as for suppressing the natural convection effect under gravity if it is needed.

2. Mathematical Formulation

2.1. Problem Description

Consider a steady, laminar flow of an incompressible ternary hybrid nanofluid in a two-dimensional rectangular enclosure exposed to a uniform magnetic field of strength B 0 at an inclination γ in the presence of uniform internal heat generation. The left vertical wall is partly subjected to constant high temperature θ h while the right vertical wall is partly subjected to cold temperature θ c as shown in Figure 1. The enclosure’s remaining boundaries are thermally insulated. Inside the enclosure, the radiation heat transfer effect is insignificant. Constant fluid properties are considered, and the variation in density that causes buoyant forces is approximated using the Boussinesq approximation.

2.2. Governing Equations

The continuity, change of linear momentum, and energy equations for steady, two-dimensional, laminar, incompressible flow in Cartesian coordinates in dimensional form are given below:
u x + v y = 0
where u and v denote the x and y components of velocity, respectively.
u u x + v u y = 1 ρ t h n f p x + μ t h n f ρ t h n f 2 u + 1 ρ t h n f σ t h n f B 0 2 ( v s i n γ c o s γ u s i n 2 γ )
u v x + v v y = 1 ρ t h n f p x + μ t h n f ρ t h n f 2 v + ( ρ β ) t h n f ρ t h n f g θ θ c + 1 ρ t h n f σ t h n f B 0 2 ( u s i n γ c o s γ v c o s 2 γ )
u θ x + v θ y = α t h n f 2 θ + Q ( ρ C p ) t h n f ( θ θ c )
where the fluid velocity and temperature are represented by v , T , and p being the pressure of fluid, respectively, where x and y -components of velocity are indicated by 𝑢 and 𝑣, respectively. ρ t h n f is the effective density and σ t h n f represents electrical conductivity, α t h n f being the thermal diffusivity of the ternary hybrid nanofluid and μ t h n f indicating the dynamic viscosity of the ternary hybrid nanofluid. Thermal expansion coefficient is represented by 𝛽 and magnetic field strength density as B 0 , respectively.
The dimensional boundary conditions for the governing equations are
u = 0 , v = 0 , θ = θ h   a t   x = 0 ,   a y b
u = 0 , v = 0 , θ = θ c   a t   x = 1.5 ,   c y d
u   =   0 ,   v   =   0 , θ x =   0   a t   y   =   0   ,   0     y     a ,   b     y   1
u   =   0 ,   v =   0 , θ x   =   0   a t   y   =   1.5   ,   0     y     c ,   d     y     1
u = 0 , v = 0 , θ y = 0   a t   y = 0   a n d   y = 1

2.3. Thermophysical Properties

Nanofluids demonstrate altered viscosity, density, specific heat capacity, and surface tension, among other properties. These modifications are highly dependent on factors like nanoparticle concentration, size, shape, and material. As a result, nanofluids offer immense potential in fields like electronics cooling, solar energy systems, automotive industries, and thermal management of electronic devices.
In this study, the thermophysical properties of nanofluid are calculated based on ref. [20]
ρ t h n f = ϕ 1 ρ p s 1 + ϕ 2 ρ p s 2 + ϕ 3 ρ p s 3 + 1 ϕ 1 ϕ 2 ϕ 3 ρ b f
β t h n f = ϕ 1 β s 1 + ϕ 2 β s 2 + ϕ 3 β s 3 + 1 ϕ 1 ϕ 2 ϕ 3 β b f
( ρ C p ) t h n f = ϕ 1 ρ C p s 1 + ϕ 2 ρ C p s 2 + ϕ 3 ρ C p s 3 + 1 ϕ 1 ϕ 2 ϕ 3 ρ C p b f
μ t h n f = ϕ 1 μ n f 1 + ϕ 2 μ n f 2 + ϕ 3 μ n f 3 ϕ ,   k t h n f = ϕ 1 k n f 1 + ϕ 2 k n f 2 + ϕ 3 k n f 3 ϕ α t h n f = k t h n f ( ρ C p ) t h n f
μ t h n f ,   k t h n f β t h n f ,   α t h n f , and ρ t h n f represent the dynamic viscosity, thermal conductivity, thermal expansion coefficient, thermal diffusivity, and density of the ternary hybrid nanofluid, respectively.
μ n f 1 = ( 1 + 2.5 ϕ 1 + ϕ 1 2 ) μ b f
μ n f 2 = ( 1 + 13.5 ϕ 2 + 904.4 ϕ 2 2 ) μ b f
μ n f 3 = ( 1 + 37.1 ϕ 3 + 612.6 ϕ 3 2 ) μ b f
k n f 1 = k b f k 1 + 2 k b f + 2 ϕ 1 k 1 k b f k 1 + 2 k b f ϕ 1 k 1 k b f
k n f 2 = k b f k 2 + 3.9 k b f + 3.9 k 2 k b f k 2 + 3.9 k b f ϕ 2 k 2 k b f
k n f 3 = k b f k 3 + 4.7 k b f + 4.7 k 3 k b f k 3 + 4.7 k b f ϕ 3 k 3 k b f
where 1, 2, and 3 represent A l 2 O 3 , CNT, and Graphene, respectively. μ n f i , k n f i where i = 1, 2, and 3 represent the viscosity and thermal conductivity of nanoparticles. ϕ i and C p i represent the volume fraction and specific heat of the nanoparticles, respectively.
For n and ψ , representing the shape factor and sphericity of nanoparticles, respectively, the relationship is as follows.
n = 3 ψ
ψ = 1 , n = 3 (Spherical-shaped nanoparticle)
ψ = 0.61 , n = 4.9 (Cylindrical-shaped nanoparticle)
ψ = 0.52 , n = 5.7 (Platelet-shaped nanoparticle)
The thermophysical properties of thermal conductivity and dynamic fluid viscosity of the ternary hybrid nanofluid also vary with particle shape. This means that the shape of particles (cylindrical, spherical, and platelet-shaped particles) can have a significant effect on these properties.

2.4. Non-Dimensional Form

The following are the dimensionless groups used in the governing equations, which are presented in their dimensionless form as:
X , Y = x , y L , U = u L α f , V = v L α f , P = p L 2 ρ n f α f 2 , T = θ θ C θ h θ c , R a = β f g L 3 θ h θ c ν f α f , P r = ν f α f , H a = B 0 L σ n f ρ n f ν f
U X + V Y = 0
U U X + V U Y = ρ f ρ t h n f P X + P r μ t h n f μ f ρ t h n f ρ f 2 U + H a 2 P r ( V s i n γ c o s γ U s i n 2 γ ) U
U V X + V V Y = ρ f ρ t h n f P Y + P r μ t h n f μ f ρ t h n f ρ f 2 V + R a P r β t h n f β f T + H a 2 Pr ( U s i n γ c o s γ V c o s 2 γ ) V
U T X + V T Y = α t h n f α f 2 T + λ T

2.5. Boundary Conditions

The dimensionless boundary conditions used to solve the governing equations are
    U = 0 , V = 0 ,                         T = T h   a t   X = 0 ,   a Y b
      U = 0 , V = 0 ,                           T = T c   a t   X = 1.5 ,   c Y d
                                                                    U   =   0 ,   V   =   0 ,             T X =   0   a t   X   =   0   ,   0     Y     a ,   b     Y     1
                                                                      U   =   0 ,   V   =   0 ,             T X   =   0   a t   X   =   1.5   ,   0     Y     c ,   d     Y     1
U = 0 , V = 0 ,                         T Y = 0   a t   Y = 0   a n d   Y = 1

2.6. Nusselt Number

Heat transfer enhancement within an enclosure is indicated by the dimensionless Nusselt number, with an increase in Nusselt number corresponding to thermal enhancement. The local and mean Nusselt number are given by the following expressions, respectively:
N u Y = k t h n f k b f T X X = 0
N u ¯ = 1 l 0 l N u ( Y ) d Y
where l represents the length of thermally active vertical walls.
The thermophysical properties for ternary hybrid nanofluid are influenced by the properties of water and nanoparticles, as outlined in Table 1.

2.7. Numerical Scheme and Code Validation

The Navier–Stokes and heat energy equations are mathematically modeled as partial differential equations using COMSOL Multiphysics v6.0, a commercial software that employs the finite element method. The computational domain is partitioned into several elements and a convergence criterion of ξ n + 1 ξ n ξ n 10 6 ¯ is set for all cases where ξ depicts the pressure, velocity, and temperature components. Mesh independence test is performed for the rectangular cavity with partially thermally sensitive vertical walls for R a = 10 5 as shown in Table 2. A grid sensitivity test is performed by taking into consideration five different elements, 2088, 3592, 8972, 22,644, and 37,430, and the grid density containing 37,430 elements is selected for all computations. The present scheme is first validated by verifying the published work of S. Dutta et al. [22]. The N u variation on the bottom heated wall for R a = 10 6 and ϕ = 0.05 (a) for numerical work [22] and (b) for the present work is shown in Figure 2. Table 3 presents the results of the validation of the method, where N u ¯ values for various H a values at a fixed ϕ = 0.05 and R a = 10 6 were compared to the findings of S. Dutta et al. [22]. The obtained values are in good agreement with that of S. Dutta et al. [22].

3. Results and Discussion

A numerical analysis was performed in this work and the results are displayed in terms of isotherms and streamlines for H a = 0 , 5 , 10 , 15   a n d   20 , A = 1.5 ,   ϕ = 0.03 , and 10 3 R a 10 6 . The flow of fluid and heat transfer enhancement within an enclosure depend on parameters such as Rayleigh number ( R a ), the strength of the magnetic field ( H a ) , and nanoparticle volume fraction ϕ .

3.1. Analysis of Flow Fields

The influence of R a and H a on flow behavior was investigated in terms of streamlines and is shown in Figure 3, Figure 4, Figure 5 and Figure 6. The Rayleigh number plays a crucial role in determining the flow pattern during natural convection. When there are temperature differences within a fluid, free convection can occur due to changes in density caused by these temperature variations. Specifically, when the temperature increases, the density of fluid decreases, leading to an upward motion of the fluid driven by buoyancy forces. For R a 10 4 , flow is not induced as the viscous forces predominate the buoyant and the Lorentz force, and the heat transfer is primarily through the conduction mechanism. It can be seen that the fluid near to the left vertical wall becomes heated and moves to the right vertical cold wall. The fluid descends after coming in contact with the cold vertical wall and then rises again along the left vertical heated wall, thus developing recirculation zones. Convection is considerably reduced near insulated boundaries. The recirculation cell becomes smaller with an increasing value of H a for R a 10 4 , as can be seen in Figure 3 and Figure 4. On the other hand, with increasing H a , the cell formed in the enclosure moves to the side walls and becomes bigger for higher values of R a 10 5 . Thus, the dominance of convective effects can be seen at higher values of R a which eventually increases the strength of vortices. It is worth noting that at lower H a , convective effects cannot be reduced by the magnetic field.

3.2. Analysis of Temperature Fields

A detailed investigation of the isotherms can help to understand the physics of heat transfer and thus is presented in Figure 3, Figure 4, Figure 5 and Figure 6. Isotherms are evenly distributed and nearly parallel to one another with a slight twisting, particularly for R a 10 4 , which shows the dominance of the conduction heat transfer mechanism in an enclosure at low R a 10 4 . An intriguing observation is made that for larger values of R a , the isotherms approach the cold vertical side wall with an increasing Hartmann number, which can be ascribed to stronger convective effects within an enclosure with the warm fluid moving to the heated wall and the corresponding low temperature fluid to the right colder wall. For R a = 10 6 , the isothermal lines are twisted for each H a , indicating the predominant convective effects within an enclosure.

3.3. Combined Effect of R a , H a , and ϕ on Heat Transfer Rate

Table 4 presents the impact of nanoparticle volume fraction ( ϕ ) , Hartmann number ( H a ) , and Rayleigh number ( R a ) on the average Nusselt number ( N u ¯ ) for the heated left vertical wall. The relationship between N u ¯ and H a is heavily influenced by R a . This can be observed from the variations in N u ¯ as R a changes. At low R a , heat transfer via conduction dominates, and the magnetic field has a minor impact. For R a 10 3 , the magnetic field has no significant effect on N u ¯ . For R a = 10 4 , the effect of increasing H a on N u ¯ is smaller, and this is true for all values of ϕ . However, at higher values of R a = 10 6 , there is a significant impact of magnetic field in reducing N u ¯ . When the Rayleigh number ( R a ) is sufficiently high, convection heat transfer becomes more dominant than the conduction mode of heat transfer, especially in the absence of a magnetic field. This is because the buoyancy force, which increases with R a , becomes stronger than viscous forces, leading to the transition from a conduction-dominated mechanism to a convection-dominated mechanism. Nonetheless, in the presence of a magnetic field within the enclosure, the drop in N u ¯ values is more significant for larger H a values at R a values of 10 5 and 10 6 compared to R a values of 10 3 and 10 4 . Therefore, the magnetic field plays a critical role in reducing N u ¯ values, and this reduction is clearly evident at larger R a values. Moreover, nanofluid presence causes a rise in N u ¯ with increasing ϕ , resulting in an increased heat transfer rate.

3.4. Influence of Magnetic Field Inclination Angle on Heat Transfer Rate

The effect of magnetic field inclination angle on average Nusselt for the given R a values, H a and ϕ , is shown in Figure 7. Fluid flow in a natural convection-dominated regime is caused by the difference in temperature between the thermally active right and left enclosure walls. When the applied magnetic field is 90° (or 270°) inclined, the magnetic field direction coincides with that of the high-temperature fluid (or low-temperature fluid).
This results in improved heat energy transportation, which also improves the efficiency of heat transfer. However, the magnetic effect prevents the high/low-temperature fluid from flowing when the applied magnetic field has an inclination angle of 0° (or 180°). Consequently, N u ¯ drops. The magnetic field’s vertical component facilitates fluid flow and heat energy transportation effects inside an enclosure for all other magnetic field orientation angles. Mean value of Nusselt number therefore exhibits a slight improvement. Hence, a magnetic field directed in a horizontal direction result in significant decrease in N u ¯ .

3.5. Effect of Hartmann Number on Heat Transfer Rate

The impact of H a on heat transfer efficiency within an enclosure filled with a ternary hybrid nanofluid, with a volume fraction of nanoparticles ϕ = 0.05 , is depicted in Figure 8. As R a rises, the heat transfer rate also increases. Conversely, an increase in H a from 0 to 20 leads to a decrease in heat transfer rate. The Hartmann number has a notable impact on heat transfer within an enclosure due to its influence on the strength and orientation of the magnetic field, which in turn affects fluid flow and heat transfer. At low H a , natural convection flow predominates, and the magnetic field is weak, while at high Hartmann numbers, the magnetic field dominates over viscous, resulting in significant suppression of fluid motion and decrease in convective heat transfer. In fact, the N u ¯ remains constant, as H a is increased for R a = 10 3 due to the predominance of the conduction heat transfer mechanism, and the magnetic field has no noticeable impact on heat transfer enhancement. When R a = 10 6 , convection heat transfer predominates and the N u ¯ decreases as H a increases.

3.6. Influence of Rayleigh Number on Heat Transfer Performance

The impact of R a on N u ¯ for distinct values of H a is shown in Figure 9. Clearly, a rise in R a corresponds to a rise in N u ¯ . R a is the crucial variable that governs how the flow is configured in natural convection. Free convection occurs when a fluid’s density changes owing to a temperature difference in the fluid. Buoyancy forces are what cause the fluid to rise when the temperature rises since the density of the fluid changes. Thus, it can be inferred that when R a rises, N u ¯ along the vertically heated section rises as a result of improved free convection heat transfer. The effect of H a on N u ¯ at R a = 10 3 , where heat transfer by conduction is dominant, is negligible. Contrarily, it may be observed that H a only significantly affects the N u ¯ at higher R a 10 6 values when convective flow is prevented by the magnetic field.

3.7. Influence of Hartmann Number on Velocity Profile

Figure 10 shows the impact of H a on vertical velocity component for different values of Rayleigh number. With an increase in H a , the velocity component’s amplitude initially increases and then decreases. This is because magnetic effects have a resistive nature, which reduces the velocity profile and, as a result, reduces the amount of heat transported by convection in the regime. In other words, buoyancy forces are negatively impacted by magnetic forces, which causes conduction to predominate the convection mechanism at higher values of H a . The magnetic field exerts Lorentz forces that resist the flow distribution, leading to a decrease in velocities and a suppression of convection as the H a increases. Increasing H a results in a decrease in velocity due to the impact of a strong magnetic field on fluid flow.

3.8. Influence of Rayleigh Number on the Velocity Profile

The effect of R a on vertical velocity component for different values of H a is shown in Figure 11. The velocity component is nearly linear for R a 10 4 . For R a = 10 5 there is a slight change in amplitude of the velocity curve. The amplitude of the velocity curve increases drastically for R a = 10 6 . In natural convection within an enclosed system, the Rayleigh number has a significant impact on the fluid’s velocity. With an increase in Rayleigh number, the buoyancy forces become more dominant than the viscous forces, resulting in faster fluid movement. For instance, when the left wall is heated and the right wall is cooled, the fluid near the hot surface expands and becomes less dense, causing it to rise due to buoyancy force. As it rises, it cools and becomes denser, resulting in it descending toward the cooler surface. This creates a convection current that propels the fluid motion within an enclosure. With an increase in the R a , the buoyancy forces become even stronger relative to the viscous forces, causing a more vigorous convection current, and therefore an increase in the fluid’s velocity.

3.9. Impact of Internal Heat Generation

The variation of N u ¯ along the heated wall is depicted in Figure 12, as a function of the heat generation coefficient ( λ ) and Rayleigh number ( R a ) for H a = 20 and ϕ = 0.05 . According to the results, the rate of heat transfer along the heated vertical wall reduces as the heat generation coefficient ( λ ) increases, regardless of the given R a values. When λ 0 , the temperature of nanofluids is more elevated than that of hot and cold vertical walls for R a = 10 3 and 10 4 , respectively. The nanofluids heat up both side walls as a result. This results in a negative Nusselt number at the heated wall. Figure 12b illustrates the influence of the heat production coefficient and the solid volume fraction on the heat transfer rate at the heated wall. With an increase in the heat generation parameter’s value, the N u ¯ decreases. For the high value of the heat-generation condition i.e., λ 0 , the addition of the nanoparticle has a significant effect. Figure 12c depicts the impact of heat generation or absorption coefficient on the N u ¯ of the hot wall for several H a . The figure shows that for any H a , the heat transfer rate at the heated wall decreases as the heat generation coefficient increases. The temperature in the cavity is higher under the heat-generating condition than it is under the non-heat generation condition ( λ = 0 ) . As a result, when the internal heat generation is present, the rate of heat transfer close to the heated walls of the cavity is lower than it would be otherwise. Thus, a decrease in H a is observed for increasing value of λ .
The impact of heat generation on temperature and flow fields, respectively, is presented in Figure 13. Initially, a single clockwise circulation cell that fills the whole cavity initially defines the flow field, and the isotherms exhibit a minor twisting toward the cold vertical wall, showing a dominance of conduction. This fact can be attributed to the external heating (temperature gradient). Internal heat generation eventually has an effect that is comparable to external heating. The formation of small secondary cells on the opposing sides of an enclosure, to the right and left of the heated and cold walls, is generated by increasing the heat generation parameter. This causes intense flow circulation inside an enclosure. An internal heating speeds up the enclosure’s rate of heat transmission in the temperature field. The temperature of the enclosure’s upper liquid layers tends to rise as a result of additional internal heat in the fluid. This alters the isothermal lines’ shape, and the thermal fields are concentrated towards an enclosure’s right half. The temperature of the nanofluid rises as the coefficient of heat generation increases, resulting in decreased heat transfer near the heated wall and promoting it towards the cold wall. The concentration of the isotherms near the cold wall for λ = 5 suggests a significant temperature gradient.

3.10. Comparison among Different Nanofluids

Figure 14 illustrates the variation of N u ¯ with solid volume fraction ϕ = 0.01, 0.03, and 0.05 at H a = 20 and R a = 10 3 for various nanoparticles ( A l 2 O 3 , CNT, Graphene) .  Table 1 lists the thermophysical characteristics of the mentioned nanoparticles. Figure 14 demonstrates that nano thermal conductivity is a key factor in heat transfer, and that water–Graphene nanofluid enhances the heat transport performance more than water– A l 2 O 3 and water–CNT nanofluid. Increasing the volume fraction of nanoparticles has a noticeable impact when natural convection is very low, or when the H a is higher, and the R a is lower.

4. Conclusions

This study uses numerical simulations to investigate natural convection in a rectangular enclosure filled with ternary hybrid nanofluid that is partially thermally active and subject to an inclined magnetic field with internal heat generation. The study demonstrates that the Rayleigh number ( R a ) , Hartmann number ( H a ) , and nanoparticle volume fraction ( ϕ ) affect the flow characteristics and heat transport rate in the enclosure. The following is a brief overview of the key findings.
The streamlines are concentrated on active sections of the vertical walls, and isotherms are distorted for higher R a values due to the dominant convection effects in the enclosure.
The rate of heat transport inside the enclosure is affected by flow parameters including R a , H a , and ϕ .
The N u ¯ decreases with increasing H a and at larger R a values. Rayleigh numbers play a significant role in controlling flow characteristics in natural convection due to stronger buoyancy effects.
At higher R a values, N u ¯ shows a decreasing trend with high ϕ due to the dominant Lorentz force.
The rate of heat transfer was found to be most significantly reduced when the magnetic field was directed horizontally.
Increasing the value of λ , the rate of heat transmission at the hot wall decreases for all R a values, while it increases close to the cold wall.

Author Contributions

Conceptualization, S.C.S.; Methodology, S.u.I., Z.Z. and S.A.; Software, M.S.; Validation, S.u.I.; Formal analysis, S.u.I.; Investigation, Z.Z., M.S. and S.A.; Resources, M.S.; Writing—original draft, S.u.I. and Z.Z.; Writing—review & editing, S.C.S., M.S. and S.A.; Supervision, S.C.S. and S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A aspect ratio ( W / H ) Greek symbols
C P Specific heat at constant pressure ( J K g 1 K 1 )
B 0 magnetic field α thermal diffusivity ( m 2 s 1 )
g acceleration due to gravity ( m s 2 ) β coefficient of thermal expansion ( K 1 )
R a Rayleigh Number μ dynamic viscosity ( N s m 2 )
H a Hartmann number ν kinematic viscosity ( m 2 s 1 )
N u local Nusselt number ρ density (kg m 3 )
N u ¯ average Nusselt number σ electrical conductivity
p pressure (N m 2 ) ϕ nanoparticle volume fraction
P r Prandtl number γ inclination angle
k thermal conductivity (W m 1 K 1 ) θ temperature (°C)
T dimensionless temperature λ heat generation coefficient
u , v velocity components (m s 1 ) S u b s c r i p t s
U , V dimensionless velocity components c cold
x , y dimensional coordinates (m)hhot
X , Y dimensionless coordinates t h n f ternary hybrid nanofluid
t dimensional time (s) b f base fluid

References

  1. Ostrach, S. Natural convection in enclosures. J. Heat Transf. 1988, 110, 1175–1190. [Google Scholar] [CrossRef]
  2. Ozoe, H.; Okada, K. The effect of the direction of the external magnetic field on the three-dimensional natural convection in a cubical enclosure. Int. J. Heat Mass Transf. 1989, 32, 1939–1954. [Google Scholar] [CrossRef]
  3. Garandet, J.P.; Alboussiere, T.; Moreau, R. Buoyancy driven convection in a rectangular enclosure with a transverse magnetic field. Int. J. Heat Mass Transf. 1992, 35, 741–748. [Google Scholar] [CrossRef]
  4. Venkatachalappa, M.; Subbaraya, C.K. Natural convection in a rectangular enclosure in the presence of a magnetic field with uniform heat flux from the side walls. Acta Mech. 1993, 96, 13–26. [Google Scholar] [CrossRef]
  5. Alchaar, S.; Vasseur, P.; Bilgen, E. Natural Convection Heat Transfer in a Rectangular Enclosure with a Transverse Magnetic Field. J. Heat Transf. 1995, 117, 668–673. [Google Scholar] [CrossRef]
  6. Ece, M.C.; Büyük, E. Natural-convection flow under a magnetic field in an inclined rectangular enclosure heated and cooled on adjacent walls. Fluid Dyn. Res. 2006, 38, 564–590. [Google Scholar] [CrossRef]
  7. Aydin, O.; Yesiloz, G. Natural Convection in a Quadrantal Cavity Heated and Cooled on Adjacent Walls. J. Heat Transf. 2011, 133, 052501. [Google Scholar] [CrossRef]
  8. Sathiyamoorthy, M.; Chamkha, A. Effect of magnetic field on natural convection flow in a liquid gallium filled square cavi-ty for linearly heated side wall(s). Int. J. Therm. Sci. 2010, 49, 1856–1865. [Google Scholar] [CrossRef]
  9. Rudraiah, N.; Barron, R.M.; Venkatachalappa, M.; Subbaraya, C.K. Effect of a magnetic field on free convection in a rectangu-lar enclosure. Int. J. Eng. Soc. 1995, 33, 1075–1084. [Google Scholar] [CrossRef]
  10. Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Argonne National Lab: Argonne, IL, USA, 1995.
  11. Nnanna, A.G.A. Experimental Model of Temperature-Driven Nanofluid. J. Heat Transf. 2006, 129, 697–704. [Google Scholar] [CrossRef]
  12. Oztop, H.F.; Abu-Nada, E. Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids. Int. J. Heat Fluid Flow 2008, 29, 1326–1336. [Google Scholar] [CrossRef]
  13. Aminossadati, S.M.; Ghasemi, B. Natural convection cooling of a localized heat source at the bottom of a nanofluid-filled enclosure. Eur. J. Mech.-B Fluids 2009, 28, 630–640. [Google Scholar] [CrossRef]
  14. Boroomandpour, A.; Toghraie, D.; Hashemian, M. A comprehensive experimental investigation of thermal con-ductivity of a ternary hybrid nanofluid containing MWCNTs-titania-zinc oxide/water-ethylene glycol (80:20) as well as bi-nary and mono nanofluids. Synth. Met. 2020, 268, 116501. [Google Scholar] [CrossRef]
  15. Rizvi, S.M.M.; El Far, B.; Shin, D. Heat capacity and viscosity of ternary carbonate nanofluids. Int. J. Energy Res. 2021, 45, 6350–6359. [Google Scholar] [CrossRef]
  16. Raisi, A.; Ghasemi, B.; Aminossadati, S.M. Magnetic field effect on natural convection in a nanofluid-filled square enclosure. Int. J. Therm. Sci. 2011, 50, 1748–1756. [Google Scholar] [CrossRef]
  17. Kahveci, K.; Öztuna, S. MHD natural convection flow and heat transfer in a laterally heated partitioned enclosure. Eur. J. Mech. E B Fluids 2009, 28, 744–752. [Google Scholar] [CrossRef]
  18. Kandaswamy, P.; Sundari, S.M.; Nithyadevi, N. Magnetoconvection in an enclosure with partially active vertical walls. Int. J. Heat Mass Transf. 2008, 51, 1946–1954. [Google Scholar] [CrossRef]
  19. Mahmoudi, A.; Mejri, I.; Abbassi, M.A.; Omri, A. Analysis of the entropy generation in a nanofluid-filled cavity in the presence of magnetic field and uniform heat generation/absorption. J. Mol. Liq. 2014, 198, 63–77. [Google Scholar] [CrossRef]
  20. Shao, Y.; Nayak, M.K.; Dogonchi, A.S.; Chamkha, A.J.; Elmasry, Y.; Galal, A.M. Ternary hybrid nanofluid natural convection within a porous prismatic enclosure with two movable hot baffles: An approach to effective cooling. Case Stud. Therm. Eng. 2022, 40, 102507. [Google Scholar] [CrossRef]
  21. Sahoo, R.R. Thermo-hydraulic characteristics of radiator with various shape nanoparticle-based ternary hybrid nanofluid. Powder Technol. 2020, 370, 19–28. [Google Scholar] [CrossRef]
  22. Dutta, S.; Pati, S.; Baranyi, L. Numerical analysis of magnetohydrodynamic natural convection in a nanofluid filled quadrantal enclosure. Case Stud. Therm. Eng. 2021, 28, 101507. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the problem.
Figure 1. Schematic diagram of the problem.
Energies 16 04462 g001
Figure 2. N u variation on bottom heated wall for different values of H a where R a = 10 6 and ϕ = 0.05 . Comparison with S. Dutta et al. [22]: (a) numerical work by S. Dutta et al. [22]; (b) present work for R a = 10 6 and ϕ = 0.05 .
Figure 2. N u variation on bottom heated wall for different values of H a where R a = 10 6 and ϕ = 0.05 . Comparison with S. Dutta et al. [22]: (a) numerical work by S. Dutta et al. [22]; (b) present work for R a = 10 6 and ϕ = 0.05 .
Energies 16 04462 g002
Figure 3. (a) denotes isotherms and (b) denotes streamlines in a rectangular enclosure corresponding to R a = 10 3 ,   A = 1.5 ,   λ = 0 , and ϕ i = 0.03 where i = 1, 2, 3 for varying H a = 0 , 5 , 10 , 15 , and 20.
Figure 3. (a) denotes isotherms and (b) denotes streamlines in a rectangular enclosure corresponding to R a = 10 3 ,   A = 1.5 ,   λ = 0 , and ϕ i = 0.03 where i = 1, 2, 3 for varying H a = 0 , 5 , 10 , 15 , and 20.
Energies 16 04462 g003
Figure 4. (a) denotes isotherms and (b) denotes streamlines in rectangular enclosure corresponding to R a = 10 4 ,   A = 1.5 ,   λ = 0 , and ϕ i = 0.03 where i = 1, 2, 3 for varying H a = 0 , 5 , 10 , 15 , and 20 .
Figure 4. (a) denotes isotherms and (b) denotes streamlines in rectangular enclosure corresponding to R a = 10 4 ,   A = 1.5 ,   λ = 0 , and ϕ i = 0.03 where i = 1, 2, 3 for varying H a = 0 , 5 , 10 , 15 , and 20 .
Energies 16 04462 g004
Figure 5. (a) denotes isotherms and (b) denotes streamlines in rectangular enclosure corresponding to R a = 10 5 ,   A = 1.5 ,   λ = 0 , and ϕ i = 0.03 where i = 1, 2, 3 for varying H a = 0 , 5 , 10 , 15 , and 20 .
Figure 5. (a) denotes isotherms and (b) denotes streamlines in rectangular enclosure corresponding to R a = 10 5 ,   A = 1.5 ,   λ = 0 , and ϕ i = 0.03 where i = 1, 2, 3 for varying H a = 0 , 5 , 10 , 15 , and 20 .
Energies 16 04462 g005
Figure 6. (a) denotes streamlines and (b) denotes isotherms in a rectangular enclosure corresponding to R a = 10 6 ,   A = 1.5 ,   λ = 0 , and ϕ i = 0.03 where i = 1, 2, 3 for varying H a = 0 , 5 , 10 , 15 , and 20.
Figure 6. (a) denotes streamlines and (b) denotes isotherms in a rectangular enclosure corresponding to R a = 10 6 ,   A = 1.5 ,   λ = 0 , and ϕ i = 0.03 where i = 1, 2, 3 for varying H a = 0 , 5 , 10 , 15 , and 20.
Energies 16 04462 g006
Figure 7. Variation of N u ¯ with H a for the values of γ labelled on graph where ϕ = 0.05 and R a = 10 5 .
Figure 7. Variation of N u ¯ with H a for the values of γ labelled on graph where ϕ = 0.05 and R a = 10 5 .
Energies 16 04462 g007
Figure 8. Variation of N u ¯ with H a for the values of R a labelled on the graph where ϕ = 0.05 , P r = 6.2 , and λ = 0 .
Figure 8. Variation of N u ¯ with H a for the values of R a labelled on the graph where ϕ = 0.05 , P r = 6.2 , and λ = 0 .
Energies 16 04462 g008
Figure 9. Variation of N u ¯   a l o n g   h o t   w a l l with R a for the values of H a labelled on the graph where ϕ = 0.05 , P r = 6.2 , and λ = 0 .
Figure 9. Variation of N u ¯   a l o n g   h o t   w a l l with R a for the values of H a labelled on the graph where ϕ = 0.05 , P r = 6.2 , and λ = 0 .
Energies 16 04462 g009
Figure 10. Velocity profiles: (a) Vertical velocity profile along horizontal centerline for the values of H a labelled on a graph where R a = 10 3 , ϕ = 0.05 , P r = 6.2 , and length of heated wall = 0.5; (b) R a = 10 4 ; (c) R a = 10 5 ; (d) R a = 10 6 .
Figure 10. Velocity profiles: (a) Vertical velocity profile along horizontal centerline for the values of H a labelled on a graph where R a = 10 3 , ϕ = 0.05 , P r = 6.2 , and length of heated wall = 0.5; (b) R a = 10 4 ; (c) R a = 10 5 ; (d) R a = 10 6 .
Energies 16 04462 g010aEnergies 16 04462 g010b
Figure 11. (a) Vertical velocity profile along horizontal centerline for the values of R a labelled on the graph, where P r = 6.2 and (a) H a = 0 ; (b) H a = 5 ; (c) H a = 10 ; (d) H a = 15 ; (e) H a = 20 .
Figure 11. (a) Vertical velocity profile along horizontal centerline for the values of R a labelled on the graph, where P r = 6.2 and (a) H a = 0 ; (b) H a = 5 ; (c) H a = 10 ; (d) H a = 15 ; (e) H a = 20 .
Energies 16 04462 g011aEnergies 16 04462 g011bEnergies 16 04462 g011c
Figure 12. N u ¯ variation for various λ and (a) R a at H a = 20 and ϕ = 0.01 ; (b) ϕ at H a = 20 and R a = 10 5 ; (c) H a at ϕ = 0.01 and R a = 10 5 .
Figure 12. N u ¯ variation for various λ and (a) R a at H a = 20 and ϕ = 0.01 ; (b) ϕ at H a = 20 and R a = 10 5 ; (c) H a at ϕ = 0.01 and R a = 10 5 .
Energies 16 04462 g012
Figure 13. (i) Streamlines and (ii) isotherms in a rectangular enclosure corresponding to (a) λ = 0 , (b) λ = 5 , (c) λ = 10 , and (d) λ = 15 .
Figure 13. (i) Streamlines and (ii) isotherms in a rectangular enclosure corresponding to (a) λ = 0 , (b) λ = 5 , (c) λ = 10 , and (d) λ = 15 .
Energies 16 04462 g013
Figure 14. N u ¯ variation with ϕ for different nanoparticles where H a = 20 and R a = 10 3 .
Figure 14. N u ¯ variation with ϕ for different nanoparticles where H a = 20 and R a = 10 3 .
Energies 16 04462 g014
Table 1. Numerical values of thermo-physical properties of base fluid (Water) and nanoparticles (CNT, A l 2 O 3 , and Graphene) [21].
Table 1. Numerical values of thermo-physical properties of base fluid (Water) and nanoparticles (CNT, A l 2 O 3 , and Graphene) [21].
PropertyWaterAl2O3GrapheneCNT
k (W/(m K))0.6134050003007.4
ρ ( k g / m 3 ) 997.1397022002100
β ( K 1 ) 2.1 × 10−45.80 × 10−6−0.0000080.00002
Cp (J/(kg K))4179765790410
μ (kg/(ms))8.91 × 10−5-
Shape-sphericalplateletcylindrical
Table 2. Mesh dependency for rectangular enclosure, R a = 10 5 .
Table 2. Mesh dependency for rectangular enclosure, R a = 10 5 .
RL#Elements #DOFAvg Nu% Error
Normal208850442.37590.88
Fine359282682.39722.20
Finer897220,2882.45101.45
Extra Fine22,64449,7922.48720.004
Extremely Fine37,43079,3642.4871-
Table 3. N u ¯ for several H a where ϕ = 0.05 and R a = 10 4 .
Table 3. N u ¯ for several H a where ϕ = 0.05 and R a = 10 4 .
Ha     N u ¯
S. Dutta et al. [22]Present WorkRelative Difference %
04.9454.90330.84
204.2734.17492.30
404.0293.97261.40
604.0003.95441.14
803.9953.95101.10
1003.9933.95011.07
Table 4. N u ¯ along heated vertical wall for varying Ra, Ha and ϕ .
Table 4. N u ¯ along heated vertical wall for varying Ra, Ha and ϕ .
N u ¯
ϕ H a R a = 10 3 R a = 10 4 R a = 10 5 R a = 10 6
0.0100.644291.39102.95755.5226
50.631171.33792.91045.4951
100.608551.20982.78235.4154
150.593131.06202.60225.2904
200.585210.929652.39895.1294
0.0300.655371.21822.72735.1881
50.651171.18232.68875.1642
100.643161.09302.58235.0947
150.636860.985832.43034.9853
200.633120.887582.25574.8439
0.0500.692151.06782.48714.8571
50.690971.04562.45624.8364
100.688460.989142.37024.7760
150.686150.919602.24504.6807
200.684560.854742.09864.5570
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Saha, S.C.; Islam, S.u.; Zia, Z.; Saleem, M.; Ahmad, S. Thermal Analysis of Magneto-Natural Convection Flows within a Partially Thermally Active Rectangular Enclosure. Energies 2023, 16, 4462. https://doi.org/10.3390/en16114462

AMA Style

Saha SC, Islam Su, Zia Z, Saleem M, Ahmad S. Thermal Analysis of Magneto-Natural Convection Flows within a Partially Thermally Active Rectangular Enclosure. Energies. 2023; 16(11):4462. https://doi.org/10.3390/en16114462

Chicago/Turabian Style

Saha, Suvash C., Shams ul Islam, Zahida Zia, M. Saleem, and Shafee Ahmad. 2023. "Thermal Analysis of Magneto-Natural Convection Flows within a Partially Thermally Active Rectangular Enclosure" Energies 16, no. 11: 4462. https://doi.org/10.3390/en16114462

APA Style

Saha, S. C., Islam, S. u., Zia, Z., Saleem, M., & Ahmad, S. (2023). Thermal Analysis of Magneto-Natural Convection Flows within a Partially Thermally Active Rectangular Enclosure. Energies, 16(11), 4462. https://doi.org/10.3390/en16114462

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop