Natural convection heat transfer in corrugated annuli with H2O-Al2O3 nanofluid.

The present work includes a numerical study of natural convection heat transfer in symmetrical and unsymmetrical corrugated annuli filled with H2O-Al2O3 nanofluid. In this study, higher and lower temperatures were kept constant at inner and outer cylinders of the annulus; respectively. Eight mathematical models with an aspect ratio of 1.5 were developed to find the best model giving the highest heat transfer rates. The stream-vorticity formulation in curvilinear coordinates was used to solve the governing equations of heat transfer and fluid motion. The influences of Rayleigh number. ( ( 10 3 ≤ Ra ≤ 10 6 ) and volume fraction of nanoparticles. ( 0 ≤ ϕ ≤ 0.25 ) on isotherms, streamlines, local and average Nusselt numbers on the inner and outer cylinder were investigated. The results show that the heat transfer rate is significantly increased with an increase in nanoparticles volume fraction and Rayleigh number. The activity of the heated surface is increased with an increase in the undulation number, but the flow motion tends to be most difficult in the spaces between active undulation walls. Moreover, the heat transfer rates in unsymmetrical annuli are relatively higher than the rates in the symmetrical annuli. There are no evident changes in isotherms with an increase in the nanofluid volume fraction. Correlations for the mean Nusselt number on the inner and outer walls of annulus were deduced as a function of Rayleigh number and nanoparticles volume fraction for eight models with an accuracy range of 8-15 %.

Introduction

Natural convective heat transfer characteristics in enclosures have a great significance in industrial applications such as solar collectors, thermal storage systems, cooling of nuclear reactors, etc. In past years, different methods were used to augment the heat transfer in enclosures such as using porous media [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15] and nanotechnology [16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43]. The corrugated annulus is a type of enhanced heat transfer techniques. Corrugation aims to raise the secondary currents by producing radial velocity components and high flow mixing. As a result, the increase in the wet perimeter with remaining the cross-section area constant leads to an increase in the convective surface area. Therefore, it was extensively used in modern heat exchangers applicated in steam generation, condensation in power and cogeneration plants, chemical and agricultural products, automobile industry, etc. [44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58].

Hasanuzzamana et al. 2007 [44], studied the natural convection in a square cavity having two V-wavy parallel vertical plates and two flat horizontal plates. It was noticed that the average Nusselt number increased with the increase in the input heat energy. Noman et al. 2008 [45], concluded that the natural convection phenomena in a square enclosure with a hot sinusoidal corrugated top surface are affected by changing the inclination angles. Hasan et al. (2012) [46], noticed that the transient phenomena of natural convection in a square enclosure are affected by Rayleigh Number, corrugation amplitude, and frequency. Salam et al. 2013 [47], studied the inclination effect of hot thin plate located concentrically in an enclosure on the behavior of thermal field and fluid flow. The results showed that the heat transfer rate increased with increase in the corrugation frequency and the inclination angle of plate. Souad and Amina 2015 [48], found that the heat transfer rate in a square cavity decreased, and the total entropy generation raised with the increase in Bijan number and the amplitude of the hot wavy right wall. Chordiya and Sharma 2019 [49], studied the effect of shape, position, amplitude, and inclination angle of corrugated diathermal partitions in a porous cavity on the thermal distribution and fluid flow. They concluded that, for higher Rayleigh numbers, the existence of corrugated partitions led to the increase in heat transfer rate by about 61%. Behrouz and Saeed 2015 [50], employed H2O, Al2O3-H2O nanofluid, and Al2O3-Cu-H2O hybrid nanofluid to examine the natural convection heat transfer in a cavity having sinusoidal corrugated vertical walls. It was observed that the enhancement of heat transfer increased at higher Rayleigh numbers and volume concentrations. Saha 2017 [51], used H2O-TiO2 nanofluid to enhance the heat transfer inside a cavity with vertical sinusoidal walls and adiabatic horizontal flat walls. The results revealed that the heat transfer rate using TiO2 nanofluid was higher compared with using base fluid only. Mitchell 2017 [52], analyzed the heat transfer by natural convection inside a wavy enclosure containing Al2O3-water nanofluid. They concluded that the heat transfer rate increased with the increase in surface waviness, aspect ratio, and nanoparticles volume fraction. Soroush et al. 2017 [53], studied the turbulent free convection, conduction, and surface thermal radiation in a rectangular cavity with different shapes of heat source located at the bottom wall. It was noticed that employing the circular obstacles for heat source gave a higher heat transfer rate than other shapes only in low emissivity values. Walid 2018 [54], analyzed the heat transfer and fluid flow in a three-dimensional triangular solar collector with a hot corrugated bottom wall. It was noticed that the flow strength and the heat transfer rate at the heated surface increased with the increase in Rayleigh number. Ammar et al. 2019) [55], used various undulation numbers of a hot cylinder located inside a cooled wavy-walled enclosure filled with Ag nanofluid layer on the right half and a fully saturated porous media with Ag nanofluid on the left half. It was found that the average Nusselt number decreased with increase in the thickness of the porous layer. Rujda and Mahapatra 2019 [56], studied the effect of magnetic field on the natural convection and mass transfer in a wavy-walled cavity filled with water-based Al2O3 nanofluid. It was observed that increasing the undulations number and buoyancy ratio led to reducing the performance of heat and mass transfer rate. Dutta et al. 2020 [57], studied the natural convection in a porous corrugated rhombic cavity. It was noticed that the heat transfer process was indicated by the interaction between Rayleigh number and Darcy number together with a number of undulations. Abeer 2020 [58], enhanced the natural convection inside a square porous cavity by 10% with an increase in nanoparticle concentration. It was concluded that the average Nusselt number of the corrugated cylinder was slightly better than the Nusselt number for the smooth cylinder at the same conditions.

