
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38040956
48704
10.1038/s41598-023-48704-2
Article
Heat generation/absorption effect on natural convective heat transfer in a wavy triangular cavity filled with nanofluid
Islam Tarikul 1
Alam Md. Nur 2
Niazai Shafiullah shafiullahniazai@lu.edu.af

3
Khan Ilyas i.said@mu.edu.sa

4
Fayz-Al-Asad Md. 5
Alqahtani Sultan 6
1 https://ror.org/043pwc612 grid.5808.5 0000 0001 1503 7226 CMUP, University of Porto, Porto, Portugal
2 https://ror.org/011xjpe74 grid.449329.1 0000 0004 4683 9733 Department of Mathematics, Bangabandhu Sheikh Mujibur Rahman Science and Technology University, Gopalganj, 8100 Bangladesh
3 Department of Mathematics, Education Faculty, Laghman University, Mehtarlam City, Laghman 2701 Afghanistan
4 https://ror.org/01mcrnj60 grid.449051.d 0000 0004 0441 5633 Department of Mathematics, College of Education, Majmmah University, 11952 Al-Majmaah, Saudi Arabia
5 grid.411512.2 0000 0001 2223 0518 Bangladesh University of Engineering and Technology (BUET), Dhaka, Bangladesh
6 https://ror.org/052kwzs30 grid.412144.6 0000 0004 1790 7100 College of Engineering Mechanical Engineering Department, King Khalid University, Abha, Saudi Arabia
1 12 2023
1 12 2023
2023
13 2117127 8 2023
29 11 2023
© The Author(s) 2023
https://creativecommons.org/licenses/by/4.0/ 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/.
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.

Subject terms

Engineering
Mathematics and computing
Nanoscience and technology
issue-copyright-statement© Springer Nature Limited 2023
==== Body
pmcIntroduction

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–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–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.

Problem formulation

