Original Article
Int. J. Chem. Sci, Volume: 15( 3)

Radiation and Chemical Reaction Effects on MHD Free Convective Heat and Mass Transfer Flow of a Viscoelastic Fluid Past a Porous Plate with Heat Generation/Absorption

*Correspondence:
Ramana Murthy MV, Department of Mathematics and Computer Science, Osmania University, Hyderabad, Tel: 040-2768 2444; E-mail: [email protected]

Received: May 02, 2017; Accepted: August 03, 2017; Published: August 07, 2017

Citation: Santhosha B, Younus S, Kamala G, et al. Radiation and Chemical Reaction Effects on MHD Free Convective Heat and Mass Transfer Flow of a Viscoelastic Fluid Past a Porous Plate with Heat Generation/Absorption. Int J Chem Sci. 2017;15(3):170

Abstract

Radiation and chemical combined effects on MHD free convective heat and mass transfer flow of viscous, incompressible, conducting elastic fluid through porous medium finite by a porous plate within the presence of heat generation/absorption is studied. The momentum, energy and mass diffusion equations are coupled non-linear partial differential equations; these are solved by applying two term perturbation methodology and mistreatment software system MATLAB. The options of the fluid flow, heat and mass transfer characteristics are analyzed by plotting graphs and mentioned qualitatively.

Keywords

MHD; Chemical reaction; Radiation; Viscoelastic; Skin-friction; Porous medium

Introduction

The heat and mass transfer in free convective through porous media has been broadly investigated in the current years due to its applications in engineering. The study of non-Newtonian fluid flow by engineers and experts has risen in recent times due to its useful application in various fields. Keeping its importance in mind an extensive range of mathematical models has been established to simulate the various hydrodynamic behaviors of these non-Newtonian fluids. Various viscoelastic fluids classical like the Rivlin-Ericksen 2nd order model, Oldroyd model and Johnson-Seagalman model. The fluid, which exhibits the elasticity property of solids and viscous property of liquids are called viscoelastic fluid. Fluid models of this type come across petrochemical, biomedical and environmental engineering including polypropylene combination sintering.

Beard and Walter [1] had presented the boundary layer treatment for unidealized viscoelastic fluid. The heat transfer in convection flow of viscoelastic fluid of Walter model was investigated by Rajagopal [2]. Ibrahim [3] discussed the effects of the radiation and chemical reaction on the unsteady MHD convection flow through semi-infinite vertical porous plate with heat source and suction.

Muthucumaraswamy and Ganesan [4] studied the unsteady flow thorugh impulsively semi-infinite vertical plate, which subjected to uniform heat flux of a chemical reaction effect. Chamkha [5] was studied the MHD flow of a numerical of uniformly stretched vertical permeable surface in the presence of heat and chemical reaction. Chamkha [6] discussed unsteady convective heat and mass transfer through semi-infinite porous plate with heat absorption. Hossain [7] were analyzed natural convection flow of a vertical surface with uniform temperature in the presence of heat generation/absorption. Hady [8] are studied the problem of free convection flow along a vertical wavy surface embedded in electrically conducting fluid saturated porous media in the presence of internal heat generation or absorption effect. Mohamed [9] has discussed double-diffusive convection radiation interaction on unsteady MHD flow over a vertical moving porous plate with heat generation and Soret effects.

Shanker et al. [10] were presented a numerical solution for radiation and mass transfer effects on unsteady MHD free convective fluid flow embedded in a porous medium with heat generation/absorption using Galerkinfinite element method. Kumar and Sivaraj [11] studied MHD Visco-elastic fluid non-Darcy flow along a moving vertical cone. Vidyasagar and Raman [12] reported a study on the radiation effect on MHD convection flow of Kuvshinshiki fluid with mass transfer past a vertical porous plate through porous medium. Bhikshu [13] investigated peristaltic flow of a conducting Williamson fluid in a vertical asymmetric channel with heat transfer through porous medium. Devika [14] investigated MHD oscillatory flow of a visco-elastic fluid in a porous channel with chemical reaction. Bhikshu [15] studied the effects of magnetohydrodynamics on the peristaltic flow of fourth grade fluid in an inclined channel with permeable walls. Sivaraj and Kumar [16] were studied an unsteady MHD dusty visco-elastic fluid Couette flow in an irregular channel with varying mass diffusion. Authors are discussed in many chemical engineering processes; the chemical reaction does occur between mass and fluid in which plate is moving [17-20].

These processes occur in varied industrial applications like compound production, producing of ceramics or glasswork and food process. The aim of this work is to investigate the influence of radiation and chemical on MHD convective heat and mass transfer flow of viscous, in-compressible, electrically conducting elastic fluid through porous medium delimited by a porous plate within the presence of heat generation and absorption. Regular perturbation method is adopted to solve governing equations are solved by the perturbation technique [21,22].

Mathematical Formulation of the Problem

