General Solution for Unsteady Natural Convection Flow with Heat and Mass in the Presence of Wall Slip and Ramped Wall Temperature

Saeed Ullah Jan,Sami Ul Haq,Syed Inayat Ali Shah,Ilyas Khan,and Muhammad Atif Khan

1Faculty of Mathematics and Statistics,Ton Duc Thang University,Ho Chi Minh City,Vietnam

2Department of Mathematics,Islamia College University,Peshawar 25000,Pakistan

3Department of Mathematics,Kohat University of Science&Technology,Kohat 26000,Pakistan

AbstractThis work is focused on the effect of heat and mass transfer with unsteady natural convection flow of viscous fluid along with ramped wall temperature under the assumption of the slip wall condition at the boundary.Analytical solutions are obtained by using Laplace transformation to the non-dimensional set of governing equations containing velocity,temperature and concentration.Moreover,the expression for skin-friction is derived by differentiating the analytical solutions of fluid velocity.Numerical tables for Skin-friction,Sherwood number and Nusselt-number are examined.For the physical aspects of the flow,we use various values of involved physical parameters such as Prandtl number(Pr),slip parameter(η),Schmidt number(Sc),buoyancy ratio parameter(N),Sherwood number(Sh),and time(t).Additionally,the general solutions are plotted graphically and a comprehensive theoretical section of numerical discussions is included.

Key words:slip boundary condition,heat and mass transfer,skin-friction

Nomenclature

1 Introduction

In the literature,the theoretical study of laminar free convection flow from a vertical plate is continuously receiving the attention of researchers due to their technological and industrial applications.In fluid mechanics,the solutions of the non-dimensional governing equations have been attempted by analytical and numerical techniques.A number of unsteady natural convection flows past a vertical plate covered by various set of thermal conditions at the bounding plate have been investigated.[1−3]Many applied problems may be demanded by considering wall conditions,which are arbitrary or non-uniform.Hayday et al.[4]added the numerical solutions for non-similar boundary layer laminar flows over a vertical plate under the effect of temperature.Fully expanded laminar flows are examined by Schetz.[5]He obtained the approximate and analytical solutions for such two dimensional flow under thermal effects.Later on,the analytical work to re fine earlier solutions for the variation of various types of wall temperature was tried by Kao,[6]Kelleher,[7]and Yovanovich.[8]Soundalgekar[9]and Ingham[10]firstly found the exact solution for the incompressible viscous fluid near vertical plates for the natural convection flow.The work of Soundalgekar[9]and Ingham[10]over the convectional flows was further extended by Singh and Kumar,[11]Soundalgekar and Gupta[12]and Raptis et al.[13]Magnetic effects on the flow of viscous fluid were also studied by Raptis and Singh.[13]These solutions provided a new streamline for the researchers to investigate different heat flow problems under the influence of magnetic field.Solution for the oscillating flow of viscous fluid with slip boundary condition is established by Khaled and Vafai.[14]Buoyancy forces that are the results of density distinction are responsible for the free convection flow of fluid.Many researchers used free convection flow due to enormous applications in industrial and natural sciences,such as liquid crystals,cooling of electron components,chemical suspensions,and dilute suspensions.[15]

In mechanics of fluids we can use two types of boundary conditions in different flow problems i.e.slip and no slip boundary conditions If there is no corresponding motion between the fluid and wall directly in connection with the wall,then it would be no-slip boundary condition.Maybe for the simpli fi cation of difficult situation,this condition is used extensively but it has also some impediment.There are few sufficiently smooth surfaces where the wall and fluid have no connections and the fluid slips against the wall.[16]For example,the no-slip condition is not able to work in the capillaries.[17]However,these impediments which are mentioned by the Navier in his beginning work,[18]are overcome by the slip condition.Slip condition is also called Navier condition.The slip condition has various applications in different fields of medical science,extrusion,lubrication,particularly in flows through porous media,polishing arti fi cial heart valves,micro and Nano fluids,studies of friction at different surfaces and biological fluids.[19−20]In the last few years,many authors used the slip condition for the flows of viscous fluids.[21−23]

