Next Article in Journal
A Review of Nanofluids as Coolants for Thermal Management Systems in Fuel Cell Vehicles
Next Article in Special Issue
Enhanced Catalytic Performance of Ag NP/0.95AgNbO3-0.05LiTaO3 Heterojunction from the Combination of Surface Plasma Resonance Effect and Piezoelectric Effect Using Facile Mechanical Milling
Previous Article in Journal
Improved Resistive and Synaptic Characteristics in Neuromorphic Systems Achieved Using the Double-Forming Process
Previous Article in Special Issue
A Bilayered Wood-Poly(3,4-ethylenedioxythiophene):Polystyrene Sulfonate Hydrogel Interfacial Evaporator for Sustainable Solar-Driven Sewage Purification and Desalination
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Free Convection in a Square Ternary Hybrid Nanoliquid Chamber with Linearly Heating Adjacent Walls

1
Department of Mathematics, GITAM (Deemed to Be University), Hyderabad Campus, Hyderabad 502329, Telangana, India
2
Laboratory on Convective Heat and Mass Transfer, Tomsk State University, 634050 Tomsk, Russia
*
Author to whom correspondence should be addressed.
Nanomaterials 2023, 13(21), 2860; https://doi.org/10.3390/nano13212860
Submission received: 8 September 2023 / Revised: 16 October 2023 / Accepted: 26 October 2023 / Published: 28 October 2023

Abstract

:
In this study, mathematical modeling of the energy transfer and flow characteristics of ternary nanoliquid in a square enclosure is performed. In the cavity considered, the left and bottom borders are warmed uniformly or non-uniformly when the rest of the borders are cooled. The robust finite element method with quads and triangles as elements is used to work out the control equations of the problem. The current study is validated against previously published works, and good agreement is shown. The isolines are investigated for various Rayleigh numbers at uniform and non-uniform thermal boundary conditions. The impact of ternary hybrid nanofluids on the mean Nusselt number at hot borders is explored in dependence on the Rayleigh number and nanoparticle concentration. A comparative study of different fluids for the mean Nusselt number at heated borders is also conducted and analyzed with appropriate graphs and tables. It has been shown that ternary nanofluids can be more effective compared to mono- and hybrid nanofluids, with a more essential growth of the energy transport rate with nanoadditives concentration.

1. Introduction