However, despite the significance of corrugation in different industrial applications, few investigations were concerned with corrugation utilization in a heat transfer augmentation study. Therefore, the present work comprises a numerical analysis of natural convection heat transfer in corrugated annuli filled with Al2O3 nanofluid. Eight mathematical models have been developed to find the best model giving the highest heat transfer rates. The ranges of Rayleigh number and volume fraction of nanoparticles were Ra and ; respectively.

Mathematical modeling

Physical domain

In the present work, eight models of two-dimensional corrugated annuli with radius ratio = 2.66, and filled with Al2O3-H2O nanofluid are studied as can be shown in Figure 1. The inner hot cylinder and the outer cold cylinder of annuli are maintained at and ; respectively. The characteristic length is taken as . The following mathematical functions are applied to have the corrugation walls of annuli [37].For symmetry simulation, changed from to Where,

Figure 1

Geometry configurations for eight models.

at inner cylinder

at outer cylinder number of undulations

amplitude.

= azimuth angle (deg)

Governing equations and formulation

As previously mentioned, the gap between the inner and outer cylinders of the annulus is filled with Al2O3-H2O nanofluid having ten values of volume fraction extended from 0 to 0.25. The thermophysical properties of the nanofluid are given in Table 1 [59]:

Table 1

Thermophysical properties of water and nanoparticles [59].

Material	C_p(Jkg.K)	ρ(kgm^3)	k(Wm.K)	β(1K)
Water	4179	997	0.613	2.1×10^−4
Al_2O_2	765	3970	40	2.4×10^−5

The thermophysical properties of H2O-Al2O3 nanofluid at different volume fractions must be estimated. These properties include the density , thermal expansion coefficient. and specific heat energy. , thermal conductivity , and dynamic viscosity , which are given as follows [59]:

Maxwell-Garanetts and Brinkman models were used in the present work to find the thermal conductivity and dynamic viscosity of H2O-Al2O3 nanofluid; respectively [50, 60, 61, 62, 63]. A steady-state single-phase incompressible flow and heat transfer by natural convection inside an annulus filled with nanofluid is considered in the present work. Additionally, the two-dimensional continuity, momentum, and energy equations are given below [35, 64]:

