Numerical Assessment of Nanofluid Natural Convection Using Local RBF Method Coupled with an Artificial Compressibility Model

Muneerah Al Nuwairanand Elmiloud Chaabelasri

1Department of Mathematics and Statistics,College of Science,King Faisal University,Al-Ahsa,31982,Saudi Arabia

2PTPME Laboratoy,Faculty of Sciences,Mohamed First University,Oujda,60000,Morocco

ABSTRACT In this paper, natural heat convection inside square and equilateral triangular cavities was studied using a meshless method based on collocation local radial basis function(RBF).The nanofluids used were Cu-water or Al2O3-water mixture with nanoparticle volume fractions range of 0 ≤ φ ≤ 0.2.A system of continuity,momentum,and energy partial differential equations was used in modeling the flow and temperature behavior of the fluids.Partial derivatives in the governing equations were approximated using the RBF method.The artificial compressibility model was implemented to overcome the pressure velocity coupling problem that occurs in such equations.The main goal of this work was to present a simple and efficient method to deal with complex geometries for a variety of problem conditions.To assess the accuracy of the proposed method,several test cases of natural convection in square and triangular cavities were selected.For Rayleigh numbers ranging from 103 to 105, a validation test of natural convection of Cu-water in a square cavity was used.The numerical investigation was then extended to Rayleigh number 106,as well as Al2O3-water nanofluid with a volume fraction range of 0 ≤φ ≤0.2.In a second investigation,the same nanofluids were used in a triangular cavity with varying volume fractions to test the proposed meshless approach on non-rectangular geometries.The numerical results appear to be in agreement with those from earlier investigations.Furthermore,the suggested meshless method was found to be stable and accurate,demonstrating that it may be a viable alternative for solving natural heat transfer equations of nanofluids in enclosures with irregular geometries.

KEYWORDS Natural heat convection;nanofluids;cavities;meshless method;radial basis function

Abbreviations

Greek Symbols

Subscripts

1 Introduction

In heat transfer, the heat convection occurs when heat is transferred from one part of the fluid to another.Natural convection has become an active research topic due to its lower cost and wide applications in engineering.Several models were developed for calculating the efficiency of thermal conductivity in solid and fluid systems [1–6].In [7], Choi showed that adding some nano metal particles to the base fluid increases its conductivity.Such mixtures,known as nanofluids,have many applications in science and engineering [8].These applications include fuel cells, hybrid-powered engines, chillers, and heat exchangers.The usefulness of a nanofluid for heat transfer applications can be verified by modeling the convective transfer in the nanofluid.

Several numerical methods, such as finite elements, finite volumes, and finite differences have been used to solve the thermal convection equations and simulate convective heat flow in various geometries.These methods are mesh based methods and are widely used.Several researchers have used these methods to model natural convection in nanofluids.Oztop et al.[9] used the method of finite volumes to investigate thermal transfer in a square cavity filled with a nanofluid.The investigation examined the effect of the size of the nanoparticles and the type of nanofluid.Mahmoodi et al.[10]used the same method with the SIMPLER algorithm to investigate heat transfer in a square cavity filled with nanofluid and having adiabatic square bodies at its center.These two studies found that increasing the volume fraction of the nanofluid increases heat transfer and that the type of nanoparticle is a key factor in heat transfer improvement.Jasim et al.[11]studied the influence of an inner adiabatic rotating cylinder on mixed convection of hybrid nanofluid in more complex geometries.Their research gave important insights for rotary heat exchanger designers.The set of conservation equations and their associated boundary conditions were discretized using the finite volume technique.

Finite elements have also been widely used.For example,Bhowmick et al.[12]used the Galerkin Finite Element Method to discretize the governing equations of natural convection heat transfer and entropy generation for a square enclosure,containing a heated circular or square cylinder subjected to non-uniform temperature distributions on the left vertical and bottom walls.The finite element method was also used by Islam et al.[13] to analyze temperature transfer within a prismatic cavity filled withCu-water nanofluid under two different temperature boundary conditions.

