Peristaltic Flow of Dusty Nanofluids in Curved Channels

Z.Z.Rashed and Sameh E.Ahmed

1Mathematics Department,Jouf University,Qurayyat,Saudi Arabia

2Department of Mathematics,King Khalid University,Abha,Saudi Arabia

Abstract:In this paper,numerical investigations for peristaltic motion of dusty nanofluids in a curved channel are performed.Two systems of partial differential equations are presented for the nanofluid and dusty phases and then the approximations of the long wave length and low Reynolds number are applied.The physical domain is transformed to a rectangular computational model using suitable grid transformations.The resulting systems are solved numerically using shooting method and mathematical forms for the pressure distributions are introduced.The controlling parameters in this study are the thermal buoyancy parameter Gr,the concentration buoyancy parameter Gc,the amplitude ratio ,the Eckert number Ec,the thermophoresis parameter Nt and the Brownian motion parameter Nb and the dusty parameters Ds,s.The obtained results revealed that an increase in the Eckert number enhances the temperature of the fluid and dusty particles while the nanoparticle volume fraction is reduced.Also,both of the temperature and nanoparticles volume fraction are supported by the growing of the Brownian motion parameter.

Keywords:Peristaltic motion;dusty particles;nanofluid;curved channels

1 Introduction

Study of the peristalsis flow has many applications such as movement of the food in the intestine tract,the urine passage from a kidney to the bladder blood flow in small veins and arteries of the blood circulation,transferring the ovum in the Fallopian tube and the movement of sperm in the channels.Also,there are many applications of peristaltic motion in the biomedical devices,such as blood pumps and heart lung machines.The physiology of the gastrointestinal tract discussed by Latham[1],Mishra et al.[2].Also,many studies interested with peristaltic flow in a curved channel[3–11].Sato et al.[3]discussed the peristaltic flow in a curved channel.Ali et al.[4]studied the peristaltic flow in a curved channel with a long wavelength approximation.Ali et al.[5,6]examined effects of the heat transfer and fluid flow of a non-Newtonian third grade fluid in a curved channel.Hayat et al.[7]discussed the Newtonian fluid peristaltic flow,heat and mass transfer in a curved channel with compliant walls.The investigation in[7]was extended by Hayat et al.[8]and Hina et al.[9]to include case of a third grade fluid.Hina et al.[9]considered case of the peristaltic motion in curved channels contain compliant walls using Johnson-Segalman fluid.Hina et al.[10,11]studied the combined heat and mass transfer effects as well as the influence of wall properties on the peristaltic flow of Johnson-Segalman and the peristaltic flow of pseudoplastic fluid,respectively.

In the recent years,the researchers focused on studying the nanofluids due to their applications in various fields.The first study was introduced by Choi[12]who studied the pure fluids with suspended nanoparticles.He discussed the substantial augmentation of the heat transported in suspensions of copper or aluminum nanoparticles in water or other liquids.Buongiorno[13]takes into his account the Brownian diffusion as well as thermophoresis in writing the transport equations.The literature survey indicates that attempts on peristalsis of the nanofluid model in a curved channel are little.Hina et al.[14]discussed numerically the peristaltic flow of a nanofluid in a curved channel.Ayub et al.[15]studied the mixed convection in the presence of a thermal radiation and a chemical reaction,analytically.The results indicated that the heat transfer rate decreases with the increase in thermophoresis parameter.Narla et al.[16]studied the peristaltic transport of a Jeffrey nanofluid in a curved channel and examined effects of various parameters on the fluid flow and the temperature distributions.Noreen et al.[17]discussed the induced magnetic field effects on the peristaltic flow in a curved channel.They found that an increase in the Brownian motion and thermophoresis parameters causes an increase in the temperature profiles.Hayat et al.[18]studied the peristaltic motion of a copper-water based nanoliquid with the thermal slip conditions.Hayat et al.[19]studied the MHD peristaltic flow of Sisko nanofluids with the Joule heating effects.They used a numerical treatment for the governing equations.They found that the increasing values of the curvature parameter results in symmetric behaviors at the centerline of the channel for the velocity,temperature and concentration distributions.Tanveer et al.[20]studied the peristaltic motion of a Sisko fluid with homogeneous-heterogeneous reaction effects.The results revealed that the lower velocity,temperature and concentration profiles are obtained in case of the higher bending.Tanveer et al.[21]studied the peristaltic flow of Eyring-Powell nanofluids in a curved channel with compliant walls.They found that the Eyring-Powell parameters tend to decrease the velocity and temperature of the nanofluid while the concentration bears a dual response.Hayat et al.[22]discussed the peristalsis of MHD Jeffery nanofluids in a curved channel with a porous medium.

