With a growing demand for improvements in heat transfer performances, the development of advanced technologies to manage thermal issues more efficiently are one of today’s hottest topics across different engineering and industrial segments. An innovative technique that recently attracted much attention consist of the utilization of HNFs to enhance NC HT. HNFs are formulated via an innovative process in which two different types of NPs are suspended in a base fluid (BF), with the intention to capture substantial heat transfer qualities from each type. For this purpose, a wide variety of applications are found in the literature for natural convection use of hybrid nanofluids. These are used in electronics cooling such as microprocessors or power electronic components, solar thermal systems and building thermal management. Recent works mainly dedicated to the use of HNFs, which is a combination of different nanoparticles, have shown their advantage in the improvement of the convective HT coefficient. Due to their enhanced properties, better stability and optimized thermal performance/cost, the applications of HNFs in convective heat transfer became the interest of several researchers.

Studies related to hybrid nanofluids showed important improvements in thermal conductivity (TC) and HT performance. Omri et al. (2022) presented a study on the performance’s enhancement of a vertical coiled heat exchanger working with CuO-Gp/water HNF. After ultrasonic dispersion (to avoid agglomeration), the authors measured the TC and viscosity of the HNF. It was found that the HT coefficient increased by 23.65% at a 0.2% mass fraction and by up to 79.68% at a 1% mass fraction. In addition, the increase of Reynolds number (Re) caused the augmentation of the HT coefficient. The HT and FF in a wavy cross-flow micro heat exchanger working with HNF have been studied by Ghachem et al. (2021). They considered a 3D configuration and studied the effects of different factors, like the fluid velocity, the waves number, and nanoparticles concentration. The HNF, which combines carbon nanotubes (CNT) and aluminum oxide (Al₂O₃) nanoparticles, showed significant improvements in both thermal conductivity and heat transfer rates, especially when the flow velocity was high. The researchers discovered that finding the optimal combinations of these parameters could reduce the size of the heat exchanger while boosting its efficiency. Dero et al. (2022) used the bvp4c solver to study the thermal stability of HNFs in a porous media using. The dual branch framework was used to highlight the effects of viscous dissipation and suction/injection effects on the thermal performances and stability of Cu-Al₂O₃/water hybrid nanofluid. Under certain conditions, dual solutions were found, this behavior indicates the complex nature of the thermal and fluid dynamics.

