Chemical Reaction and Thermal Diffusion Effects on Mass Transfer Flow through an Inclined Plate

A numerical investigation of boundary layer mass transfer flow through an inclined plate with the effect of chemical reaction and thermal diffusion is presented in this study. The governing partial differential equations (PDE) are transformed to a system of dimensionless non-similar coupled PDEs. The transformed, non-similar conservations equations (momentum balance equation, energy balance equation and concentration balance equation) are then solved using a numerical approach known as explicit finite difference method (EFDM). Basically EFDM introduced for the unsteadiness in the momentum, temperature, and concentration fluid fields is based on the time dependent fluid velocity, temperature and concentration of the boundary surface. During the course of discussion, it is found that the various parameters related to the problem influence the calculated resultant expressions. The computed numerical solution results for the velocity, temperature, and concentration distribution with the effect of various important dimensionless parameters (Grashof number, Modified Grashof number, Prandtl number, Schmidt number, Soret number, Dufour number, chemical reaction parameter and inclination parameter) entering into the problems are critically analyzed and discussed graphically. It can be seen that two physical phenomena chemical reaction and thermal diffusion can greatly effect on the boundary layer fluid flows through an inclined plate.


Introduction
The effect of thermal diffusion and chemical reaction on heat and mass transfer fluid flow has a great importance to the engineering community and to investigators dealing with the problems in many industrial processes and technological fields which experience not only temperature difference but also concentration difference.The concentration difference may sometimes produce quantitative and qualitative changes to the rate of heat transfer.Some industrial applications are found in food drying, food processing or wood and paper production.Besides, the characteristics of the flow of blood through arteries and veins are of considerable medical interest.In view of the above interests, a series of investigations have been made by different scholars.Pera and Gebhart [1] were the first authors to investigate and combine buoyancy effects of thermal and mass diffusion on natural conviction flow.Soundalgekar and Ganesan [2] presented a finite difference analysis of unsteady natural convective flow past a semi-infinite plate with mass transfer.Hydro-magnetic heat and mass transfer in flow of a viscous incompressible fluid past an infinite vertical porous plate was analyzed by Singh et al. [3].Using a finite difference technique, an unsteady free convective heat and mass transfer flow from an infinite vertical plate with constant suction have been numerically investigated by Ambethkar [4].Recently Sivaiah et al. [5] studied heat and mass transfer effects on MHD free convective flow past a vertical porous plate.
The problem of mass transfer flow through an inclined plate has generated much interest from astrophysical, renewable energy system and also hypersonic aerodynamics researchers for a number of decades.Umemura and Law [6] developed a generalized formulation for the natural convection boundary layer flow over a flat plate with arbitrary inclination.They found that the flow characteristics depend not only on the extent of inclination but also on the distance from the leading edge.Chamkha and Khaled [7] investigated the problem of coupled heat and mass transfer by hydromagnetic free convection from an inclined plate in the presence of internal heat generation or absorption, and similarity solutions were presented.Reddy and Reddy [8] performed an analysis to study the natural convection flow over a permeable inclined surface with variable temperature, momentum and concentration.Singh [9] studied heat and mass transfer in MHD boundary layer flow past an inclined plate with viscous dissipation in porous medium.Micropolar fluid behavior on MHD free convection and mass transfer with constant heat and mass fluxes is studied numerically by Ali et al. [10].Recently Islam et al. [11] have studied mass transfer flow through an inclined plate with porous medium.
In many transport processes existing in nature and in industrial applications in which heat and mass transfer is a consequence of buoyancy effects caused by diffusing of heat and chemical species.The study of such process is useful for improving a number of chemical technologies, such as polymer production, enhanced oil recovery, underground energy transport, manufacturing of ceramic and food processing.In many hydrometallurgical and chemical technology industries, heat and mass transfer is of considerable importance due to chemical reaction effects.Basically chemical reaction takes place between an external mass and the fluid in many chemical engineering processes.It can be designated as either a mixed or stable process, which be contingent on whether it occurs at a boundary or as a single-phase volume reaction.Chaudhary and Jha [12] studied the effect of chemical reaction on MHD micropolar fluid flow past a vertical plate in slip flow regime.Effects of mass transfer with a chemical reaction on unsteady flow past an accelerated isothermal vertical plate are investigated by Muthuchumaraswamy et al. [13].Rajesh et al. [14] have discussed the effects of chemical reaction and radiation effects on MHD flow past an infinite vertical plate with variable temperature.The importance of thermal-diffusion and diffusion-thermo effects for various fluid flows has been studied by Eckert and Drake [15].Olajuwon [16] examined convection heat and mass transfer in a hydromagnetic flow of a second grade fluid past a semi-infinite streching sheet in the presence of thermal diffusion and thermal radiation.Kumar et al. [17] have investigated thermal diffusion and radiation effects on unsteady MHD flow through porous medium with variable temperature and mass diffusion in the presence of heat source or sink.
There have been few studies on heat transfer boundary layer flow in the literature however mass transfer flow through an inclined plate with the effect of chemical reaction and thermal diffusion effects still need more attraction to the researchers [18] [19].
The objective of this research is to investigate the fluid behavior with the effect of chemical reaction and thermal diffusion on mass transfer flow in an inclined plate.In this research the governing equations of the problem contain a system of partial differential equations which are transformed by usual transformation into a non-dimensional system of partial coupled non-linear differential equations.The obtained non-similar partial differential equations are solved numerically by EFDM.The results of this research are discussed for the differ-ent values of the well-known dimensionless parameters and are shown graphically.