It is worth noting that,over the last two decades,the lattice Boltzmann method has seen significant advances and is now commonly used in heat transfer problems, particularly for nanofluids.In Izadi et al.[14], natural convection of multi-wall carbon nanotubes-Iron Oxide nanoparticles/water hybrid nanofluid inside a ⊥-shaped enclosure was investigated.The effects of different parameters such as the nanoparticles volume fraction,the cavity obstruction ratio,the heat source aspect ratio,the heat source position,and the Rayleigh number on the hydrodynamic and thermal characteristics were studied using an in-house code based on the lattice Boltzmann method.Naseri Nia et al.[15]analyzed,using the lattice Boltzmann method,the effect of baffle configuration on natural convection inside an L-shaped enclosure filled withCu-water nanofluid,for different parameter ranges of Rayleigh number and nanoparticle volume.Their main finding was that at high Rayleigh numbers,only the longer baffle improves convection regardless of its positioning,while increasing the volume fraction of nanoparticles always improves the thermal performance.To reduce the needed computations and to deal with more complex geometries,the meshless methods offer very good alternatives.These methods do not require a mesh,using instead a set of nodes in the space of the study domain.

Recently,Petrov-Galerkin local meshless method was used to study natural convection in a porous medium[16],Wijayanta et al.[17]analyzed a local radial basis function method to obtain numerical solutions of the conjugate natural convection heat transfer problem for square cavity.Zhang et al.[18]tested the robustness and accuracy of the variational multiscale element free Galerkin method,another meshless method,on several cases of natural convection,including a semicircular cavity,a triangular cavity with a flat wall, and a zigzag shape.The results obtained in these studies clearly showed that meshless methods are stable, accurate, and efficient.In addition to these classic methods, artificial neural networks(ANN)were used in studying convention heat transfer.An ANN based method was developed by Rostami et al.[19]to predict the thermal conductivity of multiwall carbon nanotubeswater and Cu-water nanofluids.Using experimental data,their algorithm selects ANNs with optimal performance on the prediction problem at hand.This study found that ANNs have excellent predictive ability and can make important contributions to the determination of the most desirable performance.

Other lines of investigation relevant to our study were explored.Kavusi et al.[20]studied the effect of different nanofluids, prepared using alumina, copper oxide and silver nanoparticles, at different concentrations and particle diameters on the performance of heat pipe,a heat exchange device with high efficiency and performance in heat transfer.Their study found that the use of a nanofluid in place of water resulted in an increase in thermal efficiency and a reduction in heat at the heat pipe wall.Rahmati et al.[21]simulated heat transfer in a microtube.The cooling fluid used was a combination of 1%and 1.5%volume fractions of CuO nanoparticles in a non-Newtonian pseudo-plastic fluid.Several challenges were identified in this work,namely the effect of slip velocity and the choice of the power law model for dynamic viscosity estimation.This study found that increasing the volume fraction of solid nanoparticles and the slip velocity coefficient resulted in an increase in heat transfer.

The present work aims to develop an accurate method using local radial basis function combined with the artificial compressibility model to solve the natural convective heat transfer problems in complex geometries.These methods have rarely been used in the case of nanofluids.A square cavity and an equilateral triangular cavity filled with aCu-water orAl2O3-water based nanofluid with different volume fractions were studied.The use of a combination of the local radial basis function and the artificial compressibility method to solve the natural convection problems in a nanofluid is the main contribution of this work.The numerical algorithm uses the explicit two-stage Range-Kutta scheme for the pseudo-temporal discretization.We demonstrate that the radial basis function method allows for easy treatment of non-rectangular geometries.Indeed,the generation of collocation nodes in the computational domain and the treatment of boundary conditions do not require additional effort.In addition, we demonstrate that approximating the derivatives can be performed accurately using local radial basis functions.