Natural convection corresponds to fluid motions caused by density variations due to temperature gradients. It occurs in several applications such as cooling electronic devices, solar energy, and in buildings. The conventional HT fluids such as water, ethylene glycol, and lubricants frequently lead to low thermal performances. Single nanofluids, are efficient in the enhancements of thermal conductivity and convective heat transfer, but they have limitations related to particle aggregation and sedimentation. The use of hybrid nanofluids solve these problems, in fact dispersing different nanoparticles in a base fluid allows to benefit from the distinct thermal, chemical, and mechanical characteristics of each nanoparticle type to attain better stability, higher TC, and enhanced specific heat capacity. HNFs exhibit enhanced effective thermal conductivity in comparison to conventional nanofluids, leading to enhanced heat transfer efficiency. Jino et al. (2024) introduced a new technique for cooling electric vehicle batteries. They used hybrid nanofluids in inverted right-angled triangular porous cavities. This study emphasises the significance of thermal management in the enhancement of the performances, safety, and longevity of electric vehicle batteries. The authors analyzed the impact of various parameters (Ra, Da, Ha, and φ ), on HT, FF, and irreversibility. The findings demonstrated that the HT rate was considerably increased when an inverted triangular cavity is used compared to traditional square-shaped cavities. Khan et al. (2024) examined the characteristics of a HNF in a permeable quadrantal cavity with internal heat generation. The researchers employed the FEM to solve the nondimensionalized equations. They considered the impact of heat generation on the FF and the HT. The study observed a positive correlation between the φ and Nu_a, suggesting a notable improvement in heat transport. The NC of MWCNT-Fe₃O₄/water HNF was studied by Khan N. H. et al. (2022) with a focus on the impact of changes in Ra, aspect ratio, and φ on the HT and entropy formation. The authors mentioned that Nu_a is proportional to Ra but reversely proportional to the aspect ratio. In addition, it was found that the use of MWCNT-Fe₃O₄ nanoparticles improves the HT rate. Goudarzi et al. (2020) studied the effects of NP migration caused by the thermophoresis and Brownian motion on the natural convection of Ag-MgO/water HNF in a wavy enclosure. The numerical study considered the diffusion of Ag nanoparticles in MgO-water NF and vice versa. Khan S. A. et al. (2022) investigated numerically the NC of Cu-Al₂O₃/water HNF inside an annuli. In addition, they performed a multi-response optimization to identify the optimal parameter to maximize the thermal and hydrodynamic performances. The results indicate that inclination angle and heat flux significantly enhance heat transfer, while the aspect ratio and nanoparticle concentration also play vital roles in improving the thermal performances of the HNF. Çiçek et al. (2023) presented a study on the NC of MWCNT-Fe₃O₄ hybrid nanofluid within a differentially heated square enclosure. The study consists of an analysis of the impacts of Ra and nanoparticle volume fraction on thermal and flow behaviors. The authors mentioned that the increase of φ leads to an improvement of Nu_a at low Ra values, with a 0.1% volume fraction which is identified as the optimal for the HT enhancement at Ra equal to or greater than 10⁴. Nayak et al. (2022) studied the effect of wavy baffle structures on the HT performance of double-diffusive NC within a C-shaped cavity filled with a HNF. The results indicate that streamlines intensify significantly with increased buoyancy ratio and Rayleigh number, enhancing heat and mass transfer. Furthermore, the presence of wavy baffles markedly improves thermal performance by amplifying local Nusselt numbers, especially at higher Rayleigh numbers. Aich et al. (2023) performed a three-dimensional study on the NC and entropy production in a C-shaped enclosure filled with a CNT–Al₂O₃/water HNF. The aim of the study was to examine the effects of the aspect ratio, Ra, and φ on HT, FF, and irreversibility production. It was found that the aspect ratio has a significant effect on the HT and entropy. At high Ra the entropy increases due to the enhancement of the convective HT. In addition, the inclusion of hybrid nanoparticles has a benefic effect on Nu_a, which is an indication of the improvement of heat transfer performance. Almeshaal et al. (2020) conducted a 3D study on NC within a T-shaped enclosure filled with a water-based CNT–aluminum oxide HNF. The was based on the vorticity-vector potential formalism and it was shown that the increase of the size of the enclosure, Ra, and φ cause and enhancement of the HT rate.

