Next Article in Journal
A Novel Constraint-Aware Flexible Model with Ant Colony Optimization for Symmetrical Travel Recommendation
Next Article in Special Issue
Effect of an Adiabatic Obstacle on the Symmetry of the Temperature, Flow, and Electric Charge Fields during Electrohydrodynamic Natural Convection
Previous Article in Journal
Relation-Preserving Functional Contractions Involving a Triplet of Auxiliary Functions with an Application to Integral Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Control of Three-Dimensional Natural Convection of Graphene–Water Nanofluids Using Symmetrical Tree-Shaped Obstacle and External Magnetic Field

1
Department of Mechanical Engineering, College of Engineering, University of Ha’il, Ha’il City 81451, Saudi Arabia
2
Laboratory of Metrology and Energy systems, Department of Energy Engineering, University of Monastir, Monastir 5000, Tunisia
3
Department of Computer Science, Applied College, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11432, Saudi Arabia
4
Department of Mechanical Engineering, College of Engineering, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
5
Department of Electrical Engineering, University of Ha’il, Ha’il City 81451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Submission received: 14 April 2024 / Revised: 11 May 2024 / Accepted: 23 May 2024 / Published: 4 June 2024
(This article belongs to the Special Issue Symmetry in Thermal Fluid Sciences and Energy Applications)

Abstract

:
This numerical investigation explores the enhanced control of the 3D natural convection (NC) within a cubic cavity filled with graphene–water nanofluids, utilizing a bottom-center-located tree-shaped obstacle and a horizontal magnetic field (MF). The analysis includes the effects of the Rayleigh number (Ra), the solid volume fraction of graphene ( φ ), the Hartmann number (Ha), and the fins’ length (W). The results show complex flow patterns and thermal behavior within the cavity, indicating the interactive effects of nanofluid properties, the tree-shaped obstacle, and magnetic field effects. The MHD effects reduce the convection, while the addition of graphene improves the thermal conductivity of the fluid, which enhances the heat transfer observed with increasing Rayleigh numbers. The increase in the fins’ length on the heat transfer efficiency is found to be slightly negative, which is attributed to the complex interplay between the enhanced heat transfer surface area and fluid flow disruption. This study presents an original combination of non-destructive methods (magnetic field) and a destructive method (tree-shaped obstacle) for the control of the fluid flow and heat transfer characteristics in a 3D cavity filled with graphene–water nanofluids. In addition, it provides valuable information for optimizing heat transfer control strategies, with applications in electronic cooling, renewable energy systems, and advanced thermal management solutions. The application of a magnetic field was found to reduce the maximum velocity and total entropy generation by about 82% and 76%, respectively. The addition of graphene nanoparticles was found to reduce the maximum velocity by about 5.5% without the magnetic field and to increase it by 1.12% for Ha = 100. Varying the obstacles’ length from W = 0.2 to W = 0.8 led to a reduction in velocity by about 23.6%.

1. Introduction

