The Impact of Nanofluid Flow in a Rhombus-Shaped Porous Enclosure on a Constant Heated Obstacle
##plugins.themes.bootstrap3.article.main##
In this study, the focus is on investigating the characteristic of a nanofluid consisting of water as the base fluid and Copper nanoparticles. The analysis examines the impact of this nanofluid flow on a heated obstacle placed in a rhombus-shaped enclosure filled with porous media. The driving force behind the mixed convection is assumed to be the temperature difference created by the moving walls, where the top and bottom walls have a cold-to-hot flow in opposite directions. To obtain numerical solutions, the Galerkin weighted residual method is employed. The study explores the effects of the Darcy number and Richardson number on varies aspects, including streamlines, isotherms, dimensionless temperature, velocity profiles, average Nusselt numbers, and average fluid temperature. These effects are visualized through graphical representations. The findings reveal that both the Richardson number and Darcy number significantly influence the streamlines and isotherms within the system. Furthermore, the Darcy number proves to be an important control parameter for heat transfer in fluid flow through a porous medium in the enclosure. A heat transfer correlation for the average Nusselt number is provided, considering different Darcy numbers and Richardson numbers. To ensure the validity of the research, a comparison between the obtained results and previously published findings is conducted, showing a favorable agreement between them.
Introduction
Studying the flow of nanofluids within a rhombus-shaped enclosure filled with porous media and featuring a constant heated obstacle is a topic of great relevance in today’s world. Such enclosures, when filled with electrically conducting fluids, play a central role in various engineering and geophysical systems. The outcomes of this research have the potential to be applied in a range of practical applications including escalators, lifts, ventilation systems, heating and cooling technologies, building climate control, heat exchangers, solar collectors, as well as in astrophysics and biology. In-depth discussions on these application areas can be found in the comprehensive reviews by Nield and Bejan [1] and Ingham and Pop [2]. Improving heat transfer in cavities has become a critical concern in both industrial and energy sectors. Over the past years, researchers have conducted extensive investigations using experimental, analytical, and numerical methods to enhance heat transfer within cavities. One notable study by Xuan and Li [3] focused on the utilization of nanofluids to enhance heat transfer. By introducing CuO nanoparticles into the base fluid, they discovered that nanofluids have significant potential in improving the heat transfer process. In order to meet the growing demands of modern technologies, such as chemical production, power stations, and microelectronics, there is a need to develop fluids that can effectively facilitate heat exchange. A promising class of fluids in this regard is nanofluids. These nanofluids are a type of heat transfer fluid that contains a small amount of nanoparticles suspended in a base fluid. It has been extensively studied by numerous researchers [4]–[6] that the addition of nanoparticles to the base fluid, even at low volume fractions (1%–5%), can significantly enhance the thermal conductivity by approximately 20%.
In a related study, Rahman et al. [7] investigated the effects of natural convection on heat and mass transfer within a curvilinear triangular cavity. Additionally, Chamkha and Naser [8] conducted research on double-diffusive convection occurring in an inclined porous enclosure characterized by opposing temperature and concentration gradients. Finite volume method was employed by Teamah [9] to conduct numerical simulations of double-diffusive natural convective flow in a rectangular enclosure. Their study considered the influence of a magnetic field and heat source on the flow. Basak et al. [10] examined the characteristics of lid-driven mixed convection in a porous square cavity. They found that at Reynolds number of 100, the dominant form of heat transfer was forced convection when the left wall linearly heated and the right wall was cooled. Boulahia et al. [11] investigated the combined convection of nanofluids within a square enclosure with a triangular heated block. Their findings indicated that the heat transfer rate increased with higher volume fractions of nanoparticles and lower Richardson numbers. Furthermore, Dawood and Teamah [12] conducted a study on hydro-magnetic mixed convection double-diffusive flow in a lid-driven square cavity. Oztop et al. [13] investigated MHD mixed convection in a lid-driven cavity with a heater located at one of the corners. Iwastu et al. [14] carried out numerical simulations on mixed convection within a driven cavity featuring a stable vertical temperature gradient. Munshi et al. [15] analyzed the numerical simulation of mixed convection heat transfer using nanofluids in a lid-driven square enclosure filled with a porous medium. Additionally, Munshi et al. [16] investigated hydrodynamic mixed convection in a lid-driven square cavity, which included an elliptic-shaped heated block with a corner heater. Saha et al. [17] conducted research on mixed convection heat transfer in a lid-driven cavity with a wavy bottom surface. Maxwell-Garnett [18] writes colures in metal glasses and in metallic films. Brinkman [19] analyzed the viscosity of concentrated suspensions and solution. Zienktewicz et al. [20] analyzed the finite element method.
The problem of mixed convection heat transfer in a lid-driven enclosure with various aspects has received limited attention in existing literature. Therefore, this study aims to address this research gap by proposing the impact of nanofluid flow in a rhombus-shaped porous on a constant heated obstacle. The primary objective of the present study is to simulate and analyze the aforementioned scenario.
Physical Configuration
Fig. 1 illustrates a schematic representation of the simulation setup for studying the effects of a constant heated obstacle on nanofluid flow in a rhombus-shaped enclosure filled with porous media. The top wall of the enclosure is characterized by a cold temperature filled with porous media. The top wall of the enclosure is characterized by a cold temperature and the lid moves from left to right, while the bottom wall also has a cold temperature but with the lid moving in the opposite direction, from right to left. The left and right walls of the enclosure are considered adiabatic, and within the enclosure, there exists a constant heated obstacle.
Mathematical Formulation
The given mathematical equations can be expressed as follows:
Brinkman [19] as follows:
Maxwell (by the Maxwell model [18]) as follows:
The effective dynamic viscosity of the Copper-water nanofluid is calculated according to the Brinkman model [19] (Table I).
Physical properties | Water | Copper |
---|---|---|
Cp | 4179 | 385 |
997.1 | 8933 | |
K | 0.613 | 401 |
μ | 0.855 × 10−3 | – |
Boundary Conditions (Dimensional)
The dimensional boundary conditions for the solving the governing (1)–(4) are presented in Table II.
Walls | Boundary conditions |
---|---|
Top wall | |
Left wall | |
Right wall | |
Bottom wall | |
Constant heated obstacle |
Dimensionless Analysis
Using the following dimensionless parameters, the governing equations can be converted to the dimensionless forms:
The dimensionless forms of the governing equations expressed are:
And the dimensionless governing parameters:
where Pr, Da, Re, Gr, Ri and Ha are Prandtl number, Darcy number, Reynolds number, Grashof number, Richardson number and Hartmann number, respectively.
Boundary Conditions (dimensionless)
The dimensionless boundary conditions for the solving the governing (11)–(14) are as follows (Table III):
Walls | Dimensionless boundary conditions |
---|---|
Top wall | |
Left wall | |
Right wall | |
Bottom wall | |
Constant heated obstacle |
Numerical Technique
The Galerkin Weighted residual finite-element method is used to convert the governing partial differential equations into a set of integral equations, each equation term is integrated using the Gauss quadrature method. After applying boundary conditions, the resulting nonlinear algebraic equations are modified and transformed into linear algebraic equations using newton’s method. Finally, the Triangular factorization method is utilized to solve these linear equations.
According to Table IV, the average Nusselt number was determined for various grids. The table clearly demonstrates that the grid with 6743 nodes and 12270 elements yields a desirable solution for the numerical study at hand.
Nodes | 832 | 1114 | 1725 | 6743 | 24136 |
Elements | 1601 | 2146 | 3328 | 12270 | 48468 |
8.64276 | 8.66746 | 8.67727 | 8.70372 | 8.70841 | |
Time (s) | 8 | 11 | 14 | 20 | 35 |
Program Validation
Fig. 2 displays a comparison of streamlines and isotherms in the graphical solution. The figures indicate a remarkable consistency between the two results, highlighting a strong agreement between them.
Results and Discussion
In this section, the numerical and graphical results of simulating the flow of nanofluid in a rhombus-shaped enclosure filled with porous media containing a heated obstacle are presented. The results have been obtained for the Darcy number ranging from and the Richardson number .
Variation of Darcy Number
Fig. 3 illustrates the variation of streamlines at different Darcy numbers. The figures clearly demonstrate that when symmetric boundary conditions are applied to the top and bottom sides of the enclosure, the flow exhibits symmetry along this line. As the Darcy number and buoyancy force increase, the vortices’ centers move in the upward and downward directions towards the lid, indicating a stronger influence of the moving lid.
At , the isotherms in Fig. 4 indicate that conduction dominates the heat transfer. As the Darcy number increases, the isotherms become condensed near the upper and lower side walls, indicatind an increase in heat transfer through convection. The presence of a thermal boundary layer can be observed from the isotherms at When , the isotherms show pure conduction. A similar temperature distribution is seen for
Fig. 5 depicts the variation of the local Nusselt number along the X-axis of the enclosure for different Darcy numbers, with . The figure shows district curves representing the minimum and maximum values. It can be observed that at , there is a substantial change in the local Nusselt number. However, as the Darcy number increases, the impact on the Nusselt number becomes less significant. Fig. 6 illustrates the changes in the velocity profile along the X-axis of the enclosure for various Darcy numbers, considering . The figure displays distinct curves representing the maximum and minimum values. Specifically, for X values less than 0.6, the velocity curves exhibit higher values, while for X values greater than 0.6, the velocity curves show lower values. In Fig. 7, the dimensionless temperature is plotted along the X-axis for various Darcy numbers, while considering . The figure indicates that as the Darcy numbers increase, the temperature values also increase. The parabolic shape of the temperature profile is a result of the symmetric shape of the heated block. Lower Darcy numbers exhibit a less significant change in temperature, whereas higher Darcy numbers demonstrates a more noticeable effect on temperature.
Heat Transfer Rates
This section presents a numerical simulation of nanofluid flow within a rhombus-shaped enclosure filled with porous media, where a constant heated obstacle is present. The relationship between the average Nusselt number and Darcy numbers, considering various Richardson numbers for heat transfer rates, is depicted in Fig. 8. It can be observed that the average Nusselt number increases as the Darcy number increases, Generally, higher Darcy number result in higher temperature gradients, which consequently lead to an increase in the average Nusselt number across different Richardson numbers.
Fig. 9 displays the relationship between tha average fluid temperature and the Darcy number within the enclosure, considering different values of the Richardson number ranging from 0.1 to 10, while keeping the Hartmann number at 50 and the porosity at 5%. The findings indicate that an elevated Richardson number within the porous enclosure leads to an augmentation in the average fluid temperature. Furthermore, the rate at which heat transfer increases distinctly relies on the value of the Darcy number.
Variation of Richardson Number
This section presents a numerical simulation of the flow of nanofluid in a rhombus-shaped enclosure filled with porous media, where there is a constant heated obstacle. The figures, Figs. 10 and 11, demonstrate the changes in streamlines and isotherms within the constant heated obstacle for various values of the Richardson number. The parameters used in the simulation are It is evident from the figures that symmetrical boundary conductions exist in the system. At lower values of the Richardson number, the fluid inside the enclosure raises more, and the motion of the lid becomes more vigorous. Conversely, for higher Richardson numbers, three elliptic-shaped regions, referred to as ‘eyes’ form inside the enclosure, indicating increased flow strength. The isotherms shown in Fig. 11 indicate that heat transfer is primarily dominated by conduction. As the Richardson number increases, the isotherms tend to concentrate near the constant heated obstacle, and their lines become more curved, suggesting enhanced heat transfer through convection.
Fig. 12 illustrates the local Nusselt number along the horizontal wall of the enclosure for various Richardson numbers, considering parameters The curves in the plote exhibit minimum and maximum shapes. For values of X less than 0.3, the curves take on a minimum shape, while for values of X greater than 0.3, the curves assume a maximum shape. Additionally, the absolute value of the local Nusselt number, representing the heat transfer rate, increases as the Richardson number rises. In Fig. 13, the velocity profiles along the horizontal wall of the enclosure are depicted for different Richardson numbers, taking into account parameters . It is evident from the figure that both the maximum and minimum velocities increase as the Richardson number increases. The change in velocity is relatively small for lower Richardson numbers, but it becomes more pronounced for higher Richardson numbers. Fig. 14 illustrates the temperature profiles along the horizontal wall of the enclosure for various Richardson numbers, considering parameters parameters . It can be observed from the figure that the temperature decreases as the Richardson number increases. The change in temperature is relatively less significantly for lower Richardson numbers, whereas for higher Richardson numbers, the temperature exhibits a more pronounced and significant change.
Heat Transfer Rates
This section presents a numerical simulation of nanofluid flow in a rhombus-shaped enclosure filled with porous media, featuring a constant heated obstacle. The simulation focuses on the average Nusselt number in relation to the Richardson number for different values of the Darcy number, while considering parameters
Fig. 15 displays the relationship between the average Nusselt number and the Richardson number for various Darcy numbers. It can be observed from the figure that as the Darcy number increases, the average Nusselt number also increases. This can be attributed to higher temperature gradients associated with larger Darcy numbers. Overall, the average Nusselt number exhibits an increasing trend with Darcy numbers for different values of the Richardson number.
Fig. 16 presents the relationship between the average fluid temperature and the Richardson number inside the constant heated enclosure, considering different Darcy numbers, while keeping constant. It is evient from the figure that increasing the Darcy number inside the constant heated porous enclosure leads to an increase in the average fluid temperature. Additionally, the rate of increment in heat transfer is dependent on the specific value of the Darcy number.
Conclusion
The impact of nanofluid flow in a rhombus-shaped porous on a constant heated obstacle has been conducted. This study encompasses a wide range of dimensionless groups, including the Darcy number and Richardson number. Both the Darcy number and the movement of the lids have a substantial impact on the flow and thermal fields within the enclosure. The Richardson number serves as a useful control parameter for heat transfer in fluid flow through porous media enclosures. The velocity of the fluid is enhanced as the solid volume fraction of nanoparticles increases, resulting in an increase in the heat transfer rate. The highest rate of heat transfer is observed in the vicinity of the heated obstacle within the enclosure. Furthermore, an increase in the Darcy number of the heated enclosure leads to an augmentation in the heat transfer rate.
References
-
Nield DA, Bejan A. Convection in Porous Media. 2nd ed. New York: Springer; 1999.
Google Scholar
1
-
Ingham DB, Pop I. Transport Phenomena in Porous Media. Science Direct: Pergamon; 1998.
Google Scholar
2
-
Xuan Y, Li Q. Heat transfer enhancement of nanofluids. Int J Heat Fluid Flow. 2000;21:58–64.
Google Scholar
3
-
KhanaferK, Vafai K, Lightstone M. Buoyancy-driven heat transfer enhancement in a two dimensional enclosure utilizing nanofluids. Int J Heat Mass Trans. 2003;46:3639–53.
Google Scholar
4
-
Abu-Nada A, Masoud Z, Hijazi A. Natural convection heat transfer enhancement in horizontal concentric annuli using nanofluids. Int Commn Heat Mass Trans. 2008;35:657–65.
Google Scholar
5
-
Ahmed SE, Rashad AM, Goria RSR. Natural convection in triangular enclosures filled with nanofluid saturated porous media. J Thermo Heat Trans. 2013;27:700–6.
Google Scholar
6
-
Rahman MM, OztopHF, AhsanA, Orfi J.Natural convection effects on heat and mass transfer in a curvilinear triangular cavity. Int J Heat Mass Trans. 2012;55:6250–9.
Google Scholar
7
-
Chamkha AJ, Naser HA. Double-diffusive convection in an inclined porous enclosure with opposing temperature and concentration gradients. Int J Therm Sci. 2001;40(3):227–44.
Google Scholar
8
-
Teamah MA. Numerical simulation of double diffusive natural convection in rectangular enclosure in the presences of magnetic field and heat source. Int J Therm Sci. 2008;47:237–48.
Google Scholar
9
-
Basak T, Roy S, Singh SK, Pop I. Analysis of mixed convection in a lid-driven porous square cavity with linearly heated side wall(s). Int J Heat Mass Trans. 2010;53:1819–40.
Google Scholar
10
-
Boulahia Z, Wakif A, Sehaqui R. Numerical investigation of mixed convection heat transfer of nanofluid in a lid driven square cavity with three triangular heating blocks. Int J Therm Sci. 2016;143(6):37–45.
Google Scholar
11
-
Dawood MMK, TeamahMA. Hydro-magnetic mixed convection double diffusive in a lid driven square cavity. Eur J Sci Res. 2012;85:336–55.
Google Scholar
12
-
Oztop HF, Al-Salem K, Pop I. MHD mixed convection in a lid-driven cavity with corner heater. Int J Heat Mass Trans. 2011;54:3494–504.
Google Scholar
13
-
Iwatsu R, Hyan JM, Kuwahara K. Mixed convection in a driven cavity with a stable vertical temperature gradient. Int J Heat Mass Trans. 1993;36:1601–8.
Google Scholar
14
-
Munshi MJH, Alim MA, Bhuiyan AH, Ahmed KFU. Numerical simulation of mixed convection heat transfers of nanofluid in a lid-driven porous medium square enclosure. AIP Conf. Procedia Eng. 2019;2121:1–9.
Google Scholar
15
-
Munshi MJH, Alim MA, Bhuiyan AH, Ali M. Hydrodynamic mixed convection in a lid-drivensquare cavity including elliptic shape heated block with corner heater. Published by Elsevier, Procedia Eng. 2017;194:442–9.
Google Scholar
16
-
Saha LK, Somadder MC, Uddin KMS. Mixed convection heat transfer in a lid driven cavity with wavy bottom surface. Am J Appl Math. 2013;1:92–101.
Google Scholar
17
-
Maxwell-Garnett. Colures in metal glasses and in metallic films. Philosophical Transaction of the Royal Society of London. Series A. Containing Papers of a Mathematical or Physical Character. 1904;203:385–420.
Google Scholar
18
-
Brinkman HC. The viscosity of concentrated suspensions and solution. J Chem Phys. 1952;20:571–81.
Google Scholar
19
-
Zienktewicz OC, Taylor RL. The Finite Element Method. 4th ed. McGraw-Hill: American publishing company; 1991.
Google Scholar
20