Physical model

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(π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 (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 B=(B0,0,0) can be written as follows:Figure 1 Schematic view for sinusoidal triangular enclosure.

Table 1 Thermophysical properties of the base fluid and solid particles (see57,58).

Base fluid/nanoparticles	cp [Jkg−1K−1]	ρ [kgm−3]	k [Wm−1K−1]	μ [kgm−1s−1]	β × 10−5 [K−1]	σ [Sm−1]	Pr	
H2O	4179	997.1	0.613	0.001003	21	5.5 × 10−6	6.8377	
Cu	385	8933	400	–	1.67	5.96 × 107	–	

1 σnf(u×B)×B=Fm=(0,-σnfB02v,0)

The feature of the vertical wavy wall is employed to mimic the following model: 2 X=A(1-cos(2πλx))

where A represents wave amplitude and λ represents undulation number.

Mathematical model

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:3 ∂u∂x+∂v∂y=0

4 u∂u∂x+v∂u∂y=-1ρnf∂p∂x+μnfρnf∂2u∂x2+∂2u∂y2

5 ρnfu∂v∂x+v∂v∂y=-∂p∂y+μnf∂2v∂x2+∂2v∂y2+ρβnfg(T-Tc)-σB02ν

6 u∂T∂x+v∂T∂y=αnf∂2T∂x2+∂2T∂y2+q(ρcp)nfT-Tc

Boundary conditions (dimensional)

7 At the horizontal sidewall:u=0,v=0,T=sin(πx/L)

8 At the vertical wavy sidewall:u=0,v=0,T=Tc

9 At the inclined wall:u=0,v=0,∂T∂x=0

Thermophysical properties of nanofluids

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:10 αnf=knfρcpnf

11 ρcpnf=1-ϕρcpbf+ϕρcpsp

12 ρβnf=1-ϕρβbf+ϕρβsp

13 μnf=μbf1-ϕ-2.5

14 ρnf=1-ϕρbf+ϕρsp

15 σnf=σsp+2σbf-2ϕσbf-σspσsp+2σbf+2ϕσbf-σspσbf

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:16 knf=kstatic+kBrownian

where,kstatic=ksp+(n-1)kbf-(n-1)(kbf-ksp)ϕksp+(n-1)kbf+(kbf-ksp)ϕkbf,kBrownian=ϕρpcp,sp22KBTref3πdμnf

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:17 X=xL,Y=yL,U=uLαbf,V=vLαbf,P=pL2ρnfαbf2,andθ=T-TcTh-Tc

The non-dimensionless governing equations by employing dimensionless variables are as follows57,58:18 ∂U∂X+∂U∂Y=0

19 U∂U∂X+V∂U∂Y=-ρbfρnf∂P∂X+Prρbfρnf∂2U∂X2+∂2U∂Y2

20 U∂V∂X+V∂V∂Y=-ρbfρnf∂P∂Y+Prρbfρnf∂2V∂X2+∂2V∂Y2+ρβnfρnfβbfRaPrθ-ρbfρnfσnfσbfHa2PrV

21 U∂θ∂X+V∂θ∂Y=αnfαbf∂2θ∂X2+∂2θ∂Y2+αnfαbfQθ

where, Pr=υbfαbf is the Prandtl number, Ra=gβbfTh-TcL3υbfαbf is the Rayleigh number, Ha=BoLσbf/μbf is the Hartmann number, and Q=qL2αnfρcpnf is the heat generation /absorption parameter.

Non-dimensional boundary conditions are as follows:22 Atthehorizontalwall:U=0,V=0,θ=sin(πX)

22 Attheverticalwavywall:U=0,V=0,θ=0

23 Attheinclinedwall:U=0,V=0,∂θ∂X=0

Average Nusselt number

For this model, average Nusselt number on horizontal warmed wall has been calculated from the following expression57,58:24 Nuav=-knfkbf∫01∂θ∂YdX

Computational procedure

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 Υm+1-Υm≤10-5, where Υ represents subsidiary variables U,V,θ and m represent iteration number.

Grid independence test

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.Figure 2 Average Nusselt number (Nuav) of different element numbers.

Code validation

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.Table 2 Code authentication with Ghasemi et al.60 using of average Nusselt number when n = 3, Pr = 6.8377, and Ha = 30 for Al2O3–H2O nanofluid.

	Average Nusselt number (Nuav)	
ϕ = 0	ϕ = 0.02	ϕ = 0.04	
Ghasemi et al.60	Present study	Relative error (%)	Ghasemi et al.60	Present study	Relative error (%)	Ghasemi et al.60	Present study	Relative error (%)	
103	1.002	1.002	0	1.060	1.072	1.13	1.121	1.121	0	
104	1.183	1.182	0.08	1.212	1.235	1.90	1.249	1.269	1.60	
105	3.150	3.138	0.38	3.138	3.201	2.01	3.124	3.216	2.94	
106	7.907	7.815	1.16	7.979	8.018	0.49	8.042	8.087	0.56	

Results and discussion

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.

Impact of Rayleigh number

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.Figure 3 Outcomes of Rayleigh number (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.Figure 4 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.

Impact of Hartmann number

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.Figure 5 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 (λ).Figure 6 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.

Impact of heat generation/absorption

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.Figure 7 Outcomes of Q (< 0) for various undulation parameter (λ) using isotherms when ϕ = 0.04, Ha = 15, Ra = 105, n = 3, and d = 10 nm.

Figure 8 Outcomes of Q (> 0) for various undulation parameter (λ) using isotherms when ϕ = 0.04, Ha = 15, Ra = 105, n = 3, and d = 10 nm.

Effect of average Nusselt number

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.Figure 9 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.

Effect of nanoparticles shape, size and Brownian motion

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%.Table 3 Variation of average Nusselt number (Nuav) along heated wall for different nanoparticles volume (ϕ), diameter (d) and shape (n) of nanoparticles considering with/without Brownian with Ra = 105, Ha = 15, and λ = 3.

Brownian effect (BM)	Nanoparticles volume (ϕ)	Nanoparticles shape (n)	Nanoparticles diameter (d)	Average Nusselt number (Nuav)	Increases (%)	
with BM	0	3.0	10 nm	4.267974931830414	–	
0.01	5.140180145823340	20.43	
0.02	5.930500451043041	38.95	
0.05	7.971343224154655	86.77	
0.04	3.0	10 nm	7.254760603941890	–	
3.7	7.334404774458736	1.10	
4.9	7.470160718151543	2.97	
5.7	7.560131307489343	4.21	
8.6	7.882830676803941	8.66	
0.04	3.0	100 nm	5.398074838809371	–	
50 nm	5.779106109891497	7.06	
20 nm	6.519641772718621	20.78	
10 nm	7.334404744750127	35.87	
1 nm	13.02908943936327	141.37	
Without BM	0	3.0	–	4.267974931830413	–	
0.01	–	4.309365483763966	0.97	
0.02	–	4.353923846604179	2.01	
0.05	–	4.503711767742641	5.52	
0.04	3.0	–	4.362500816778986	–	
3.7	–	4.451326343833879	2.03	
4.9	–	4.602177643914426	5.49	
5.7	–	4.701783415723639	7.78	
8.6	–	5.056793437004171	15.92	

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.

Conclusions

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.

List of symbols

A Wave amplitude [m]

B0 Magnitude of magnetic field [kgs−2A−1]

cp Specific heat at constant pressure [Jkg−1K−1]

g Gravitational acceleration [ms−2]

Ha Hartmann number

k Thermal conductivity [Wm−1K−1]

L Length of the enclosure at bottom wall

Nuav Average Nusselt number

p Dimensional pressure [kgm−1s−2]

P Dimensionless pressure

Pr Prandtl number

Ra Rayleigh number

T Fluid temperature [K]

u, v Dimensional velocity components [ms−1]

U, V Dimensionless velocity components

x, y Dimensional coordinates [m]

X, Y Dimensionless coordinates

Greek letters

α Thermal diffusivity [m2s−1]

β Thermal expansion coefficient [K−1]

ϕ Volume fraction of nanoparticles

μ Dynamic viscosity [kgm−1s−1]

υ Kinematic viscosity [m2s−1]

θ Non-dimensional temperature

ρ Density [kgm−3]

σ Electric conductivity

λ Undulation number

Subscripts

h Hot surface

c Cold surface

nf Nanofluid

sp Solid particle

bf Base fluid

Acknowledgements

The authors extend their appreciation to the Ministry of Education in KSA for funding this work through Project Number KKU-IFP2-H-21.

Author contributions

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.

Data availability

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Choi SUS Enhancing thermal conductivity of fluids with nanoparticles ASME Fluids Eng. 1995 231 99 105
2. Wong KV Leon OD Applications of nanofluids: current and future Adv. Mech. Eng. 2010 2 519 659 10.1155/2010/519659
3. Izadi M Mohebbi R Karimi D Sheremet MA 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. 2018 125 56 66 10.1016/j.cep.2018.01.004
4. Alshuraiaan B Shahrestani AB Izadi M Numerical studies on passive parameters of a fluid-solid interaction problem derived by natural convection in a circular enclosure Alex. Eng. J. 2023 63 415 426 10.1016/j.aej.2022.07.057
5. Izadi M Alshehri HM Hosseinzadeh F Rad MS Hamida MBB 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. 2023 150 612 623 10.1016/j.enganabound.2023.02.046
6. Xiong Q Hajjar A Alshuraiaan B Izadi M Altnji S Shehzad SA State-of-the-art review of nanofluids in solar collectors: A review based on the type of the dispersed nanoparticles J. Cleaner Prod. 2021 310 127528 10.1016/j.jclepro.2021.127528
7. Izadi M Shahmardan MM Behzadmehr A Rashidi AM Amrollahi A Modeling of effective thermal conductivity and viscosity of carbon structured nanofluid Trans. Phenom. Nano Micro Scales 2015 3 1 1 13
8. 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 1999 42 2345 2355 10.1016/S0017-9310(98)00319-6
9. Holtzman GA Hill RW Ball KS Laminar natural convection in isosceles triangular enclosures heated from below and symmetrically cooled from above J. Heat Transfer 2000 122 485 491 10.1115/1.1288707
10. Buongiorno J Convective transport in nanofluids J. Heat Transfer 2006 128 240 250 10.1115/1.2150834
11. Das SK Choi SUS Yu W Pradeep T Nanofluids: Science and technology 2007 Wiley
12. Varol Y Oztop HF Yilmaz T Natural convection in triangular enclosures with protruding isothermal heater Int. J. Heat Mass Transfer 2007 50 13–14 2451 2462 10.1016/j.ijheatmasstransfer.2006.12.027
13. Basak T Roy S Babu SK Balakrishnan AR 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 2008 51 17–18 4496 4506 10.1016/j.ijheatmasstransfer.2007.12.018
14. Kent EF Numerical analysis of laminar natural convection in isosceles triangular enclosures for cold base and hot inclined walls Mech. Res. Commun. 2009 36 497 508 10.1016/j.mechrescom.2008.11.002
15. Kamiyo OM Angeli D Barozzi GS Collins MW Olunloyo VOS Talabi SOA Comprehensive review of natural convection in triangular enclosures Appl. Mech. Rev. 2010 63 060801 60813 10.1115/1.4004290
16. Bhuvaneswari M Sivasankaran S Kim YJ Magneto-convection in a square enclosure with sinusoidal temperature distributions on both side walls Numer. Heat Transfer A 2011 59 3 167 184 10.1080/10407782.2011.541219
17. 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 2013 2 2 120 126 10.1166/jon.2013.1044
18. Tayebi T Chamkha AJ MHD buoyancy-driven flow in a nanoliquid filled-square enclosure divided by a solid conductive wall Math. Meth. Appl. Sci. 2020 10.1002/mma.6598
19. Tayebi T Öztop HF Chamkha AJ MHD natural convection of a CNT-based nanofluid-filled annular circular enclosure with inner heat-generating solid cylinder Eur. Phys. J. Plus. 2021 136 150 10.1140/epjp/s13360-021-01106-7
20. Alsabery AI Tayebi T Roslan R Chamkha AJ Hashim I Entropy generation and mixed convection flow inside a wavy-walled enclosure containing a rotating solid cylinder and a heat source Entropy 2020 22 606 10.3390/e22060606 33286378
21. Alsabery AI Tayebi T Chamkha AJ 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 2018 11 2507 10.3390/en11102507
22. Hussain S Tayebi T Armaghani T Rashad AM Nabwey HA Conjugate natural convection of non-Newtonian hybrid nanofluid in wavy-shaped enclosure Appl. Math. Mech. 2022 43 447 466 10.1007/s10483-022-2837-6
23. Nasrin R Joule heating effect on MHD combined convection in a wavy chamber having conducting square cylinder J. Sci. Res. 2012 4 1 39 49 10.3329/jsr.v4i1.8014
24. Bhardwaj S Dalal A Analysis of natural convection heat transfer and entropy generation inside right-angled triangular enclosure Int. J. Heat Mass Transfer 2013 65 500 513 10.1016/j.ijheatmasstransfer.2013.06.020
25. Hamida MBB 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 2015 1 1 1 16
26. Suvash CS Gu YT Natural convection in a triangular enclosure heated from below and non-uniformly cooled from top Int. J. Heat Mass Transfer 2015 80 529 538 10.1016/j.ijheatmasstransfer.2014.09.047
27. Rahman MM Alam MS Al-Salti N Eltayeb IA Hydromagnetic natural convective heat transfer flow in an isosceles triangular cavity filled with nanofluid using two-component nonhomogeneous model Int. J. Therm. Sci. 2016 107 272 288 10.1016/j.ijthermalsci.2016.04.009
28. Mirabedin SM CFD modeling of natural convection in right-angled triangular enclosures Int. J. Heat Technol. 2016 34 3 503 506 10.18280/ijht.340322
29. Alam MS Islam T Uddin MJ Mathematical modelling for heat transfer of a micropolar fluid along a permeable stretching/shrinking wedge with heat generation/absorption Math. Model. Eng. Probl. 2016 3 1 1 9 10.18280/mmep.030101
30. Triveni M Panua R Sen D Effects of variant positions of cold walls on natural convection in a triangular cavity J. Appl. Fluid Mech. 2016 9 1 185 193 10.18869/acadpub.jafm.68.224.24002
31. 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. 2016 143 6 0975 8887
32. Uddin MJ Hoque AKMF Rahman MM 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. 2019 6 4 19060101
33. Alam EMM Sumona H Uddin MS Rahman MM Effect of sinusoidal thermal boundary condition on unsteady magnetohydrodynamics convection in a square enclosure filled with Fe3O4-water Ferrofluid Int. J. Stat. Appl. Math. 2019 4 6 111 127
34. Al-Asad MF Munshi MJH Bhowmik RK Sarker MMA MHD free convection heat transfer having vertical fin in a square wavy cavity Int. J. Stat. Appl. Math. 2019 4 3 32 38
35. Zahan I Nasrin R Alim MA Mixed convection hybrid Nanofluid flow in lid-driven undulated cavity: Effect of MHD and Joule heating J. Naval Arch. Mar. Eng. 2019 16 109 126 10.3329/jname.v16i2.40585
36. 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 2022 13 182 10.3390/mi13020182 35208306
37. Abderrahmane A Jamshed W Abed AM Smaisim GF Guedri K Devi SU Akbari OA Younis O Baghae S Heat and mass transfer analysis of non-Newtonian power-law nanofluid confined within annulus enclosure using Darcy-Brinkman-Forchheimer model Case Stud. Therm. Eng. 2022 40 102569 10.1016/j.csite.2022.102569
38. Alshare A Abderrahmane A Guedri K Younis O Al-Asad MF Ali HM Al-Kouz W 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 2022 12 9 1392 10.3390/nano12091392 35564102
39. Izadi M Tayebi T Alshehri HM Hajjar A Hamida MBB Galal AM Transient magneto-buoyant convection of a magnetizable nanofluid inside a circle sensible storage subjected to double time-dependent thermal sources J. Therm. Anal. Calorim. 2023 148 8511 8531 10.1007/s10973-023-12242-w
40. Izadi M Shahmardan MM Maghrebi MJ Behzadmehr A The developed laminated mixed convolution of Al2O3/water nanofluid in an annulus was studied numerically Chem. Eng. Commun. 2013 200 7 878 894 10.1080/00986445.2012.723077
41. Izadi M Effects of porous material on transient natural convection heat transfer of nanofluids inside a triangular chamber Chin. J. Chem. Eng. 2018 28 5 1203 1213 10.1016/j.cjche.2020.01.021
42. Izadi M Hajjar A Alshehri HM Sheremet M Galal AM 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 2022 138 106349 10.1016/j.icheatmasstransfer.2022.106349
43. Izadi M Behzadmehr A Shahmardan MM Effects of inclination angle on mixed convection heat transfer of a nanofluid in a square cavity Int. J. Comput. Methods Eng. Sci. Mech. 2015 16 1 11 21 10.1080/15502287.2014.976674
44. Izadi M Ghalambaz M Mehryan SAM Location impact of a pair of magnetic sources on melting of a magneto-Ferro phase change substance Chin. J. Phys. 2020 65 377 388 10.1016/j.cjph.2020.03.002
45. Izadi M Bastani B Sheremet MA Numerical simulation of thermogravitational energy transport of a hybrid nanoliquid within a porous triangular chamber using the two-phase mixture approach Adv. Powder Technol. 2020 31 6 2493 2504 10.1016/j.apt.2020.04.011
46. Sajjadi H Delouei AA 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. 2021 29 292 318 10.4208/cicp.OA-2020-0001
47. Izadi M Sheremet M Hajjar A Galal AM Mahariq I Jarad F Hamida MBB Numerical investigation of magneto-thermal-convection impact on phase change phenomenon of Nano-PCM within a hexagonal shaped thermal energy storage Appl. Therm. Eng. 2023 223 119984 10.1016/j.applthermaleng.2023.119984
48. Izadi M Mohebbi R Delouei AA Sajjadi H Natural convection of a magnetizable hybrid nanofluid inside a porous enclosure subjected to two variable magnetic fields Int. J. Mech. Sci. 2019 151 154 169 10.1016/j.ijmecsci.2018.11.019
49. Izadi M Alshuraiaan B Hajjar A Sheremet MA Hamida MBB 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. 2023 192 108404 10.1016/j.ijthermalsci.2023.108404
50. Islam T Alam MN Asjad MI Parveen N Chu YM Heatline visualization of MHD natural convection heat transfer of nanofluid in a prismatic enclosure Sci. Reports 2021 11 10972
51. 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 2022 8 3 e09015 10.1016/j.heliyon.2022.e09015 35265763
52. Islam T Yavuz M Parveen N Al-Asad MF Impact of non-uniform periodic magnetic field on unsteady natural convection flow of nanofluids in square enclosure Fractal Fract. 2022 6 2 101 10.3390/fractalfract6020101
53. Saha T Islam T Yeasmin S Parveen N Thermal influence of heated fin on MHD natural convection flow of nanofluids inside a wavy square cavity Int. J. Thermofluids 2023 18 100338 10.1016/j.ijft.2023.100338
54. Islam S Bairagi T Islam T Rana BMJ Reza-E-Rabbi SK Rahman MM Heatline visualization in hydromagnetic natural convection flow inside a prismatic heat exchanger using nanofluid Int. J. Thermofluids 2022 16 2022 100248 10.1016/j.ijft.2022.100248
55. Islam, T. & Nasrin, R. Thermal operation by nanofluids with various aspects: A comprehensive numerical appraisal. Waves in Random and Complex Media (2022).
56. Khatun MA 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 2022 16 100204 10.1016/j.ijft.2022.100204
57. Al-Kalbani KS Rahman MM Convective heat transfer in the flow of a nanofluid in an inclined square enclosure J. Eng. Phys. Thermophys. 2019 92 5 1150 1170 10.1007/s10891-019-02030-5
58. Al-Weheibi SM Rahman MM Alam MS Vajravelu K Numerical simulation of natural convection heat transfer in a trapezoidal enclosure filled with nanoparticles Int. J. Mech. Sci. 2017 131–132 599 612 10.1016/j.ijmecsci.2017.08.005
59. Zienkiewicz, O.C. & Taylor, R.L. The finite element method, 4th edn. (McGraw-Hill, 1991).
60. Ghasemi B Aminossadati SM Raisi A Magnetic field effect on natural convection in a nanofluid filled square enclosure Int. J. Therm. Sci. 2011 50 9 1748 1756 10.1016/j.ijthermalsci.2011.04.010