On the other hand,the practical applications of the dusty fluid flow can be found in atmospheric,engineering and physiological fields,for example,conveying of powdered materials,purification of crude oil,environmental pollutants,dust in gas cooling systems,petroleum industry.Many other applications are included in the valuable book written by Rudinger[23].Farbar et al.[24]studied the heat transfer by flow of the gas-solid mixtures in a circular tube.Saffman[25]presented a dusty fluid model in the laminar flow.Many researchers[26–33]extended the dusty fluid topic with different physical circumstance.Recently,this topic is generalized to case of the dusty nanofluid by many researchers.Siddiqa et al.[34]conducted an analysis of a two-phase natural convection flow of dusty nanofluid along a vertical wavy surface.They found that presence of the dust particles have a notable influence on the temperature distribution as the isotherms get stronger for the dusty water.Begum et al.[35]studied the gyrotactic bioconvection of the dusty nanofluid along an isothermally heated vertical wall.They applied a numerical treatment for the mathematical model using the two-point implicit finite difference method.Gireesha et al.[36],Gireesha et al.[37]studied the Hall effect on a two-phase transient flow with stretching sheet using KVL model and irregular heat generation/consumption,respectively.Good recently studies in the nanofluid topics are found in[38–43].

The authors in all the mentioned papers disregarded case of the peristaltic motion of the dusty nanofluids in complex shapes.Therefore the main objective of this paper is to study the peristaltic flow in a curved channel using the dusty nanofluids.The two-phase nanofluid model is used to simulate this case and approximations of the long wave length and low Reynolds number are applied.Also,one of the objectives of this study is to express the pressure distributions in the flow domain and examining effects of the dusty and geometry parameters on the nanofluid flow,heat and nanoparticles distributions.

2 Discerption of the Problem and Mathematical Formulation

Consider an unsteady two dimensional peristaltic motion of a dusty nanofluid inside a curved channel.Fig.1 shows conditions of the problem and coordinates system.In this figure,theR-axis is taken normal to surface of the channel and theX-axis is taken along walls of the channel.Width of the channel is 2d1 surrounded in a circle of a radiusR*and a centerO.Boundaries of the channel walls are determined as:wherec,aand λ are the speed,amplitude and length of the wave.In addition,the nanofluid is modeled using the two-phase model in which effects of the Brownian motion and thermophoresis are included.The base fluid,dusty particles and nanoparticles are in a thermal equilibrium model.Viscous dissipation effect is considered and a linear Boussinesq approximation is taken into account.A uniform size of the dusty particles is assumed and they distribute equally in the mixture.

Figure 1:The physical model of the channel

Under all the mentioned assumptions,the governing equations of the problem are introduced as,see[14,26,27]

2.1 Nanofluid Phase

2.2 Dusty Particles

In Eqs.(1)–(9),U,V,Us,Vsare the velocity components of the fluid and dust phase,in the laboratory frame(R,X);P,Psis the pressure,of the fluid and dust phase,ρf,ρsis density of the fluid and dust phase,μ is the dynamic Viscosity,ν is the kinematic viscosity,κ is the thermal conductivity,cpis the specific heat at constant pressure,Cis the concentration andTis the temperature of the fluid,α is the coefficient of linear thermal expansion of the fluid,β is the coefficient of expansion with concentration,andgis acceleration due to gravity,Ifr,(x)andu,(v),us,(vs)are the coordinates and velocity components in the wave frame then

Also,the following non-dimensional quantities are introduced:

In Eq.(11),the subscriptsrefers to the dusty phase and 0 refers to the conditions at the channels walls.Moreover,definitions the stream function for the nanofluid phase ψ and the dusty particles phase ˉψ are expressed as:

Substituting Eqs.(10)–(12)in the systems of Eqs.(1)–(5)and(6)–(9)and applying the approximations of the low Reynolds number and long of the wave length,the governing equations become:

2.3 Nanofluid Phase

2.4 Dusty Particles

where

In Eq.(20),Pris the Prandtl number,Reis the Reynolds number,Ecis the Ekert number,Gris the thermal buoyancy parameter,Gcis the concentration buoyancy parameter,Nbis the dimensionless Brownian motion,Ntis the dimensionless thermophores parameter,γ is the specific heat ratio of the mixture,Dsis the mass concentration of particle phase and αdis the dust parameter.

Eqs.(13),(14),(17)and(18)after eliminating the pressure terms are written as:The corresponding boundary conditions are given by:

On the other hand,rate of the volume flows in the laboratory frame and in the wave frame are,respectively,given by:

The dimensionless mean flows in the laboratory and in the wave frameFare defined as:

From the previous equations,the following relations are obtained:

3 Method of Solution

To solve the governing Eqs.(15),(16),(19),(21)and(22),it is needed to map the wavy boundaries into a rectangular computational domain.Therefore,the following new independent variables are introduced:

The partial derivatives for the dependent variables are obtained as follows:

Using Eqs.(31),(32),the computational domain is transformed to -1 ≤η′≤1 which makes the applying of the numerical method is available.Here the Runge-Kutta method with shooting technique is used to solve the resulting system of the equations.The number of the grid points are taken to be equal 401 and the convergence criteria is 10-6.In addition,a validation test consisting of comparisons with previously published results is performed and presented in Fig.2.It is found very good agreements are observed between the presented study(in special cases)and those obtained by Hina et al.[14].

Figure 2:Validation tests at αd =0 and Ds =0

4 Results and Discussion

Discussion of the obtained results is notified in this section.A set of graphical results in terms of the velocity profiles for the dusty particles,temperature distributions and nanoparticle volume fraction are presented in Figs.3–15.During these computations,the governing parameters are considered in wide ranges,i.e.,range of the thermal buoyancy parameter 1 ≤Gr≤4,range of the concentration buoyancy parameterGcis 1 ≤Gc≤4,range of amplitude ratio ϵ is 1 ≤ϵ ≤5,range the Eckert number is 0.1 ≤Ec≤0.5,range of the thermophoresis parameterNtis 0.1 ≤Nt≤0.5 and the Brownian motion parameterNbis varying from 0.1 to 0.5.Here it should be mentioned that from Eq.(19),the profiles of the temperature of the dusty particles are the same of nanofluid temperature.Also,values of the mass concentration of the dusty particlesDsand the dusty parameters αdare assumed 0.1 and 10,respectively.

Figure 3:Profiles of the fluid temperature for variations of Gr at Gc=1,k =100,ϵ=0.2,Ec=1,Nt =0.5,Nb=0.5,αd =10 and Ds =0.1

Figure 4:Profiles of the nanoparticles volume fraction for variations of Gr at Gc=1,k =100,ϵ=0.2,Ec=1,Nt =0.5,Nb=0.5,αd = 10 and Ds =0.1

Figs.3 and 4 display profiles of the temperature distributions for both fluid and dusty particles and nanoparticles volume fraction for different values of the thermal buoyancy parametersGratGc=1,k=100,ϵ=0.2,Ec=1,Nt=0.5,Nb=0.5,αd=10 andDs=0.1.It is found that,in the curvature domain,the temperature distributions are enhanced as the thermal buoyancy parameter increases.The interpretation of this behavior,physically,related to the temperature differences inside the flow domain that is enhanced asGrincreases and hence the fluid and dusty particles temperature are supported.In addition,it is noted that profiles of the nanoparticles volume fraction are reduced asGris growing due to the fact that the increase in the temperature differences diminishes the concentration distributions.

Figure 5:Profiles of the fluid temperature for variations of Gc at Gr=1,k =100,ϵ=0.2,Ec=1,Nt =0.5,Nb=0.5,αd =10 and Ds =0.1

Figure 6:Profiles of the fluid temperature for variations of Gc at Gr=1,k =100,ϵ=0.2,Ec=1,Nt =0.5,Nb=0.5,αd =10 and Ds =0.1

In Figs.5 and 6,variations of the concentration buoyancy parameterGcand their effects on the dusty particles and fluid temperature and nanoparticles volume fraction are presented.The other parameters are fixed atGr=1,k=100,ϵ=0.2,Ec=1,Nt=0.5,Nb=0.5,αd=10 andDs=0.1.The results revealed that the increase inGccauses a clear reduction in the temperature distributions.The physical explanation of this behavior is due to the buoyancy forces due to the concentration differences that minimize profiles of the temperature.However,an obvious enhancement in the nanoparticles volume fraction is seen asGcincreases due to the support in concentration differences in the curvature domain.

Profiles of velocity of the dusty particles,the fluid and dusty particles temperature and the nanoparticles volume fraction for different values of the amplitude ratio ϵ atGr=Gc=1,k=100,Ec=1,Nt=0.5,Nb=0.5,αd=10 andDs=0.1 are shown in Figs.7,8 and 9,respectively.It is observed that an obvious reduction in the velocity profiles is seen as ϵ is increased due to an increase in area of the channel.Also,based on the fact that the decrease in the fluid velocity enhances the fluid temperature,the fluid temperature is supported as ϵ increases.Additionally,the nanoparticles volume fraction distributions are wasted as the amplitude ratio increases due to the increase in the temperature difference inside the flow domain.