The stream function and vorticity can be defined as follows:

Variables Definitions’ in the dimensionless form are given below:

,

By substituting the dimensionless variables in Eqs. (8), (9), (10), and (11) and using Eq. (12), the governing equations in streamline-vorticity form become as given below [35, 64]:

The boundary conditions are:

, at

, at ,

at .

The local and average Nusselt numbers for the inner hot cylinder and outer cold cylinder are given in the following equations:Where, is the normal direction to the wall, and is the wall length.

The major advantages of applied scheme (stream-vorticity formulation) in two-dimensional natural convection heat transfer and incompressible flow are:

continuity equation is automatically satisfied

one vorticity equation will be solved only

the streamlines of the flow are given by contour lines of the stream function.

the vorticity is a conserved quantity.

The major scheme limitation is difficulty in extending this formulation to 3D and also due to some difficulties in the vorticity boundary treatment.

Grid generation

Grid Generation is conclusive for saving of time and cost, and to obtain good results quality for all cases studied in the present work. So, the physical domain should be transformed into a computational domain. The most applicable, programmable, and usage method is the elliptical partial differential equations method which produces grids with quietly changing cell size and slope of the grid lines. Moreover, the Poisson equation is used to control the orthogonality and the required spacing close to the boundaries of the generating system. The transformation functions are defined as where & are arbitrary functions required to modify the local density of the grids.

FORTRAN code has been developed to solve the governing equations by iterative scheme named finite difference discretization. In the first step, initial grids in the interior domain are generated using a simple algebraic grid generation method. Then, the solution and iteration process are starting with the initial grid. The final grid does not depend on the initial grid. Finally, the Poisson equations are solved iteratively using the iteration method. The boundary grid points are defined as Dirichlet boundary conditions. Generally, a converged grid is sufficiently obtained after more than 100 iterations. The final grid generation shapes for eight models obtained from FORTRAN code are shown in Figure 2.

Figure 2

Grid used in the present study for eight models of corrugated annuli.

Then, the set of Eqs. (13), (14), and (15) are written in transformation form as given below [65]:Where

, , , , ,

, ,

The dimensionless boundary conditions are given as follows:

at ,

at ,

at ,

, at .

Eqs. (25) and (26) were calculated using Simpson's rule 1/3 method.

Selection of grids and validation

To obtain grid-independence, some runs are carried out to make sure the present results are disposed of the mesh calibration, as shown in Table 2. The results included the average Nusselt number on the inner hot cylinder , and maximum/minimum stream functions ,. It is noticed that the grid size of 41 41 is soft enough to achieve accurate results. So, it was adopted.

Table 2

Grid independency test at = 106, = 0.1.

Model	Grids
	21×21		31×31		41×41		51×51
	Nu_m	ψ_max,ψ_min	Nu_m	ψ_max,ψ_min	Nu_m	ψ_max,ψ_min	Nu_m	ψ_max,ψ_min
1	5.202	0.883,−20.54	5.504	0.933,−22.66	5.604	0.991,−23.71	5.604	0.990,−23.72
2	5.207	0.01,−23.44	5.400	0.002,−25.88	5.527	0,−26.21	5.527	0,−26.20
3	5.098	0.02,−26.01	5.175	0.001,−27.23	5.200	0,−28.19	5.200	0,−28.21
4	4.977	0.01,−21.11	5.149	0.004,−26.21	5.271	0,−25.77	5.271	0,−25.76
5	5.276	0.02,−25.32	6.034	0.01,−29.01	6.154	0,−30.15	6.154	0,−30.14
6	5.254	0.01,−23.08	6.006	0.01,−24.55	6.177	0,−26.61	6.177	0,−26.62
7	7.165	0.01,−31.13	6.472	0.01,−29.14	6.105	0,−28.10	6.105	0,−28.11
8	6.987	0.01,−30.31	6.353	0.01,−26.77	6.163	0,−25.12	6.163	0,−25.11