Natural convection flow of a fluid which is electrically conducting and flowing at a vertically moving plate having ramped temperature,is examined by Ghara et al.[24]Seth et al.[25]discussed the impulsive plate driven flows and examined the Soret as well as Hall effects over such MHD flows,where the suction phenomenon in such flows through a porous plate is considered by Reddy.[26]Samiulhaq et al.[27]investigated a second-grade fluid for Magneto hydrodynamic natural convection flow. The unsteady MHD natural convection flows near an accelerated vertical plate along temperature source and constant heat flux was investigated by Narahari and Debnath.[28]Micropolar fluid with slip condition and periodic temperature passed through a porous medium in Magneto hydrodynamic natural convection flow was studied by Mansour et al.[29]

The purpose of present work is the investigation of the impact of slip wall conditions on unsteady natural convection flows with the phenomena of heat transfer and the ramped wall temperature.We consider the slippage as well as ramped wall temperature,at the boundary.Analytical solutions of the non-dimensional governing equation i.e.temperature,concentration and velocity are obtained by integral transform method(Laplace transform).Numerical tables for Skin-friction,Sherwood number and Nusselt-number are examined.At the end,graphical explanation for velocity,temperature,concentration,Skinfriction,Sherwood number and Nusselt-number are presented and the impacts of involved physical parameters are observed.

2 Analysis of the Governing Problem

Let us consider the unsteady flow of incompressible viscous fluid with the influence of slippage over the boundary.The pressure variations are assumed to be negligible.Motion of the fluid is induced subject to slippage of fluid at the surface of the plate.The plate is considered along x′-axis in the upward direction and y′-axis is perpendicular to the plate.At the initial stage for t′≤ 0,both the fluid and the plate is assumed to be stationary at the same temperature T∞′and concentration C∞′.When 0t0the temperature of the plate is kept at uniform temperature Tw′and the level concentration of the fluid at plate is elevated to Cw′or concentrations is accumulated to the plate at constant rate.By Boussinesq’s approximation the unsteady incompressible viscous fluid,which flows near infinite vertical plate is investigated by.[30−32]

The corresponding initial and boundary conditions were under as:

where λ is the dimensional slip parameter.The imposed condition on the velocity is identical to those obtained by Samiulhaq et al.[33]while on the temperature and concentration is identical to those obtained by Narahari et al.[34]and Nehad et al.[35]To convert the governing equation into non-dimensional form,we introduce a set of nondimensional variables

As the flow is fully developed over a plate of infinite length,so the physical variables are functions of y′and t′,where u′denotes the dimensional velocity in x′-direction.Using these non-dimensional variables(5)in Eqs.(1),(2),and(3),the governing equations(1),(2),and(3)take the simplest forms

where υ denotes the velocity in non-dimensional form,ζ is coordinate axis which is normal to the plates in nondimensional form,τ is time in non-dimensional form,Θ is temperature in non-dimensional form,Ψ is concentration in non-dimensional form,Gm is mass Grashof number,Gr is thermal Grashof number,Sc is Schmidt number,Pr is Prandtl number,N is the mass to thermal buoyancy ratio parameter.Moreover,characteristics time t0can be defined as

is identical to that of Samiulhaq et al.[33]and Narahari et al.[36]Corresponding initial and boundary condition(4)in the non-dimensional form which are:

where η is the dimensionless slip parameter.

3 Analytical Laplace Transform Solutions

3.1 Solution of Ramped Temperature

Applying Laplace transforms to Eqs.(6),(7),(8)and using initial and boundary condition(10),we obtain the following problem for temperature,concentration and velocity for Pr=1,Sc=1,and Pr=1,Sc=1.