The use of internal obstacles in cavities filled with hybrid nanofluids corresponds to an effective strategy for the enhancement of the convective HT. When strategically placed, the obstacles cause the disruption of the FF, and create secondary vortices that lead to a more uniform temperature distribution and increased heat transfer rates. Hybrid nanofluids, consist of the combination of the superior thermal properties of different nanoparticles, for the enhancement of the thermal conductivity and the convective heat transfer. This method is highly beneficial for applications such as electronic cooling, solar collectors, and energy-efficient building designs. Some recent studies have examined the impact of the internal obstacles, like heat-generating bodies and fins, on the HT performance in cavities filled with HNF. In the study of Hassen et al. (2022) the thermocapillary and buoyancy-driven convection in a cavity containing a heated obstacle filled with a MWCNT-Fe₃O₄/oil HNF was numerically studied. The study was based on 2D simulations for the examination of the effects of the obstacle’s position, size, Marangoni number, Ra, and φ on the flow and thermal behavior. The findings indicated that repositioning the obstacle from the bottom to the top of the cavity leads to the optimization of the heat transfer enhancement by over 170%. HT and magnetohydrodynamic (MHD) FF characteristics were studied by Mokaddes Ali et al. (2021) in wavy-walled cavities that held a heat-generating cylinder and a Cu-Al₂O₃/water hybrid nanofluid. The study was based of the finite element approach for the investigation of the effects of several parameters, such as the cylinder radius, heater length, hybrid nanoparticle volume fraction, Ra, and Ha, on the flow and thermal behavior. The results showed that high Hartmann number causes the decrease of the fluid motion because of the generated Lorentz force, but high Rayleigh number causes the increase of the flow intensity and the rate of HT. Furthermore, the study shows that at lower cylinder radii, shorter heater lengths, and greater hybrid nanoparticle volume fractions all lead to the enhancement of the rate of HT. Dey and Rathish Kumar (2024) conducted a finite element analysis to investigate the NC in a complex porous enclosure filled with a HNF under thermal non-equilibrium assumption. The study corresponds to an investigation of the effects of Ra, Pr, Da, porosity, and φ , on the HT and FF. The results indicated that the dispersion of the Cu-Al₂O₃ hybrid nanoparticles, in the base fluid has a benefic effect on the improvement if the HT and FF in the porous media. The LBM was used by Cheng et al. (2023) to investigate the HNF MHD NC and irreversibilities production. The aim of the research was the examination of the impact of Ha, baffle sizes, and magnetic field inclination on the generation of thermal and frictional entropy in a square cavity. The results indicated that the velocity components and temperature gradient are reduced with the increase of Hartmann number, which results in a reduction in all forms of irreversibility. In contrast, the thermal, frictional, and total entropy are augmented by an increase in the inclination angle. Roy (2022) presented a study on the MHD NC HT and FF of a HNF in a cavity equipped with multiple heat sources. The authors used the finite difference method to solve the dimensionless system of equations. The numerical model was validated by comparing with experimental and numerical results. It was found that Rayleigh number, the number and size of heat sources have a substantial influence on the flow pattern. The NC of a Cu-Al₂O₃-water HNF in a square cavity equipped with two heat-generating structures was investigated by Hidki et al. (2023). The authors investigated the impact of Ra, φ , and the presence of fins on the blocks, on the temperature field and flow structures using the finite volume method. The results showed that the maximal temperature reduces with the increase of Ras, which is a result of the improvement of the heat exchange. Furthermore, the temperature decreases by up to 18% when the hybrid nanoparticles are incorporated in the base fluid. It was also indicated that the maximal temperature decreases by up to 12% for Ra ≥ 105 when the fins are installed on the surfaces of the blocks. The NC of a non-isothermal HNF has been computationally analysed by Iftikhar et al. (2023). The flow was influenced by a hot fin and non-linear radiation. The FEM was used for the investigation of the effect of the fin height, nanoparticle volume fraction (Ag and Cu), and non-linear thermal radiation on the HT and entropy production. It was found that the NC flow is improved when the height of the hot fins is increased, which results in a considerable increase of the flow intensity and HT rate of up to 42% and 57%, respectively. Furthermore, the addition of NPs leads to a further enhancement of the heat transfer, while the entropy generation exhibits an inverse trend as the nanoparticle concentration is increased. Islam et al. (2023) studied the magnetohydrodynamic (MHD) NC in a prismatic-shaped enclosure containing Cu-TiO₂/water HNFs. The authors examined the impact of Ra, Ha, and φ , on HT and FF. It was indicated that an enhancement of the HT performances occurs using hybrid NPs especially at high Ra values, whereas the reverse effect is observed when the magnetic influence is increased. In addition to the numerical study, a sensitivity analysis was conducted using response surface methodology.

Based on the above presented literature review it can be concluded that the studies related to 3D natural convection of hybrid nanofluids in complex geometries are very limited. The 3D convective flow in hybrid nanofluid-filled bi-truncated-pyramid cavities equipped with adiabatic cylinders represents a critical advancement in the application of hybrid nanofluids for the enhancement of the thermal management. The aim of this paper is to investigate the complex interactions of natural convection within such geometrically intricate enclosures, taking in account the role of adiabatic cylinders in disrupting the FF and the creation of secondary vortices that cause the augmentation of the heat transfer. This numerical investigation consists of analyzing the effects of varying Rayleigh numbers, nanoparticle volume fractions, and cylinder size on the thermal and fluid dynamic behaviors within the cavity. This approach highlights the enhancements of the thermal properties of the hybrid nanofluids and provides a detailed understanding of the effects of geometric modifications on the optimization of the heat transfer in practical applications.

The configuration considered in this investigation (Figure 1) consist of a bi-truncated pyramid filled with a hybrid nanofluid and containing two adiabatic cylindrical obstacles (with a diameter of “D”). The left walls are maintained at a higher temperature (T_h) and the right walls are maintained at a lower temperature (T_c). The heat transmission across these boundaries is prevented by the inclusion of the adiabatic cylinders (with a diameter of “D”) and the adiabatic pyramid’s top and bottom walls.

The dimensional PDEs Equations 1–3 governing the 3D NC are as follows (Almeshaal et al., 2020): ∇ . V ′ → = 0 (1) ∂ V ′ → ∂ t ′ + V ′ .   → ∇ → V ′ → = − 1 ρ nf ∇ → P ′ + v ∆ V ′ → − β t T ′ − T ′ c g → (2) ∂ T ′ ∂ t ′ + V ′ →   . ∇ T ′ = α nf ∆ T ′ (3)