This paper is structured as follows: The problem and the mathematical implementation of the governing equations are described in Section 2.The numerical formulas of the meshless approach are detailed in Section 3.In Section 4,the numerical method is applied and the numerical results are presented, along with a discussion of the accuracy and efficiency.Section 5 summarizes the results,adding some remarks and conclusions.

2 Problems Definition and Governing Equations

2.1 Problem Definition

To demonstrate the simplicity and efficiency of the proposed method to deal with complex geometries, we choose two domains.The proposed configurations are illustrated in Fig.1, they consist of:

· A two-dimensional square cavity.The left and right walls of the square cavity are the hot wall(Th)and the cold wall(Tc),respectively.The upper and lower horizontal walls of the cavity are kept adiabatic.

· An equilateral triangular cavity where the bottom wall is the hot wall while the other two upper walls are kept adiabatic.

Figure 1: Schematic geometry of physical problems: (a) Differentially heated square cavity and (b)Uniformly heated triangular cavity from bottom

The cavities considered are filled with a nanofluid composed of a nanoparticle-water mixture.The nanofluid is considered Newtonian,and the flow is laminar.The thermo-physical properties of the base fluid and the nanoparticle are taken as[22]and listed in Table 1.They are assumed constant except for the density which is given by the Boussinesq approximation.

Table 1: Thermophysical properties of water and nanoparticles

2.2 Governing Equations

The equations of the continuity,momentum,and energy for buoyancy-driven laminar fluid flow and heat transfer of the nanofluid inside a cavity are

whereuandvare the components of the velocity in thexandydirections,Tis the temperature,pis the pressure,the meaning of the other variables are found in the abbreviations at the end of the paper.Using the dimensionless variables:

L is the characteristic length as well as the length of the cavity edge.For this work,we set it to 1.The above equations can be written in the dimensionless form:

For the cases studied, we have used three types of walls; hot, cold or adiabatic wall.The corresponding boundary conditions are

·U(X,Y)=V(X,Y)=0 for all walls.

·θ(X,Y)=0 for cold walls.

·θ(X,Y)=1 for hot walls.

For calculating the local and average Nusselt numbers, the following formulas were used,respectively

In this work, the values of the Nusselt number were calculated on horizontal or vertical walls,where for a horizontal wallξ=yand for a vertical wallξ=x.

2.3 Nanofluid Thermophysical Properties

The thermophysical parameters of the base fluid,such as density,viscosity,and heat conductivity,are affected by the addition of nanoparticles.In this study,the models that describe the thermophysical properties were obtained from the literature and are as follows:

The effective density of the nanofluid is given as

whereφis the volume fraction of the nanoparticles,ρfis the density of base fluid andρpis the density of nanoparticles.Thermal diffusivity of the nanofluid is given by

whereknfis the effective thermal conductivity.Following Maxwell’s model[23]for spherical particles,we used the following expression:

Herekfis the thermal conductivity of the base fluid,andkpthe thermal conductivity of dispersed particles.The heat capacity of the nanofluid introduced in Eq.(12)is given by

The thermal expansion coefficient of the nanofluid is given by

Finally,the effective viscosity of the nanofluid,based on Brinkman model[24],is given by:

whereμfis the viscosity of the base fluid.

3 Numerical Method

3.1 Radial Basis Function(RBF)Method

We employ the radial basis functions method,which is a meshless method originally proposed by Kansa [25], for the numerical solution of natural convection equations.In the following, we outline the main principles of this method.

Radial basis function interpolation approximates the solution function by an expansion.According to the Kansa method, the solution is approximated on a set ofNcollocation nodes by a linear combination of local radial basis function as follows[26]:

whereλF,jare the unknown coefficients to be calculated,φ(x)is a radial basis function,‖x－xj‖is the distance between the nodes x andxj,andsi,nis a stencil that contains the indexiand indices of thennearest neighboring of the nodexi.The stencil consists of the indices of node points falling within three times the distance to the nearest node.Thus,the number of nodes indexed by the stencil is variable,as it depends on the distribution of nodes in the computational domain.