where

The inverse Laplace transform of Eqs.(11),(12),and(13)are

where

H(τ−1)is the unit step function defined as

where c is an arbitrary constant.

(ii)Case 2 Pr=1,Sc=1,

where

The inverse Laplace transform of Eqs.(18),(19),and(20)are

where

(iii)Nusselt Number The corresponding Nusselt number in case of ramped temperature,which is the measure of the rate of heat transfer at the plate can be determined by differentiating Eq.(14)partially with respect to“ζ”

using ζ=0,we get

Also the corresponding Nusselt number can be determined by differentiating Eq.(21)partially with respect to “ζ”

and using ζ=0,we get

(iv)Sherwood Number The corresponding Sherwood number in case of ramped temperature which the measure of the ratio of convective and diffusive mass transfer can be determined by differentiating Eq.(15)partially with respect to “ζ”

using ζ=0,we get

Also the corresponding Sherwood number can be determined by differentiating Eq.(22)partially with respect to“ζ” and using ζ=0,we get

(v)Skin-friction The corresponding Skin-friction in case of ramped temperature,which is the measure of the rate of shear stress at the plate can be determined by differentiating Eq.(16)partially with respect to “ζ”

using ζ=0,we get

where

Also,corresponding Skin-friction can be determined by differentiating Eq.(23)partially with respect to “ζ” and using ζ=0,we get

3.2 Solution of Stepped Temperature

To investigate the ramped temperature effect of the fluid flow on the boundary,it is essential to compare the flows with the one nearby the isothermal plate.Under this assumption,it can be displayed that the non-dimensional temperature,concentration,and velocity variable can be expressed by applying Laplace transforms to Eqs.(6),(7)and(8)and using initial and boundary condition(10).

(i)Case 1 Pr=1,Sc=1,

where

The inverse Laplace transform of Eqs.(30),(31),and(32)are

where

(ii)Case 2 Pr=1,Sc=1,

The inverse Laplace transform of Eqs.(36),(37),and(38)are

(iii)Nusselt Number The corresponding Nusselt number in case of stepped temperature,which is the measure of the rate of heat transfer at the plate can be determined by differentiating Eq.(33)partially with respect to“ζ”

using ζ=0,we get

Also the corresponding Nusselt number can be determined by differentiating Eq.(39)partially with respect to “ζ”and using ζ=0,we get

(iv)Skin-friction The corresponding Skin-friction in case of stepped temperature which is the measure of the rate of shear stress at the plate can be determined by differentiating Eq.(35)partially with respect to “ζ”

using ζ=0,we get

Also the corresponding Skin-friction can be determined by differentiating Eq.(41)partially with with respect to “ζ”and using ζ=0,we get

4 Special Cases

(i) Ψ(ζ,τ)=0,υ(0,τ)=0

We take Ψ(ζ,τ)=0 and use the boundary condition υ(0,τ)=0 into Eq.(6).After lengthy,but straightforward computation,the fluid velocity takes the form:

The result is identical to those obtained by Chandran et al.[1]for ζ=y,τ=t,and

(ii) υ(0,τ)=0

Introducing υ(0,τ)=0 into Eq.(6),we find that fluid velocity takes the form:

The result is identical to those obtained by Narahari et al.[36]for

(iii) Ψ(ζ,τ)=0

Consider Ψ(ζ,τ)=0 into Eq.(6),we find that fluid velocity takes the form:

As it was expected the corresponding result is identical to those obtained by Samiluhaq et al.[33]for ζ=y, τ=t,and

5 Result and Discussion