The phenomenon of buoyancy-induced energy transport in enclosures has been the subject of a huge amount of study in the past, and it is an essential component for a diverse range of technological purposes. This is the case, for example, with the cooling of electronic systems, where a huge number of high-power dissipation units are frequently packaged in cabinets to minimize the amount of space required as well as the number of external cooling sources. In addition, heat transfer engineers often choose not to make use of mechanical devices for the circulation of the coolant because of the high amount of power it consumes, the excessive working noise it creates, or worries regarding the system’s reliability. On the other hand, the efficiency of thermal energy removal from equipment whose temperature management depends on thermal convection is typically severely limited by the poor thermal conductivity of standard thermal fluids like water, ethylene glycol, and mineral oils. Nanofluids, or liquid suspensions containing tiny solid particles, offer a potential solution to this issue because of their effective thermal conductivity, which is considerably greater when compared to the pure host fluid.
In their crucial work, Mahmoudi et al. [1] performed a meticulous simulation of the intricate phenomenon of free convection cooling. Specifically, they focused on the captivating scenario of two heat sources that were affixed in a vertical orientation to the horizontal walls of a confined enclosure. The study’s findings revealed a significant correlation between the Rayleigh numbers and the heater’s location with respect to the flow field and thermal patterns in the chamber. The mathematical modeling performed by Cho et al. [2] explored the thermal convection parameters of Al2O3–H2O nanoliquid within an enclosure. The cavity was defined by vertical isothermal borders featuring a complicated wavy sheet, while the horizontal borders remained straight and subjected to adiabatic conditions. The findings elucidated in this investigation offer valuable perspectives on potential methodologies to augment the convective thermal energy transport efficiency inside the enclosure featuring intricate, undulating border surfaces. Abu-Nada and Chamkha [3] conducted a computational investigation of the steady mixed convection phenomenon in a lid-driven chamber. The working fluid considered was a nanosuspension with CuO nanoadditives dispersed in water. The bottom surface of the chamber was characterized by a wavy pattern. The enclosure was thermally isolated along its left and right boundaries, while the undulating horizontal boundaries were held at fixed temperatures. The uppermost border was identified as the thermally active boundary, functioning as the heated enclosure boundary while simultaneously exhibiting a uniform and unchanging velocity. As a result, it has been observed that the inclusion of nanoparticles induces a notable enhancement in thermal energy transport for various Richardson numbers and ratios of bottom wall geometry. The numerical solution presented by Sheremet et al. [4] pertains to the MHD thermal convection of a Cu–H2O nanoliquid within an irregular, porous, semi-open chamber. The chamber contains a heat source that maintains a constant temperature. The solution was obtained by employing the finite difference technique. The researchers have observed an enhancement in energy transfer with a rise in the Rayleigh number and a decrease in thermal energy transport with a rise in the Hartmann number. The numerical investigation conducted by Mondal and Mahapatra [5] focuses on magnetohydrodynamics (MHD) mixed convection flow and mass transport of Al2O3–H2O nanoliquid within a trapezoidal chamber. The study specifically examines the presence of discrete energy and mass sources placed at the lower border of the chamber. The present investigation has observed that the combined convective flow under a weak magnetic effect and within a chamber of low aspect ratio and rectangular shape consistently yields favorable results in terms of minimizing total entropy production. Furthermore, it has been determined that the central region of the bottom border of the chamber is particularly well-suited for the introduction of thermal energy and mass, ensuring a uniform distribution throughout the entirety of the cavity. In their numerical investigation, Reddy and Sreedevi [6] examined the impact of radiation on the transport of thermal energy and mass in a nanoliquid within a square cavity. This analysis was conducted by incorporating Buongiorno’s mathematical model. Thermophoresis and Brownian diffusion are employed for analysis. In this investigation, it has been observed that the strength of the thermal energy transfer rises with an increase in the Rayleigh number. A few other significant contributions using nanofluids to enhance heat transfer can be found in Xue et al. [7], Wang et al. [8], Qi et al. [9], Wang et al. [10], and Ge et al. [11]. Additionally, important investigations related to the cavity flow can be found in Chen et al. [12] and Luo and Yang [13].
Hybrid nanoliquids represent a novel class of working liquids that have been meticulously designed to possess augmented physical characteristics. The hybrid nanoliquids leverage the physical characteristics inherent in multiple varieties of nanoparticles. Recently, researchers have employed a hybrid nanofluid, consisting of a combination of disparate nanoparticles dispersed within a host fluid, as the operational medium to further augment the heat transfer properties. Henceforth, Mehryan et al. [14] undertook an investigation on the phenomenon of free convective thermal energy transport pertaining to the hybrid nanoliquid composed of Al2O3–Cu–H2O within a chamber that has been occupied by a porous medium. The experimental findings demonstrated a decrease in thermal conductivity through the implementation of nanoparticles within a permeable medium. The observed diminution of the thermal energy transport strength exhibits a significantly greater magnitude for the hybrid nanoliquid compared to the mono-nanoliquid. Alsabery et al. [15] conducted an investigation on the conjugate free convection and entropy production of Cu–Al2O3/H2O hybrid nanoliquids. The study was carried out inside a partially divided chamber featuring wavy and cold sidewalls. The experimental results have demonstrated that the utilization of hybrid nanofluids results in a growth in the Nusselt number in comparison with mono-nanofluids. Tayebi and Chamkha [16] employed a hybrid nanofluid to investigate the entropy production resulting from MHD thermal convection within a square chamber. The study considered the inclusion of an irregular circular conductive cylinder, accounting for the presence of irreversibilities arising from the magnetic effect. The model proposed by Khan et al. [17] examines heat management and thermal enhancement within a split lid-driven square chamber. In this investigation, a Y-shaped obstruction is positioned within the enclosure, which is saturated with hybrid nanoliquids composed of Al2O3, Cu, and water. The observation is made that a Y-shaped obstruction results in an augmentation in the strength of thermal energy transport within the enclosed space. Belhaj and Ben-Beya [18] conducted a study with the objective of analyzing the entropy production and MHD thermal convection of a hybrid nanoliquid in a chamber. The enclosure contained a warmed elliptical body positioned at its center, and a periodic-variable magnetic influence was present during the investigation. Several notable investigations on mono- and hybrid nanoliquids have been performed by Rajesh et al. [19], Sarkar et al. [20], Rajesh et al. [21], Sidik et al. [22], Dhinesh Kumar and Valan Arasu [23], Vemula et al. [24], and Sheremet et al. [25].
The investigation of thermal energy transfer in tri-hybrid nanofluids, wherein three distinct nanoparticles are dispersed within a base fluid, has garnered significant attention in recent studies. The process of synthesizing nanofluids is undertaken with the objective of augmenting the physical attributes inherent to said nanofluids. This methodology facilitates the generation of a heat carrier exhibiting augmented thermal and rheological properties. The experimental outcomes have shown that the heat conductivity of tri-hybrid nanofluid surpasses that of both mono-phased and double-phased nanofluids. In the investigation conducted by Sahoo [26], a novel coolant for radiators was examined. This coolant consisted of a ternary hybrid nanofluid containing nanoparticles of different shapes, specifically spherical (Al2O3), cylindrical (CNT), and platelet (graphene). The influence on heat transfer performance combined with the exergetic study, resulting from variations in coolant flow strength, the volume fraction of nanoadditives, and air velocity, has led to the conclusion that alterations in nanoadditives concentrations have an essential influence on heat transport performance, primarily owing to the morphology of the nanoadditives. Zahan et al. [27] conducted research on a ternary nanosuspension contained in distilled water and ethylene glycol, incorporating Co, Ag, and Zn nanoadditives at varying solid volume fractions. The objective was to examine the liquid dynamics and heat behavior within a CD nozzle. In this investigation, the enhanced strength of heat transport is determined by employing laminar-shaped nanoadditives. The investigation conducted by Shao et al. [28] involved the study of thermal convection in a prismatic chamber containing two mobile hot obstacles and porous material. The enclosure was saturated with a ternary hybrid nanoliquid. The findings suggest that a rise in the Darcy number characterizes a significant enhancement of both vertical and horizontal velocities, with improvements of 88.51% and 83.77%, respectively. Additionally, there is a substantial increase of 88.12% in the streamline patterns, as well as notable enhancements in Nusselt numbers and isotherms. The article by Sneha et al. [29] elucidates the dynamics of an unsteady flow involving a couple of stresses through the implementation of a ternary hybrid nanosuspension on a stretching sheet accompanied by porous material. Few other significant works using ternary hybrid nanofluids are presented in papers by Hafeez et al. [30], Adnan and Ashraf [31], Al Oweidi et al. [32], Puneeth et al. [33], and Mahmood and Khan [34].
Based on the comprehensive analysis of the existing literature, it becomes apparent that numerous investigations have been conducted thus far pertaining to the phenomenon of natural convection occurring within an enclosed space, employing both conventional fluids, nanofluids, and hybrid nanofluids. However, research on the thermal convection phenomena occurring within a confined space utilizing a ternary hybrid nanofluid has been relatively limited. In light of heat transfer applications in the realm of engineering and technology, the current investigation seeks to unveil the fundamental influence of a ternary hybrid nanofluid on the circulation and thermal energy transport properties inside a square chamber. Specifically, the left and bottom borders of the cavity are subjected to uniform or linear heating, while the rest are subjected to cooling. The analysis of the Nusselt number at the left and lower walls, combined with the variations in velocity and temperature within the cavity for different fluids, is conducted using FEM in COMSOL Multiphysics software 6.1. This investigation includes the influence of important governing factors such as the Rayleigh number and the nanoadditives concentration. The results are shown through the utilization of streamlines, isotherms, and other appropriate graphical representations and tabulated data.