The heat transfer and fluid flow destructive and non-destructive control techniques have drawn the attention of researchers and engineers with the principal goal to optimize the energy systems. These techniques involve mixing nanoparticles (NPs) with conventional base fluids, the use of fins and obstacles, and the application of magnetic fields. Liu et al. [1] explored the impacts of NP concentration, MF intensity and direction, and cylindrical grooves on nanofluid convective flow 2D cavities. The results show that stepped grooves and specific magnetic configurations considerably improve heat transfer. Tayebi et al. [2] performed a study on the convective nanofluid flow and irreversibility production in an annular enclosure under the application of a MF and using fins having different lengths. They examined the heat transfer and entropy production, with an economic analysis on the use of the nanofluid, to balance between performance improvement and cost. Saha et al. [3] conducted the study of the MHD convective flow of an Al2O3-H2O nanoliquid within an enclosure having a corrugated top wall and a hot fin located at the bottom wall. It was revealed that blade-shaped nanoparticles improve the thermal transport by 7.65%, compared to 2.86% for spherical nanoparticles. Iftikhar et al. [4] explored the influences of hot fins and radiation on the NC flow of Ag and Cu nanoparticles dispersed in ethylene glycol. The study demonstrates an improvement of the buoyancy-driven flow with an increase in fin dimensions and NP concentration; the enhancement reached a 42% increase in fluid velocity and 57% of the heat transfer rate. Nayak et al. [5] studied the unsteady nanofluid MHD–NC flow in a hexagonal cavity incorporating a cold diamond-shaped obstacle. The findings show an important improvement in heat transfer efficiency: a 76.16% improvement in the Nuav as the Ra increases, and the optimal obstacle position significantly enhances the heat transfer. Islam et al. [6] investigated the hybrid nanoliquid MHD NC in a prismatic cavity. The study highlighted an improvement of the thermal performance with an increase in the Ra and nanoparticle concentrations, while magnetic influence shows opposite effects. In the investigation of Vinodhini and Prasad [7], the MHD NC of nanofluids in a differentially heated 2D cavity was investigated. The heat generation/absorption and an oblique magnetic field were considered, with a focus on the effects of Brownian motion and thermophoresis. Al-Farhany et al. [8] presented a study focusing on the effect of fin lengths and location on the nanofluid MHD free convective flow in a porous cavity. The findings showed an improved heat transfer with increasing fin length, fin spacing, and modified Rayleigh number, while an increase in the Ha decreases the heat transfer. Cheng et al. [9] analyzed water/copper nanofluid convective flow in a confined cavity. The study finds that the magnetic forces reduce the velocity and temperature gradients, decreasing entropy generation, while varying the magnetic field angle and baffle dimensions significantly impacts the generation of entropy. Abderrahmane et al. [10] numerically explored the MHD 3D free convection in an undulated cubic cavity. The results indicate the sensitivity of entropy generation to Darcy and Rayleigh numbers, the Hartmann number, and the undulation number. Acharya et al. [11] performed a study on the MHD flow of a hybrid nanofluid flow in a discreetly heated, fin-equipped circular cavity. The study reveals that increasing Rayleigh numbers accelerate velocity, while magnetic fields and higher concentrations of nanoparticles tend to slow the flow. Waqas et al. [12] studied the impact of the fin number, NP fraction, and Ra on laminar and steady NC flow in a horizontal ring. It was found that the insertion of fins and the NPs improve the heat transfer efficiency, especially at higher Ra values. Armaghani et al. [13] considered a two-phase model to investigate the Al2O3-nanofluid convective flow in an E-shaped cavity. A novel correlation was used to get accurate values of the thermal conductivity and it was found that the entropy generation was reduced by adding nanoparticles. Khan et al. [14] conducted a computational analysis on the 2D nanoliquid NC in a closed cavity having partially active walls. They demonstrated that the addition of nanoparticles and the increase in the Ra enhance velocity and heat transfer performance. The ternary nanofluid convective in a tilted porous cavity under various conditions was investigated by Thirumalaisamy et al. [15]. It was found that the specific combination of the different kinds of nanoparticles can lead to an optimized heat transfer. The conjugated MHD nanofluid free convection and irreversibility production considering heat-generating blocks was studied by Tasnim et al. [16]. The authors mentioned that higher heat transfer occurs with pure water than with nanofluids and that the position of the heat-generating blocks and the enclosure’s tilt angle play a significant role on the thermal performance. The MHD–NC flow of a micropolar nanoliquid in a semi-annular cavity has been considered by Seyyedi et al. [17]. The heat transfer efficiency was found to decline with the increase in Ha, and it was improved at a higher Ra and by adding nanoparticles. The effect of using a porous fin nanofluid natural convective flow was examined by Siavashi et al. [18]. The heat transfer was improved at high Darcy number values, while lower values reduced it. Additionally, low concentrations of nanoparticles enhanced the heat transfer more effectively compared to higher concentrations. Çiçek et al. [19] numerically analyzed a hybrid nanofluid natural convective flow, considering the migration and deposition nanoparticles phenomena. They found that adding nanoparticles and increasing Ra led to an improved heat transfer. Acharya [20] studied the hydrothermal effects of laminar NC of a hybrid nanofluid in a circular cavity equipped with I-shaped fins. The study reveals that fin length significantly influences thermal control and flow behavior, and variations in Ra, Ha, and nanoparticle concentration affect the temperature field and flow dynamics. Ali et al. [21] explored the MHD NC of a non-Newtonian nanofluid in a U-shaped enclosure, focusing on the impacts of Ra, Ha, NP volume fraction, ratio appearance of the deflectors, the angle of inclination, and the power index. They noted a notable impact of the power index on heat transfer at higher Ra. Qi et al. [22] studied the effects of rotation angles, metal foam, and adjustable magnetic fields on the nanoliquids convective heat transfer. Their study reveals that the optimal nanoparticle mass fraction for maximum heat transfer is 0.3% and that vertical magnetic fields improve thermal performance, particularly at a 90° cavity rotation angle. Liao and Li [23] performed a study on the heat transfer performance of nanoliquids under an inclined MF. The authors mentioned that there are three distinct heat transfer structures based on the magnetic field magnitude and analyzed the critical Hartmann numbers, allowing the transitions between these structures. In addition, the authors proposed an empirical expression, allowing the prediction of the critical Ha based on the magnetic field orientation. Vahedi et al. [24] presented a study on the thermal performance of a hybrid nanofluid under MHD effects. Based on sensitivity analysis and optimization, the authors mentioned that although the increase in the number of fins led to an enhancement of the heat transfer by about 30%, longer fins led to similar results even with fewer fins, reducing the viscous entropy generation. The hybrid nanoliquid MHD hydrothermal behavior and entropy production in an octagonal enclosure has been explored by Acharya [25], with a special focus on the effect of fin lengths. The authors concluded that the fin length and the addition of nanoparticles have a significant influence on the irreversibilities and heat transfer. Cao et al. [26] discussed the 2D MHD NC of a nanoliquid. It was mentioned that the magnetic fields and NPs have an important effect on the heat transfer, the NPs transport area, and on the temperature field. Selimefendigil et al. [27] evaluated the influences of a non-uniform MF on the mixed convection in a 3D vented enclosure filled with encapsulated PCM through which a nanoliquid flows. They concluded that the phase change is accelerated with higher values of Re and concentration of nanoparticles. Alazzam et al. [28] explored the heat transfer dynamics of nano-encapsulated PCMs in a cubic corrugated enclosure under a MF. They found that the Ra and Ha significantly influence heat transfer, where a higher Ra is improving and a higher Ha is decreasing the average Nusselt number. Additional investigations into the influence of combining nanofluids with magnetic fields and fins in complex cavity designs on heat transfer characteristics can be found in the references [29,30].
Based on the above-mentioned literature survey, it is evident that there is a lack of investigations on the three-dimensional aspect of nanofluids’ flows in confined spaces. In fact, most of the studies dealing with the effects of nanofluids and magnetic fields focused on two-dimensional configurations. Researchers generally avoid the consideration of 3D MHD flows due to the complexity of the governing equation and the need to resolve the current density and electrical potential equation. In addition, two-dimensional studies often do not give accurate results, compared to the real convective phenomena; this is mainly due to the neglection of the transversal flow. Thus, more research is necessary to explore the more in-depth MHD convective flows, especially in the presence of internal obstacles. Our research aims to address this gap through a computational study focused on the enhanced control of the 3D NC in a cubic cavity filled with graphene–water nanofluid by means of a uniform MF and a tree-shaped obstacle, which can potentially contribute to the development of more effective and energy-efficient cooling solutions, ensuring the reliable and sustained operation of advanced electronic devices under high thermal loads.

2. Problem Statement and Mathematical Modeling