The main theme of this paper is the investigation of thermal and mass diffusion in the unsteady,slip flow of a viscous fluid having the property of incompressibility.Various numerical values of the skin-friction and velocity were calculated as function of time τ for different values of physical parameters which are Schmidt number Sc,Prandtl number Pr,slip parameter η and buoyancy ratio parameter N.The influence of these parameters on ramped and constant profiles of velocity υ,temperature Θ,Concentration Ψ,skin-friction Cf,Nusselt number Nu and Sherwood number Sh are shown graphically.From the views of Figs.1 and 2,it is clear that growing values of slip parameter accelerates the velocity of fluid in case of both ramped and stepped wall temperature.In Fig.3,we compare the impact of ramped and constant temperature corresponding to the plate on the fluids velocity.It has been noted that the velocity of fluid(free stream)for constant temperature is higher than the velocity for ramped temperature.These predictions have been illustrated in Fig.3.In order to compare the effect of ramped and constant temperature on the fluid velocity for various value of the buoyancy ratio parameter N,we have plotted Fig.4.We observe from Fig.4 that the motion of fluid rises for large value of buoyancy ratio parameter N(i.e,N>0).The reason is that the mass buoyancy force acts in the same direction as the thermal buoyancy force.It is also perceived that velocity in case of constant temperature is greater than velocity in case of ramped temperature.Physically this is true because in ramped temperature case heating takes place gradually when compared to isothermal temperature case.

Fig.1 Non-dimensional velocity outline along ζ for various values of slip parameter η corresponding to ramped temperature of plate with τ=0.8,Pr=0.71,and Sc=0.16.

In Fig.5,we discuss the behavior of velocity for different value of time parameter τ,for both the cases,that is ramped and constant temperature.When we increase time,the velocity of fluid is increased.Also in this Fig.5,we compare the effect of ramped and constant temperature.Furthermore,it has been observed that velocity of the fluid for constant temperature is comparatively higher than ramped temperature.Further analysis suggests that the fluid velocity is maximum near the plate and subsequently decreases as we move away from the plates.The motion of fluid becomes zero in particular case,where the distance possibly enlarges.Moreover the velocity is increasing function of time.Figure 6 shows that an increasing Sc implies that viscous forces dominate over the diffusional effects.Schmidt number in free convection flow regimes represents the relative effectiveness momentum and mass transfer by diffusion in the velocity and concentration.Smaller values correspond to lower molecular weight species diffusing i.e.Hydrogen in air(Sc∼0.16)and higher values to denser hydrocarbons diffusing in air i.e.Ethyl benzene in air(Sc∼2.0).effectively therefore if Sc is increased then momentum diffusion will counteract,since viscosity effects will increase and molecular diffusivity will be reduced.The flow will therefore be decelerated with a rise in in both cases,namely ramped temperature of the plate and constant temperature of the plate as testified to by Fig.6.It also discovered that the rate of velocity in case of constant temperature is greater than in case of ramped temperature.This is expected since in the case of ramped wall temperature the heating of the fluid takes place more gradually than in the isothermal plate case.This feature is important in for example achieving better flow control in nuclear engineering applications,since ramping of the enclosing channel walls can help to decrease velocities.

Figures 7 and 8 show the velocity of fluid in case of ramped temperature and constant temperature for different value of τ and two various values of Pr i.e.(Pr=0.71)and(Pr=7).

Fig.2 Non-dimensional velocity outline for different type values of slip parameter η corresponding to stepped temperature of the plate with τ=1.2,Pr=0.71,and Sc=0.16.

Fig.3 Non-dimensional velocity outline for different values of slip parameter η corresponding to stepped and ramped temperature of the plate with τ=3,Pr=0.71,and Sc=2.01.

Fig.4 Non-dimensional velocity outline for various value of N corresponding to ramped and stepped temperature of the plate with τ=3,Pr=0.71,and Sc=2.01.

Fig.5 Non-dimensional velocity outline for various value of τ corresponding to ramped and stepped temperature of the plate with η=0.9,Pr=0.71,and Sc=2.01.

Fig.6 Non-dimensional velocity outline for various value of Sc corresponding to ramped and stepped temperature of the plate with η =0.9,Pr=0.71,and τ=3.

