Mathematical Problems in Engineering

Volume 2015, Article ID 724596, 10 pages

http://dx.doi.org/10.1155/2015/724596

## Diffusion of Chemically Reactive Species in Casson Fluid Flow over an Unsteady Stretching Surface in Porous Medium in the Presence of a Magnetic Field

School of Mathematics, Statistics and Computer Science, University of KwaZulu-Natal, Private Bag X01, Scottsville, Pietermaritzburg 3209, South Africa

Received 31 July 2014; Accepted 14 September 2014

Academic Editor: Stanford Shateyi

Copyright © 2015 Gilbert Makanda et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A study is performed on two-dimensional flow and diffusion of chemically reactive species of Casson fluid from an unsteady stretching surface in porous medium in the presence of a magnetic field. The boundary layer velocity, temperature, and concentration profiles are numerically computed for different governing parameters. The paper intends to show unique results of a combination of heat transfer and chemical reaction in Casson fluid flow. The resulting partial differential equations are converted to a system of ordinary differential equations using the appropriate similarity transformation, which are solved by using the Runge-Kutta-Fehlberg numerical scheme. The results in this work are validated by the comparison with other authors.

#### 1. Introduction

The study of Casson fluid has attracted attention to many researchers due to its application in the field of metallurgy, food processing, drilling operations, and bioengineering operations. Its application extends to the manufacturing of pharmaceutical products, coal in water, china clay, paints, synthetic lubricants, and biological fluids such as synovial fluids, sewage sludge, jelly, tomato sauce, honey, soup, and blood due to its contents such as plasma, fibrinogen, and protein, making the study of Casson fluid important in fluid dynamics. Casson fluid is classified as a non-Newtonian fluid due to its rheological characteristics. These characteristics show shear stress-strain relationships that are significantly different from Newtonian fluids. The study of non-Newtonian fluids has not been thoroughly covered due to the complex representation of their constitutive equations. It is therefore important to undertake this study of Casson fluid. Most studies have concentrated on viscoelastic fluids in which different constitutive equations have been suggested. This work can be applied to chemical processing equipment in which some fluids react chemically with some ingredients present in them.

The driving force for mass transfer is a combination of temperature and concentration gradients. In this study the effect of chemical reaction on the fluid is considered as in Mukhopadhyay and Vajravelu [1]. The study of boundary layer flow over a stretching sheet has been studied by Mukhopadhyay et al. [2], among others, who investigated Casson fluid flow over an unsteady stretching surface. In their work they did not consider mass transfer and they considered a different wall temperature expression. Abd El-Aziz [3] studied mixed convection flow of a micropolar fluid from an unsteady stretching surface with viscous dissipation. In this work he considered a similar stretching velocity, wall temperature, and wall concentration distribution. We extended the work of Grubka and Bobba [4] who investigated heat transfer characteristics of a continuous stretching surface with variable temperature in which we introduced the MHD and porous medium source terms and chemical reaction effects. Sharidan et al. [5] studied similarity solutions for the unsteady boundary layer flow and heat transfer due to a stretching sheet; Nadeem et al. [6, 7] and Ahmed and Nazar [8] also studied Casson fluid over a stretching sheet and in their paper they assumed that the velocity of the stretching surface is linearly proportional to the distance from fixed origin.

In the study of non-Newtonian fluids many authors have studied the flow of blood as Casson fluid. The studies were carried out by Rohlf and Tenti [9], among others, who investigated the role of Womersley number in pulsatile blood flow, a theoretical study of the Casson model; Sankar and Lee [10, 11] investigated two-fluid nonlinear model for flow in catheterized blood vessels and two-fluid Casson model for pulsatile blood flow through stenosed arteries, respectively. Shaw et al. [12] studied pulsatile Casson fluid flow through stenosed bifurcated artery. In relation to blood flow there are other research works that were done in different geometries such as flows in microslit channels, slightly curved channels, and peristaltic transport as in [13–15].

The study of Casson fluid in porous media was also studied by Nadeem et al. [16] who considered MHD three-dimensional Casson fluid flow past a porous linearly stretching sheet. Dash et al. [17] studied Casson fluid flow in a pipe filled with homogeneous porous medium. Tripathi [18] investigated the transient peristaltic heat flow through a finite porous channel. Pramanik [19] studied Casson fluid flow and heat transfer past an exponentially porous stretching surface in the presence of thermal radiation. Ramachandra et al. [20] investigated flow and heat transfer of Casson fluid from a horizontal circular cylinder with partial slip in a non-Darcy porous medium. In their work they considered slip conditions on the wall.

