Unsteady Hydro-Magnetic Heat and Mass Transfer Flow of a Non-Newtonian Power-Law Fluid past a Flat Plate in the Presence of Homogeneous Chemical Reaction ()
1. Introduction
The fluids which are encountered in chemical and allied processing applications are known as non-Newtonian- fluids. The study of non-Newtonian fluid flows has considerable interest for their numerous engineering appli- cations. During the past four decades the study of non-Newtonian fluids has gained interest because of their nu- merous technological applications, including manufacturing of the plastic sheets, performance of lubricants and
movement of biological fluids. To explain the behavior of non-Newtonian fluid different models have been proposed. Among these the power law fluid has gained importance. The order of chemical reactions depends on several factors. One of the simplest chemical reactions is the first order reaction in which the rate of reaction is directly proportional to the species concentration. Now a days, due to the growing use of these non-Newtonian substances in various manufacturing and processing industries, considerable efforts have been directed towards understanding their friction and heat transfer characteristics. By the application of a magnetic field hydromagnetic techniques are used for the purification of molten metal. The problem of steady flow and heat transfer in power law fluid by free convection along a vertical plate has been investigated by many researchers. Vujanovic et al. [1] investigated a variational solution of the Rayleigh problem for a power law non-Newtonian conducting fluid. Padhy & Pattnayak [2] studied the mass transfer and free convective effects of a power law fluid past an impulsively started vertical plate. Murthy [3] investigated effects of double dispersion on mixed convection heat and mass transfer in a non-Darcy porous medium. Muthucumaraswamy et al. [4] discussed on diffusion and first order chemical reaction on impulsively started infinite vertical plate with variable temperature. Naseer [5] investigated the problem of unsteady free convection with heat & mass transfer from an isothermal vertical flat plate to a non-Newtonian power law fluid immersed in a saturated porous medium. Chamkha et al. [6] investigated unsteady natural convective power law fluid flow past a vertical plate embedded in a non-Darcian porous medium in the presence of a homogeneous chemical reaction. Khan et al. [7] discussed non-Newtonian MHD mixed convective power law fluid flow over a vertical stretching sheet with thermal radiation, heat generation and chemical reaction effects. Olajuwon et al. [8] studied convection heat mass transfer in a power law fluid with non-constant relaxation time past a vertical porous plate in the presence of thermo and thermal diffusion. Olajuwon [9] examined effects of thermo diffusion and chemical reaction on heat and mass transfer in a power law fluid over a flat plate with heat generation. Then in 2014 many researchers had shown interest to examine on the subject heat and mass transfer in a non-Newtonian fluid. Uwanta et al. [10] investigated heat and mass transfer flow past an infinite vertical plate with variable thermal conductivity heat source and chemical reaction. Jothimani and Vidhya [11] studied non-Newtonian fluid flow and heat transfer over a non-linearly stretching surface along with porous plate in porous medium. Recently Madhu et al. [12] studied effect of viscous dis- sipation and thermal stratification on chemical reacting fluid flow over a vertical stretching surface with heat source.
The aim of the present work is to investigate the unsteady hydro magnetic non-Newtonian power law fluid past a flat plate with heat and mass transfer effect. The governing equations, describing the model are highly nonlinear coupled partial differential equations in nature. Hence closed form solutions are not possible. Suitable implicit finite difference scheme has been used to get the solution of the problem. Graphs have been plotted against various flow parameters to study the characteristics of velocity, temperature and concentration of the fluid.
2. Formulation of the Problem
There exist different types of non-Newtonian fluids but the simplest and most common type is the power-law fluid for which the rheological equation of the state between stress components and strain rate components defined by Vujanovic is
(1)
where, P is the pressure, is the Kronecker delta, K and n are the consistency and flow behavior indices of the fluid respectively. When n > 1 the fluid is described as dilatant, n < 1 as pseudo-plastic and when n = 1 it is known as the Newtonian fluid.
Consider the unsteady free convection heat and mass transfer flow of a two-dimensional, viscous, incompressible, electrically conducting and chemically reactive non-Newtonian power-law fluid along an infinite non- conducting vertical flat plate in the presence of a uniform magnetic field B0 applied in a transverse direction to fluid flow. Let x′-axis be along the plate in upward direction, y′-axis is normal to it & z′-axis is normal to x′y′-plane. Initially, at time, the fluid and plate are at rest and at a uniform temperature. When the plate is maintained at constant temperature and constant species concentration. Since the plate is of infinite extent in x′ direction and is electrically non-conducting, except pressure all other physical quantities are functions of y′ and t′ only. The governing equations describing the model are
(2)
(3)
(4)
(5)
where g is acceleration due to gravity, α represents the thermal diffusivity, βT is coefficient of thermal expansion of fluid, βC is volumetric coefficient of expansion or contraction, k is thermal conductivity of the fluid, ρ is fluid density, n is power law index, u′ & v′ are stream wise and transverse velocity respectively.
Similarly x′ and y′ are stream wise and transverse co-ordinate. T' is temperature of the fluid and t′ is time, D is the coefficient of mass diffusivity, kc is the rate of chemical reaction, and are the free stream temperature and concentration of the fluid respectively, and and are the temperature and concentration at the wall respectively.
The initial and boundary conditions are
(6)
The dimensionless variables are defined as follows:
where, and l is the suitable length scale. Substituting the above non-dimensional variables into
Equations (2)-(5) yield the following dimensionless equations
(7)
(8)
(9)
(10)
where is the Reynold number, is the Grashof number, is the modified Grashof number, is the Prandtl number, is the Schmidt number, is the magnetic parameter, is the chemical reaction parameter, is the specific heat at constant pressure, is the kinematic viscosity and µ is the constant viscosity of the fluid in boundary layer region.
Accordingly, the initial and boundary conditions will be reduced to
(11)
The special significance of this type of flow with heat and mass transfer situation are the skin-friction coefficient Cf, the local Nusselt number Nu and Sherwood number Sh. These physical quantities are defined in non- dimensional form, respectively, as follows:
(12)
(13)
(14)
3. Solution of the Problem
The Equations (7)-(10) are solved by implicit finite difference method. For discretization in space and time a uniform mesh of step and along x & y direction respectively and time are employed so that the grid points are for, &. The discretized form of Equation (7), (8), (9) and (10) are obtained respectively as,
(15)
(16)
(17)
(18)
The above discretized Equations (15)-(18) are solved iteratively using the following algorithm.
Step I
Initialize
Step II
For
For
For
Step III
Step IV
Step V
The Steps (II)-(V) are repeated until the relative errors of two consecutive values of are less than a given tolerance.
4. Results & Discussion
The non-linear governing Equations (7)-(10) with the boundary conditions (11) are solved using finite difference method. The velocity, temperature, and concentration of the fluid for different Reynold numbers are shown in Figures 1(a)-(c). The velocity, temperature and concentration decrease as Re increases.
The velocity and temperature of the fluid for different Prandtl numbers are shown in Figure 2(a) and Figure 2(b). Prandtl number increases the viscous diffusivity of the fluid at the surface which enhanced the velocity of the fluid near the surface as depicted in Figure 2(a), increase in Pr implies flow of liquid with low thermal diffusivity and high viscous stress, which increases thermal boundary layer thickness near the surface as shown in Figure 2(b).
In Figure 3(a) and Figure 3(b) we have seen as chemical reaction parameter increases the velocity and species concentration decrease, but it is reverse in the case of Schmidt number. The velocity and species concentration increase as Sc increases it is reflected through Figure 4(a) and Figure 4(b).
It is evident from Figure 5 that the presence of transverse magnetic field has a retarding effect on velocity field. But from Figures 6-8 it is observed that with an increase in Gr, Gm or N the velocity increases.
For the physical interest in view we found the influence of power law index, magnetic parameter, Prandtl number, Reynold number, Schmidt number, chemical reaction parameter, thermal Grashof number and modified
(a) (b) (c)
Figure 1. Variation of Re on (a) fluid velocity, (b) temperature and (c) species concentration when N = 1; Gr = 5; Gm = 5; M = 1; Pr = 1; Sc = 0.5; Kr = 1.
(a) (b)
Figure 2. Variation of Pr on (a) fluid velocity and (b) temperature when N = 1; Gr = 5; Gm = 5; M = 1; Re = 1; Sc = 0.5; Kr = 1.
(a) (b)
Figure 3. Variation of Kr on (a) Velocity and (b) species concentration when N = 1; Gr = 5; Gm = 5; M = 1; Re = 1; Pr = 1; Sc = 2.
Grashof number on the skin friction Cf, local Nusselt number Nu and Sherwood number Sh is shown in the Table 1. It is interesting to note that the increase in magnetic parameter decreases the velocity of the fluid that helps to reduce the skin friction at the surface, but the local heat transfer and mass transfer are not influenced by the magnetic parameter. When the fluid is dilatant the skin friction at the surface decreases, but local Nusselt number does not change. As Prandtl number increases the skin friction increases, local Nusselt number decreases but
(a) (b)
Figure 4. Variation of Sc on (a) fluid velocity and (b) species concentration when N = 1; Gr = 5; Gm = 5; M = 1; Re = 1; Pr = 1; Kr = 0.4.
Figure 5. Variation of M on velocity when N = 1; Gr = 5; Gm = 5; Pr = 1; Sc = 0.6; Kr = 1.
Figure 6. Variation of N on velocity when Gr = 5, Gm = 5, M = 1, Re = 2, Pr = 1, Sc = 0.6, Kr = 1.
Figure 7. Variation of Gr on velocity when N = 1, Gm = 5, M = 1, Re = 1, Pr = 1, Sc = 0.6, Kr = 0.4.
Figure 8. Variation of Gm on velocity when N = 1; Gr = 5; M = 1; Re = 1; Pr = 1; Sc = 0.6; Kr = 0.4.
Table 1. Skin friction, Nusselt number and Sherwood number.
Cf = Skin Friction :; Nu = Local Nusselt Number:; Sh = Local Nusselt Number:.
concentration does not change with it. It is evident from the table that with an increase in Re the skin friction decreases while a reverse effect is seen in case of Nusselt number Nu and Sherwood number Sh is very much affected by Reynold number.
5. Conclusions
Unsteady free convective heat and mass transfer in the flow of a two dimensional viscous incompressible electrically conducting and chemically reactive non-Newtonian power-law fluid along an infinite non-conducting vertical flat plate in the presence of uniform magnetic field are studied. It is found that,
・ With an increase in Re velocity, temperature & concentration of the fluid decrease.
・ With the increasing value of chemical reaction parameter fluid velocity & concentration decrease near the plate, but the species concentration shows reverse characteristics as depicted in the skin friction table.
・ Magnetic field has a retarding effect on the fluid flow while the thermal radiation has a reverse effect on it.
The solutions obtained are well agreed with the Newtonian case and they give improved results, taking into consideration of the behaviour of the magnetic field. This method well suits for other non-Newtonian fluid flow problems.