The schematic view depicted in Figure 1 involves a cubic cavity filled with graphene–water nanofluid, where natural convection is controlled using an external magnetic field and a tree-shaped obstacle positioned at the bottom center, with the associate coordinate system. The left and right walls are considered as hot(Th) and cold(Tc), respectively, while all the remaining lateral walls are thermally insulated.
Based on the aforementioned assumptions and the Boussinesq approximation for the buoyancy force, the governing equations are presented as follows [31]:
· V = 0
V t + ( V · ) V = 1 ρ n f P + 1 ρ n f ( J × B ) + v V β t   ( T T c )   g
T t + V · T = α n f T
J = σ n f   ( ф + V × B )
· J = 0
The selection of dimensionless variables is as follows:
T = t L 2 α ; (x,y,z) = ( x , y , z ) L ; ( V x , V y ,   V z ) = ( V x , V y , V z ) α L ; P = P ρ ( α L ) ; J = J σ · υ 0 · B 0 ; B = B B 0 ; ф = ф l · υ 0 · B 0 and T = (T′ − T c )/( T h T c ).
Therefore, the governing equations are as follows:
· V = 0
V t + ( V · ) V = P + σ n f ρ f σ n ρ n f   Ha 2   Pr   ( J × B ) + Pr   V n f V f V + β n f β f   Ra   Pr   T g
T t + V · T = α n f α f 2 T
j = ф + V × e B
2 ф = · ( V × B )
with:
Pr = υ f α f ;   Ra = g · β f · T · l 3 υ f · α f   and   Ha = B 0 L σ ρ f V f
Boundary conditions are defined as:
T = 1 ,   a t   x = 0
T = 0 ,   a t   x = 1
T n = 0   a t   y = 0 , z = 0   a n d   z = 1
V x = V y = V z = 0 ,   a t   a l l   w a l l s
ф n = 0   a t   a l l   w a l l s
J · n = 0   o n   a l l   w a l l s
Taking into account the magnetohydrodynamic effects, the generation of entropy is described as follows [31]:
S g e n = k n f T 0 2 T x 2 + T y 2 + T z 2 + µ n f T 0 2 V x x 2 + V y y 2 + V z z 2 + V y x + V x y 2 + V z y + V y z 2 + V x z + V z x 2 + 1 T 0   1 σ n f   ( J x 2 + J y 2 + J z 2 )
Based on the previously mentioned dimensionless variable, the locally generated entropy can be expressed as follows:
N s = k n f k f T x 2 + T y 2 + T z 2 + φ s μ n f μ f 2 V x x 2 + V y y 2 + V z z 2 + V y x + V x y 2 + V z y + V y z 2 + V x z + V z x 2 + I s μ f μ n f σ n f σ f   Ha 2   J x 2 + J y 2 + J z 2
I s = α f l T 2   T 0 is the irreversibility coefficient.
I s is fixed at 10−4 in this study.
The total generated entropy is:
S t o t = v N S d v = v N S , t h + N S , f r + N S , m a g d v = S t h + S f r + S m a g
The local and average Nusselt numbers are expressed as:
N u l o c = k n f k f T x x = 0
and
N u a v = 0 1 0 1 N u l o c · d y · d z
The density, specific heat, and electrical conductivity of the nanofluid are evaluated using the following expressions:
ρ n f = ( 1 φ )   ρ f + ф   ρ n p
c p , n f = 1 φ   ρ f   c p , f + φ   ( ρ n p   c p , n p )   ρ n f
σ n f = σ f 1 + 3 · φ · σ s σ f 1 σ n p σ f + 2 φ · σ n p σ f 1
It should be noted that achieving the uniformity of NP dispersion in the base fluid is essential to enhance the thermal conductivity and heat transfer. But practically, this uniform dispersion can hardly be achieved. Sedimentation, aggregation, clogging, etc., can become the common problems that would deteriorate the efficiency of nanofluid-based systems. Furthermore, it is often assumed that the properties of NPs and the base fluid are constant. But their properties depend on the conditions of the system, such as temperature, pressure, NP concentration, etc. Hence, the performance of the system can be impacted in practice. The properties of water and graphene are presented in Table 1.
In the present study, the thermal conductivity and dynamic viscosity of the nanofluid are measured using a KD2-Pro thermal properties analyzer and a Malvern Kinexus pro rheometer, respectively. More details on the measurement procedure can be found in Ref. [32]. Based on these measurements, the following expressions are established. The power law regression model has been used for the thermal conductivity and the enhanced power law regression model, with the added interaction and polynomial terms, for the dynamic viscosity.
k   ( W m · K ) = 0.509 · T 0.013 · φ 0.311 ( R 2 = 0.9666 ) 0 % φ 2 %   and   20   C T 60   C
µ   ( c P ) = 0.00622 · e 3.67 · l o g T + 0.52 · l o g φ 0.59 · l o g T 2 + 0.088 · l o g φ 2 0.028 · l o g ( T ) · l o g ( φ ) ( R 2 = 0.978 ) 0 % φ 2 %   and   20   C T 60   C
It is to be mentioned that the thermal conductivity and dynamic viscosity are taken at an average temperature of 30 °C for all the performed numerical simulations.
The finite element method, utilizing Galerkin’s weighted residual technique, is applied for the resolution of the system of Equations (6)–(10). The Lagrange polynomials are represented as:
H = j = 1 N s ζ j s H j
  • ζ : shape function
  • H: value of the variable at the nodes of the elements.
  • The time-dependent components are handled using the backward differentiation formula.
The convergence criterion is: f n + 1 f n f n + 1 10 6 .
Model Verification and grid dependency test.
The present model’s verification (Table 2) is based on a comparison with the results of Ozoe and Okada [33]. The analysis demonstrates a notable correspondence between the results obtained, affirming the accuracy of the current model.
A grid sensitivity check was performed with four different grid sizes (G1, G2, G3, G4) for Ra = 105, φ = 0.02, W = 0.6, and Ha = 25, as Table 3. The Nuav is calculated for each grid size, and the results show that as the mesh gets finer, the Nuav also increases. The percentage increase in the Nuav is highest when moving from G1 to G2 and then starts to decrease with subsequent refinements, which is a sign that the solution is approaching grid independence. The incremental increase reduces for finer grids, with an incremental increase of about 0.6% when moving from G3 to G4. Compared to G4, the G3 mesh provides accurate results with less computational cost. It is important to select a grid size that balances precision with computational time. Thus, G3 is chosen over G4 for all numerical executions to save computation time while still providing accurate results.

3. Results and Discussion

