Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 13, Article number: 21171 (2023 ) Cite this article Ebara Absorption Heat Pump
This study is numerically executed to investigate the influence of heat generation or absorption on free convective flow and temperature transport within a wavy triangular enclosure filled by the nanofluid taking the Brownian effect of nanoparticles. The water (H2O) is employed as base fluid and copper (Cu) as nanoparticles for making effective Cu–H2O nanofluids. The perpendicular sinusoidally wavy wall is cooled at low temperature while the horizontal bottom sidewall is heated non-uniformly (sinusoidal). The inclined wall of the enclosure is insulated. The governing dimensionless non-linear PDEs are executed numerically with the help of the Galerkin weighted residual type finite element technique. The numerically simulated results are displayed through average Nusselt number, isothermal contours, and streamlines for the various model parameters such as Hartmann number, Rayleigh number, heat generation or absorption parameter, nanoparticles volume fraction, and undulation parameter. The outcomes illustrate that the temperature transport rate augments significantly for the enhancement of Rayleigh number as well as nanoparticles volume fraction whereas reduces for the increment of Hartman number. The heat transfer is significantly influenced by the size, shape, and Brownian motion of the nanoparticles. The rate of heat transport increases by 20.43% considering the Brownian effect for 1% nanoparticle volume. The thermal performance increases by 8.66% for the blade shape instead of the spherical shape of nanoparticles. In addition, heat transfer is impacted by the small size of nanoparticles. The thermal transport rate increases by 35.87% when the size of the nanoparticles reduces from 100 to 10 nm. Moreover, the rate of heat transmission increases efficiently as the undulation parameter rises. It is also seen that a crucial factor in the flow of nanofluids and heat transmission is the heat generation/absorption parameter that influences temperature distribution, heat transfer rates, and overall thermal performance.
The investigation of natural convection using nanofluids remains a promising field for raising heat transfer efficiency and thermal management in a variety of applications. Natural convective heat transfer within a triangular cavity has widespread applications in numerous industrial, and engineering systems like heat exchangers, ventilation systems, solar collectors, fire prevention, refrigeration, geothermal reservoirs, and so on. A number of researchers1,2,3,4,5,6,7 investigated the nanofluid under different circumstances. Aydin et al.8 investigated free convective flow into a rectangular enclosure warmed by one side and consoled from the roof. Holtzman et al.9 performed free convective laminar flow into an isosceles triangle cavity warmed from the bottom side. Buongiorno10 investigated convective transport in nanofluids. Das et al.11 observed about nanofluids and their applications for science and technology. Free convective within a triangle-shaped cavity by an extending isothermal cooker was observed by Varol et al.12. Basak et al.13 investigated free convection simulation for both uniform and non-uniform systems within an isosceles triangle enclosure. Kent14 observed a free convective numerical solution in a triangle chamber where the base sidewall is cooled and the inclined wall is heated. Kamiyo et al.15 performed an extensive study on natural convective flow within a triangular cavity. Bhuvaneswari et al.16 performed hydromagnetic flow within a square-shaped cavity using a non-uniform heated sidewall.
The natural convection in a square chamber filled with a water-based nanofluid (water containing Cu particles) was explored by Tayebi et al.17. Tayebi and Ali J. Chamkha18 researched to comprehend the properties of natural convection heat transfer and flow using Cu/H2O in the square domain subjected to a horizontal magnetic field and different conductive boundary configurations. The natural convection flow mechanism and heat exchange under a magnetic field within a concentric circular annulus between a heat-generating conductive internal cylinder and an isothermally cold exterior cylinder filled with a CNTs-water-based nano-liquid were examined by Tayebi et al.19. The study conducted by Alsabery et al.20 examined the creation of entropy and mixed convection flow within a wavy-walled enclosure that contains a heat source and a rotating solid cylinder. The influence of natural convection within a square cavity, in addition to the presence of a conducting solid block and a corner heater, was numerically examined by Alsabery et al.21. The heat and flow exchange properties in a new configuration saturated with a Ag–MgO hybrid nanofluid were explored numerically by Hussain et al.22. Nasrin23 investigated MHD combined convective flow into the wavy cavity by account of Joule heating. Bhardwaj et al.24 performed an exploration of natural convective temperature disposal with entropy generation into a triangle-shaped cavity. Hamida and Charrada25 studied temperature generation or absorption effects upon free convective temperature flow into a square cavity using nanofluids including the magneto-hydrodynamic effect. Suvash et al.26 investigated free convective flow within a triangle-shaped cavity using uniform thermal boundary conditions.
Several common types of nanoparticles are obtainable for commercial purposes like Al2O3, Fe3O4, Cu, TiO2, and so on. Some of these nanoparticles are also effectively employed to enhance thermal performance. Hydromagnetic natural convective temperature flow within an isosceles triangle shape cavity using nanofluid was investigated by Rahman et al.27. Mirabedin28 performed computational simulations on natural convective flow into a right-angled triangle-shaped cavity. Alam et al.29 studied numerically convective temperature flowing of micropolar fluid including heat generation/absorption. Triveni et al.30 studied free convective flowing into the triangle cavity using the different locations of the cold wall. Boulahia et al.31 performed combined convection temperature disposal within a square enclosure charged with nanofluid using triangle-shaped warmed-up obstacles. Uddin et al.32 studied computational simulation numerically for convective heat flow into nanofluid. Unsteady convective MHD numerical impact of sinusoidal warmed boundary condition with Fe3O4-water ferrofluid was performed by Alam et al.33. Al-Asad et al.34 studied hydrodynamics impacts of convective flow using a perpendicular fin within a wavy square cavity. Zahan et al.35 studied the combined magneto-hydrodynamic convective flow in a wavy cavity considering Joule heating. The magnetohydrodynamics natural convection and entropy generation of a Cu–water nanofluid constrained in a porous annulus was numerically discovered by Mourad et al.36. In order to examine the heat transport processes of a non-Newtonian hybrid nanofluid in an annular enclosure, Abderrahmane et al.37 explored the MHD convective transport. Alshare et al.38 investigated convection in a nanofluid-filled, lid-driven hollow exposed to a magnetic field.
An emerging magnetic field is also employed to control the mechanism in material manufacturing processes. Izadi et al.39 examined the thermal and flow characteristics by considering both uniform and nonuniform magnetic field effects as well as time-dependent heat sources. Izadi et al.40 also investigated the convective flow of Al2O3/water nanofluid an annulus. They reported the heat transfer coefficient is directly influenced by the Richardson number. Convection heat transfer rate improved when the single nano-fluid (water/aluminum oxide) inside the chamber's nano-particle volume fraction increased41. The behavior of phase change material (PCM) was also marginally affected by the volume fraction42. As the concentration of nanoparticles increases, the bottom heated wall's mean Nusselt number falls43. Izadi et al.44 also showed how the placement of the heat sources affects how quickly the melting front advances. The increase in the average Nusselt number and energy transport intensity is seen for the Rayleigh number45. Additionally, it has been noted that when the flow motion gets stronger due to convection, there is an increase in instability in the temperature and flow fields46. When the Rayleigh number rises, the rate of heat transfer increases47. The average Nusselt number rises to 61% for Ra = 105 and 12% for Ra = 10348. When the magnetic field is ignored (Ha = 0), the dimensionless melting time increases by 266% when the Ha is imposed up to 50049. In comparison to α = 60°, the average Nusselt number for the position of the heat source and sink (α = − 60°) is up to 124.5% greater. The heat transfer coefficient rises by 11.45% at Re 4000 by adding 1% nanoparticle volume to the base fluid49. Tarikul et al.50,51,52,53,54,55,56 recently studied numerically MHD influence on free convective temperature distributions into the different nanofluids with different thermal conditions.
It appears from this study that there is a lot of scope for research on nanofluids under different thermal conditions for checking thermal performance which has specific applications such as solar collectors, cooling/heating systems, and so on. The main intention of the current research is to explore convective flow and thermal transport into a wavy triangle-shaped cavity utilizing the nanofluids with the presence of heat generation/absorption and sinusoidal thermal boundary conditions. The influence of various physical parameters such as undulation parameter, heat generation/absorption, Rayleigh number, nanoparticle volume, Hartmann number, nanoparticles shape and size, and Brownian motion and other flow behaviors of nanofluids on the fluid flow and heat transfer performance are numerically studied and discussed from the physical standpoint.
A schematic view of the wavy triangular enclosure considered for the present research, boundary conditions, and the coordinates are displayed in Fig. 1. In this enclosure, the vertical sinusoidal wall is nursed at a low-temperature T = Tc, whereas its horizontal wall is warmed sinusoidally \(T = \sin \,(\pi x/L)\) . The hypotenuse of the enclosure is insulated. The copper–water nanofluid (Cu–H2O) is employed in the enclosure where water (H2O) is taken as base fluid and copper (Cu) as nanoparticles. The strength (B0) of the uniform magnetic field is employed towards the horizontal direction and normal to the vertical wavy wall. g = (0, g) is the gravitational acceleration that works downward direction towards the y-axis. The nanofluid is considered to be laminar, viscous, and incompressible. The two-dimensional time-independent flow is also assumed. The horizontal wall measures the x-axis and the vertical wall measures the y-axis. Every rigid boundary is accepted like a no-slip wall. The characteristics of base fluid and nanoparticles have been given in Table 1. In the absence of the electric field (\({\varvec{E}}\) = 0), the interaction of the fluid flow with the magnetic field applied to it gives rise to the Lorentz force and the uniform magnetic force with \({\varvec{B}}=( {B}_{0}, 0, 0)\) can be written as follows:
Schematic view for sinusoidal triangular enclosure.
The feature of the vertical wavy wall is employed to mimic the following model:
where A represents wave amplitude and λ represents undulation number.
The governing non-linear PDEs of conservation of mass, momentum, and energy for nanofluid for the aforementioned assumptions are displayed in the Cartesian coordinate system in dimensional form57,58:
The physical properties of nanofluids such as thermal diffusivity, specific heat, thermal expansion coefficient, viscosity, density, and electrical conductivity of nanofluids that appear in governing Eqs. (3)–(6) are given by the following formulations57,58:
The Brownian motion of nanoparticles may be added to the Maxwell classical correlation in order to determine the effective and efficient heat conductivity of a nanofluid taking the account of the random movement of the nanoparticles in it as follows58:
where ψ denotes sphericity, which is defined as the ratio of spherical surface area to real particle surface area with the same volumes, and n is the form factor of nanoparticles, expressed as n = 3/ψ. The reference temperature, Tref; the Boltzmann constant, KB = 1.38064852 × 10−23 JK−1; the blade, platelet, brick, cylinder, and spherical-shaped nanoparticles, and their diameters, d, are represented by the integers n = 8.6, 5.7, 4.9, 3.7, and 3.0, respectively.
The following transformation of variables (Eq. 7) is used to make the Eqs. (3)–(6) including boundary conditions (7)–(9) in the dimensionless form57,58:
The non-dimensionless governing equations by employing dimensionless variables are as follows57,58:
where, \(\Pr = \frac{{\upsilon_{bf} }}{{\alpha_{bf} }}\) is the Prandtl number, \({\text{Ra}} = \frac{{g\beta_{bf} \left( {T_{h} - T_{c} } \right)L^{3} }}{{\upsilon_{bf} \,\alpha_{bf} }}\) is the Rayleigh number, \({\text{Ha}} = B_{o} L\sqrt {\sigma_{bf} /\mu_{bf} }\) is the Hartmann number, and \(Q = \frac{{qL^{2} }}{{\alpha_{nf} \left( {\rho c_{p} } \right)_{nf} }}\) is the heat generation /absorption parameter.
Non-dimensional boundary conditions are as follows:
For this model, average Nusselt number on horizontal warmed wall has been calculated from the following expression57,58:
The non-dimensional highly non-linear PDEs (18)–(21) including boundary conditions (22)–(24) are numerically solved using Galerkin weighted residual finite element technique that is a powerful tool to handle these kinds of non-linear equations. The information about the numerical technique is narrated by Zienkiewicz and Taylor59. In this numerical technique, the solution domain is discretized firstly within finite grid numbers that are composed of multiform triangle elements. Six node triangle shape elements are employed to improve the finite element equations. Then, the Galerkin weighted residual scheme is executed into a system of integral equations from non-linear PDEs. To solve these integral equations, Gauss’s quadrature method is also employed. Thereafter, the boundary necessary boundary conditions are employed to modify the non-linear algebraic. Newton–Raphson iteration method is employed for solving the global non-linear algebraic equation in terms of Matrix. The standard of convergence including error inference for the solution is worked as \(\left| {\Upsilon^{m + 1} - \Upsilon^{m} } \right| \le 10^{ - 5}\) , where \(\Upsilon\) represents subsidiary variables \(U,\,\,V,\,\,\theta\) and m represent iteration number.
To guarantee the numerical scheme, a comprehensive mesh testing action is worked for the present problem of wavy triangular enclosure when Ra = 105 with fixed Ha = 15, ϕ = 0.04, Q = 4, λ = 0, and n = 3. We observed the following five different grid systems designated as normal, fine, finer, extra fine, and extremely fine including elements number into the resolution field: 1161, 1739, 4995, 13,315, and 16,753. The numerical design is executed employing average Nusselt number (Nuav) for checking the grid fineness for the above-mentioned elements is displayed by Fig. 2. The elements 13,315 depicts less difference in the outcomes than the others elements for present study. Therefore, for the numerical simulation, elements 13,315 is founded for fulfil the requirements of the grid independent solution.
Average Nusselt number (Nuav) of different element numbers.
Through the data regarding average Nusselt number (Nuav), the current results are equated with Ghasemi et al.60 in order to verify the validity of the current numerical scheme using several Rayleigh numbers (Ra) and nanoparticle volume (ϕ) using Al2O3–H2O nanofluid with Ha = 30, Pr = 6.8377, and n = 3 (Table 2). In contrast, the same code was used for the present nanofluid model. It seems that using the current numerical system fulfills with requirements.
This portion, we illustrate the outcomes of free convectional temperature transport and fluid flow of nanofluids into a wavy triangle shape enclosure by employing finite element method of Galerkin weighted residual type. The simulated outcomes form the numerical calculations of current problem is discussed for different model parameters designated as temperature generation/absorption parameter (Q), Rayleigh number (Ra), number of undulations (λ), volume fraction of nanoparticles (ϕ), Hartmann number (Ha), nanoparticles shape (n), size (d) and Brownian motion upon flow characteristics in isotherm contours, streamline patterns and average Nusselt number. The following domains of non-dimensional parameters for numerical computation are used: 104 ≤ Ra ≤ 106, 0 ≤ Ha ≤ 50, 0 ≤ ϕ ≤ 0.1, − 12 ≤ Q ≤ 12, 0 ≤ λ ≤ 3 and Pr = 6.8377, 1 nm ≤ d ≤ 100 nm, and n = 3, 3.7, 4.9, 5.7, and 8.6.
The fluid's buoyancy forces are measured by the Rayleigh number. It measures how significant natural convection is in relation to other heat transfer mechanisms like forced convection and conduction. Figure 3 demonstrates outcomes of the Rayleigh number (Ra = 105, 106, 107) using streamlines with fixed ϕ = 0.04, Ha = 15, n = 3, d = 10 nm and Q = 4. It is seen from this figure, for Ra = 105, the streamline contours are not forcefully influenced by nanofluid particles while the sinusoidally warmed-up boundary is imposed on the horizontal wall which is a good syndrome of potential conduction. The form of streamline contours is nearly uniform for all undulation parameters (λ). As increases of buoyancy-driven parameter Ra (> 105), fluid within the cavity is more debauched, and streamline contours capacity increases which indicates the convection mood of temperature transport. The streamlines also grab the maximum area within the enclosure for a higher Rayleigh number. In addition, for the increases of Ra, the denseness and vortex strength within the cavity are enhanced as well as turn familiar concerning the free convection. Also, for raising the value of the undulation parameter (λ) the extension of flow circulation turns out to impair with all Rayleigh numbers (Ra). This is due to the fact that the mounting value of λ produces a friction force between the wavy surface and the liquid that generates temperature which diminishes liquid rapidity within the enclosure.
Outcomes of Rayleigh number (\({\text{Ra}}\) ) for various undulation parameter (λ) using streamlines when ϕ = 0.04, Ha = 15, n = 3, d = 10 nm, and Q = 4.
Isotherms are useful to visualize the temperature flow style whether it is convection or conduction. It is also extremely efficient to measure the proficiency of temperature transference into nanofluid. Figure 4 depicts Rayleigh number influence (Ra = 105, 106, 107) on isotherms with fixed ϕ = 0.04, Ha = 15, n = 3, d = 10 nm and Q = 4. These figures illustrate that for low Ra, streamlines are superfluous compact near the bottom hot and perpendicular wavy cold walls into the enclosure for conduction type of temperature disposal. It is also seen that for the Rayleigh number, the isotherm's contours are parabolic in shape near horizontal hot walls on account of conduction.
Outcomes of Rayleigh number (Ra) for various undulation parameter (λ) using isotherms when ϕ = 0.04, Ha = 15, n = 3, d = 10 nm, and Q = 4.
The pattern of the isotherm contours as well as heated field become non-persistent with the variation of Ra for a fixed undulation parameter (λ). As the buoyancy driven parameter augments Ra (> 105), the crowded of isotherm contours into the cavity enhances and intimates neighbor horizontal hot and wavy cold walls which indicates that the formation of heated boundary layer on horizontal wall for conduction style of temperature changing. It is also seen that a dilute heated boundary is elaborated near perpendicular wavy wall. The wavy wall of the triangular domain can improved by the turbulence and vortices that accelerates the heat transfer rate in the nanofluid. It is also seen that better thermal performance is observed from the wavy surface's acceleration of nanofluid mixing as compared to a triangular cavity with plain walls. Because of the higher flow resistance caused by wavy walls, the pressure drop in the flow is also enhanced. Furthermore, isotherm contours are more compact and concentrate when undulation number increases. The physical significance is that a high Rayleigh number implies that buoyancy-driven convection dominates.
The outcomes of Hartmann number Ha (= 0, 10, 50) on streamlines with fixed ϕ = 0.04, Ra = 105, n = 3, d = 10 nm, and Q = 4 (Fig. 5). It is seen from the figure for the absence of the Hartmann number (Ha = 0) a clockwise circulation is formed nearly the whole enclosure. A tiny revolving chamber with a vortex is also generated at the center of the enclosure. As increases Ha, the pattern of the streamline has been changed and flow strength diminishes within the cavity. Because an external strong magnetic field generates a Lorentz force that weakens the streams within the enclosure. This figure also depicts that a small secondary counterclockwise vortex is created on the right corner of the enclosure for a higher Hartmann number. This is due to the fact that the flow enlarges in the middle and this leads to create a secondary tiny circulation cell inside the enclosure. These results also narrate that adjacent fluid near the bottom hot wall after getting more heated, and attempts to go overhead by virtue of buoyant force whilst comparatively refrigerant liquid neighbor cooled wavy wall moves to the horizontal hot wall which creates a secondary small vortex. In addition, for the increasing values of undulation parameter λ (= 1, 2, 3), a tiny circulation grows neighbor to the horizontal hot wall. As increases λ, the enclosure kernel fills with more streamlines. It is also seen that the alternation in λ causes an ordinary change in streamlines for a particular Hartmann number.
Outcomes of Hartmann number (Ha) for various undulation parameter (λ) using streamlines when ϕ = 0.04, Ra = 105, n = 3, d = 10 nm, and Q = 4.
The impact of Ha (= 0, 10, 50) using distributions of isotherms is demonstrated in Fig. 6 with fixed ϕ = 0.04, Ra = 105, n = 3, d = 10 nm, and Q = 4. The figures depict that isotherm lines are compact neighbor horizontal hot wall and vertical wavy cold wall inside the enclosure because Ha creates interior incident into the fluid that generates more temperature into the fluid. In the absence of an external magnetic effect (Ha = 0), isotherm contour patterns are more conspicuous which indicates that strong convection occurs for the nonappearance of the magnetic field. For higher Hartmann numbers, isotherm contour lines walk towards the vertical wavy surface and are almost parallel to each other. This is due to the fact that the strong Lorentz force from the Hartmann number reduces the temperature distribution. Besides, for a particular Hartmann number, a higher temperature transport is observed at λ = 3 than other values of the undulation parameter (λ).
Outcomes of Hartmann number (Ha) for various undulation parameter (λ) using isotherms when ϕ = 0.04, Ra = 105, n = 3, d = 10 nm, and Q = 4.
Therefore, for upper values of λ and lower Hartmann numbers, the isothermal contours are slightly denser near the wavy surface. A magnetic field's strength is an important factor that have varying effects on heat transfer rates, flow patterns, nanoparticle behavior, and energy efficiency in the context of nanofluid in a wavy triangular cavity. The velocity of charged nanoparticles within the nanofluid can be affected by a magnetic field. The convective motion is often suppressed by the Lorentz force acting on the charged nanoparticles. Because nanofluids have higher thermal conductivity when subjected to magnetic fields, these properties can change heat transfer in nanofluids. The alignment and motion of nanoparticles are influenced by the magnetic field, which optimizes heat transfer rates. The dispersion of nanoparticles in the nanofluid can also be impacted by the magnetic field. By preventing or controlling agglomeration, it can assist keep nanoparticles evenly distributed throughout the fluid. This is essential to preserving the nanofluid's essential characteristics.
The effects of heat generation Q (> 0) or heat absorption Q (< 0) on isotherm contours for variation of wave number λ (= 0, 1, 2, 3) when ϕ = 0.04, Ha = 15, Ra = 105, n = 3, and d = 10 nm are presented in Figs. 7 and 8, respectively. The compactness of the isotherms is denser near the hot and cold wall for the heat sink than to the heat source parameter. We observed that the denseness of the isotherms within the enclosure for the heat sink condition Q = − 12. But, the crowded isotherm lines are diminished slowly for non-dimensional temperature enhanced up to the higher values of heat generation case Q = 12. We also observed that for temperature absorption situations, enclosure heat is lesser compared to the temperature in fluid including heat generation. The temperature gradient near the enclosure wavy wall with the heat source case is inferior compared to the heat sink case. This phenomenon is observed for all undulation parameters (λ). The heat source increases fluid temperature and reduces temperature differences inside the enclosure. In addition, heat sink enhances, isotherm contours compact neighbor the wavy surface within the cavity. Also, isotherm contours are parabolic in appearance neighboring the warmed-up wall and a dilute heated boundary layer is developed neighbor the wall within the cavity.
Outcomes of Q (< 0) for various undulation parameter (λ) using isotherms when ϕ = 0.04, Ha = 15, Ra = 105, n = 3, and d = 10 nm.
Outcomes of Q (> 0) for various undulation parameter (λ) using isotherms when ϕ = 0.04, Ha = 15, Ra = 105, n = 3, and d = 10 nm.
Figure 9, respectively, demonstrates the influence of mean Nusselt number (Nuav) along the bottom hot wall varying (a) Rayleigh number (Ra), (b) Hartmann number (Ha), (c) Nanoparticles volume (ϕ) and (d) heat generation/absorption (Q) against the wave number, λ (= 1, 2, 3) taking as default values other parameters. Figure 9a, depicts that the average rate of temperature transfer increases for all values of buoyancy-driven parameter Ra. It is also seen that the temperature transfer rate is enhanced effectively with the augmentation of wave number λ. Therefore, the average rate of temperature transfer is more efficient for ascending values of the undulation parameter (λ). In addition, the enhancement of Ra promotes the buoyancy effect, therefore a large quantity of temperature is shifted into the enclosure from the horizontal warmed wall. Furthermore, the average rate of temperature transfer increases gradually by virtue of the corrugated surface of the enclosure. To get the required heat production or absorption properties, one can adjust the nanofluid design, nanoparticle type, concentration, and size. In order to maximize heat transfer and energy efficiency in a variety of industries, including chemical processes, microelectronics cooling, the oil and gas sector, solar thermal collectors, heat exchangers, and cooling systems, researchers are still investigating and creating novel nanofluids with particular qualities.
Average Nusselt number for (a) Rayleigh number, (b) Hartman number, (c) nanoparticles volume fraction and (d) Heat generation or absorption parameter.
Figure 9b depicts the mean Nusselt number (Nuav) reduced for the augmentation of Ha. For Ha, Nuav constantly turns higher within the enclosure. The Hartmann number (Ha) produces a strong Lorentz force which acts in opposition to the velocity fields and reduces the circulation strength into fluid within the cavity. Concerning the ascending as well as the collapsing model for wavy surface, the outcomes represent Nuav increases sufficiently for the improvement of the undulation parameter (λ). Figure 9c depicts that as nanoparticle volume (ϕ) rises, Nuav also rises monotonically. Therefore, the composition of nanoparticles in base fluid enhances Nuav. Figure 9d depicts heat generation or heat absorption rises for ascending values of λ. This figure also depicts that higher Nuav is observed for the heat sink. The temperature transport rate reduces for the rising heat source. In contrast, the temperature transfer rate is enhanced with the increase of heat sink. In total, the undulated surface enhances the mean Nusselt number (Nuav). In comparison to a cavity with a flat wall, a triangular cavity with a wavy wall can improve heat transmission, but it also creates issues with pressure drop, nanofluid stability, and flow pattern development. For nanofluid flow and heat transfer, a triangular hollow with a wavy or corrugated wall adds complexity and can significantly affect the flow and heat transfer properties. Heat transfer may be improved by the turbulence and vortices that the wavy wall geometry may produce in the nanofluid flow. Better thermal performance results from the wavy surface's acceleration of nanofluid mixing as compared to a triangular cavity with plain walls. Furthermore, because of the higher flow resistance caused by wavy walls, the pressure drop in the flow is also enhanced. The flow velocity and the particular size and waviness of the walls determine the effect on pressure drop. Moreover, the stability of the nanofluid may be impacted by the triangular enclosure with a wavy wall. Because of their propensity to aggregate or settle out of suspension, nanoparticles can affect the flow phenomena caused by irregular flow patterns. The flow domain's boundary layer growth can also be affected by the wavy wall geometry.
In order to investigate how Brownian motion influences the rate of heat transfer, the heat transport rate in regard to the average Nusselt number along the heated diameter is computed. Table 3 examines the effects of the Brownian effect, volume (ϕ), shape (n), and size (d) of the nanoparticles, with fixed values of λ = 3, Ha = 15, and Ra = 105. The result suggests that Brownian motion is a significant factor in the enhancement of the temperature transfer rate. The Brownian motion of the nanoparticles adds to the transmission of additional heat in the nanofluids and to the micro-convection of the fluid surrounding individual nanoparticles as a result of their movement into the base fluid. The outcome demonstrates that, for both scenarios without/with the Brownian motion of the nanoparticles, the effect of the volume fraction of the nanoparticles on the average Nusselt number is more effective for considering the Brownian effect. Furthermore, when taking into account the Brownian motion of the nanoparticles, the average heat transfer increases by 86.77% for 5% nanoparticle volume fractions while without the Brownian effect, it increases by only 5.52%.
This Table also shows that the average rate of heat transfer decreases with the increase of the diameter of nanoparticles. The small size of nanoparticles is more effective for the enhancement of thermal performance. The physical meaning of it, by decreasing the nanoparticle diameter, the specific area increases, which helps to enhance nanofluid thermal conductivity and consequently increases the average Nusselt number. The heat transport rate increases by 20.78% when the nanoparticle diameter reduces 100 nm to 20 nm, it increases by 35.87% when the nanoparticle size reduces 100 nm to 10 nm, and it rises by 141.37% when the nanoparticle diameter reduces 100 nm to 1 nm. The blade shape of nanoparticles generates the highest thermal performance. The heat transport rate increases by 8.66% for the blade shape of nanoparticles compared to the spherical shape of nanoparticles. In many different industries, including heat exchanger design, solar energy systems, chemical engineering, renewable energy systems, microfluidics, and so on, the study of heat generation/absorption effects on natural convective heat transfer in wavy triangular cavities filled with nanofluids has a wide range of practical applications. This research helps to improve energy efficiency, system performance, and the development of novel technologies.
The two-dimensional time independent natural convection laminar flow and temperature transport in a sinusoidally wavy triangle loaded by cupper-water nanofluid in the presence of heat generation/absorption parameter is investigated numerically. Galerkin weighted residual finite element technique is executed for numerical calculations. The impacts of several model parameters such as heat generation/absorption parameter (Q), Rayleigh number (Ra), number of undulations (λ), volume fraction of nanoparticles (ϕ), Hartmann number (Ha), nanoparticles shape (n), size (d) and Brownian motion on the flow and thermal field are presented using isotherm lines, streamline contours and average Nusselt number and interpreted them. The following major findings are listed:
The wavy surface has a significant impact on the flow domain as well as the temperature field. The ascending number of wave surfaces (λ = 3) augments temperature transfer sufficiently than the flat surface (λ = 0).
The strength of isotherms reduces successively for the rise of heat generation up to high heat generation (Q = 12). Consequently, higher temperature into nanofluid within the enclosure that reduces heat transfer rate.
By the increases of Rayleigh number and nanoparticles volume, fluid flow as well as temperature flow change remarkably. The higher Rayleigh number and lower Hartmann number correspond to convection heat transfer.
The small size of nanoparticles has a key role in thermal transport. The thermal transport rate increases by 35.87% when the nanoparticle diameter reduces from 100 to 10 nm.
A major contributing component to the improvement of heat transmission is Brownian motion. The heat transport rate increases by 20.43% for considering the Brownian effect for only 1% of nanoparticle volume.
The blade shape of nanoparticles confirms better thermal performance. The thermal transport increases by 8.66% compared to the default spherical shape of nanoparticles.
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Magnitude of magnetic field [kgs−2A−1]
Specific heat at constant pressure [Jkg−1K−1]
Length of the enclosure at bottom wall
Choi, S. U. S. Enhancing thermal conductivity of fluids with nanoparticles. ASME Fluids Eng. 231, 99–105 (1995).
Wong, K. V. & Leon, O. D. Applications of nanofluids: current and future. Adv. Mech. Eng. 2, 519–659 (2010).
Izadi, M., Mohebbi, R., Karimi, D. & Sheremet, M. A. Numerical simulation of natural convection heat transfer inside a ┴ shaped cavity filled by a MWCNT-Fe3O4/water hybrid nanofluids using LBM. Chem. Eng. Process. Process Intensif. 125, 56–66 (2018).
Alshuraiaan, B., Shahrestani, A. B. & Izadi, M. Numerical studies on passive parameters of a fluid-solid interaction problem derived by natural convection in a circular enclosure. Alex. Eng. J. 63, 415–426 (2023).
Izadi, M., Alshehri, H. M., Hosseinzadeh, F., Rad, M. S. & Hamida, M. B. B. Numerical study on forced convection heat transfer of TiO2/water nanofluid flow inside a double-pipe heat exchanger with spindle-shaped turbulators. Eng. Anal. Bound. Elem. 150, 612–623 (2023).
Xiong, Q. et al. State-of-the-art review of nanofluids in solar collectors: A review based on the type of the dispersed nanoparticles. J. Cleaner Prod. 310, 127528 (2021).
Izadi, M., Shahmardan, M. M., Behzadmehr, A., Rashidi, A. M. & Amrollahi, A. Modeling of effective thermal conductivity and viscosity of carbon structured nanofluid. Trans. Phenom. Nano Micro Scales 3(1), 1–13 (2015).
Aydin, O., Ünal, A. & Ayhan, T. Natural convection in rectangular enclosures heated from one side and cooled from the ceiling. Int. J. Heat Mass Transfer 42, 2345–2355 (1999).
Article CAS MATH Google Scholar
Holtzman, G. A., Hill, R. W. & Ball, K. S. Laminar natural convection in isosceles triangular enclosures heated from below and symmetrically cooled from above. J. Heat Transfer 122, 485–491 (2000).
Buongiorno, J. Convective transport in nanofluids. J. Heat Transfer 128, 240–250 (2006).
Das, S. K., Choi, S. U. S., Yu, W. & Pradeep, T. Nanofluids: Science and technology (Wiley, 2007).
Varol, Y., Oztop, H. F. & Yilmaz, T. Natural convection in triangular enclosures with protruding isothermal heater. Int. J. Heat Mass Transfer 50(13–14), 2451–2462 (2007).
Basak, T., Roy, S., Babu, S. K. & Balakrishnan, A. R. Finite element analysis of natural convection flow in an isosceles triangular enclosure due to uniform and non-uniform heating at the side walls. Int. J. Heat Mass Transfer 51(17–18), 4496–4506 (2008).
Article CAS MATH Google Scholar
Kent, E. F. Numerical analysis of laminar natural convection in isosceles triangular enclosures for cold base and hot inclined walls. Mech. Res. Commun. 36, 497–508 (2009).
Kamiyo, O. M. et al. Comprehensive review of natural convection in triangular enclosures. Appl. Mech. Rev. 63, 060801–060813 (2010).
Bhuvaneswari, M., Sivasankaran, S. & Kim, Y. J. Magneto-convection in a square enclosure with sinusoidal temperature distributions on both side walls. Numer. Heat Transfer A 59(3), 167–184 (2011).
Article ADS CAS Google Scholar
Tayebi, T., Djezzar, M. & Kheir, S. Effect of sinusoidal thermal boundary condition on natural convection in a cavity filled with Cu-water nanofluid. J. Nanofluids 2(2), 120–126 (2013).
Tayebi, T. & Chamkha, A. J. MHD buoyancy-driven flow in a nanoliquid filled-square enclosure divided by a solid conductive wall. Math. Meth. Appl. Sci. https://doi.org/10.1002/mma.6598 (2020).
Tayebi, T., Öztop, H. F. & Chamkha, A. J. MHD natural convection of a CNT-based nanofluid-filled annular circular enclosure with inner heat-generating solid cylinder. Eur. Phys. J. Plus. 136, 150 (2021).
Alsabery, A. I., Tayebi, T., Roslan, R., Chamkha, A. J. & Hashim, I. Entropy generation and mixed convection flow inside a wavy-walled enclosure containing a rotating solid cylinder and a heat source. Entropy 22, 606 (2020).
Article ADS MathSciNet PubMed PubMed Central Google Scholar
Alsabery, A. I., Tayebi, T., Chamkha, A. J. & Hashim, I. Effects of non-homogeneous nanofluid model on natural convection in a square cavity in the presence of conducting solid block and corner heater. Energies 11, 2507 (2018).
Hussain, S., Tayebi, T., Armaghani, T., Rashad, A. M. & Nabwey, H. A. Conjugate natural convection of non-Newtonian hybrid nanofluid in wavy-shaped enclosure. Appl. Math. Mech. 43, 447–466 (2022).
Article MathSciNet MATH Google Scholar
Nasrin, R. Joule heating effect on MHD combined convection in a wavy chamber having conducting square cylinder. J. Sci. Res. 4(1), 39–49 (2012).
Bhardwaj, S. & Dalal, A. Analysis of natural convection heat transfer and entropy generation inside right-angled triangular enclosure. Int. J. Heat Mass Transfer 65, 500–513 (2013).
Hamida, M. B. B. & Charrada, K. Heat generation/absorption effect on natural convection heat transfer in a square enclosure filled with a ethylene glycol-copper nanofluid under magnetic field. Am. J. Mod. Energy 1(1), 1–16 (2015).
Suvash, C. S. & Gu, Y. T. Natural convection in a triangular enclosure heated from below and non-uniformly cooled from top. Int. J. Heat Mass Transfer 80, 529–538 (2015).
Rahman, M. M., Alam, M. S., Al-Salti, N. & Eltayeb, I. A. Hydromagnetic natural convective heat transfer flow in an isosceles triangular cavity filled with nanofluid using two-component nonhomogeneous model. Int. J. Therm. Sci. 107, 272–288 (2016).
Mirabedin, S. M. CFD modeling of natural convection in right-angled triangular enclosures. Int. J. Heat Technol. 34(3), 503–506 (2016).
Alam, M. S., Islam, T. & Uddin, M. J. Mathematical modelling for heat transfer of a micropolar fluid along a permeable stretching/shrinking wedge with heat generation/absorption. Math. Model. Eng. Probl. 3(1), 1–9 (2016).
Triveni, M., Panua, R. & Sen, D. Effects of variant positions of cold walls on natural convection in a triangular cavity. J. Appl. Fluid Mech. 9(1), 185–193 (2016).
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. Comput. Appl. 143(6), 0975–8887 (2016).
Uddin, M. J., Hoque, A. K. M. F., Rahman, M. M. & Vajravelu, K. Numerical simulation of convective heat transport within the nanofluid filled vertical tube of plain and uneven side walls. Int. J. Thermofluid Sci. Technol. 6(4), 19060101 (2019).
Alam, E. M. M., Sumona, H., Uddin, M. S. & Rahman, M. M. Effect of sinusoidal thermal boundary condition on unsteady magnetohydrodynamics convection in a square enclosure filled with Fe3O4-water Ferrofluid. Int. J. Stat. Appl. Math. 4(6), 111–127 (2019).
Al-Asad, M. F., Munshi, M. J. H., Bhowmik, R. K. & Sarker, M. M. A. MHD free convection heat transfer having vertical fin in a square wavy cavity. Int. J. Stat. Appl. Math. 4(3), 32–38 (2019).
Zahan, I., Nasrin, R. & Alim, M. A. Mixed convection hybrid Nanofluid flow in lid-driven undulated cavity: Effect of MHD and Joule heating. J. Naval Arch. Mar. Eng. 16, 109–126 (2019).
Mourad, A., Abderrahmane, A., Younis, O., Marzouki, R. & Alazzam, A. Numerical simulations of magnetohydrodynamics natural convection and entropy production in a porous annulus bounded by wavy cylinder and koch snowflake loaded with Cu–water nanofluid. Micromachines 13, 182 (2022).
Article PubMed PubMed Central Google Scholar
Abderrahmane, A. et al. Heat and mass transfer analysis of non-Newtonian power-law nanofluid confined within annulus enclosure using Darcy-Brinkman-Forchheimer model. Case Stud. Therm. Eng. 40, 102569 (2022).
Alshare, A. et al. Hydrothermal and entropy investigation of nanofluid natural convection in a lid-driven cavity concentric with an elliptical cavity with a wavy boundary heated from below. Nanomaterials 12(9), 1392 (2022).
Article CAS PubMed PubMed Central Google Scholar
Izadi, M. et al. Transient magneto-buoyant convection of a magnetizable nanofluid inside a circle sensible storage subjected to double time-dependent thermal sources. J. Therm. Anal. Calorim. 148, 8511–8531 (2023).
Izadi, M., Shahmardan, M. M., Maghrebi, M. J. & Behzadmehr, A. The developed laminated mixed convolution of Al2O3/water nanofluid in an annulus was studied numerically. Chem. Eng. Commun. 200(7), 878–894 (2013).
Izadi, M. Effects of porous material on transient natural convection heat transfer of nanofluids inside a triangular chamber. Chin. J. Chem. Eng. 28(5), 1203–1213 (2018).
Izadi, M., Hajjar, A., Alshehri, H. M., Sheremet, M. & Galal, A. M. Charging process of a partially heated trapezoidal thermal energy storage filled by nano-enhanced PCM using controllable uniform magnetic field. Int. Commun. Heat Mass Transfer 138, 106349 (2022).
Izadi, M., Behzadmehr, A. & Shahmardan, M. M. Effects of inclination angle on mixed convection heat transfer of a nanofluid in a square cavity. Int. J. Comput. Methods Eng. Sci. Mech. 16(1), 11–21 (2015).
Article MathSciNet CAS Google Scholar
Izadi, M., Ghalambaz, M. & Mehryan, S. A. M. Location impact of a pair of magnetic sources on melting of a magneto-Ferro phase change substance. Chin. J. Phys. 65, 377–388 (2020).
Izadi, M., Bastani, B. & Sheremet, M. A. Numerical simulation of thermogravitational energy transport of a hybrid nanoliquid within a porous triangular chamber using the two-phase mixture approach. Adv. Powder Technol. 31(6), 2493–2504 (2020).
Sajjadi, H., Delouei, A. A., Mohebbi, R., Izadi, M. & Succi, S. Natural convection heat transfer in a porous cavity with sinusoidal temperature distribution using Cu/Water nanofluid: Double MRT lattice boltzmann method. Commun. Comput. Phys. 29, 292–318 (2021).
Article MathSciNet MATH Google Scholar
Izadi, M. et al. Numerical investigation of magneto-thermal-convection impact on phase change phenomenon of Nano-PCM within a hexagonal shaped thermal energy storage. Appl. Therm. Eng. 223, 119984 (2023).
Izadi, M., Mohebbi, R., Delouei, A. A. & Sajjadi, H. Natural convection of a magnetizable hybrid nanofluid inside a porous enclosure subjected to two variable magnetic fields. Int. J. Mech. Sci. 151, 154–169 (2019).
Izadi, M., Alshuraiaan, B., Hajjar, A., Sheremet, M. A. & Hamida, M. B. B. Free convection of nanofluids in a porous sensible heat storage unit: Combined effect of time-periodic heating and external magnetic field. Int. J. Therm. Sci. 192, 108404 (2023).
Islam, T., Alam, M. N., Asjad, M. I., Parveen, N. & Chu, Y. M. Heatline visualization of MHD natural convection heat transfer of nanofluid in a prismatic enclosure. Sci. Reports 11, 10972 (2021).
Islam, T., Parveen, N. & Nasrin, R. Mathematical modeling of unsteady flow with uniform/non-uniform temperature and magnetic intensity in a half-moon shaped domain. Heliyon 8(3), e09015 (2022).
Article PubMed PubMed Central Google Scholar
Islam, T., Yavuz, M., Parveen, N. & Al-Asad, M. F. Impact of non-uniform periodic magnetic field on unsteady natural convection flow of nanofluids in square enclosure. Fractal Fract. 6(2), 101 (2022).
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 18, 100338 (2023).
Islam, S. et al. Heatline visualization in hydromagnetic natural convection flow inside a prismatic heat exchanger using nanofluid. Int. J. Thermofluids 16(2022), 100248 (2022).
Islam, T. & Nasrin, R. Thermal operation by nanofluids with various aspects: A comprehensive numerical appraisal. Waves in Random and Complex Media (2022).
Khatun, M. A. & Islam, T. Influence of magnetic field and heat generation/absorption on unsteady MHD convective flow along a permeable stretching/shrinking wedge with thermophoresis and variable fluid properties. Int. J. Thermofluids 16, 100204 (2022).
Al-Kalbani, K. S. & Rahman, M. M. Convective heat transfer in the flow of a nanofluid in an inclined square enclosure. J. Eng. Phys. Thermophys. 92(5), 1150–1170 (2019).
Al-Weheibi, S. M., Rahman, M. M., Alam, M. S. & Vajravelu, K. Numerical simulation of natural convection heat transfer in a trapezoidal enclosure filled with nanoparticles. Int. J. Mech. Sci. 131–132, 599–612 (2017).
Zienkiewicz, O.C. & Taylor, R.L. The finite element method, 4th edn. (McGraw-Hill, 1991).
Ghasemi, B., Aminossadati, S. M. & Raisi, A. Magnetic field effect on natural convection in a nanofluid filled square enclosure. Int. J. Therm. Sci. 50(9), 1748–1756 (2011).
The authors extend their appreciation to the Ministry of Education in KSA for funding this work through Project Number KKU-IFP2-H-21.
CMUP, University of Porto, Porto, Portugal
Department of Mathematics, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj, 8100, Bangladesh
Department of Mathematics, Education Faculty, Laghman University, Mehtarlam City, Laghman, 2701, Afghanistan
Department of Mathematics, College of Education, Majmmah University, 11952, Al-Majmaah, Saudi Arabia
Bangladesh University of Engineering and Technology (BUET), Dhaka, Bangladesh
College of Engineering Mechanical Engineering Department, King Khalid University, Abha, Saudi Arabia
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
T.I.: Model and analysis; M.N.A.; Problem formulation; S.N.: results and discussion; I.K.: plotting and numerical simulations; M.D.F.: problem solution; S.A. Computed results of average Nusselt number, Results Discussion, conclusion, revision.
Correspondence to Shafiullah Niazai or Ilyas Khan.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Islam, T., Alam, M.N., Niazai, S. et al. Heat generation/absorption effect on natural convective heat transfer in a wavy triangular cavity filled with nanofluid. Sci Rep 13, 21171 (2023). https://doi.org/10.1038/s41598-023-48704-2
DOI: https://doi.org/10.1038/s41598-023-48704-2
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
Journal of Thermal Analysis and Calorimetry (2024)
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
LiBr Absorption Heat Pump Sign up for the Nature Briefing: AI and Robotics newsletter — what matters in AI and robotics research, free to your inbox weekly.