Conjugate Effects of Heat and Mass Transfer on Natural Convection Flow along an Isothermal Sphere with Radiation Heat Loss

Conjugate effects heat and mass transfer on free convection flow across an isothermal sphere immersed in a viscous incompressible fluid in the presence of species concentration with radiation heat loss has been investigated in this paper. The governing boundary layer equations are first transformed into a non-dimensional form and the resulting nonlinear systems of partial differential equations are then solved numerically using Finite—difference method with Keller-box Scheme. We have focused our attention on the evaluation velocity profiles, temperature profiles and species concentration profiles of the fluid as well as the local skin friction coefficient, local heat transfer rate and local species concentration transfer rate for a wide range of radiation parameter Rd, Schmidt number Sc and Prandlt number Pr.


Introduction
Convection is the mode of energy transfer between a solid surface and the adjacent liquid or gas that is in motion and it involves the combined effects of conduction and fluid motion.The faster the fluid motion, the greater the heat transfers due to convection.Free convection flow is often encountered in cooling of nuclear reactors or in the study of the structure of stars and planets.The study of temperature and heat transfer is of great importance to the engineers because of its almost universal occurrence in many branches of science and engineering.The application of boundary layer techniques to mass transfer has been of considerable assistance in developing the theory of separation processes and chemical kinetics.Some of the interesting problems that have been studied are mass transfer from droplets, free convection on electrolysis in non-isothermal boundary layer.Radiation contributes substantially to energy transfer in furnaces, combustion chambers, fires, and to the energy emission from a nuclear explosion.Radiation must be considered in calculating thermal effects in rocket no-zzles, power plants, engines, and high temperature heat exchangers.Radiation can sometimes be important even though the temperature level is not elevated and other modes of heat transfer are present.Radiation has a great effect in the energy equation which leads to a highly non-linear partial differential equation.Heat, mass and momentum transfer on a continuously moving or a stretching sheet has several applications in electrochemistry and polymer processing [1][2][3].
Gebhart and Pera [4] investigated the natural convection flow resulting from the combined buoyancy effects of thermal and mass diffusion.Diffusion and chemical reaction in an isothermal laminar flow along a soluble flat plate was studied and an appropriate mass-transfer analogue to the flow along a flat plate that contains a species say, A slightly soluble in the fluid say, B has been discussed by Fairbanks and Wick [5].Hossain and Rees [6] have investigated the combined effect of thermal and mass diffusion in natural convection flow along a vertical wavy surface.The effects of chemical reaction; heat and mass transfer on laminar flow along a semi-infinite horizontal plate have been studied by Anjalidevi and Kandasamy [7].By taking advantage of the mathematical equivalence of the thermal boundary layer problem with the concentration analogue, results obtained for heat transfer characteristics can be carried directly over to the case of mass transfer by replacing the Prandtl number Pr by the Schimdt number Sc.However, the presence of a chemical reaction term in the mass diffusion equation generally destroys the formal equivalence with the thermal energy problem and moreover, generally prohibits the construction of the otherwise attractive similarity solutions.Takhar et al. [8] for example, considered diffusion of chemically reactive species from a stretching sheet.The application of the boundary layer theory with chemical reaction has been applied to some problems of free and mixed convection flow from the surface of simple geometry by the above authors.Chiang et al. [9] investigated the laminar free convection from a horizontal cylinder.Sparrow and Lee [10] looked at the problem of vertical stream over a heated horizontal circular cylinder.They have obtained a solution by expanding velocity and temperature profiles in powers of x, the co-ordinate measuring distance from the lowest point on the cylinder.The exact solution is still out of reach due to the non-linearity in the Navier-Stokes equations.By taking advantage of the mathematical equivalence of the thermal boundary layer problem with the concentration analogue, results obtained for heat transfer characteristics can be carried directly over to the case of mass transfer by replacing the Prandtl number Pr, the Schimdt number Sc and the radiation parameter Rd.
It appears that Merkin [11,12] was the first who presented a complete solution of this problem using Blasius and Görtler series expansion method along with an integral method and a finite-difference scheme.Ingham [13] investigated the free convection boundary layer flow on an isothermal horizontal cylinder.Recently, Nazar et al. [14], Huang and Chen [15] have considered the problem of natural convection flow from lower stagnation point to upper stagnation point of a horizontal circular cylinder immersed in a micropolar fluid.Effects of pressure work and viscous dissipation on natural convection flow around a sphere with radiation heat loss have been discussed in Akther [16] and Alim et al. [17] respectively.Hye et al. [18] have considered conjugate effects of heat and mass transfer on natural convection flow across an isothermal circular cylinder in presence of chemical reaction.
In the present study the focus is given on the effects of heat and mass transfer on the natural convection boundary layer flow across a sphere with radiation heat loss when the sphere is at a uniform temperature and a uniform mass diffusion.The basic equations are transformed to local non-similarity boundary layer equations, which are solved numerically using a very efficient finitedifference scheme together with the Keller-box method, described by Keller [19] and later by Cebeci and Bradshaw [20].To the best of our knowledge, this problem has not been considered before.
Numerical results have been shown in terms of local skin friction, rate of heat transfer, rate of species concentration, velocity profiles as well as temperature profiles for a selection of relevant physical parameters consisting of heat radiation parameter Rd, Prandtl number Pr and the Schmidt number Sc are shown graphically.Some results for Skin friction coefficient, rate of heat transfer and rate of species concentration for different values of Schmidt number Sc has been presented in tabular form as well.