2. Mathematical Analysis

In our analysis, we shall consider a square chamber with a characteristic size of L, as depicted in Figure 1. When the right and top borders are subjected to cooling, the left and bottom walls experience either uniform or non-uniform heating.
The working liquid employed in this research is the ternary hybrid nanofluid, specifically composed of Cu–CuO–Al2O3 nanoparticles dispersed in water. It is an incompressible Newtonian medium with constant characteristics. Moreover, the present investigation neglects the influence of radiation and viscous dissipation while assuming the absence of internal heat generation/absorption. The problem has been analyzed using the two-dimensional Navier-Stokes equations and energy equation, considering steady and laminar flow conditions. The density variations in the nanosuspension are studied via the Boussinesq approach. It is postulated that the fundamental liquid, namely, water, and the nanoadditives are in a state of heat equilibrium. The control equations, given the aforementioned assumptions, can be expressed as
v 1 x + v 2 y = 0
ρ t h n f v 1 v 1 x + v 2 v 1 y = p x + μ t h n f 2 v 1 x 2 + 2 v 1 y 2
ρ t h n f v 1 v 2 x + v 2 v 2 y = p y + μ t h n f 2 v 2 x 2 + 2 v 2 y 2 + ρ β t h n f g T * T c
v 1 T * x + v 2 T * y = κ t h n f ρ c p t h n f 2 T * x 2 + 2 T * y 2
with boundary conditions
v 1 0 , y = 0 , v 1 L , y = 0 , v 1 x , 0 = 0 , v 1 x , L = 0 , v 2 0 , y = 0 , v 2 L , y = 0 , v 2 x , 0 = 0 , v 2 x , L = 0 , T * x , 0 = T h T h T c m x L , T * x , L = T c , T * 0 , y = T h T h T c m y L , T * L , y = T c
The physical features of host liquid, nanoadditives, and ternary hybrid nanosuspension are provided in Table 1 and Table 2.
Employing the mentioned transformations
X = x L , Y = y L , V 1 = v 1 L α f , V 2 = v 2 L α f , P = p L 2 ρ f α f 2 , T = T * T c T h T c
into Equations (1)–(4) above, we obtain
V 1 X + V 2 Y = 0
ρ t h n f ρ f V 1 V 1 X + V 2 V 1 Y = P X + P r μ t h n f μ f 2 V 1 X 2 + 2 V 1 Y 2
ρ t h n f ρ f V 1 V 2 X + V 2 V 2 Y = P Y + P r μ t h n f μ f 2 V 2 X 2 + 2 V 2 Y 2 + R a   P r ρ β t h n f ρ β f T
ρ c p t h n f ρ c p f V 1 T X + V 2 T Y = κ t h n f κ f 2 T X 2 + 2 T Y 2
with conditions
V 1 0 , Y = 0 , V 1 1 , Y = 0 , V 1 X , 0 = 0 , V 1 X , 1 = 0 , V 2 0 , Y = 0 , V 2 1 , Y = 0 , V 2 X , 0 = 0 , V 2 X , 1 = 0 , T X , 0 = 1 m X , T X , 1 = 0 , T 0 , Y = 1 m Y , T 1 , Y = 0
Here R a = g β f T h T c L 3 υ f α f is the Rayleigh number and P r = υ f α f is the Prandtl number.
The paramount parameter of technical significance in this investigation is the mean Nusselt number, as expressed by
N u L = κ t h n f κ f 2 2 m 0 1 T X X = 0 d Y is the Nusselt number at the left border and
N u B = κ t h n f κ f 2 2 m 0 1 T Y Y = 0 d X is the Nusselt number at the bottom border.