In this section, the effect of the governing parameters on the control of three-dimensional natural convection of graphene–water nanofluids, through the use of a symmetrical tree-shaped obstacle and an external magnetic field, are presented. The goal is to understand how these control methods influence the thermal and velocity fields and to optimize the heat transfer and entropy generation.
Figure 2 shows the effects of a MF and the addition of NPs on the velocity magnitude along particle trajectories for Ra = 105 and fin’s length W = 0.2. The particle trajectories in pure water without a MF show a complex flow with high velocity magnitudes, indicating intense convection, especially near the hot surfaces. The addition of NPs to the fluid enhances the thermal conductivity and alters the flow characteristics. The particle trajectories still show a 3D complex flow pattern similar to the pure water case, but the velocity magnitudes are slightly more intense due to enhanced heat transfer and, therefore, increased buoyancy forces, leading to more vigorous convection. Introducing a MF (Ha = 50) to pure water significantly changes the flow structure. The particle trajectories are more organized and less chaotic, indicating that the Lorentz force is damping the fluid motion. Accordingly, the velocity magnitudes are reduced, which can be attributed to the fact that the MF has a stabilizing effect on the convective flow. When a MF is applied, the particle trajectories show that the fluid’s velocity magnitudes are higher than in the pure water case with a MF, but still less chaotic than without a MF. When the MF is applied at Ha = 50, a 50% reduction in the maximum velocity magnitude, indicated by the shift in color scales, is noticed. Due to the viscous effect, the use of nanoparticles slightly decreases the maximum velocity by around 16.7%.
The iso-surfaces of temperature are plotted in Figure 3, comparing the effects of a MF on the thermal field with and without graphene nanoparticles. In the absence of a MF, the iso-surfaces for both cases are curved, showing the dynamic nature of the convective flow. The complexity of the surfaces is due to the presence of the conduction and convection regimes, with the dominance of the convection for high Ra values. The presence of NPs in the nanofluid enhances the thermal conductivity and affects the flow behavior compared to pure water by enhancing the convective heat transfer. The tree obstacle induces flow disruptions, visible by the deformation of the iso-surfaces around it. Applying a MF, the gray iso-surfaces for pure water become more stretched with a horizontal stratification and become similar to a conductive regime due to the Lorentz force. The Lorentz forces arising from the interaction between the MF and the moving conductive fluid act to dampen the flow intensity, which reduces the convection. The nanofluid iso-surfaces under the influence of a MF show that the field’s damping effect is also present but to a lesser degree than in pure water. The enhanced thermal conductivity of the nanofluid counteracts some of the magnetic damping, maintaining a higher degree of thermal mixing and convective activity than in the pure water case. It is also to be mentioned that the application of the MF leads to a vertical symmetry of the flow structure.
The mid-plane velocity vector projections and thermal field are illustrated in Figure 4 for Ra = 105, W = 0.4 with and without graphene nanoparticles. In a non-magnetic environment, the pure water flow structure shows a classic convective pattern with strong recirculation zones indicated by the closed-loop velocity vectors. The tree obstacle has an important effect on the flow structure by creating a local vortex between the tree branches, especially for Ha = 0. When nanoparticles are added, the flow becomes more intense due to the increased thermal conductivity that favorizes the convective heat transfer. The temperature field for pure water is distorted at the top region due to the dominance of the convection and is quasi-vertical in the bottom region due to the presence of the obstacle that opposes the flow. For pure water under magnetic influence, the flow structure is noticeably more organized and symmetric, with reduced recirculation zones. This is an indication of the damping effect of the Lorentz forces that opposes the buoyancy force. The graphene–water nanofluid’s flow still exhibits greater motion compared to pure water under the same magnetic conditions. With a magnetic field, the temperature field in pure water is more uniform due to the reduction in the convective flow, leading to a more conduction-like heat transfer regime. When adding nanoparticles, even with the magnetic field, the compressed isotherms near the hot wall remain but are more stretched, showing that the nanofluid retains higher heat transfer rates than pure water but is also affected by the magnetic field. It can be observed that a rise in magnetic field strength from Ha = 0 to Ha = 100, in the case of both pure water and nanofluid, resulted in an important decrease in the maximum velocity values by about 82%. This implies that the magnetic field has a very strong damping effect on fluid flow intensity. The presence of nanoparticles has different effects on the fluid’s maximum velocity at different magnetic field conditions. With no magnetic field, the maximum velocity in the presence of nanoparticles is 5.45% less than for the pure water case. On the other hand, with a strong magnetic field, nanoparticles increase fluid maximum velocity by 1.12%, meaning a complex interaction between the magnetic field and adding nanoparticles.
Figure 5 displays the effects of the Rayleigh number on the mid-plane velocity vectors projection and the thermal field for fin’s length W = 0.6, φ = 0.02, and Ha = 50. At a relatively low Rayleigh number (Ra = 103), the flow structure exhibits a modest amount of convection with weak velocity vectors intensity. The isotherms are symmetric, parallel, and equidistant due to the dominance of the thermal conduction over convection, and the flow field is relatively stable and organized. Increasing the Rayleigh number to 104 enhances the buoyancy forces, leading to a more pronounced convective pattern in the flow structure. The velocity vectors are stronger, and the presence of multiple convection cells can be observed. The isotherms become more distorted, particularly close to the active walls, indicating stronger convection, which improves the heat transfer. At Ra = 105, the flow is dominated by the convection and has a complex structure due to the intense flow. The isotherms are highly distorted on the heated side due to the intense thermal convection. The influence of the tree obstacle on the flow structure and temperature field is more pronounced at this level. The Ra causes a strong increase in maximum velocity values as it increases. In fact, when the Ra is increased from 103 to 104, the increase in velocity is very large and is about 900%. Similarly, when the Ra is increased from 104 to 105, the increase is about 500%, illustrating the dominance of the convective heat transfer mode.
Figure 6 showcases the mid-plane velocity vectors projection and the thermal field for graphene–water nanofluid under a magnetic field (Ha = 25). Various fins’ lengths are presented, ranging from W = 0.2 to W = 0.8, while keeping the Ra constant at Ra = 105, which indicates that the flow is considerably influenced by buoyancy-driven convection. The fins’ lengths are a crucial aspect of the heat transfer enhancement in the cavity, and their effects on the thermal and flow fields are markedly significant. With shorter fins (W = 0.2), the velocity vectors show a well-defined large central convective roll, which is indicative of strong buoyancy forces. The thermal field shows a relatively smooth gradient with evenly spaced isotherms in the middle portion. Increasing the fin length to W = 0.4, we observe the development of more complex flow patterns, with the large central vortex being disrupted by secondary circulations. The thermal field shows tighter isotherms on the heated side, indicating enhanced heat transfer. Further increasing the fins’ length, the flow becomes more intense, as seen by the increased number and complexity of circulation patterns. The thermal field reflects this with more distorted isotherms, especially around the obstacle, indicating that the mixing effect of the flow is intensifying with the longer fins. The reduction in velocity is about 23.61% as the W increases from 0.2 to 0.8. In fact, when the fins become longer, they create more resistance to the flow, which slows down the velocity.
Figure 7 displays the total entropy generation for Ra = 105 and fin’s length W = 0.2 within the central plane of the differentially heated cavity with and without graphene nanoparticles. Without a magnetic field and without the addition of nanoparticles, there are high levels of entropy generation near the heated and cooled walls, where temperature gradients, and thus thermal irreversibilities, are greatest. The addition of NPs enhances the thermal conductivity, and thus improves the heat transfer. However, the entropy generation remains high close to the active walls but with different distributions. The application of a MF has a noticeable effect on the entropy generation. The patterns are altered, and we observe a localized decrease in the entropy generation where magnetic damping of the flow is strongest. The entropy generation levels are lower due to the considerable suppression of the NC. Moreover, the MF reduces the entropy generation in the nanofluid; however, the reduction is not as pronounced as in pure water. The addition of nanoparticles, which enhance thermal properties, leads to a resistance against the magnetic stabilizing effect. Nevertheless, the distribution of entropy generation magnitudes is lower compared to the case of absence of a MF. Due to the increased viscous effects, the use of nanoparticles increases the entropy generation by 74.44% and 108.86%, for Ha = 0 and 50, respectively. Due to the magnetic damping effect, the increase of the Ha from 0 to 50 leads to a decrease in the total entropy generation by 76.24% and 71.55% for φ = 0 and φ = 0.02, respectively.
The variation of the Nuav as a function of the Ha and nanoparticle concentration for shorter fins (W = 0.2) and longer fins (W = 0.8) is illustrated in Figure 8 for Ra = 105. For W = 0.2, as the Ha increases, there is a general decrease in the Nuav, indicating that the MF has a damping effect on the convective heat transfer due to the Lorentz force suppressing the convection and therefore reducing the efficiency of convective heat transfer. For a given Ha, increasing φ increases the Nuav, which is attributed to the enhanced thermal properties with adding graphene nanoparticles that improve the efficiency of convective heat transfer despite magnetic damping. Similarly for W = 0.8, an increase in the Ha generally leads to a reduction in the Nuav. This effect appears to be more pronounced for larger fins’ length, indicating an even stronger damping effect of the MF on convection. As with the shorter fins, increasing φ for a given Ha enhances the Nuav; however, the overall Nusselt numbers appear to be slightly lower than for W = 0.2. While the nanoparticles improve heat transfer, the larger fins could create more flow disturbance and resistance, somewhat counteracting the benefits of the enhanced thermal conductivity.
Figure 9 presents the surface plot showing the change in the Nuav against the Ra and fins’ length for Ha = 50 and φ = 0.02. As the Ra increases, the buoyancy forces become more significant, and the convection becomes more efficient, which is reflected in an increase in the Nuav. The steeper slope of the surface at a lower W indicates that the Nuav is more sensitive to changes in the Ra when the fins are shorter. The plot shows that as the W increases, there is a very slight decrease in the Nuav, which is attributed to the fact that wider fins impede fluid flow. In addition, it is clear that the increase in the Ra increases the Nuav for all the values of W. The effect of the W on the Nuav is less pronounced than that of Ra, but there is a general trend of decreasing Nuav with increasing W, especially at higher Ra values.
Figure 10 illustrates how the different types of entropy generation, thermal (Sth), frictional (Sfr), magnetic (Smag), and total (Stot), vary with the Ha and φ at Ra = 105 and W = 0.4. An increase in the Ha results in increased Smag but also opposes the flow, which potentially reduces thermal and frictional entropy generation. The stabilization becomes more effective so that the overall system becomes more ordered, thus reducing Smag. The addition of NPs enhances the thermal conductivity and leads to stronger temperature gradients and higher thermal entropy generation (Sth). For Ra = 105, the temperature gradients are more important, and higher Sth occurs. The frictional entropy (Sfr) is generated due to fluid viscosity and flow resistance. The flow around the tree obstacle creates frictional losses, and the fluid’s viscosity is increased by the addition of nanoparticles. Thus, Sfr is considerably affected by the increase in φ. The Stot is the sum of all forms of entropy generations. The figure shows a decrease in the Stot with the increase in the MF magnitude and an increase with the increase in the NPs volume fraction.
Figure 11 represents a series of surface plots showing the variations of thermal (Sth), frictional (Sfr), magnetic (Smag), and total (Stot) entropy generation as functions of the Ra and the fin length for Ha = 50 and a NPs volume fraction φ = 0.02. The variation of Sth shows that as the Ra increases, indicating stronger convective flows due to higher temperature differences, Sth increases. This is because larger temperature gradients, which are more prevalent at a higher Ra, lead to greater thermal irreversibility. Similarly, varying the W may show that longer fins initially increase Sth due to greater surface area affecting heat transfer but could potentially lead to a decrease in Sth at some point where the flow becomes too disrupted by the fins, leading to less efficient heat exchange. Sfr generally arises from viscous dissipation in the fluid flow. As the Ra increases, the flow becomes more intense, which increases frictional losses, thus increasing Sfr. However, the increase in the fin length W could show a more complex trend, as it introduces more surface for friction but also changes the flow pattern, which could either increase or decrease Sfr. With a consistent Ha of 50, Smag is influenced primarily by the flow velocity and electrical properties of the fluid, both affected by φ. The plot for Smag shows a less pronounced variation with Ra, which indicates that the magnetic field effectively dampens the convection, leading to a more orderly flow despite the increase in buoyancy forces. However, the increase in the fin’s length for higher Rayleigh numbers reduces Smag because as the W increases, the flow becomes more obstructed in such a way that it has less interaction with the magnetic field, which decreases the magnetic entropy generation. Overall, at a higher Ra and decreasing W, the Stot shows an increase due to the combined effects of thermal, frictional, and magnetic entropy generations, indicating that there is an optimal fin’s length where the Stot is minimized for a more efficient system design.
Regression analysis:
The aim of this section is to establish correlations allowing the direct calculation of the Nuav and the Stot. Due to a number of parameters and their mutual interactive effects, the Lasso regression model is chosen to reach a high accuracy of predictivity. The use of the Lasso regression model with interaction terms and squared variables allows for simultaneous feature selection and regularization. In fact, Lasso regression helps in reducing the complexity of the model by penalizing the absolute size of the regression coefficients. It also allows for the shrinking of some coefficients to zero and effectively selects a simple model to avoid overfitting.
The Lasso regression model can be expressed as minimizing the following objective function:
m i n β 1 2 n i = 1 p y i β 0 j = 1 p X i j β j 2 + λ j = 1 p β j
where:
  • n is the number of observations,
  • yi is the observed value of the dependent variable for the ith observation,
  • Xij is the value of the jth predictor for the ith observation,
  • βj are the coefficients (including the intercept β0 and those for all predictors and their interactions),
  • p is the total number of predictors (including interaction terms),
  • λ is the regularization parameter that controls the strength of the penalty (the higher the value, the more significant the amount of shrinkage).