CFD code written in the FORTRAN was used to solve numerically the governing equations. The present numerical results are compared with the results obtained by kim et al. [66] and Moukalled [67] for circular hot cylinder enclosed by a square cylinder filled with air as shown in Table 3. It is noticed that the average Nusselt numbers at a certain Rayleigh number for all works are identified with a small difference ranging from 0.2 % to 1.7 %. Moreover, the results of streamlines and isotherms for the present code according to case study of Kim et al. [66] for Pr = 0.7 are shown in Figure 3. It is noticed that the present results have given a good agreement with the results of kim et al. [66]. Additional validation shown in Figure 4 is carried out by comparison the present results with the published work of Eiyad and Hakan [68] which is the natural convection in two wavy enclosures filled with Al2O3–H2O nanofluid at Ra = 105 and = 0.01. The figure shows that the patterns of streamlines and isotherms for present work are identical to results obtained by [68]. The values of minimum stream function for both works shown in Table 4 are close together with a very little difference ranging from 0.58 % to 4.2 %. As a result of these comparisons, it is concluded that the present numerical analysis is reliable and efficient.

Table 3

Code validation.

Ra	Mean Nusselt number for hot inner cylinder (air, Pr = 0.7)
	Present study	Kim et al. [66]	Error %	Moukalled [67]	Error %
104	3.3874	3.4140	-0.7830	3.331	1.6672
105	5.1643	5.1385	0.4998	5.080	1.6326
106	9.3554	9.3900	-0.3691	9.373	-0.1981
107	15.7012	15.6650	0.2309	15.790	-0.5651

Figure 3

Streamlines and isotherms in a square enclosure for present numerical results according to case study of [66], Pr = 0.7.

Figure 4

Streamlines and isotherms in an enclosure with different geometry configurations for present numerical results according to case study of [68] for Ra = 105, water-Al2O3 nanofluid with = 0.1 (a) AR = 0.85, (b) AR = 1.1.

Table 4

Comparison of present workstream function with Eiyad and Hakan work [68].

Case	EiyadandHakan [68]		Presentresults		Error %
	a	b	a	b	a	b
Ψ_min_pureFluid	−8.58	−11.86	−8.89	−11.93	-0.034	-0.0058
Ψ_min_nanofluid	−11.35	−15.8	−11.85	−15.62	-0.042	0.011

Numerical solution

The governing equations associated with boundary conditions were solved numerically by the finite difference schemes. Central difference and forward and backward upwind difference approximations were utilized for partial differential derivatives and convective terms; respectively. While the explicit method was used for velocity and temperature fields. The stream function calculations were performed by using the successive over-relaxation (SOR) method with tolerance 10−6.

Results and discussion

Streamlines and isotherms