3. Numerical Solution

The PDEs (7)–(10) with conditions (11) modeling the problem have been unraveled employing the finite element method with COMSOL Multiphysics. The extra-fine mesh containing 16,946 elements is selected as the optimum mesh based on a grid-independent test, the details of which are provided in Table 3a,b. The mesh statistics of the selected effective mesh type (extra fine) are provided in Table 4. The computational mesh chosen for the present research is shown in Figure 2. Further, to validate the current model, the mean Nusselt numbers of this study for different Ra when Pr = 0.71, b1 = 0, b2 = 0, and b3 = 0 were compared to the findings in [36,37] with the same boundary conditions in Figure 3 and found excellent agreement between the results. Figure 4 shows a good agreement with data of Sathiyamoorthy and Chamkha [38]. This further confirms the accuracy of the current numerical results.

4. Results and Discussion

In the current investigation, we maintained Ra between 103 and 105, Pr at a value of 6.2, and the volume fractions of Al2O3, CuO, and Cu nanoadditives (b1, b2, and b3, respectively) at 0.04 each. Additionally, we considered two cases: m = 0 for constant/uniform heating and m = 1 for linear/non-uniform heating. These values were held constant throughout the study, unless explicitly stated otherwise for the specific parameter under examination.
Figure 5 and Figure 6 present isolines for various Rayleigh numbers at uniform and non-uniform heating of the left and lower borders. For the non-uniform heating (see Figure 5), when the temperature on the left and lower borders is changed owing to the linear function of the coordinate, one can see an appearance of global clockwise flow and heating of the chamber from these walls. For a low Rayleigh number (Ra = 103), a weak global eddy appeared inside the chamber, with upstream flows near the heated left surface and downstream near the cold right surface. The temperature field illustrates the heating of the cavity from the bottom left corner, where the maximum temperature is fixed. Heating occurs like a heat conduction mechanism due to the low value of the Rayleigh number. A raise of Ra (=104) results in a strengthening of convective motion and thermal energy transport. As a result, a weak deformation of the convective cell core can be seen with the appearance of the thermal plume along the left vertical surface. Further rise of Ra (=105) characterizes a reduction in the heat plume thickness with more significant warming of the chamber, while the major convective cell is degraded with a formation of small recirculation within the upper left corner.
For the uniform heating (see Figure 6), when the temperature on the left and bottom walls is high and constant, one can also reveal an appearance of global clockwise motion and heating of the chamber from these walls. For Ra = 103, the major thermal energy transport mechanism is heat conduction, and as a result, the isotherms do not reflect the appearance of a heat plume. The growth of the temperature difference with a rise of Ra (=104) illustrates significant liquid motion, and isotherms show the development of ascending and descending thermal plumes along vertical walls. For Ra = 105, a global circulation has an elongated convective cell core, and isotherms within the cavity reflect the appearance of the thermal stratification core.
An impact of Ra and working fluid structure on NuL and NuB is shown in Figure 7 for non-uniform heating. These distributions characterize the thermal energy transport augmentation with a transition between the host liquid and ternary hybrid nanofluid through alumina/water nanoliquid and copper oxide–alumina/water hybrid nanoliquid. Essential intensification of the convective energy transfer occurs at the bottom surface due to an interaction between the cold descending fluid flow along the right cold border and the heated bottom border. In the case of the left heated border, the mean Nusselt numbers are not as significant as for the bottom heated wall because there is no essential temperature difference here compared to the lower heated border.
For uniform heating (Figure 8), the integral thermal flux from the hot walls is greater than this heat flux for non-uniform wall heating, and as a result, the average Nusselt numbers for m = 0 are greater than for the case of m = 1. The behavior of NuB with the working fluid is the same as in the case of m = 1, but for the left border, some differences are obtained. The reason for these differences is the formation of a cold zero temperature at the top left corner with a small recirculation (see Figure 5c) in the case of non-uniform heating.
More detailed outcomes for the average Nusselt numbers for various hybrid nanofluids are shown in Table 5 and Table 6, where one can find that the maximum thermal energy transport strength can be found for Cu–CuO/water.
Figure 9 demonstrates the impact of the copper nanoadditives concentration on the mean Nusselt number for the non-uniform and uniform heating of ternary hybrid nanofluid at Ra = 1000. An increase in the copper nanoadditives concentration reflects a rise in the thermal energy transfer strength of the ternary hybrid nanoliquid, regardless of the considered cases.