The first term represents the residual sum of squares (RSS), and the second term is the Lasso penalty.
As a result of the Lasso regression model, the following expressions are obtained for the Nuav and the Stot:
N u a v = 1.5106 0.0805 · W 0.2091 · H a + 0.1883 · φ + 0.6700 · R a + 0.0497 · W 2 + 0.1035 · W · H a 0.0323 · W · φ 0.0981 · W · R a + 0.2088 · H a 2 0.1078 · H a · φ 0.5896 · H a · R a + 0.0319 · φ 2 + 0.1680 · φ · R a + 0.1485 · R a 2 + ϵ M S E = 0.02 R 2 = 0.970 λ = 3.74 × 10 4
S t o t = 1.7635 0.1608 · W 0.3276 · H a + 0.1979 · φ + 1.0870 · R a + 0.0828 · W 2 + 0.2396 · W · H a 0.0558 · W · φ 0.3548 · W · R a + 0.2496   H a 2 0.1163 · H a · φ 0.8366 · H a · R a + 0.0405 · φ 2 + 0.1662 · φ · R a + 0.4805 · R a 2 + ϵ M S E = 0.0463 R 2 = 0.959 λ = 7.56 × 10 4
Equations (22) and (23) allow the prediction of the Nuav and the Stot as functions of the variables Ra, Ha, W, φ, and their significant interactions and squared terms. ϵ represents the model error term (deviations of the predictions from the actual values).
For the Nuav and the Stot, the Lasso model has effectively captured the complexity of the relationships while also applying regularization to prevent overfitting. In fact, the R-squared values of 0.970 and 0.959 mean that approximately 97% of the variance in the Nuav and 95.9% of the variance in the Stot are explained by the models, which means the models fit the data well. The MSE values are relatively low, indicating that the model’s predictions are, on average, close to the actual values. The values of λ are obtained using a cross-validation process between fitting the data well and keeping the model complexity under control.
Figure 12 presents the plots of residuals versus predicted values for the Nuav and the Stot. The plot of residuals allows the assessment of the homoscedasticity of the residuals and the identification of non-linearities or outliers that can affect the models. From the plots of residuals versus predicted values, it is clear that the residuals are scattered around the horizontal line at zero, corresponding to a well-fitting model. The residuals are distributed across the range of the predicted Nuav and Stot values. Thus, it can be mentioned that the obtained expressions predict the functions with accuracy and without obvious biases or non-linearities. The variance of the residuals is consistent across the range of predicted values, which indicates the homoscedasticity of the model. This property is desirable because it indicates that the model is stable. Finally, it can be concluded that the absence of clear patterns or systematic bias in these plots suggests that the Lasso regression model is performing well.