Figure 7:Dusty velocity for different values of ϵ at Gc=Gr=1,k =100,Ec=1,Nt =0.5,Nb=0.5,αd =10 and Ds =0.1

Figure 8:Profiles of the temperature for variations of ϵ at Gr=Gc=1,k =100,Ec=1,Nt =0.5,Nb=0.5,αd =10 and Ds =0.1

Figure 9:Profiles of the nanoparticles volume fraction for variations of ϵ at Gr=Gc=1,k =100,Ec=1,Nt =0.5,Nb=0.5,αd =10 and Ds =0.1

Figure 10:Profiles of the temperature for variations of Ec at Gr=Gc=1,k =100,ϵ=0.2,Nt =0.5,Nb=0.5,αd =10 and Ds =0.1

Figure 11:Profiles of the nanoparticles volume fraction for variations of Ec at Gr=Gc=1,k =100,ϵ=0.2,Nt =0.5,Nb=0.5,αd =10 and Ds =0.1

Figure 12:Profiles of the temperature for variations of Nt at Gr=Gc=1,k =100,ϵ=0.2,Ec=1,Nb=0.5,αd =10 and Ds =0.1

Figure 13:Profiles of the nanoparticles volume fraction for variations of Nt at Gr=1,Gc=1,k =100,ϵ=0.2,Ec=1,Nb=0.5,αd =10 and Ds =0.1

Effects of the viscous dissipation represented by variation of the Eckert numberEcon the temperature distributions and nanoparticles volume fraction are examined with help of Figs.10 and 11.These figures are plotted atGr=Gc=1,k=100,ϵ=0.2,Nt=0.5,Nb=0.5,αd=10 andDs=0.1.As it expected,a clear enhancement in profiles of the temperature is obtained asEcincreases.Physically,the fluid is heated up asEcincreases due to the kinetic energy in the flow area.On the other hand,the temperature differences are increased and hence the nanoparticles volume fraction is detracted asEcis growing.

Figs.12 and 13 analyzed effects of the thermophoresis parameterNton the temperature distributions and nanoparticles volume fraction atGr=Gc=1,k=100,ϵ=0.2,Ec=1,Nb=0.5,αd=10 andDs=0.1 The results show that effects ofNton the temperature are clearer than the nanoparticles volume fraction.In addition,both the temperature and concentration profiles are reduced asNtincreases.In fact,this outcome agrees with the previous study reported by Hina et al.[14].

Figure 14:Profiles of the temperature for variations of Nb at Gr=Gc=1,k =100,ϵ=0.2,Ec=1,Nt =0.5,αd =10 and Ds =0.1

Figure 15:Profiles of the nanoparticles volume fraction for variations of Nb at Gr=1,Gc=1,k =100,ϵ=0.2,Ec=1,Nt =0.5,αd =10 and Ds =0.1

Effects of the Brownian motion parameterNbon the temperature distributions and profiles of the nanoparticles volume fraction are displayed in Figs.14 and 15.It is found that the increase inNbresults in an increase in both of the temperature and nanoparticles volume fraction.Also,the results revealed that the concentration distribution is higher at the walls in comparison with the central part of the curved channel.All these effects are examined the a referenced caseGr=1,Gc=1,k=100,ϵ=0.2,Ec=1,Nt=0.5,αd=10 andDs=0.1.

5 Conclusions

In this investigation,a peristaltic flow of a dusty nanofluid in a curved channel was,numerical studied.Approximations of the low Reynolds number and the long wave length are considered.Two systems of the equations are presented for the nanofluid phase and the dusty particles phase.A mathematical form for the pressure distributions is introduced and validation tests with previously published results are performed.The following findings can be summarized:

Values of the dusty temperature are equal values of the nanofluid temperature regardless variations of the governing parameter.

An increase in the thermal buoyancy parameter causes an enhancement in the temperature distributions while a clear reduction in profiles of the nanoparticles volume fraction is noted.

The nanoparticles volume fractions are supported by the increase in the concentration buoyancy parameter but the both the nanofluid and the dusty particles are decreased.

The viscous dissipation helps in the increase of temperature of the nanofluid and the dusty particles while it decreases the nanoparticles volume fraction.

Velocity of the dusty particles is reduced as the amplitude ratio increases(particularly at the midsection of the channel.

Funding Statement:The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through research groups program under Grant Number(R.G.P2/72/41).

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