Figures 5 and 6 show the influences of Rayleigh number and the nanoparticles volume fraction on the streamlines and isotherms; respectively, for eight models of annuli. In general, increasing of nanoparticles volume fraction causes an increase in the intensity of streamlines. The lines of flow and temperature are symmetric about the mid vertical line of annulus, so it can be taken one half of the annulus only. Heating the corrugated inner wall leads to movement of the nanofluid in the enclosure. In model one symmetric annulus) at Ra = 104, the streamlines create pair of rotating cells bounded the hot corrugated wall on each half of the annular gap, one is located on the upper undulation of the heated wall, while the other one is located close to the lower undulation of the inner cylinder. The lower and upper cells emerge together with the increase in Rayleigh number. When the cells circulated clockwise, the maximum strength of the streamlines is adopted by . The left cells have a positive stream functions, while the right cells have negative stream function. The hot fluid is displaced upward towards the hot wall. Then it moves downward near to the outer cylinder wall. This leads to generating two-fluid flow directions, clockwise on the right side, and counterclockwise on the left side. On the other side, addition of nanoparticles to the base fluid causes an increase in the thermal conductivity of nanofluid. Increasing the volume fraction of nanoparticles leads to increasing the stream function intensity. So, the streamlines become strong. While, the isotherms seem to be more smoothing as the nanoparticles volume fraction increases. Actually, the existence of nanoparticles (Al2O3) has given an accumulation of the isotherms near the hot wall. This causes improving the heat transfer rates which is expressed by the increase in the Nusselt number. The figure shows also that the behavior of temperature and fluid fields changes with changing of undulation number for inner and outer cylinders (). The activity of heated surface increases by increasing the undulations number, but the flow motion tends to be most difficult in the spaces between active undulation walls. Moreover, there are no vortices at the lower region of the enclosure in which the flow is inert and stably stratified. It is noticed that the strength of the flow movement is decreased with the increase in corrugation number. Additionally, the streamlines cells extend toward upward because the reduction in the hot area causes changing the pattern of vortices. The multiple cellular motions are still continual at model two (, symmetric annulus), with weaker flow circulation than model one ). In this model, three major vortices are found on each side of the enclosure, located close to the upper, side, and lower undulations; respectively. These vortices tend to emerge together with the increase in Rayleigh number. The intensity of vorticity (stream function) increases with the increase in number of undulations for inner and outer cylinders to six as shown in model-3 ). In this model, three cells on each side have generated at the top, middle, and bottom of the annulus. In each half of model-4 ), one major vortex is found in the middle of the annulus. It moves towards the top of the annulus as the Rayleigh number is increased.

Figure 5

Streamlines for eight models of corrugated annuli at nanoparticle volume fractions .

Figure 6

Isotherms for eight models of corrugated annuli at nanoparticle volume fractions .

At a small Rayleigh number, the temperature distribution seems to appear with a shape similar to the shape of the corrugated inner cylinder because of weak fluid flow. At high Rayleigh numbers, the hot corrugated cylinder tends to transmit the heat to the nanofluid because of its highest temperature compared to temperature of the nanofluid inside the annular gap. The streamlines cells nearby the left and right surfaces begin to extend horizontally. Moreover, the streamlines move vertically far from the lower portion of the enclosure. The isotherms pattern distorts and the flow is almost not impacted adjacent to the bottom wall. Increasing Rayleigh number means fast flow circulation and thicker thermal boundary layer in the upper region. This leads to moving the isotherms pattern vertically, which means increasing the thermal currents at the upper part of the annulus. At a high Rayleigh number (Ra = 106), the thermal currents are powerful enough to push heat toward the left and right corrugated walls. Generally, the patterns of isothermals for the symmetrical models (1–4) are nearly similar. The isotherms extend horizontally due to the flow strongly impacts the top cold wall boundary.

One main major vortex located in the first quarter of annulus is found in model-5 (). Additionally, minor vortex generates down the lower part of the enclosure and emerges with the main vortex as Rayleigh number increases. In this model, the large active volume causes increasing the temperature gradient in the upper portion of the enclosure. The strength of the fluid circulation enlarges the significantly. Two major vortices are generated in model-6 ). The first one is located at the upper part of the annulus, while the other is located on the bottom side. They emerge together with the increase in Rayleigh number. In model-7 ), the mini vortex is penetrated in the middle horizontal line of the annulus, while two major vortices are generated above and below the middle horizontal line of annulus on each side of the enclosure. These cells tend to reverse their nature if the number of undulations is reversed between the inner and outer cylinders as shown in model-8 ). It is noticed that the major vortex permeates the horizontal axis of the enclosure and the other two minor vortices are generated above and below the horizontal axis of the annulus. The major-minor vortices emerge together with the increase in Rayleigh number. Increasing the corrugation number enlarges the hot surface area and increasing the temperature gradient. As a result, the flow is accelerated. The intensity of stream function increases significantly with the increase in the Rayleigh number. This means that the buoyancy force increases and overcomes the viscous force. As a result, the natural convection is the dominating factor in the heat transfer process. The conduction heat transfer is dominated at low Rayleigh number. In this case, the isotherms arrange in paths parallel to the shape of annulus walls. The behavior of thermal fields changes with increase in Rayleigh number to 106. The dominating convective heat transfer makes the isotherms to be undulant. There are no significant changes in isotherms contour with the changing of nanoparticles volume fraction.