5. Conclusions

Mathematical modeling of thermal convection of mono-, hybrid-, and ternary nanofluids in a differentially heated chamber with uniform and non-uniform heating from the left vertical and lower horizontal walls has been performed for a wide range of the Rayleigh number and nanoparticle volume fraction. It has been shown that non-uniform heating for high values of the Rayleigh number reflects a formation of secondary recirculation in the upper left corner where the cold temperature is kept. A rise in the Rayleigh number and nanoparticle concentration illustrates the thermal energy transfer enhancement for the ternary nanoliquid, while for the considered hybrid nanofluids, the more effective is Cu–CuO/water. Thus, a growth of the Rayleigh number between 103 and 105 allows increasing the average Nusselt number at the bottom wall by about 76% for the ternary nanoliquid in the case of uniform heating, while in the case of non-uniform heating, one can find a growth of about 4.5 times. A rise in the nanoparticle volume fraction for the ternary nanoliquid for b3 ranges between 0.0 and 0.04, leading to a growth of the average Nusselt number at the bottom wall at about 11% in the case of uniform heating and at about 10% in the case of non-uniform heating. The performed analysis shows that ternary nanofluids can be effective “smart” working fluids that allow intensification of the convective energy transference within the closed chambers.

Author Contributions

Conceptualization, V.R. and M.S.; software, V.R.; validation, V.R.; investigation, V.R. and M.S.; writing—original draft preparation, V.R. and M.S.; writing—review and editing, V.R. and M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mahmoudi, A.H.; Shahi, M.; Shahedin, A.M.; Hemati, N. Numerical modeling of natural convection in an open cavity with two vertical thin heat sources subjected to a nanofluid. Int. Commun. Heat Mass Transf. 2011, 38, 110–118. [Google Scholar] [CrossRef]
  2. Cho, C.C.; Chen, C.L.; Chen, C.K. Natural convection heat transfer performance in complex-wavy-wall enclosed cavity filled with nanofluid. Int. J. Therm. Sci. 2012, 60, 255–263. [Google Scholar] [CrossRef]
  3. Abu-Nada, E.; Chamkha, A.J. Mixed convection flow of a nanofluid in a lid-driven cavity with a wavy wall. Int. Commun. Heat Mass Transf. 2014, 57, 36–47. [Google Scholar] [CrossRef]
  4. Sheremet, M.A.; Oztop, H.F.; Pop, I.; Al-Salem, K. MHD free convection in a wavy open porous tall cavity filled with nanofluids under an effect of corner heater. Int. J. Heat Mass Transf. 2016, 103, 955–964. [Google Scholar] [CrossRef]
  5. Mondal, P.; Mahapatra, T.R. MHD double-diffusive mixed convection and entropy generation of nanofluid in a trapezoidal cavity. Int. J. Mech. Sci. 2021, 15, 208. [Google Scholar] [CrossRef]
  6. Reddy, P.S.; Sreedevi, P. Buongiorno’s model nanofluid natural convection inside a square cavity with thermal radiation. Chin. J. Phys. 2021, 72, 327–344. [Google Scholar] [CrossRef]
  7. Xue, W.; Wang, Y.; Chen, Z.; Liu, H. An integrated model with stable numerical methods for fractured underground gas storage. J. Clean. Prod. 2023, 393, 136268. [Google Scholar] [CrossRef]
  8. Wang, Y.; Shi, H.; Fang, B.; Zakin, J.; Yu, B. Heat transfer enhancement for drag-reducing surfactant fluid using photo-rheological counterion. Exp. Heat Transf. 2012, 25, 139–150. [Google Scholar] [CrossRef]
  9. Qi, Y.; Weavers, L.K.; Zakin, J.L. Enhancing heat-transfer ability of drag reducing surfactant solutions with ultrasonic energy. J. Non-Newton. Fluid Mech. 2003, 116, 71–93. [Google Scholar] [CrossRef]
  10. Wang, Y.; Wang, Y.; Cheng, Z. Direct numerical simulation of gas-liquid drag-reducing cavity flow by the VOSET method. Polymers 2019, 11, 596. [Google Scholar] [CrossRef]
  11. Ge, W.; Zhang, Y.; Zakin, J.L. Surfactant turbulent drag reduction in an enclosed rotating disk apparatus. Exp. Fluids 2007, 42, 459–469. [Google Scholar] [CrossRef]
  12. Chen, K.T.; Tsai, C.C.; Luo, W.J.; Lu, C.W.; Chen, C.H. Aspect ratio effect on multiple flow solutions in a two-sided parallel motion lid-driven cavity. J. Mech. 2015, 31, 153–160. [Google Scholar] [CrossRef]
  13. Luo, W.J.; Yang, R.J. Multiple fluid flow and heat transfer solutions in a two-sided lid-driven cavity. Int. J. Heat Mass Transf. 2007, 50, 2394–2405. [Google Scholar] [CrossRef]
  14. Mehryan, S.A.M.; Kashkooli, F.M.; Ghalambaz, M.; Chamkha, A.J. Free convection of hybrid Al2O3-Cu water nanofluid in a differentially heated porous cavity. Adv. Powder Technol. 2017, 28, 2295–2305. [Google Scholar] [CrossRef]
  15. Alsabery, A.I.; Hashim, I.; Hajjar, A.; Pour, M.S. Entropy generation and natural convection flow of hybrid nanofluids in a partially divided wavy cavity including solid blocks. Energies 2020, 13, 2942. [Google Scholar] [CrossRef]
  16. Tayebi, T.; Chamkha, A.J. Entropy generation analysis during MHD natural convection flow of hybrid nanofluid in a square cavity containing a corrugated conducting block. Int. J. Numer. Methods Heat Fluid Flow 2020, 30, 1115–1136. [Google Scholar] [CrossRef]
  17. Khan, Z.H.; Khan, W.A.; Hamid, M.; Liu, H. Finite element analysis of hybrid nanofluid flow and heat transfer in a split lid-driven square cavity with Y-shaped obstacle. Phys. Fluids 2020, 32, 093609. [Google Scholar] [CrossRef]
  18. Belhaj, S.; Ben-Beya, B. Thermal performance analysis of hybrid nanofluid natural convection in a square cavity containing an elliptical obstacle under variable magnetic field. Int. J. Numer. Methods Heat Fluid Flow 2022, 32, 1825–1860. [Google Scholar] [CrossRef]
  19. Rajesh, V.; Chamkha, A.J.; Sridevi, C.; Al-Mudhaf, A.F. A numerical investigation of transient MHD free convective flow of a nanofluid over a moving semi-infinite vertical cylinder. Eng. Comput. 2017, 34, 1393–1412. [Google Scholar] [CrossRef]
  20. Sarkar, J.; Ghosh, P.; Adil, A. A review on hybrid nanofluids: Recent research, development and applications. Renew. Sustain. Energy Rev. 2015, 43, 164–177. [Google Scholar] [CrossRef]
  21. Rajesh, V.; Sheremet, M.A.; Öztop, H.F. Impact of hybrid nanofluids on MHD flow and heat transfer near a vertical plate with ramped wall temperature. Case Stud. Therm. Eng. 2021, 28, 101557. [Google Scholar] [CrossRef]
  22. Sidik, N.A.C.; Adamu, I.M.; Jamil, M.M.; Kefayati, G.H.R.; Mamat, R.; Najafi, G. Recent progress on hybrid nanofluids in heat transfer applications: A comprehensive review. Int. Commun. Heat Mass Transf. 2016, 78, 68–79. [Google Scholar] [CrossRef]
  23. Dhinesh Kumar, D.; Valan Arasu, A. A comprehensive review of preparation, characterization, properties and stability of hybrid nanofluids. Renew. Sustain. Energy Rev. 2018, 81, 1669–1689. [Google Scholar] [CrossRef]
  24. Vemula, R.; Chamkha, A.; Mallesh, M.P. Nanofluid flow past an impulsively started vertical plate with variable surface temperature. Int. J. Numer. Methods Heat Fluid Flow 2016, 26, 328–347. [Google Scholar] [CrossRef]
  25. Sheremet, M.A.; Groşan, T.; Pop, I. Free convection in shallow and slender porous cavities filled by a nanofluid using Buongiorno’s model. J. Heat Transf. 2014, 136, 082501. [Google Scholar] [CrossRef]
  26. 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]
  27. Zahan, I.; Nasrin, R.; Khatun, S. Thermal performance of ternary-hybrid nanofluids through a convergent-divergent nozzle using distilled water–ethylene glycol mixtures. Int. Commun. Heat Mass Transf. 2022, 137, 106254. [Google Scholar] [CrossRef]
  28. 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]
  29. Sneha, K.N.; Vanitha, G.P.; Mahabaleshwar, U.S.; Laroze, D. Effect of couple stress and mass transpiration on ternary hybrid nanoliquid over a stretching/shrinking sheet with heat transfer. Micromachines 2022, 13, 1694. [Google Scholar] [CrossRef]
  30. Hafeez, M.B.; Krawczuk, M.; Nisar, K.S.; Jamshed, W.; Pasha, A.A. A finite element analysis of thermal energy inclination based on ternary hybrid nanoparticles influenced by induced magnetic field. Int. Commun. Heat Mass Transf. 2022, 135, 106074. [Google Scholar] [CrossRef]
  31. Adnan; Ashraf, W. Thermal efficiency in hybrid (Al2O3-CuO/H2O) and ternary hybrid nanofluids (Al2O3-CuO-Cu/H2O) by considering the novel effects of imposed magnetic field and convective heat condition. Waves Random Complex Media 2022. [Google Scholar] [CrossRef]
  32. Al Oweidi, K.F.; Shahzad, F.; Jamshed, W.; Usman; Ibrahim, R.W.; El Din, E.S.M.T. Partial differential equations of entropy analysis on ternary hybridity nanofluid flow model via rotating disk with hall current and electromagnetic radiative influences. Sci. Rep. 2022, 12, 20692. [Google Scholar] [CrossRef] [PubMed]
  33. Puneeth, V.; Anandika, R.; Manjunatha, S.; Ijaz Khan, M.; Imran Khan, M.; Althobaiti, A. Implementation of modified Buongiorno’s model for the investigation of chemically reacting rGO-Fe3O4-TiO2-H2O ternary nanofluid jet flow in the presence of bio-active mixers. Chem. Phys. Lett. 2022, 786, 139194. [Google Scholar] [CrossRef]
  34. Mahmood, Z.; Khan, U. Unsteady three-dimensional nodal stagnation point flow of polymer-based ternary-hybrid nanofluid past a stretching surface with suction and heat source. Sci. Prog. 2023, 106, 003685042311527. [Google Scholar] [CrossRef]
  35. Turkyilmazoglu, M. Unsteady convection flow of some nanofluids past a moving vertical flat plate with heat transfer. J. Heat Transf. 2014, 136, 031704. [Google Scholar] [CrossRef]
  36. Khanafer, K.; Vafai, K.; Lightstone, M. Buoyancy-driven heat transfer enhancement in a two-dimensional enclosure utilizing nanofluids. Int. J. Heat Mass Transf. 2003, 46, 3639–3653. [Google Scholar] [CrossRef]
  37. De Vahl Davis, G. Natural convection of air in a square cavity: A bench mark numerical solution. Int. J. Numer. Methods Fluids 1983, 3, 249–264. [Google Scholar] [CrossRef]
  38. Sathiyamoorthy, M.; Chamkha, A.J. Natural convection flow under magnetic field in a square cavity for uniformly (or) linearly heated adjacent walls. Int. J. Numer. Methods Heat Fluid Flow 2012, 22, 677–698. [Google Scholar] [CrossRef]