Many radial basis functions were proposed In the literature.In the current study, we use the infinitely smooth multiquadric radial basis function defined asφ(‖x－xj‖)=whereε0 is a shape parameter that plays a substantial role in the resulting accuracy.Generally, the determination of the range of shape value can be obtained by using numerical tests.In this study,we,in addition,use the procedure proposed by Sarra in[27].Therefore,an optimal shape parameter has been identified.In all simulations,ε=in whichdminis the minimum distance between two nodes andε0=10－2for square cavity andε0=0.15 for triangle cavity.

The numerical solution of an equation involving partial derivatives can be approximated by using linear approximation of the partial derivatives at the x nodes,based on local radial basis functions as follows:

where(Φ[i])－1is the inverse of the matrixΦ[i].By removing the dependence on the expansion coefficients,the equation in(19)becomes:

Note that in the governing equations the first and second derivatives are calculated for each node by using several matrix operations onn×nmatrices.This slowdown the calculation process.In order to overcome the problem of successive inversion of manyn×nmatrices,Yao et al.[26]reformulated the local interpolation and differentiation matrices in terms of a global sparse matrix by appending zeros to the proper entries.This reformulation was also used by[28]for the space-time derivatives.The entries of the global differentiation matrixDLare then defined as follows:

Using this approach,the first and second spatial derivatives of a functionF(x,t)can be expressed as follows:

whereDx,Dy,DxxandDyyare first and second derivative matrices with respect toxandy.

3.2 Numerical Procedure

In attempting a numerical solution of the Eqs.(6)–(9), one is faces with several difficulties.The momentum equations contain a pressure gradient,moreover,these equations can not be solved if the pressure term is not specified.The method of artificial compressibility proposed by Chorin [29] is adopted in this work for the treatment of the pressure-velocity coupling.This method adds a pseudotemporal derivative to the continuity equation to couple pressure with velocity,as follows:

The little mermaid could not help thinking of her first rising out of the sea, when she had seen similar festivities and joys; and she joined in the dance, poised116 herself in the air as a swallow when he pursues his prey117, and all present cheered her with wonder

Hereγis an artificial compressibility parameter,a free but important parameter that affects the convergence of the solution.In this study,to obtain the convergence of the solution,a value of 10－5was employed.

For simplicity,the Eqs.(25),(7)–(9)are rewritten in a compact form as follows:

where,

with the superscriptTdenoting the transpose and

On a set ofNnodes,and using the differentiation matrix defined above,the right-hand side of the partial differential problem(25)can be approximated as

The problem in Eq.(26)is integrated explicitly using the second Runge-Kutta method.Hence,we divide the pseudo-time into intervals[ξn,ξn+1],and denote Wnthe value of the function W at timeξn.The recursive computation is given by

where

andL(Wn)is given in Eq.(26).

It should be noted that radial basis function approaches for time-dependent partial differential equations with diffusive terms can be stably advanced in time with a suitable choice of time step size[27].Although,the equations used in our study are time dependent,the aim of the study is to examine the steady-state solution.Thus, the method is used to reach the steady-state as the limit when the variations of variables approaches zero.The convergence criteria for steady state is

and

Herenandn+1 denote two adjacent time steps andεerris taken equal to 10－5in all cases.

3.3 Grid Independence and Validation Tests

Several preliminary tests were performed to evaluate the sensitivity of the results to the number of nodes.These tests were used to determine the optimal distribution of nodes that ensures good accuracy with reasonable computational cost,Three distributions are used for this investigation and are shown in Fig.2.For these testsCu-water nanofluid with volume fractionφ=0.2 is considered as the working fluid,and the Rayleigh number is taken as the constant 105.A horizontal cross section aty= 0.5 of the horizontal velocityUwas made and shown in Fig.3.The same figure shows the variation of the local Nusselt number at the left wall.From these figures,it can be observed that the solutions become independent of the node distribution from 31×31.Therefore,to ensure good accuracy of the results,the distribution 41×41 was selected for the numerical investigations.