We think about two-dimensional unsteady convective flow of Viscous-elastic fluid (Kuvshinski fluid), electrically conducting, an incompressible, past a semi-infinite vertical periodical porous plate in the presence of magnetic field, chemical reaction, heat generation and radiation absorption. Let x* -axis is taken along the porous plate in the rising direction and y* -axis is standard to it. It is assumed that there is no applied voltage of which implies the absence of an electric field. The squarely applied magnetic field and magnetic Reynolds number are very small and hence the made magnetic field is very small. The MHD term is defined from an order of magnitude analysis of the Naiver-stokes equation. We take the porous medium as a collection of minor identical spherical particles fixed in space. The chemical reactions are occurring in the flow and physical properties are taken constant. The magnetic and viscous dissipations are very small. Due to the assumption that the plate in x* -direction is of unlimited length, all the flow variables except pressure are functions of x* and t* only.

The governing boundary layer equation of the flow field considering the above assumptions are given by:

equation (1)

equation (2)

equation (3)

equation (4)

Where u,v are the velocity components in x, y directions respectively, t -the time, ρ-the fluid density , v -the kinematic viscosity, Cp -the specific heat at constant pressure, g -the acceleration due to gravity, β and β*-the thermal and concentration expansion coefficient, respectively, B0-the magnetic induction, D-Coefficient of chemical molecular diffusivity, K1 -Reaction rate constant, R1 -Radiation absorption coefficient, s1 -Heat source/sink constant, K -Thermal conductivity, K0 -Permeability of the porous medium, equation-the dimensional visco-elastic parameter, T -the dimensional temperature,C -the dimensional concentration.

Here the boundary conditions are:

equation (5)

Where U0 is the plate velocity, TW and CW are the wall dimensional temperature and concentration, respectively, T and C are the free stream dimensional temperature and concentration, respectively, ω -the constant. From equation (1) we get,

equation (6)

Where V0 is the suction velocity normal to the plate.

Introducing the following non-dimensional quantities:

equation (7)

Using the above non-dimensionless quantities and equation (6) the equations (2)-(4) will become,

equation (8)

equation (9)

equation (10)

Where equation

The corresponding boundary conditions are,

equation (11)

Solution of the Problem

The equations 8-10 with boundary conditions (11) are solved, by using the perturbation technique. Let the velocity, temperature and concentration fields as:

equation (12)

Where ∈ is a small quantity equation are the thermal Grash of number, Solutal Grashof number, Magnetic parameter, viscoelastic parameter, Prandtl number, heat absorption parameter, radiation absorption parameter, Schmidt number and chemical reaction parameter, respectively.

Substituting equation (12) into the equations (8)-(10) and equating the like terms, neglecting the terms of ∈2 , then we get,

equation (13)

equation (14)

equation (15)

equation (16)

equation (17)

equation (18)

Where equation

Here the primes denote the differentiation with respect to y.

The boundary conditions are:

equation (19)

In view of (19), solving the equation (13) – (18)

We obtain,

equation (20)

equation (21)

equation (22)

equation (23)

equation (24)

equation (25)

On substituting the expressions of equation in the equations (12) the

expressions for velocity, temperature and concentration are:

equation (26)

equation (27)

equation (28)

Skin-friction at the plate in dimensionless form is given by:

equation (29)

Nusselt number

Knowing the rate of heat transfer on the well in terms of Nusselt number Nu is given by:

equation (30)

Sherwood number

Knowing the rate of mass transfer on the wall is given by:

equation (31)

Result and Discussion

The physical characteristic of the problem is obtained for the above physical quantities by computing numerically for different values of the governing parameters viz. the Grashof number (Gr), Modified Grashof number (Gm), Prandtl number (Pr), Schmidt number (Sc), magnetic parameter, Permeability parameter (K), Eckert number (Ec) and Heat absorption parameter (Q), heat source parameter (S) radiation parameter (F) and chemical reaction parameter (Kr).

Figure 1 shows the impact of field on magnetic profiles within the boundary layer. Here we have an observed to determined that the velocity profiles are decreases with increasing the magnetic parameter (M). The transverse magnetic field opposes the motion of the fluid and the velocity of transport is significantly reduced. It’s due to with the increase within the magnetic parameter (M), force will increase and it produces a lot of resistance to the flow. Because it will increase, the thermal boundary layer thickness will increase however the momentum boundary layers get belittled. The contribution of the porosity parameter of the fluid bed on the speed field is illustrated in Figure 2. In general, it's noticed that, because the porousness will increase, the speed profiles is additionally will increase for different parameters area unit mounted constants.

international-journal-chemical-sciences-velocity-profiles

FIG. 1. Effects of velocity profiles for different values of magnetic parameter (M).

international-journal-chemical-sciences-permeability-parameter

FIG. 2. Effects of velocity profiles for different values of permeability parameter (K).

Figure 3 shows that the impact of solutal Grashof range (Gm) on velocity profiles. It's detected that the dimensionless velocity will increase with associate degree increasing the Grashof range. Figure 4 illustrates distinct velocity profiles within the boundary layer for numerous values of solutal Grashof range (Gm). The solutal Grashof range (Gm) defines the quantitative relation of the species buoyancy force to the viscous hydro dynamic force. It’s detected that the velocity will increase with associate degree increasing values of the solutal Grashof range. Additionally, it's found that an increase in metric weight unit ends up in acceleration within the velocity profiles. Figures 5 to 11 illustrate velocity and temperature profiles for various values of Prandtl range (Pr). The numerical results show that the impact of increase values of Prandtl range leads to a decreasing velocity. From Figures 11 to 14, it's detected that an increase within the Prandtl range results a decrease of the thermal physical boundary layer and normally lower average temperature among boundary layer. In the case of smaller Prandtl numbers as the boundary layer is thicker and the rate of heat transfer reduced.