Fig.7 Non-dimensional velocity outline for various value of Pr and τ corresponding to ramped temperature of the plate with η=0.9 and N=0.8.

Fig.8 Non-dimensional velocity outline for various value of Pr and τ corresponding to stepped temperature of the plate with η=0.9 and N=0.8.

The two most frequently encountered fluids in engineering are air and water and these values of Pr correspond to these two cases,respectively.This approach was established by Ostrach and Simon[37]at NASA in the early 1950s.It is obviously follows from Figs.7 and 8 that if we increase time then the velocity of fluid is increasing but decreasing for large value of Pr.Since higher Pr fluids will possess greater viscosities and this will serve to reduce velocities.Furthermore from Figs.7 and 8 that the velocity of fluid for the air(Pr=0.71)is larger than that of water(Pr=7)in both cases.These Figs.7 and 8 clearly show that when time is zero,the velocity satisfies the initial condition given in Eq.(10).Figures 9 and 10 show the variations of temperature outline for various value of τ and for two different values of Pr such as for air(Pr=0.71)and for water is(Pr=7).These analyses predict that from Figs.9 and 10 that temperature profiles decrease if we increase the value of Pr for the both cases of ramped and constant temperature of the plate.In Fig.11,we compare the variation of temperature profile in both cases,namely ramped and constant temperature of the plate.From Fig.11 it is clear that constant temperature profiles are found larger than ramped temperature profiles.In addition,it is is observed from Fig.11 that for large values of time,the temperature of the fluid is increasing in previously both cases.All the observation and result obtained from these figures are identical to the boundary condition given in Eq.(10)at temperature.Figure 12 shows the variations of concentration profile for different value of τ and for two different values of Sc.It is observed that concentration profiles decreases if we increase the value of Sc,while increases for large value of τ.

Fig.9 Non-dimensional temperature outline for different values of Pr and τ corresponding to ramped temperature of the plate.

Fig.10 Non-dimensional temperature outline for different values of Pr and τ corresponding to stepped temperature of the plate.

Fig.11 Non-dimensional temperature outline for different values of τ and Pr=0.71 corresponding to ramped and stepped temperature of the plate.

Fig.12 Non-dimensional concentration outline for different values of Sc and τ corresponding to ramped temperature of the plate.

Fig.13 Modification of skin-friction Cffor various value of slip parameter η corresponding to the ramped temperature of the plate with N=0.8,Pr=7,and Sc=2.01.

Fig.14 Modification of skin-friction Cffor various value of slip parameter η corresponding to the stepped temperature of the plate with N=0.8,Pr=7,and Sc=2.01.

Figures 13 and 14 display the modifications of skin-friction along with time for various value of slip parameter η.Figures 13 and 14 show that if we increase the value of slip parameter η then skin-friction is decreasing in both possible cases,discussed in previous section i.e.ramped and constant temperature of the plate.In Fig.15 we compare the modification of skin-friction in both cases namely,ramped and constant temperature plate.It is noted from Fig.15 that Modification of skin-friction in case of constant temperature plate is found smaller than ramped temperature of the plate.Figure 16 shows the modification of Sherwood number along with time τ for various value of Sc.It is clearly observed that Sherwood number is increasing for the growing value of Sc.Figures 17 and 18 show the modification of Nusselt number along with time τ for various value of Pr.Figures 17 and 18 inference that if we increase the value of Pr both in cases,namely ramped temperature of the plate as well as in constant temperature of the plate then Nusselt number is increasing.It is also obviously noted that Nusselt number for mercury(Pr=0.015)is smaller than air(Pr=0.71),electrolytic solutions(Pr=1)and for water(Pr=7).These behaviors in current study have a proper agreement with the observations made in Ref.[33].

Fig.15 Modification of skin-friction Cffor various value of slip parameter η corresponding to the ramped and stepped temperature of the plate with N=0.6,Pr=7,and Sc=2.01.