Local Nusselt number

Figures 7 and 8 show the local Nusselt number variation around the inner outer cylinders from = (top of wall) to (bottom of wall) for model 1 and model 6; respectively, at two Rayleigh numbers (Ra = and and three-volume fractions (). As can be shown in Figure 7 that the number of peaks for Nusselt number is similar to the number of undulations. The maximum Nusselt numbers on the inner cylinder after which the curve begins to descend are located at for model one, and at ( for model six. While, the minimum values after which the curve begins to rise are located at ( for model one, and at ( for model six, as shown in Figure 7. As we have said previously, the activity of heated surface increases by increasing the number of undulations. But, the flow motion tends to be most difficult in the space between active undulation walls. So, these regions are almost sluggish leading to reduction in the heat transfer rate. Moreover, the heat transfer rate decreases at the top and bottom of the inner cylinder. Additionally, the maximum and minimum Nusselt numbers on the outer cylinder are located at the top region ( and close to the bottom region (nearly at ; respectively, as shown in Figure 8. Generally, increasing the local Nusselt number with Rayleigh number and nanoparticle volume fraction is related to increase vortex intensity and the heat transfer activity between the nanoparticles and water. Furthermore, increasing the thermal conductivity of the nanofluid gives more heat transport to the cold wall compared with the hot wall. Additionally, the induced buoyancy force raises the hot fluid towards the inner hot wall while the colder fluid is displaced downward towards the outer cold wall. More importantly, the heat transfer behavior along the outer cylinder is slightly affected by increasing the number of undulations for the inner hot cylinder which causes increasing the heat transfer rate on the inner cylinder.

Figure 7

Local Nusselt number for the inner cylinder at different Rayleigh numbers and nanoparticle volume fractions.

Figure 8

Local Nusselt number for the outer cylinder at different Rayleigh numbers and nanoparticles volume fractions.

Average Nusselt number

Figure 9 shows the effects of nanoparticles volume fraction on the average Nusselt number for the inner and outer cylinders, for eight models of corrugated annuli at different Rayleigh numbers. Increasing the volume fraction causes an increase in the average of Nusselt number. The reason behind that the heat capacity, density, and thermal conductivity of the nanofluid increase as the nanoparticles volume fractions increase. Moreover, the heat transfer rates from the inner hot cylinder to the nanofluid are higher than similar ones from the nanofluid to the outer cold cylinder. There is no doubt that the heat transfer process enhances with the increase in Rayleigh number because the stronger buoyancy force increases the movement of fluid particles and transport heat to the boundary walls. This leads to accelerating the Nusselt number on the inner hot wall and outer cold wall of the corrugated annulus.

Figure 9

Variation of average Nusselt number with nanoparticles volume fractions for inner cylinder (I. C.) and outer cylinder (O. C.) for eight models of corrugated annuli at different Rayleigh numbers.

One of the important objectives of the present study is finding the best model giving the highest heat transfer rates. The variation of on the inner and outer walls of corrugated annulus versus Ra for eight models is shown in Figure 10, at three-volume fractions ( = 0.01, 0.1, and 0.25); respectively. The figure shows that there is a slight influence for increasing Ra from 104 to 105 on the average Nusselt number. While, the average Nusselt number increases with an increase in Rayleigh number to 106. The data extracted from the current study reveals that the highest heat transfer rates can be sequenced for the eight models as follows: (5, 6, 7, 8, 2, 1, 3, 4) for the inner hot cylinder, and (6, 8, 7, 5, 2, 1, 3, 4) for the outer cold cylinder. It is noticed that the unsymmetrical models (5, 6, 7, 8) produce the higher heat transfer rate than symmetric models (2, 1, 3, 4). The model-5 produces the highest heat transfer rate than other models. This due to that, one great main vortex is generated in this model witch its intensity increases with the increase in nanoparticle volume fraction. The wide gap of model-5 relatively compared with other models causes growing the large eddies easily which then merge to generate a great main vortex at the top region of the annulus gap. It is seen that increasing the number of undulations of unsymmetrical annulus leads to decreasing the heat transfer rate. Model-4 produces the lowest heat transfer rate than other models. This due to the weak wavy vortex resulted from a high number of undulations for inner and outer corrugated walls (. As a result, the convection currents are moving along the undulation paths leads to producing lower heat transfer rates.

Figure 10

Mean Nusselt number for the hot and cold cylinders of the corrugated annuli against Rayleigh number.

Average Nusselt number correlations

The logarithmic mean Nusselt number for the inner and outer corrugated cylinders against the logarithmic for inner and outer cylinders are shown in Figure 11-a & b; respectively. The ranges of Rayleigh number and volume fraction of nanoparticles are Ra and (; respectively. The general correlation for the average Nusselt number as a function of Rayleigh number and nanoparticles volume fraction has been deduced as follows:Where c, n, and m are constants given in Table 5.

Figure 11

Logarithmic mean Nusselt number versus logarithmic for the inner and outer cylinders with different models.

Table 5

Constants of Eq. (27) for eight models.

Inner Cylinder
Model	c	n	m
1	0.7888	0.1571	1.1150
2	0.7451	0.1649	1.1211
3	0.6687	0.1783	1.1315
4	0.6611	0.1798	1.1327
5	0.7152	0.1818	1.1343
6	0.6855	0.1886	1.1397
7	0.6822	0.1878	1.1390
8	0.6400	0.1972	1.1465
Outer Cylinder
1	0.7123	0.1235	1.0894
2	0.6897	0.1332	1.0967
3	0.6488	0.1408	1.1025
4	0.6591	0.1386	1.1008
5	0.7314	0.1198	1.0867
6	0.7932	0.1081	1.0778
7	0.6869	0.1347	1.0978
8	0.7385	0.1224	1.0885

Conclusions

The present work displays a numerical study of natural convection heat transfer of H2O-Al2O3 nanofluid inside symmetrical and unsymmetrical corrugated annuli at different nanoparticles volume fractions and Rayleigh numbers. It can be concluded from the current study that:

Adding the nanoparticles (Al2O3) to the base fluid leads to more accumulation of the isotherms near the hot wall and an increasing the flow strength and heat transfer rate.

Increasing the nanoparticles volume fraction and Rayleigh number leads to increasing the stream function intensity and heat transfer rate.

The activity of heated surface increases by increasing the undulations number and making the flow motion tends to be most difficult in the spaces between active undulation walls. This leads to a decrease in the heat transfer rate.

The heat transfer rate in unsymmetrical annuli is higher than that in the symmetrical annuli.

Model 5 = 3 and = 6) produces the highest heat transfer rate than other models.

Model 4 ( = = 8) produces the lowest heat transfer rates than other models.

The conduction heat transfer is dominated at a low Rayleigh number. In this case, the isotherms arrange in paths parallel to the shape of annulus walls.

There are no significant changes in isotherms contour with the changing of nanoparticles volume fraction.

The values and trends of the local Nusselt number for the outer cold cylinder are slightly affected by increasing the undulations number of the inner cylinder (from . While this increase has given a behavior change and increased the rate of heat transfer in the inner cylinder.

General correlations of mean Nusselt number for the inner and outer corrugated cylinders as a function of Rayleigh number and nanoparticles volume fraction have been deduced for eight models.

Declarations

Author contribution statement

Sattar Aljabair: Conceived and designed the experiments; Performed the experiments; Analyzed and interpreted the data.

Akeel Abdullah Mohammed: Analyzed and interpreted the data; Contributed reagents, materials, analysis tools or data; Wrote the paper.

Israa Alesbe: Performed the experiments; Analyzed and interpreted the data; Contributed reagents, materials, analysis tools or data.

Funding statement

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Competing interest statement

The authors declare no conflict of interest.

Additional information

No additional information is available for this paper.