There are many studies that investigated fluid flow with chemical reactions. Kameswaran et al. [21] investigated homogeneous-heterogeneous reactions in a nanofluid flow due to a porous stretching sheet, Shaw et al. [22] studied homogeneous-heterogeneous reactions in a nanofluid flow due to a porous stretching sheet, and Chamkha et al. [23] investigated similarity solutions for unsteady heat and mass transfer from a stretching surface embedded in a porous medium with suction/injection and chemical reaction effects.

Although there are many applications and use of non-Newtonian fluids in industry, the study of Casson fluid has not been thoroughly investigated for heat and mass transfer past a stretching surface. In this work we extend the work of Mukhopadhyay and Vajravelu [1] and Grubka and Bobba [4] in which the energy equation, the source terms for porous medium, and magnetic field are introduced. Similarity transformations are used to convert the partial differential equations into ordinary differential equations which are then solved by using Runge-Kutta-Fehlberg integration scheme and the successive linearization method described by Makanda et al. [24]. In this work we investigate the effect of varying unsteadiness parameter, Casson, Schmidt, and Prandtl numbers, and the reaction rate parameter on the velocity, temperature, and concentration profiles with the depiction of graphical illustrations.

#### 2. Mathematical Formulation

Consider two-dimensional laminar boundary layer flow, temperature, and mass transfer of an incompressible Casson fluid flow over an unsteady stretching sheet. The flow of heat and mass transfer starts at . The sheet is pulled out of the slit at the origin and moves with velocity , , are constants, and is the initial stretching rate. The rheological equation of state for an isotropic and incompressible flow of a Casson fluid is given as in [1, 2, 20] bywhere and is the th component of the deformation rate, is the product of the deformation rate with itself, is a critical value of this product based on the non-Newtonian model, is the plastic dynamic viscosity of the non-Newtonian fluid, and is the yield stress of the fluid. Given that and are, respectively, the temperature and concentration at the sheet and and are, respectively, the ambient conditions, the positive coordinate is measured along the stretching sheet and the positive coordinate is measured perpendicular to the sheet. It is assumed that both temperature and concentration at the surface vary with distance from the origin and time. The temperature and concentration at the surface are therefore given bywhere and are constants. The surface temperature and surface concentration increase if and are positive and reduce if they are negative from and at the origin to and the temperature and concentration increase/decrease along the sheet. It is assumed that radiation effects and viscous dissipation are negligible. The expressions , , and are only valid for but not when . Under these assumptions the governing equations in this flow are given aswhere is kinematic viscosity of Casson fluid, is the non-Newtonian Casson parameter, is the permeability of the porous medium, is the electrical conductivity, is the strength of the magnetic field, is the density of the Casson fluid, is the diffusion coefficient of species in the fluid, is the thermal diffusivity, and is the time dependent reaction rate, where represents destructive reaction, represents constructive reaction, and is a constant. The boundary conditions are given aswhere the subscript refers to the ambient condition.

We introduce the nondimensional variableswhere is the stream function which satisfies the continuity equation (3). The velocity components are defined asThe governing equations reduce towith boundary conditionswhere is the unsteadiness parameter, is the Prandtl number, is the Schmidt number, and is the reaction parameter. is the permeability parameter; is the magnetic parameter. The nondimensional temperature and concentration are, respectively, given by and

The parameters of engineering interests are the local skin friction and the Nusselt and Sherwood numbers which are defined asThe local Nusselt and Sherwood numbers are defined asUsing expressions (11) and (6), where is the Reynolds number defined as It is important at this stage to mention that the system of equations (8)-(9) reduce to those of Grubka and Bobba [4] when , . The present problem reduces to that of Grubka and Bobba [4], denote steady flow, and in their paper they obtained an exact solution in terms of Kummer’s functions written in terms of the confluent hypergeometric functions.

The solution of the boundary value problem (8)-(9) was solved using the Runge-Kutta-Fehlberg integration scheme. In the method we choose finite values of . This value is the boundary layer thickness given by . We begin by choosing an initial guess of to obtain the values , , and . The solution is repeated with new values until two consecutive values differ by .

#### 3. Results and Discussion

To obtain a clear understanding of the flow of Casson fluid, we discuss the physics of the problem by studying the effects of the unsteadiness (), permeability (), magnetic (), Prandtl (Pr), Schmidt (Sc), and reaction rate () numbers on velocity, temperature, and concentration profiles. We also study the variation of skin friction, the local Nusselt number, and the Sherwood number with unsteadiness parameter. For validation of the numerical method used in this study, results for the Nusselt number for the Newtonian fluid were compared to those of Abd El-Aziz [3] and Grubka and Bobba [4] for the unsteadiness parameter . The comparison is shown in Table 1 and it is found to be in agreement with at least four decimal places. To further verify the accuracy of the numerical scheme used, the successive linearization method (SLM) was used and there was agreement with Runge-Kutta-Fehlberg integration scheme.