Fig.16 Modification of Sherwood number Sh for various value of Pr corresponding to stepped temperature of the plate.

Fig.17 Modification of Nusselt-number Nu for various value of Pr corresponding to the ramped temperature of the plate.

Fig.18 Modification of Nusselt-number Nu for various value of Pr corresponding to the stepped temperature of the plate.

Table 1 Skin-friction for ramped temperature.

Table 2 Skin-friction for stepped temperature.

Numerical results of skin friction,Nusselt number,and Sherwood number at the plate(ζ=0)are presented in Tables 1 to 5 for various values of slip parameter η,Prandtl number Pr,Schmidt number Sc,mass to thermal buoyancy ratio parameter N and time τ.It is observed from Table 1 that skin friction Cfincreases with an increase in Prandtl number Pr,mass to thermal buoyancy ratio parameter N and time τ while the result is reversed with increase in slip parameter η and Schmidt number Sc for ramped wall temperature as well as for stepped wall temperature in Table 2.Numerical results of Nusselt number at the plate(ζ=0)are expressed in Tables 3 and 4 for different values Prandtl number Pr and time τ.

Table 3 shows that the Nusselt number Nu which determines the rate of heat transfer at the plate increases as Prandtl number Pr and time τ progresses for ramped wall temperature.It is perceived from Table 4 that the Nusselt number Nu increases for large values of Prandtl number Pr while the result is reversed with increase in time τ for stepped wall temperature.The rate of concentration at the surface of the wall is calculated in term of Sherwood in Table 5.It is seen that Sherwood number Sh increases with an increase in Schmidt number while the result is reversed with increase in time.

Table 3 Nusselt number for ramped temperature.

Table 4 Nusselt number for stepped temperature.

Table 5 Sherwood number.

6 Conclusion

Exact solutions are obtained for the wall driven flow of a viscous fluid under the influence of ramped wall temperature.Slippage over the vertical plate is not negligible.The dimensionless mathematical model for the fluid is operated with the help of Laplace transforms.Furthermore,these obtained solutions are in terms of simpler transcendental and complementary error functions.It is verified that solutions for temperature,concentration,and velocity satisfy the imposed initial and boundary conditions.The impact of involved physical parameters on the velocity field,temperature,Nusselt number and Skin friction has been examined.The entire discussion can be reduced to the following concluding remarks:

•Velocity increases for the large values of slip parameter.

•Skin friction is decreasing function of slip parameter.

•Increasing values of Prandtl number Pr decrease the speed of fluid.Since higher Pr fluids will possess greater viscosities and this will serve to reduce velocities,thereby lowering the skin friction.

•For the large values of Prandtl number Pr the thermal diffusivity,take small space.

•Increasing values of Schmidt number Sc,stops the motion of fluid.An increasing implies that viscous forces dominate over the diffusional effects.Schmidt number in free convection flow regimes represents the relative effectiveness momentum and mass transfer by diffusion in the velocity and concentration.effectively therefore if it is increased then momentum diffusion will counteract,since viscosity effects will increase and molecular diffusivity will be reduced.The flow will therefore be decelerated with a rise in Sc.

•Relation of Schmidt number Sc and skin friction are inverse.

•Velocity and skin friction are increasing functions of buoyancy ratio parameter.The reason is that the mass buoyancy force acts in the same direction as the thermal buoyancy force.

•For large values of time,the temperature of fluid rises and it moves faster.We compare the effect of ramped and stepped temperature corresponding to the plate on the fluid velocity.For this situation our result indicates that velocity in case of constant temperature is greater than in case of ramped temperature.The present results are useful in further explaining the important class of flows in which the driving force is induced by a combination of the thermal and chemical diffusion effects.Such results have immediate relevance in industrial thermo fluid dynamics,transient energy systems and also buoyancy-driven geophysical and atmospheric vertical flows.