Figure 3:Grid independence tests for variations of the U component of the velocity vs.X coordinate at Y=0.5(left)and local Nusselt number vs.Y coordinate at left wall(right)

The meshless method proposed in this work was then validated in the case of natural convection on two problems.The first was the natural convection in a cavity filled with pure water.In the second,we considered theCu-water nanofluid with a volume fractionφ=0.2.The results are shown in Figs.4 and 5.The results presented in Fig.4 in terms of local Nusselt number at left wall for various Rayleigh number values,show that our code satisfactorily reproduces the results reported in[30].A comparison of the graphs in Fig.5 of average Nusselt for different Rayleigh numbers shows that our results are in excellent agreement with those obtained in[30].

Figure 4:Validation cases in Basak et al.[30]and present study.local Nusselt number vs.Y coordinate at left wall for Ra=103,104 and 105 and for Cu-water nanofluid with φ=0.2 and pure water

Figure 5:Validation cases in Basak et al.[30]and present study.Average Nusselt number vs.Rayleigh for φ=0,φ=0.1 and φ=0.2 of Cu-water nanofluid

4 Results and Discussion

4.1 Differentially Heated Square Cavity

In the first investigation, we ran a simulation of the effects of differentially heated boundary conditions on natural convection in a square cavity filled with a nanofluid.We used the cavity previously for the validation of meshless method,but extend the range of the Raylieh number to 106using nanofluids composed of different volume fractionsφ= 0.1, 0.15 and 0.20 using two types of nanoparticles,namelyCuandAl2O3.Their characteristics are listed in Table 1.

Figs.6 and 7 show isotherms distributions for diluteCu-water andAl2O3-water nanofluids.Isotherms are also shown for pure water for comparisons.Results are shown for all ranges of 103≤Ra≤ 106and volume fractionsφ= 0.1, 0.15, and 0.20.Overall, the effect of adding nanoparticles is seen in all ranges of the Rayleigh number,as well as for the various proposed volume fractions.We also note that the convective heat flux for water is stronger near the center in all cases,particularly for high Rayleigh numbers and that the isotherms for water are more twisted at the center than those for theCu-water orAl2O3-water nanofluid.Near the wall,the isotherms forCu-water andAl2O3-water are different from those of water.These observations can be explained by the suppression of convection at low Rayleigh numbers.This also explains why the iso-lines are nearly vertical lines.In contrast,convection improves rapidly at the center of the cavity with increasing Rieylegh number reaching high ratios forRa=106.Finally,as Rayleigh number and volume fraction grow,the intensity of heat transfer increases,and the behavior of isotherms forRa=105and 106are identical.

Figure 6:Isotherms of differentially heated square cavity with a hot wall on the left,cold wall on the right,of pure water(dotted line)and a Cu-water nanofluid(continued),for volume fraction φ=0.1,0.15 and 0.2,and the numbers of Ra=103,104,105 and 106

Figure 7:Isotherms of differentially heated square cavity with a hot wall to the left,cold wall to the right,of pure water(dotted line)and a Al2O3-water nanofluid(continued),for volume fraction φ=0.1,0.15 and 0.2,and the numbers of Ra=103,104,105 and 106