4. Conclusions

The achieved study investigates the control of three-dimensional NC of graphene–water nanofluid within a 3D cubic differentially heated cavity influenced by both a tree-shaped obstacle and an external MF. It highlights the interplay between the effects of a MF, the addition of graphene nanoparticles, and varying fins’ length on the flow structure and convective heat transfer efficiency. The main findings reveal that:
  • The MF has suppressive effects on the NC flow of graphene–water nanofluid; about an 82% decrease in the maximum velocity values occurs for both pure water and nanofluid when the Ha = 100.
  • In the absence of a MF, the addition of graphene NPs decreases the maximum velocity by 5.45%, while for Ha = 100, it increases it by 1.12%.
  • The enhanced thermal conductivity due to the addition of graphene NPs improves the convective heat transfer as Rayleigh numbers increase.
  • The incorporation of the tree-shaped obstacle allows a considerable control of the flow and affect of the temperature field. In fact, a reduction in the maximum velocity by about 23.61% occurs when the fins’ length is varied from W = 0.2 to W = 0.8.
  • The increase in the length of the branches considerably reduces the heat transfer, especially at higher Ra values.
  • The use of nanoparticles increases the entropy generation by 74.44% and 108.86% for Ha = 0 and 50, respectively.
Overall, the present study focuses on the potential of using the combination of a non-destructive technique (magnetic field) and a destructive technique (tree-shaped obstacle) on the fluid flow, heat transfer, and entropy production during graphene–water nanofluid 3D natural convection. The study can find its applications in thermal management by enhancing the operational efficiency of electronic cooling systems and other applications requiring precise thermal control.

Author Contributions

Conceptualization, I.H.-J., W.A. and L.K.; methodology, C.M., B.M.A. and A.A.; software, I.H.-J., W.A. and L.K.; validation, C.M., B.M.A. and A.A.; formal analysis, A.A. and C.M.; investigation, B.M.A. and L.K.; writing—original draft preparation, I.H.-J., W.A., L.K., C.M., B.M.A. and A.A.; writing—review and editing, I.H.-J., W.A., L.K., C.M., B.M.A. and A.A.; project administration, I.H.-J.; funding acquisition, I.H.-J. All authors have read and agreed to the published version of the manuscript.

Funding

This research project was funded by the Deanship of Scientific Research and Libraries, Princess Nourah bint Abdulrahman University, through the Program of Research Project Funding After Publication, grant No (RPFAP-24-1445).

Data Availability Statement

Data are contained with the article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

B dimensionless magnetic field = B / B 0
BeBejan number
Cpspecific heat capacity, kJ·kg−1·K−1
e B direction of the magnetic field
E dimensionless electric field
F dimensionless Lorentz force
g gravitational acceleration, m·s−2
HaHartmann number
ISirreversibility coefficient
J dimensionless density of electrical current
Kthermal conductivity, kW·m−1·K−1
NuNusselt number
NSdimensionless local generated entropy
Pdimensionless pressure
Pr Prandtl number
Ra Rayleigh number
S′gen entropy generation, J·K−1
Stotdimensionless total generated entropy
Tdimensionless temperature
t dimensionless time
V dimensionless velocity vector
Wfins’ length
x,y,zdimensionless Cartesian coordinates
Greek Characters
α thermal diffusivity, m2·s−1
β thermal expansion coefficient, K−1
μdynamic viscosity, Pa·s
νkinematic viscosity, m2·s−1
ρdensity of the fluid, kg·m−3
σelectric conductivity, Ω−1·m−1
Subscripts
avaverage
c cold
frfriction
h hot
m magnetic
nf nanofluid
p solid particle
ththermal