Using the following dimensionless variables: t = t ′ L ′ ² α   ;   x , y , z = x , y , z ′ L ′   ;   V x , V y , V z = V x , V y , V z ′ α L ′   ;   P = P ′ ρ   α L ′   ;  and  T = T ′ − T c ′     /   T h ′ − T c ′

The governing equations are converted to the dimensionless form: ∇ . V → = 0 (4) ∂ V → / ∂ t + V → .   ∇ →   V → = − ∇ → P + Pr   V nf   / V f   ∆ V → + β nf   / β f    Ra  Pr   T g → (5) ∂ T / ∂ t + V → . ∇ T = α nf α f ∇ ² T (6)

With: Pr = υ f α f   ;  and Ra = g .   β f   . ∆ T   . l 3 υ f . α f

The dimensionless Boundary conditions are: T = 1 , at the left walls T = 0 , at the right walls ∂ T ∂ n = 0  at remaining walls V x = V y = V z = 0 , at all walls

The local and average Nusselt numbers are Equations 7, 8: N u l = k nf k f   ∂ T ∂ n (7)

And Nu a = ∫ 0 1 ∫ 0 1 Nu l   . dy . dz (8)

The following expressions Equations 9–15 are used to evaluate the thermophysical properties (Almeshaal et al., 2020; Gal et al., 2022):

Density: ρ n f = 1 − φ ρ f + φ ρ s (9)

Heat capacitance: ρ C p n f = 1 − φ ρ C p f + φ ρ C p s (10)

Thermal conductivity.

- For CNT-nanofluid: k n f k f = 1 − ϕ + 2 φ k s k s − k f ln k s + k f 2 k f 1 − φ + 2 φ k f k s − k f ln k s + k f 2 k f (11)

Dynamic viscosity:

- For CNT-nanofluid: µ nf = µ f 1 + 13.5 . φ + 904.4 . φ 2 (13)

The properties of HNF are evaluated as follows: y h n f = y C N T . f r + y A l 2 O 3 . 1 − f r (15) where y is any thermophysical property and f r is the mass fraction of CNT fixed at 0.15.

The properties of water and NPs are presented in Table 1.

The FEM, based on the Galerkin’s Weighted Residual technique, is used to solve the set of Equations 4–6 along with the applied boundary conditions.

The Lagrange polynomials are expressed as: H = ∑ j = 1 N s δ j s K j where δ is the shape function and K represent the values of the variables at the nodes of the elements. The convergence criterion is fixed at f n + 1 − f n f n + 1 ≤ 1 0 − 6 .

This study aims to explore the enhancement of NC HT using HNF within a bi-truncated pyramid equipped with adiabatic cylinders. The purpose is to understand how varying Ra, φ , and D influence the thermal and fluid dynamics in complex geometries. The methodology involves a numerical simulation approach using the FEM to solve the governing equations of FF and HT. The study systematically varies the Ra from 10³ to 10⁶, φ from 0 to 0.045, and the cylinder diameter (D) from 0.05 to 0.24.

Figure 3 shows the flow patterns of the 3D convective flow in the hybrid nanofluid within the bi-truncated pyramid enclosure equipped with two internal adiabatic cylinders, for various Ra values (10³–10⁶). The left walls of the bi-truncated pyramid are hot, and the right walls are cold. This temperature differential is the driving force for the convective flow within the pyramid. As can be remarked from the figure’s legends, higher Ra numbers correspond to higher velocities within the fluid due to the increase of the buoyancy force with Ra. At Ra = 10³, the flow is more orderly with lower velocities. As Ra increases, higher velocities and more complex flow patterns, occur with the apparition of secondary flows between the cylinders. The most complex and intense flow structure occurs at Ra = 10⁶, where the highest velocities are reached, and additional secondary flow appear close to the cylinders. The figure illustrates the transition from orderly to complex convective flow as Ra increases from 10³ to 10⁶. Initially, at Ra = 10³, the HT is domination by conduction mode, and is characterized by smooth and coherent flow structures that indicate less intense convection due to the lower thermal driving force. As Ra increases, the impact of buoyancy forces overcomes the viscous damping and conduction, and the flow becomes more complex, signifying enhanced heat transfer. The presence of the two adiabatic cylinders, adds to the complexity by disrupting the flow, creating additional vortices and flow patterns around them.