Mathematical Model of Flow
By introducing Cartesian co-ordinate system, the X-axis is chosen along the plate in the direction of the flow and the Y-axis is normal to it.The physical configuration has been shown in Figure 1.Initially it has been considered that the plate as well as the fluid is at the same temperature ( ) T T ∞ and concentration level ( ) C C ∞ everywhere in the fluid is same.Also it is considered that the fluid and the plate is at rest after that the plate is to be moving with a constant velocity 0 U in its own plane and instantaneously at time 0 t > the species concen- tration and the temperature of the plate are raised to , which are there after maintained constant, where , are the concentration of the species far away from the plate.Within the framework of the above stated assumptions with reference to the generalized equations described before the equation relevant to the transient two dimensional problems are governed by the following system of coupled non-linear differential equations using the Boussinesq and boundary layer approximations.

Continuity Equation
Energy Equation 22

Concentration Equation
( ) where , x y are Cartesian co-ordinate, , u v are the components of flow velocity along , x y directions respec- tively.Here g is the local acceleration due to gravity; υ is the kinetic viscosity; ρ is the density of the fluid, m D is the coefficient of mass diffusivity, T D is the coefficient of Thermal diffusivity and C K be the chemical reaction parameter.

Mathematical Formulation
Since the solution of the governing equations under the initial and boundary conditions will be based on a finite difference method it is required to make the said equations dimensionless.
For this purpose it has been now introduced the following dimensionless variables; where τ represents the dimensionless time, X & Y be the dimensionless Cartesian coordinates, where U and V be the dimensionless velocity components, T and C is the dimensionless temperature and concentration respectively.
Using the above relation, we obtain the following non-dimensional coupled partial differential equation, where, ( ) Also the associate initial and boundary condition become At

Numerical Solution
The explicit finite difference method (EFDM) has been used to solve the governed second order nonlinear coupled dimensionless partial differential equations with the corresponding initial and boundary conditions.To obtain a system of finite difference equations, the flow region is divided into a grid or meshes of lines parallel to X and Y axes where X-axis is taken along the plate and Y-axis is normal to the plate.Here it is considered that max X (=400) i.e.X varies from 0 to 400 and regarding max Y (=25) as corresponding to Y → ∞ i.e.Y varies from 0 to 25.It is also considered that 60 m = and 60 n = grid spacing in X and Y directions respectively.We have the constant mesh size along X direction, ( ) and the constant mesh size along Y direction, ( ) with the smaller time-step 0.05 τ ∆ = Let , and denote the values of , and U C T at the end of a time-step respectively.Using the finite difference approximations we obtain the following set of finite difference equations, , cos And the initial and boundary conditions with finite difference scheme are; , , , , 0, 0, 0, 0, Here the subscripts and i j designate the grid points and x y coordinates respectively and the superscript n represents a value of time, n τ τ = ∆ where, 0,1, 2,3, n =  .From the initial condition, the values of U is known at 0 τ = .Then at the end of the anytime-step τ ∆ , the new temperature T ′ and concentration C′ , the new velocity U ′ at all interior nodal points may be obtained by successive application of temperature, concen- tration and momentum equation respectively.This process is repeated in time and provided the time-step is sufficiently small, hence , , and U V T C should eventually convergent values which approximate the steady-state solution of the problem.