international-journal-chemical-sciences-Grashof-number

FIG. 3. Effects of velocity profiles for different values of Grashof number (Gr).

international-journal-chemical-sciences-modified-Grashof

FIG. 4. Effects of velocity profiles for different values of modified Grashof number (Gm).

international-journal-chemical-sciences-Prandtl-number

FIG. 5. Effects of velocity profiles for different values of Prandtl number (Pr).

international-journal-chemical-sciences-Schmidt-number

FIG. 6. Effects of velocity profiles for different values of Schmidt number (Sc).

international-journal-chemical-sciences-chemical-reaction

FIG. 7. Effects of velocity profiles for different values of chemical reaction parameter (SC).

international-journal-chemical-sciences-radiation-parameter

FIG. 8. Effects of velocity profiles for different values of radiation parameter (F).

international-journal-chemical-sciences-absorption-parameter

FIG. 9. Effects of velocity profiles for different values of heat absorption parameter (Q).

international-journal-chemical-sciences-visco-elastic

FIG. 10. Effects of velocity profiles for different values of visco-elastic parameter (λ).

international-journal-chemical-sciences-temperature-profiles

FIG. 11. Effects of temperature profiles for different values of Prandtl number (Pr).

international-journal-chemical-sciences-absorption-parameter

FIG. 12. Effects of temperature profiles for different values of heat absorption parameter (R).

international-journal-chemical-sciences-radiation-parameter

FIG. 13. Effects of temperature profiles for different values of radiation parameter (F).

international-journal-chemical-sciences-heat-source

FIG. 14. Effects of temperature profiles for different values of heat source (Q).

From Figure 6 and Figure 15 we observe that for different values of the Schmidt number (Sc), the velocity and concentration profiles are plotted and that the reductions in the velocity and concentration profiles are accompanied by simultaneous reductions in the velocity and concentration boundary layers, which is evident.

international-journal-chemical-sciences-Schmidt-number

FIG. 15. Effects of temperature profiles for different values of Schmidt number (Sc).

Figure 7 and Figure 16 are shows that the effects of the chemical reaction (Kr) on the velocity and concentration profiles, respectively. It is noticed that as chemical reaction increases, the considerable reduction in the velocity profiles and the presence of the peak indicates that the maximum value of the velocity occurs in the body of the fluid close to the surface but not at the surface. Also, with an increase in the chemical reaction parameter, the concentration decreases. It is evident that the increase in the chemical reaction (Kr) significantly alters the concentration boundary layer thickness but does not alter the momentum boundary layers.

international-journal-chemical-sciences-concentration

FIG. 16. Effects of concentration profiles for different values of Schmidt number (Sc).

Figure 8 and Figure 13 are shows that the influence of the thermal radiation parameter (F) on the velocity and temperature (Figures 10 and 11). It is observed that an increase in the radiation parameter results in decreasing velocity and temperature within the boundary layer. The effect of increasing the value of the Heat absorption parameter (Q) is to decrease the boundary layer as shown in Figure 9 and Figure 14, which is as expected since when heat is, absorbed the buoyancy force decreases which retard the flow rate and thereby giving rise to decrease in the velocity profiles.

Figure 12 has been plotted to depict the variation of temperature profiles against for different values of heat source parameter (R) by fixing other physical parameters (Figures 13-17). From this graph, we observe that temperature decrease with increase in the Heat source parameter (R).

international-journal-chemical-sciences-chemical-reaction

FIG. 17. Effects of concentration profiles for different values of chemical reaction parameter (Kr).

Conclusion

In this study, an analytical solution of unsteady hydromagnetic free convective heat and mass transfer flow of viscous, in-compressible, electrically conducting viscoelastic fluid through porous medium bounded by porous plate in the presence of chemical reaction, radiation absorption and heat source has been investigated.

The conclusions of the study are follows:

(i) The Hartmann number has the effect of decreasing the flow field both primary and secondary velocity at all the points due to the magnetic pull of the Lorentz force acting on the flow field. So, magnetic field can effectively be used to control the flow.

(ii) The primary and secondary velocity along main flow increase with increasing thermal Grash of number.

(iii) There is an accelerating effect of Sc on the primary and the secondary velocity in presence of heavier diffusing species.

(iv) The chemical reaction parameter has the effect of decreasing the both primary and secondary velocity profiles.

(v) The Prandtl number has a retarding effect on the primary and secondary velocity of the flow field.

(vi) The primary velocity decreases with an increase in heat source parameter and radiation absorption parameter, whereas secondary velocity increases with an increase in heat source parameter and radiation absorption parameter.

(vii) The viscoelastic parameter has the influence of decreasing the primary and secondary velocity.

References