Figure 4 displays the streamlines at the central plane of the bi-truncated pyramid, comparing the flow patterns without nanoparticles ( φ = 0) and with nanoparticles ( φ = 0.045), at various Ra (10⁴, 10⁵, and 10⁶). Streamlines are presented for a better understanding of the flow of a hybrid nanofluid in the central cross-section of the pyramid. The presence of nanoparticles enhances the fluid’s thermal conductivity and influence the flow dynamics. At Ra = 10⁴, the flow is relatively moderate. Without nanoparticles ( φ = 0), the streamlines are orderly and symmetrical, showing a stable convective pattern. With nanoparticles ( φ = 0.045), the streamlines appear slightly more concentrated and exhibit stronger convection near the hot and cold walls, indicating that the nanoparticles are enhancing the convective motions. As Ra increases (10⁵), the effect of buoyancy becomes stronger, and we observe more pronounced convective vortexes. The presence of nanoparticles again amplifies the velocity, suggesting stronger convective flow and more efficient heat transfer. At the highest Ra (Ra = 10⁶), the flow pattern shows highly concentrated areas of streamlines, especially near the walls, indicating more intense convective motions. The addition of nanoparticles causes even tighter streamline patterns and larger gradients in the velocity field, suggesting an even stronger convective flow, which is indicative of enhanced heat transfer efficiency. The velocity fields for φ = 0.045 are more intense than that for φ = 0, this due to the modified thermal and flow properties of the fluid with nanoparticles. From this figure it can be concluded that, Without nanoparticles ( φ = 0), the flow is smooth, and the streamlines are well-ordered, suggesting a regular pattern of convection with moderate velocities. When nanoparticles are introduced ( φ = 0.045), the streamlines tighten and the flow velocities increase, particularly near the thermal boundaries, highlighting the nanoparticles’ role in enhancing thermal conductivity and the resultant convective HT. As the Rayleigh number escalates from 10⁴ to 10⁶, the streamlines in both cases become more compressed, indicating the strengthening in the buoyancy-driven flow. The nanoparticles amplify this effect, as the flow with φ = 0.045 consistently shows more intense streamlines and greater velocities. This effect is more pronounced at higher Ra, where the thermal driving force is greater, suggesting that nanoparticles could be instrumental in optimizing HT in engineering applications where controlling the flow intensity and heat distribution is crucial.

Figure 5 presents the 3D temperature field (left) and isotherms (right) at the central plane of the bi-truncated pyramid equipped with two internal adiabatic cylinders, for various Ra, and with a contrast between the pure fluid ( φ = 0) and the HNF ( φ = 0.045). At the low Ra (Ra = 10⁴), the temperature field is relatively undisturbed, indicating a weaker convective flow. The isotherms on the right-side show that with nanoparticles (colored iso-lines), the temperature gradient is more pronounced near the cylinders, suggesting that the nanoparticles enhance thermal mixing. At intermediate Ra (Ra = 10⁵), the temperature field becomes more varied, and the convective patterns are more evident, particularly in the presence of nanoparticles. The isotherms for φ = 0.045 are more tightly packed near the hot and cold boundaries, which points to a stronger temperature gradient due to the improved HT properties of the HNF. At the highest Ra (Ra = 10⁶), the temperature gradients are extremely sharp, and the influence of convection is dominant. For φ = 0.045, the isotherms depict a very concentrated thermal field near the hot surface, implying an intense heat transfer regime, likely resulting from the combined effects of increased thermal conductivity and convective enhancement provided by the nanoparticles. In addition, a vertical stratification of the isotherm appears and is more pronounced for the NF compared to the pure water. This figure put in evidence the critical role that nanoparticles play in manipulating the thermal characteristics of a fluid. The presence of nanoparticles tends to sharpen temperature gradients and enhance heat transfer, especially as the convective forces become stronger with higher Ra. The internal adiabatic cylinders serve to complicate the flow and temperature fields, acting as obstacles that induce local variations in heat transfer and fluid movement. From this figure, it can be concluded that at the lower Ra, the temperature fields and isotherms are relatively uniform, indicating moderate convective heat transfer. As Ra increases, the introduction of nanoparticles ( φ = 0.045) considerably alters the temperature distribution, as seen by the compressed isotherms, particularly around the hot and cold boundaries, signifying intensified thermal gradients and stronger convective currents. At the highest Ra, the impact of nanoparticles is the most pronounced, with the thermal field exhibiting tight gradients and complex patterns, indicative of a robust convection-driven heat transfer process. This thermal behavior suggests that nanoparticles effectively promote the heat transfer, which is a desired outcome in applications aiming for efficient thermal management and enhanced heat dissipation.