References

  1. Liu, Y.; Tian, Z.; Qi, C.; Li, R. Natural convection and entropy generation of Fe3O4-H2O nanofluids in square cavities with cylindrical grooves under magnetic field. Colloids Surf. A Physicochem. Eng. Asp. 2023, 679, 132564. [Google Scholar] [CrossRef]
  2. Tayebi, T.; Sattar Dogonchi, A.; Karimi, N.; Ge-JiLe, H.; Chamkha, A.J.; Elmasry, Y. Thermo-economic and entropy generation analyses of magnetic natural convective flow in a nanofluid-filled annular enclosure fitted with fins. Sustain. Energy Technol. Assess. 2021, 46, 101274. [Google Scholar] [CrossRef]
  3. Saha, T.; Islam, T.; Yeasmin, S.; Parveen, N. Thermal influence of heated fin on MHD natural convection flow of nanofluids inside a wavy square cavity. Int. J. Thermofluids 2023, 18, 100338. [Google Scholar] [CrossRef]
  4. Iftikhar, B.; Javed, T.; Siddiqui, M.A. Natural convective flow of non-isothermal hybrid nanofluid inside the cavity with the influence of a heated fin and non-linear thermal radiation: Second law analysis. Mater. Today Commun. 2023, 34, 105341. [Google Scholar] [CrossRef]
  5. Nayak, M.K.; Dogonchi, A.S.; Rahbari, A. Free convection of Al2O3-water nanofluid inside a hexagonal-shaped enclosure with cold diamond-shaped obstacles and periodic magnetic field. Case Stud. Therm. Eng. 2023, 50, 103429. [Google Scholar] [CrossRef]
  6. Islam, M.S.; Islam, S.; Siddiki, M.N.A.A. Numerical simulation with sensitivity analysis of MHD natural convection using Cu-TiO2-H2O hybrid nanofluids. Int. J. Thermofluids 2023, 20, 100509. [Google Scholar] [CrossRef]
  7. Vinodhini, N.; Prasad, V.R. NUMERICAL study of MAGNETO convective Buongiorno nanofluid flow in a rectangular enclosure under oblique magnetic field with heat generation/absorption and complex wall conditions. Heliyon 2023, 9, e17669. [Google Scholar] [CrossRef] [PubMed]
  8. Al-Farhany, K.; Al-Chlaihawi, K.K.; Al-dawody, M.F.; Biswas, N.; Chamkha, A.J. Effects of fins on magnetohydrodynamic conjugate natural convection in a nanofluid-saturated porous inclined enclosure. Int. Commun. Heat Mass Transf. 2021, 126, 105413. [Google Scholar] [CrossRef]
  9. Cheng, H.; Sajadi, S.M.; Jouzdani, M.Z. Control of hybrid nanofluid natural convection with entropy generation: A LBM analysis based on the irreversibility of thermodynamic laws. Eng. Anal. Bound. Elem. 2023, 155, 108–120. [Google Scholar] [CrossRef]
  10. Abderrahmane, A.; Manoongam, A.; Alizadeh, A.; Younis, O.; Zekri, H.; Isa, S.S.P.M.; Baghaei, S.; Jamshed, W.; Guedri, K. Investigation of the free convection of nanofluid flow in a wavy porous enclosure subjected to a magnetic field using the Galerkin finite element method. J. Magn. Magn. Mater. 2023, 569, 170446. [Google Scholar] [CrossRef]
  11. Acharya, N. Hydrothermal scenario of buoyancy-driven magnetized multi-walled carbon nanotube-Fe3O4-water hybrid nanofluid flow within a discretely heated circular chamber fitted with fins. J. Magn. Magn. Mater. 2024, 589, 1–25. [Google Scholar] [CrossRef]
  12. Waqas, H.; Khan, S.A.; Yasmin, S.; Liu, D.; Imran, M.; Muhammad, T.; Alhushaybari, A.; Farooq, U. Galerkin finite element analysis for buoyancy driven copper-water nanofluid flow and heat transfer through fins enclosed inside a horizontal annulus: Applications to thermal engineering. Case Stud. Therm. Eng. 2022, 40, 102540. [Google Scholar] [CrossRef]
  13. Armaghani, T.; Sepehrnia, M.; Molana, M.; Nour, M.M.; Safari, A. Studying Alumina–Water Nanofluid Two-Phase Heat Transfer in a Novel E-Shaped Porous Cavity via Introducing New Thermal Conductivity Correlation. Symmetry 2023, 15, 2057. [Google Scholar] [CrossRef]
  14. Khan, S.A.; Yasmin, S.; Imran, M.; Muhammad, T.; Alhushaybari, A.; Farooq, U.; Waqas, H. Computational analysis of natural convection with water based nanofluid in a square cavity with partially active side walls: Applications to thermal storage. J. Mol. Liq. 2023, 382, 122003. [Google Scholar] [CrossRef]
  15. Thirumalaisamy, K.; Sivaraj, R.; Reddy, A.S. Fluid flow and heat transfer analysis of a ternary aqueous Fe3O4 + MWCNT + Cu/H2O magnetic nanofluid in an inclined rectangular porous cavity. J. Magn. Magn. Mater. 2024, 589, 171503. [Google Scholar] [CrossRef]
  16. Tasnim, S.; Mitra, A.; Saha, H.; Islam, M.Q.; Saha, S. MHD conjugate natural convection and entropy generation of a nanofluid filled square enclosure with multiple heat-generating elements in the presence of Joule heating. Results Eng. 2023, 17, 100993. [Google Scholar] [CrossRef]
  17. Seyyedi, S.M.; Hashemi-Tilehnoee, M.; Palomo del Barrio, E.; Dogonchi, A.S.; Sharifpur, M. Analysis of magneto-natural-convection flow in a semi-annulus enclosure filled with a micropolar-nanofluid; a computational framework using CVFEM and FVM. J. Magn. Magn. Mater. 2023, 568, 170407. [Google Scholar] [CrossRef]
  18. Siavashi, M.; Yousofvand, R.; Rezanejad, S. Nanofluid and porous fins effect on natural convection and entropy generation of flow inside a cavity. Adv. Powder Technol. 2018, 29, 142–156. [Google Scholar] [CrossRef]
  19. Çiçek, O.; Sheremet, M.A.; Baytaş, A.C. Effect of natural convection hybrid nanofluid flow on the migration and deposition of MWCNT-Fe3O4 in a square enclosure. Int. J. Therm. Sci. 2023, 190, 108318. [Google Scholar] [CrossRef]
  20. Acharya, N. Buoyancy driven magnetohydrodynamic hybrid nanofluid flow within a circular enclosure fitted with fins. Int. Commun. Heat Mass Transf. 2022, 133, 105980. [Google Scholar] [CrossRef]
  21. Ali, F.H.; Hamzah, H.K.; Egab, K.; Arıcı, M.; Shahsavar, A. Non-Newtonian nanofluid natural convection in a U-shaped cavity under magnetic field. Int. J. Mech. Sci. 2020, 186, 105887. [Google Scholar] [CrossRef]
  22. Qi, C.; Tang, J.; Ding, Z.; Yan, Y.; Guo, L.; Ma, Y. Effects of rotation angle and metal foam on natural convection of nanofluids in a cavity under an adjustable magnetic field. Int. Commun. Heat Mass Transf. 2019, 109, 104349. [Google Scholar] [CrossRef]
  23. Liao, C.C.; Li, W.K. Effect of different magnetic field angles on the relationship between nanofluid concentration and heat transfer. Int. Commun. Heat Mass Transf. 2022, 135, 106137. [Google Scholar] [CrossRef]
  24. Vahedi, S.M.; Pordanjani, A.H.; Aghakhani, S.; Afrand, M.; Daneh-Dezfuli, A.; Azaiez, J. Investigation of the impact of fin number and length attached to an enclosure subjected to a magnetic field. Eng. Anal. Bound. Elem. 2023, 152, 678–687. [Google Scholar] [CrossRef]
  25. Acharya, N. On the hydrothermal behavior and entropy analysis of buoyancy driven magnetohydrodynamic hybrid nanofluid flow within an octagonal enclosure fitted with fins: Application to thermal energy storage. J. Energy Storage 2022, 53, 105198. [Google Scholar] [CrossRef]
  26. Cao, Y.; Ayed, H.; Jarad, F.; Togun, H.; Alias, H.; Issakhov, A.; Dahari, M.; Wae-hayee, M.; El Ouni, M.H. MHD natural convection nanofluid flow in a heat exchanger: Effects of Brownian motion and thermophoresis for nanoparticles distribution. Case Stud. Therm. Eng. 2021, 28, 101394. [Google Scholar] [CrossRef]
  27. Selimefendigil, F.; Öztop, H.F.; Izadi, F. Non-uniform magnetic field effects on the phase transition dynamics for PCM-installed 3D conic cavity having ventilation ports under hybrid nanofluid convection. J. Build. Eng. 2022, 49, 104074. [Google Scholar] [CrossRef]
  28. Alazzam, A.; Qasem, N.A.A.; Aissa, A.; Abid, M.S.; Guedri, K.; Younis, O. Natural convection characteristics of nano-encapsulated phase change materials in a rectangular wavy enclosure with heating element and under an external magnetic field. J. Energy Storage 2023, 57, 106213. [Google Scholar] [CrossRef]
  29. Gürdal, M.; Arslan, K.; Gedik, E.; Minea, A.A. Effects of using nanofluid, applying a magnetic field, and placing turbulators in channels on the convective heat transfer: A comprehensive review. In Renewable and Sustainable Energy Reviews; Elsevier Ltd.: Amsterdam, The Netherlands, 2022; Volume 162. [Google Scholar] [CrossRef]
  30. Narankhishig, Z.; Ham, J.; Lee, H.; Cho, H. Convective heat transfer characteristics of nanofluids including the magnetic effect on heat transfer enhancement—A review. Appl. Therm. Eng. 2021, 193, 116987. [Google Scholar] [CrossRef]
  31. Gal, S.; Kolsi, L.; Hassen, W.; Ali, N.B.; Khedher, N.B.; Chamkha, A.J. Three-dimensional study of magnetohydrodynamic natural convection, entropy generation and electromagnetic variables in a nanofluid filled enclosure equipped with inclined fins. ACS Omega 2022, 7, 12365–12373. [Google Scholar] [CrossRef]
  32. Omri, M.; Aich, W.; Rmili, H.; Kolsi, L. Experimental analysis of the thermal performance enhancement of a vertical helical coil heat exchanger using copper oxide-graphene (80–20%) hybrid nanofluid. Appl. Sci. 2022, 12, 11614. [Google Scholar] [CrossRef]
  33. 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 Transfer. 1989, 32, 1939–1954. [Google Scholar] [CrossRef]