Figure 1. Problem model and coordinate system.
Figure 1. Problem model and coordinate system.
Nanomaterials 13 02860 g001
Figure 2. Computational mesh was selected for the study.
Figure 2. Computational mesh was selected for the study.
Nanomaterials 13 02860 g002
Figure 3. Validation with average Nusselt number [36,37].
Figure 3. Validation with average Nusselt number [36,37].
Nanomaterials 13 02860 g003
Figure 4. Comparison of isotherms for Pr = 0.054 and Ra = 105. (a) Isotherms for m = 1; (b) Isotherms for m = 0 [38].
Figure 4. Comparison of isotherms for Pr = 0.054 and Ra = 105. (a) Isotherms for m = 1; (b) Isotherms for m = 0 [38].
Nanomaterials 13 02860 g004
Figure 5. Isolines at b1 = b2 = b3 = 0.04, when m = 1 (linear heating). (a) Ra = 103; (b) Ra = 104; (c) Ra = 105.
Figure 5. Isolines at b1 = b2 = b3 = 0.04, when m = 1 (linear heating). (a) Ra = 103; (b) Ra = 104; (c) Ra = 105.
Nanomaterials 13 02860 g005
Figure 6. Isolines at b1 = b2 = b3 = 0.04, m = 0 (uniform heating). (a) Ra = 103; (b) Ra = 104; (c) Ra = 105.
Figure 6. Isolines at b1 = b2 = b3 = 0.04, m = 0 (uniform heating). (a) Ra = 103; (b) Ra = 104; (c) Ra = 105.
Nanomaterials 13 02860 g006
Figure 7. Average Nusselt numbers with Ra for different fluids when m = 1 (linear heating).
Figure 7. Average Nusselt numbers with Ra for different fluids when m = 1 (linear heating).
Nanomaterials 13 02860 g007
Figure 8. Mean Nusselt numbers with Ra for various fluids when m = 0 (uniform heating).
Figure 8. Mean Nusselt numbers with Ra for various fluids when m = 0 (uniform heating).
Nanomaterials 13 02860 g008
Figure 9. Mean Nusselt number with b3 for ternary hybrid nanofluid at Ra = 1000 when m = 1 (linear heating) and m = 0 (constant heating).
Figure 9. Mean Nusselt number with b3 for ternary hybrid nanofluid at Ra = 1000 when m = 1 (linear heating) and m = 0 (constant heating).
Nanomaterials 13 02860 g009
Table 1. Physical attributes of Cu–CuO–Al2O3/H2O.
Table 1. Physical attributes of Cu–CuO–Al2O3/H2O.
PropertyTernary Nanosuspension (Cu–CuO–Al2O3/H2O)
Density (ρ) ρ t h n f = 1 b 3 1 b 2 1 b 1 ρ f + b 1 ρ s 1 + b 2 ρ s 2 + b 3 ρ s 3
Dynamic viscosity (μ) μ t h n f = μ f 1 b 1 2.5 1 b 2 2.5 1 b 3 2.5
Heat capacity (ρcp) ρ c p t h n f = 1 b 3 × × 1 b 2 1 b 1 ρ c p f + b 1 ρ c p s 1 + b 2 ρ c p s 2 + b 3 ρ c p s 3
Thermal expansion coefficient (ρβ) ρ β t h n f = 1 b 3 1 b 2 1 b 1 ρ β f + b 1 ρ β s 1 + b 2 ρ β s 2 + + b 3 ρ β s 3
Thermal conductivity (κ) κ t h n f = κ h n f κ s 3 + 2 κ h n f 2 b 3 κ h n f κ s 3 κ s 3 + 2 κ h n f + b 3 κ h n f κ s 3 , where
κ h n f = κ n f κ s 2 + 2 κ n f 2 b 2 κ n f κ s 2 κ s 2 + 2 κ n f + b 2 κ n f κ s 2 and
κ n f = κ f κ s 1 + 2 κ f 2 b 1 κ f κ s 1 κ s 1 + 2 κ f + b 1 κ f κ s 1
Table 2. Physical characteristics of host liquids and nanoadditives [35].
Table 2. Physical characteristics of host liquids and nanoadditives [35].
PropertiesH2O (f)CuO (s1)Cu (s2)Al2O3 (s3)
ρ (kg/m3)997.1632089333970
cp (J/kgK)4179531.8385765
β (1/K)21 × 10–51.8 × 10–51.67 × 10–50.85 × 10–5
κ (W/mK)0.61376.540140
Table 3. (a). Mesh refinement check with NuL and NuB for Ra = 103, m = 1, Pr = 6.2, b1 = b2 = b3 = 0.04. (b). Mesh refinement check with NuL and NuB for Ra = 103, m = 0, Pr = 6.2, b1 = b2 = b3 = 0.04.
Table 3. (a). Mesh refinement check with NuL and NuB for Ra = 103, m = 1, Pr = 6.2, b1 = b2 = b3 = 0.04. (b). Mesh refinement check with NuL and NuB for Ra = 103, m = 0, Pr = 6.2, b1 = b2 = b3 = 0.04.
(a)
Elements NumberNodes NumberNuLNuB
Extremely Coarse1921451.24651.5958
Extra Coarse3542501.246781.59534
Coarser5403611.246321.59534
Coarse10126451.246261.5953
Normal15109301.246361.59518
Fine251614751.246281.59514
Finer653637311.245941.59548
Extra fine16,94693741.24581.59564
Extremely fine26,35214,0771.245781.59566
(b)
Elements NumberNodes NumberNuLNuB
Extremely Coarse1921452.39382.707
Extra Coarse3542502.76833.0627
Coarser5403612.96653.2585
Coarse10126453.35233.6427
Normal15109303.56013.8578
Fine251614753.74834.048
Finer653637314.43094.731
Extra fine16,94693745.03015.3296
Extremely fine26,35214,0775.03115.3294
Table 4. Statistics of effective mesh type (extra fine).
Table 4. Statistics of effective mesh type (extra fine).
Grid InformationGrid Statistics
Number of elements
Quads1200
Edge elements600
Triangles15,746
Average element quality0.8028
Number of Nodes9374
Vertex elements4
Table 5. Mean Nusselt numbers with Ra for different nanosuspension when m = 1 (linear heating).
Table 5. Mean Nusselt numbers with Ra for different nanosuspension when m = 1 (linear heating).
CuO–Al2O3/WaterCu–Al2O3/WaterCu–CuO/Water
RaNuLNuBNuLNuBNuLNuB
1031.06891.466281.070541.471141.070141.47836
1040.243.39620.239863.41160.237723.4418
1050.462086.80340.46196.8180.465466.8516
Table 6. Mean Nusselt number with Ra for various hybrid nanofluids at m = 0 (uniform heating).
Table 6. Mean Nusselt number with Ra for various hybrid nanofluids at m = 0 (uniform heating).
CuO–Al2O3/WaterCu–Al2O3/WaterCu–CuO/Water
RaNuLNuBNuLNuBNuLNuB
1034.44424.78384.45444.79684.4634.8119
1044.3716.23614.38336.25594.39846.2847
1055.42998.6125.44298.64255.46358.6863
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

Rajesh, V.; Sheremet, M. Free Convection in a Square Ternary Hybrid Nanoliquid Chamber with Linearly Heating Adjacent Walls. Nanomaterials 2023, 13, 2860. https://doi.org/10.3390/nano13212860

AMA Style

Rajesh V, Sheremet M. Free Convection in a Square Ternary Hybrid Nanoliquid Chamber with Linearly Heating Adjacent Walls. Nanomaterials. 2023; 13(21):2860. https://doi.org/10.3390/nano13212860

Chicago/Turabian Style

Rajesh, Vemula, and Mikhail Sheremet. 2023. "Free Convection in a Square Ternary Hybrid Nanoliquid Chamber with Linearly Heating Adjacent Walls" Nanomaterials 13, no. 21: 2860. https://doi.org/10.3390/nano13212860

APA Style

Rajesh, V., & Sheremet, M. (2023). Free Convection in a Square Ternary Hybrid Nanoliquid Chamber with Linearly Heating Adjacent Walls. Nanomaterials, 13(21), 2860. https://doi.org/10.3390/nano13212860

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