Figure 6 presents the flow structure and temperature field at the central plane of the bi-truncated pyramid at φ = 0.045, Ra = 10⁶ and various diameters (D) of the internal adiabatic cylinders. The streamlines represent the flow patterns within the central plane of the pyramid. As the diameter of the adiabatic cylinders increases from D = 0.05 to D = 0.24, the streamlines indicate changes in flow behavior around the cylinders. With smaller diameters, the flow can move around the cylinders more freely, which leads to a more uniform flow field. When the cylinder size is increased, the obstruction to the flow becomes more significant, which leads to a more disturbed flow patterns and to the apparition of secondary vortexes. The temperature distribution is presented to show how the heat is transferred from the hot to the cold walls of the pyramid. For small cylinder diameters (D = 0.05), a smooth transitions of the temperature gradients across the entire domain, occurs, this is an indication of the improvement of the efficient HT from the hot to the cold side. When the cylinder diameter is increased, the temperature field becomes more uneven, particularly around the cylinders. This is an indication that the larger cylinders imped the heat flow, creates cooler regions behind them, and potentially reduces the overall heat transfer efficiency. The figure gives an illustration of how the variation of the diameter of internal adiabatic cylinders affects both the convective flow patterns and the thermal fields. In fact, for the smallest cylinders (D = 0.05), the flow is less disrupted, and leads to a more homogeneous pattern, corresponding to the smoother temperature gradients, indicating that the HT is more efficient. As cylinder diameter increases, the flow encounters greater resistance, generating complex patterns that suggest enhanced fluid velocity, yet also revealing more pronounced thermal stratification with cooler regions behind the cylinders. This indicates a trade-off: larger cylinders augment fluid intensity but also create areas of thermal stagnation, potentially impacting the system’s overall HT efficacy.

Figure 7illustrates the effect of φ on Nu_a at the left hot walls, for D = 0.1, across a range of Ra. When Ra increases, there is a consistent rise in the Nu_a, indicating an enhancement in convective HT caused by the intensification of the buoyancy forces. This effect is furthermore intense for higher nanoparticles volume fractions; particularly at Ra = 10⁶; in fact, the slope of the curve is steeper, indicating better heat transfer improvement capabilities. At low Ra values, Nu_a slightly varies with varying φ , this is due to the dominance of the conductive HT regime where nanoparticles addition has minimal impact. This behavior changes at higher Ra values, where the introduction of NPs disrupts the thermal boundary layer more effectively, enhancing the fluid’s thermal conductivity and convective HT efficiency.

Figure 8 illustrates the impact of increasing the diameter of internal adiabatic obstacles (D) on Nu_a at φ = 0.045, and different Ra values. As the diameter of the obstacles increases, there is a general decrease in Nu_a for all Ra values. This trend is more pronounced at higher Ra values. Specifically, at Ra = 10⁶, Nu_a exhibits a sharp decline as D increases, indicating that larger obstacles significantly reduce the fluid flow, thereby reducing heat transfer efficiency. In contrast, at lower Ra values, the decrease in Nu_a is less steep, thus it can be concluded that the influence of the adiabatic obstacles on the HT is less critical under weaker convection. This pattern highlights the nuanced relationship between obstacle size and convective heat transfer: larger obstacles disrupt the fluid path more extensively, creating regions of reduced fluid movement and increased thermal resistance. Conversely, smaller obstacles allow for more effective fluid navigation around them, maintaining better convective heat transfer.

The current study investigates the effectiveness of HNFs in enhancing natural convection HT within a bi-truncated pyramid. The simulations demonstrate that increasing Ra and φ significantly improve the HT efficiency, while the adiabatic cylinders size has an opposite effect by damping the flow and reducing convection. These findings highlight the potential of hybrid nanofluids to optimize thermal management systems, particularly in applications requiring precise control of heat transfer processes. Future research should explore the influence of different nanoparticle combinations and geometric configurations to further refine these enhancements and broaden their applicability in various engineering sectors. The main Results can be summarized as follows:

• The Nusselt number increased by up to 23% at Ra = 10⁶ and φ = 0.045.

Future research will focus on optimizing heat transfer in 3D geometries by exploring the effects of varying cylinder shapes, positions, and materials, as well as incorporating advanced hybrid nanofluid formulations with enhanced thermophysical properties. Additionally, the integration of external influences such as magnetic fields, vibration effects, or transient thermal boundary conditions.