Results and Discussion
To discuss the effects of the associated non-dimensional parameters on the flow variables, the numerical solutions of the problem are obtained by the EFDM numerical methodology with the help of a computer programming language Compaq Visual Fortran.In and Grashof number as ratio of the buoyancy to viscous force acting on a fluid was considered as 0.1, 0.2 and 0.5.An increase in Grashof number clearly enhances profiles of velocity, therefore induces a robust speeding up in the flow.A distinct velocity shoot arises for all profiles near the boundary surface (Y = 0) and this is emphasized with increasing Grashof number.With increasing Grashof number the thermal buoyancy force is augmented which aids in momentum change in the boundary layer.Also the boundary layer thickness of velocity is therefore increased as well.With superior intervene of time, the velocity was found to be considerably enhanced.
The same effect on the velocity curve is found in Figure 4 that is the velocity decreases with increase of Prandtl number r P (=0.71, 1.0 and 7.0).This is because Prandtl number implies the relative influence of momentum diffusion to thermal diffusion in the periphery layer regime.When P r = 1, the momentum diffusion rate surpasses thermal diffusion rate.As a result the velocity in the fluid regime will be decreased with a rise in Pr.The dimensionless velocity distribution of fluid was found to be strongly reduced for the decreasing effect of     chemical reaction parameter r K (=0.6, 1.0 and 1.2) found on Figure 5. Figure 6 shows that the steady state temperature of fluid decreases with the increase of Prandtl number r P (=0.71, 1.0 and 7.0).This is because the thermal boundary layer thickness decreases due to increase in this parameter.With increasing period, temperatures are yet again detected to be strongly enhanced through the boundary surface.
In Figure 7 it was observed that the fluid temperature remain unchanged with the increasing value of chemical reaction parameter r K (=0.6, 1.0 and 1.2).A decreasing effect of dimensionless mass distribution curves are observed for increasing Schmidt number C S in Figure 8 that is fluid concentration decreases with increase of Schmidt number C S (=0.96, 1.22 and 1.60).Therefore the ratio of diffusivity of momentum i.e., viscosity and diffusivity of mass has a remarkably effects on fluid concentration.And the Schmidt number can be used to describe fluid flows in which there are instantaneous processes of momentum and mass diffusion convection.The unchanged effect of Soret number o S on concentration profiles are observed in Figure 9 where the mass distribution of fluid is remain same with increasing value of Soret number o S .The ratio of thermos-diffusion and diffusion coefficient was assumed to be 0.5, 0.8 and 1.0.In Figure 10 it can be seen that the fluid concentration remains unchanged with the increasing value of chemical reaction parameter r K (=0.6, 1.0 and 1.2).From the above figures we also observed that fluid behavior is varies with dimensionless time τ .With the rise of dimensionless time τ the fluid velocity, temperature and concentration are gradually rises.

Conclusions
Mixed (combined) convection boundary layer flows are of great attention because of their innumerable industrial, engineering and scientific applications in heat and mass transfer.In this paper, a boundary layer analysis for mixed convection heat and mass transfer fluid flow through an inclined plate in the presence of chemical reaction and thermal diffusion is considered.The governing of boundary layer flow such as momentum, energy and mass balance is then numerically analyzed by explicit finite difference method (EFDM).The results are presented graphically with the effect of various dimensionless parameters.The significant findings obtained from the graphical observation are listed below:  The velocity distribution decreases with the increase of inclined angle ( ) α .
 The velocity distribution increases with the increase of Grashof number ( ) are species concentration and temperature at the wall of the plate and ,

Figure 1 .
Figure 1.The Physical model and coordinate system.O

2 .
order to analyze the physical solution of the model, we have computed the steady numerical values of the non-dimensional Velocity U , Temperature T , and Concentration C within the boundary layer for different values of Inclination ( ) α , Grash of number ( ) Chemical reaction parameter ( ) r K .To obtain the steadystate solutions of the computation, the calculations have been carried out up to non-dimensional time τ = 10 to 60. Hence the velocity, temperature and concentration profile are drawn for dimensionless time τ = 10, 30 and 60.In the Figures 2-10 x-axis taken as dimensionless co-ordinate variable (Y) which is varies 1 to 25 and y-axis taken as fluid velocity, temperature and concentration respectively.The effect of the inclination on the velocity field is presented in Figure It is observed that the dimensionless velocity distribution of fluid decreases with the rise of inclination α (=45˚, 50˚ and 60˚) for different time steps.Therefore angle of inclination has a great impact on the velocity of fluids.Figure 3 represents the evolution of dimensionless velocity with Y for different values of dimensionless times ( ) 10,30 and 60 τ =