Figs.8 and 9 illustrate the variations in the local Nusselt number along the heated wall for different values of the Rayleigh number and different volume fractions for theCu-water andAl2O3-water nanofluid,respectively.In general,the Rayleigh number has a significant impact on the variation of the local Nusselt number,the heat transfer rate increases with rising Rayleigh number due to higher buoyancy forces.ForRa= 103,the variation in local Nusselt number is small,but it grows for other values ofRa.The local Nusselt number has its highest values near the bottom of the wall, and it decreases gradually as it goes up the wall.This is mainly due to the decrease in the buoyancy of nanoparticles due to the gravitational force.This increases the conductivity of the nanofluid at the bottom of wall.It should be noted that increasing the volume fraction of nanoparticles inCuorAl2O3water, increases the conductivity of the nanofluid and consequently the increases the local Nusselt number.Thus,the natural heat convection is better enhanced with copper nanoparticles.

Figure 8:Local Nusselt number at the left wall for Cu-water nanofluid with the volume fractions φ=0,0.1,0.15 and 0.2 at Rayleigh numbers 103,104,105 and 106

Figure 9:Local Nusselt number at the left wall for Al2O3-water nanofluid with the volume fractions φ=0,0.1,0.15 and 0.2 at Rayleigh numbers 103,104,105 and 106

Fig.10 shows the variation of average Nusselt number against nanoparticle volume fraction forRa= 103, 104, 105and 106.The average Nusselt number is sensitive to nanoparticle type and volume fraction, as shown in the Figures.It is worth noting that theCu-water andAl2O3-water nanofluids followed similar trend.At low Rayleigh numbers, the effect of nanoparticles and the growth in their volume fractions is more relevant than at high Rayleigh numbers.As an example,forCu-water nanofluid,the slope of the Rayleigh number against volume fraction is nearly at 20%forRa=103and 10%atRa=106.TheAl2O3-water exhibits the same behavior.This can be explained by the fact that at low Rayleigh numbers, conduction dominates heat transfer and viscous effects are less important.As a result, the inclusion of high thermal conductivity nanoparticles will improve conduction.The viscous effects, on the other hand, become more pronounced with large Rayleigh numbers.Finally,comparisons with previous results of averaged Nusselt number in Basak et al.[30]shows that our results are in good agreements with theirs.Therefore,the meshless method based on radial basis functions coupled with the artificial compressibility,described in this work,can be used with confidence for the simulation in such problems.

Figure 10: Comparison and variation of average Nusselt number with volume fractions of Cu and Al2O3 nanoparticles at different Rayleigh number

4.2 Uniform Heated Triangular Cavity

In the second investigation, we looked at the natural convection in an equilateral triangular enclosure.The two nanofluids employed in the previous case were used.The two inclined walls of the cavity were kept cold,while the bottom wall is heated.Pranowo et al.[31]used a similar configuration filled with air to test the performance of their meshless method.The nodes distribution used for the numerical computation(Fig.11)contains 1043 non-uniformly distributed nodes.Figs.12 and 13 show the impact of Rayleigh number and volume fraction of nanoparticles on isotherms.As expected,for low Rayleigh numbers thermal conduction dominates.For the case of the triangular cavity this dominance occurs at Rayleigh numbers lower than 104.Thereafter, the behavior of the isotherms remains unchanged with the same pattern around the vertical plane atx= 0.5.Indeed, the high viscosity of the nanofluid negatively affects the flow,leading to a slower convection and a decrease in heat transfer.As the Rayleigh number increases,so does the buoyancy force,resulting in an increase in natural convection,this results in the distortion towards the top of the isotherm.

Figure 11: Node distribution used for numerical studies for triangular cavity, 1043 non-uniformly distributed nodes

Figure 12:Isotherms of uniformly heated triangular cavity with a hot wall on the bottom,cold lateral walls, of Cu-water nanofluid, for volume fraction φ=0.1, 0.15 and 0.2, and Rayleigh numbers of Ra=104,105 and 3.105

Figure 13:Isotherms of uniformly heated triangular cavity with a hot wall on the bottom,cold lateral walls, of Al2O3-water nanofluid (continued) and Cu-water (doted line), for volume fraction φ=0.1,0.15 and 0.2,and Rayleigh numbers of Ra=104,105 and 3.105