Figure 1. Studied configuration; (a) schematic view of 3D cubic cavity fitted with horizontal fins and coordinate system; (b) central plan. (c) used mesh.
Figure 1. Studied configuration; (a) schematic view of 3D cubic cavity fitted with horizontal fins and coordinate system; (b) central plan. (c) used mesh.
Symmetry 16 00692 g001
Figure 2. Particle trajectories (colored via velocity magnitude) for Ra = 105 and W = 0.2 for pure water and nanofluid with and without magnetic field.
Figure 2. Particle trajectories (colored via velocity magnitude) for Ra = 105 and W = 0.2 for pure water and nanofluid with and without magnetic field.
Symmetry 16 00692 g002aSymmetry 16 00692 g002b
Figure 3. Iso-surfaces of temperature with and without magnetic field for fins’ length. W = 0.2 and Ra = 105, φ = 0.02 (colored) and φ = 0 (gray).
Figure 3. Iso-surfaces of temperature with and without magnetic field for fins’ length. W = 0.2 and Ra = 105, φ = 0.02 (colored) and φ = 0 (gray).
Symmetry 16 00692 g003
Figure 4. Mid-plan velocity vectors projection (left) and thermal field (right) for Ra = 105, W = 0.4 with and without nanoparticles.
Figure 4. Mid-plan velocity vectors projection (left) and thermal field (right) for Ra = 105, W = 0.4 with and without nanoparticles.
Symmetry 16 00692 g004
Figure 5. Mid-plan velocity vectors projection (left) and thermal field (right) for Ha = 50, W = 0.6, φ = 0.02, and various Ra.
Figure 5. Mid-plan velocity vectors projection (left) and thermal field (right) for Ha = 50, W = 0.6, φ = 0.02, and various Ra.
Symmetry 16 00692 g005
Figure 6. Mid-plan velocity vectors projection (left) and thermal field (right) for Ra = 105, Ha = 25, φ = 0.02 and various fins’ lengths.
Figure 6. Mid-plan velocity vectors projection (left) and thermal field (right) for Ra = 105, Ha = 25, φ = 0.02 and various fins’ lengths.
Symmetry 16 00692 g006
Figure 7. Total entropy generation contours at the central plan for Ra = 105, W = 0.2 with and without nanoparticles.
Figure 7. Total entropy generation contours at the central plan for Ra = 105, W = 0.2 with and without nanoparticles.
Symmetry 16 00692 g007
Figure 8. Variation of Nuav versus Ha and φ for Ra = 105 (a) W = 0.2 and (b) W = 0.8.
Figure 8. Variation of Nuav versus Ha and φ for Ra = 105 (a) W = 0.2 and (b) W = 0.8.
Symmetry 16 00692 g008
Figure 9. Variation of Nuav versus Ra and W for Ha = 50 and φ = 0.02.
Figure 9. Variation of Nuav versus Ra and W for Ha = 50 and φ = 0.02.
Symmetry 16 00692 g009
Figure 10. Variation of thermal (Sth), frictional (Sfr), magnetic (Smag), and total entropy generation (Stot) as functions of Ha and φ, for Ra = 105 and W = 0.4.
Figure 10. Variation of thermal (Sth), frictional (Sfr), magnetic (Smag), and total entropy generation (Stot) as functions of Ha and φ, for Ra = 105 and W = 0.4.
Symmetry 16 00692 g010
Figure 11. Variation of thermal (Sth), frictional (Sfr), magnetic (Smag), and total entropy generation (Stot) as functions of Ra and W, for Ha = 50 and φ = 0.02.
Figure 11. Variation of thermal (Sth), frictional (Sfr), magnetic (Smag), and total entropy generation (Stot) as functions of Ra and W, for Ha = 50 and φ = 0.02.
Symmetry 16 00692 g011
Figure 12. Plot of residuals versus predicted values of (a) Nuav and (b) Stot for the Lasso regression model.
Figure 12. Plot of residuals versus predicted values of (a) Nuav and (b) Stot for the Lasso regression model.
Symmetry 16 00692 g012
Table 1. Properties of water and graphene at 30 °C.
Table 1. Properties of water and graphene at 30 °C.
WaterCNT
C p (J·Kg−1·K−1)4180425
ρ (kg·m−3)995.872600
β   ( K 1 ) 3.03 × 10−416 × 10−7
σ (Ω−1·m−1)0.5 × 10−34.8 × 10−9
Table 2. Comparison of N u a v with the results of Ozoe and Okada [33].
Table 2. Comparison of N u a v with the results of Ozoe and Okada [33].
Ha = 100Present Model4.4481
[33]4.457
Ha = 200Present Model2.9192
[33]2.917
Table 3. Grid sensitivity check for Ra = 105, φ = 0.02, W = 0.6, and Ha = 25.
Table 3. Grid sensitivity check for Ra = 105, φ = 0.02, W = 0.6, and Ha = 25.
NuavPercentage IncreaseIncremental Increase
G1: 3329842.6193--
G2: 8533262.84098.460276-
G3: 24722472.971313.43874.590095
G4: 99737002.989214.122090.60243
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

Aich, W.; Hilali-Jaghdam, I.; Alshahrani, A.; Maatki, C.; Alshammari, B.M.; Kolsi, L. Control of Three-Dimensional Natural Convection of Graphene–Water Nanofluids Using Symmetrical Tree-Shaped Obstacle and External Magnetic Field. Symmetry 2024, 16, 692. https://0-doi-org.brum.beds.ac.uk/10.3390/sym16060692

AMA Style

Aich W, Hilali-Jaghdam I, Alshahrani A, Maatki C, Alshammari BM, Kolsi L. Control of Three-Dimensional Natural Convection of Graphene–Water Nanofluids Using Symmetrical Tree-Shaped Obstacle and External Magnetic Field. Symmetry. 2024; 16(6):692. https://0-doi-org.brum.beds.ac.uk/10.3390/sym16060692

Chicago/Turabian Style

Aich, Walid, Inès Hilali-Jaghdam, Amnah Alshahrani, Chemseddine Maatki, Badr M. Alshammari, and Lioua Kolsi. 2024. "Control of Three-Dimensional Natural Convection of Graphene–Water Nanofluids Using Symmetrical Tree-Shaped Obstacle and External Magnetic Field" Symmetry 16, no. 6: 692. https://0-doi-org.brum.beds.ac.uk/10.3390/sym16060692

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