Formulation of the Problem
Let us consider a steady two-dimensional laminar free convective flow that flows across a uniformly heated sphere of radius a, which is immersed in a viscous and incompressible fluid.It is assumed that the surface temperature of the sphere is w , where w .Here T T T   T  is the ambient temperature of the fluid, T is the temperature of the fluid in the boundary layer, g is the acceleration due to gravity,   r x is the radial distance from the symmetrical axis to the surface of the sphere and   ˆ, u v are velocity components along the   ˆ, x y axis.The coordinates system and the configuration are shown in Figure 1.
Under the usual Boussinesq and boundary layer approximation, the governing equations of the flow are: The boundary conditions for the Equations ( 1)-( 4) are ˆˆ0, , at 0 where    ˆŝin r x a x a   is the radial distance from the centre of the sphere, g is the acceleration due to gravity,  is the coefficient of concentration expansion, µ is the viscosity of the fluid, p C is the specific heat at constant pressure , ρ is the density, D is the molecular diffusivity of the species concentration and q r is the radiative heat flux in the y direction.In order to reduce the complexity of the problem and to provide a means of comparison with future studies that will employ a more detail representation for the radiative heat flux; we will consider the optically dense radiation limit.Thus the Rosseland diffusion approximation proposed by Siegel and Howell [21] and is given by simplified radiation heat flux term as: We now introduce the following non-dimensional variables: where and Gr is the Grashof number,  is the non-dimensional temperature and  is the non-dimensional species concentration function.Substituting the variables (7) into Equations ( 1)-(4) leads to the following non-dimensional equations The boundary conditions (5) become where Rd is the radiation-conduction parameter, Pr is the Prandtl number and N is the ratio of the buoyancy forces due to the temperature & concentration are defined respectively as and To solve Equations ( 8)-( 11), subject to the boundary conditions (12), we assume the following variables , , sin where,  is the non-dimensional stream function, which is related to the velocity components in the usual way as Substituting (14) into Equations ( 9)-( 11) we get, after some algebra the following transformed equations   along with boundary conditions It has been seen that the lower stagnation point of the sphere or , Equations ( 16)-( 18) reduce to the following ordinary differential equations: Subject to the boundary conditions at In the above equations primes denote the differentiation with respect to η.
In practical applications, the physical quantities of main interest are the shearin stress, the rate heat transfer and the rate of species concentration transfer in terms of the skin-friction coefficients C f , Nusselt number Nu and Sherwood number Sh respectively, which can be written as , and where Using the variables ( 7), ( 14) and the boundary condition (19a) into (24)-(25), we get and The values of the velocity, temperature and concentration distribution are calculated respectively from the following relations:

Results and Discussion
Equations ( 16)-( 18) subject to the boundary conditions were solved numerically using a very efficient implicit finite difference method together with the Keller-box scheme, which is described by Cebeci and Bradshow [20].The numerical solutions start at the lower stagnation point of the sphere i.e. at with initial profiles given by the Equations ( 20 16)-( 18 and , obtained in the present work and obtained earlier by Hye et al. [18] and Merkin [12] has been made in Table 1.It is clearly seen that there is an excellent agreement among the re-spective results.local skin-friction coefficient and the local Sherwood number and an increase in the local Nusselt number.This may be attributed to the fact that the increase in the values of implies more interaction of species concentration with the thermal boundary layer and less interaction with the momentum boundary layer.
Variation in the Schimdt number Sc (Sc = 0.73, 0.80, 1.00, 1.10, 1.30) are considered in Figures 3(a)-(c  Here it is found that the velocity profiles increases sharply and the non-dimensional temperature and concentration profile decrease significantly with the increase of radiation parameter Rd.
In Table 2

Conclusions
Conjugate effect of heat and mass transfer in natural convection flow from an isothermal sphere with radiation heat loss has been investigated numerically.New variables to transform the complex geometry into a simple shape have been introduced and the resulting non-similarity boundary layer equations are solved by a very efficient implicit finite difference method together with the Keller-box scheme [19].From the present investigation the following conclusions may be drawn:  Both the local skin-friction coefficient fx C and the local Nusselt number

Figure 1 .
Figure 1.Physical model and coordinate system.
)-(22) along with boundary conditions (23) and proceed round the sphere up to the upper ), numerical values of the local skin-friction coefficient fx , the local Nusselt number C x Nu and the local Sherwood number x Sh , are calculated respectively from the Equations (26)-(28) for the surface of the sphere from lower stagnation point to upper stagnation point.Numerical values of fx x and C Nu x Sh are depicted in Tables 1, 2 and Figures 2-5.A comparison of the local Nusselt number x Nu and fx C for different values of curvature parameter x while 1.0 Pr  The effect of the Prandtl number Pr, on the reduced

Figure 2 .Figure 3 .Figure 4 .
Figure 2. (a) Skin friction coefficient; (b) Rate of heat transfer and (c) Rate of species concentration for different values of Pr while Rd = 0.30 and Sc = 0.73.local skin-friction coefficient C fx , local Nusselt number Nu x and the local Sherwood number Sh x is shown in

Figure 5 .
Figure 5. (a) Velocity; (b) Temperature and (c) Concentration distribution for different values of Rd while while Pr = 0.72 and Sc = 0.73.ficient fx and the local Nusselt number C

x
Nu increase and the local Sherwood number x Sh decrease when the value of the radiation parameter increases. As Rd increases both the concentration distribution and the temperature profile decrease significantly and the velocity profile increases slightly at π 2 x  of the surface. An increase in the values of Pr leads to an increase in the values of the local skinfriction coefficients fx C , local Nusselt number x Nu and the local Sherwood number x Sh . The skin-friction coefficient fx C and local Nusselt number x Nu decrease and the rate of species con-

Table 2 . Skin friction coefficient, rate of heat transfer and rate of species concentration against x for different values of the Schmidt number Sc with other controlling parame- ters Pr = 0.72 and radiation parameter Rd = 0.30.
Rd  ,it is the concentration field which induce the boundary flow.It can be stated that an increase in the values of Rd leads to an increase in the values of the local skin-friction coefficients fx , local Nusselt number C we have entered the numerical values fx C Rd  .It is observed that for increasing axial distance parameter x, the values of skinfriction increase and the values of rate of heat transfer and that of the rate of species concentration decrease for the values of Schmidt number Sc.On the other hand both the fx and C dingly as Schmidt number Sc increases.And these changes in , Sh due to an increase in Sc are consistent with a free convection boundary layer.
x Nu decrease and x Sh increases accorx