Figs.12 and 13, forRa= 106, indicate the dominance of convection.It should be noted that in spite of the high viscosity, the flow intensity is significantly improved for large Rayleigh values.Consequently,heat is clearly transferred from the hot bottom wall to the cold side walls.Furthermore,the heat transfer of the two nanofluids fits the same pattern seen in the case of the square cavity,with the heat transfer being better for theCu-water nanofluid thanAl2O3-water.

In order to quantify the heat exchange within the cavity filled withCu-water andAl2O3-water,the effects of the volume fraction of both nanoparticles on the local Nusselt number along bottom wall are shown in Figs.14 and 15 forRa=104,105and 3.105.The smallest values of the local Nusselt number are observed in the middle of the heated part,where the temperature gradient is the smallest according to the isotherms shown in Figs.12 and 13.It should be noted that the addition of nanoparticles increases the thermal conductivity of the nanofluid, and thus the heat transfer, which results in an increase in the local Nusselt number.Fig.16 illustrates a comparison of the local Nusselt number for the two nanofluids.TheCu-water mixture remains a more conductive nanofluid,promoting heat transfer from the hot to cold walls.

Figure 14:Local Nusselt number at the bottom wall for Cu-water nanofluid with the volume fractions φ=0,0.1,0.15 and 0.2 at Rayleigh numbers 104,105 and 3.105

Figure 15: Local Nusselt number at the bottom wall for Al2O3-water nanofluid with the volume fractions φ=0,0.1,0.15 and 0.2 at Rayleigh numbers 104,105 and 3.105

Figure 16: Comparison of local Nusselt number at the bottom wall for Cu-water and Al2O3-water nanofluid with the volume fractions φ=0.10 and 0.2 at Rayleigh numbers 104 and 3.105

Fig.17 shows the variation of the average Nusselt number as a function of the volume fraction for various values of the Rayleigh number.As can be seen,increasing the volume fraction results in higher average Nusselt number, this is due to the increased conductivity achieved by adding nanoparticles.Furthermore, these results confirm thatCu-water nanofluid is a better heat transfer fluid thanAl2O3-water nanofluid when used in a triangular cavity, and this appears only for high Rayleigh numbers.

Figure 17:Variation of average Nusselt number with volume fractions of Cu and Al2O3 nanoparticles at different Rayleigh number

5 Conclusions

A two-dimensional numerical investigation of natural convection in enclose cavities was carried out in this study.Two geometries were used, namely those of a square and a triangular cavity.The laminar heat transfer equations of nanofluid were solved using the radial basis function method coupled with an artificial compressibility, which is classified as a meshless method.Cu-water andAl2O3-water nanofluids with various volume fractions of nanoparticles at different Rayleigh numbers,starting at 103and going up to 106,were studied.The following are the main findings of the study:

· The implemented meshless method offers high flexibility in dealing with complex geometries due to the simplicity of the numerical evaluations of space derivative.This method provides a valuable and efficient way to analyze natural convection in nonrectangular geometries.

· To overcome the pressure velocity coupling problem that occurs in equation systems like the ones examined,the artificial compressibility model was implemented.The quality of the results was found to be of the same accuracy order as the other classical methods.

· Results clearly indicate that the addition of nanoparticles has produced a substantial enhancement of heat transfer as compared to that of the pure fluid.As mentioned above, convective transfer dominates conduction when Rayleigh numbers are increased.As a consequence,heat transfer from the hot walls to the cool walls is improved.This remains true for both of the cases investigated.Furthermore,we conclude that the heat transfer enhancement is not clear in low Rayleigh numbers or when the enclosures have geometric singularities,as was the case in our study,for the case of the triangular cavity.

Acknowledgement:The authors would like to acknowledge the financial support from King Faisal University,Saudi Arabia,Project No.AN000675.

Funding Statement:This work was supported through the Annual Funding Track by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University,Saudi Arabia[Project No.AN000675].

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.