A Mixed Convection Boundary Layer Flow over a Vertical Wall in a Porous Medium , with Exponentially Varying Fluid Viscosity

This study investigates a mixed convection boundary layer flow over a vertical wall embedded in a highly porous medium. The fluid viscosity is assumed to decrease exponentially with temperature. The boundary layer equations are transformed into a non-similar form using an appropriate nonsimilar variable ξ and a pseudo-similar variable η . The non-similar equations are solved using an efficient local non-similarity method. The effect of viscosity variation parameter on the heat transfer, skin friction and the velocity and temperature distribution within the boundary layer is investigated. The viscosity variation parameter, the viscous dissipation parameter and non-similarity variable are shown to have a significant effect on velocity and thermal boundary layer and also on the skin friction coefficient and heat transfer at the wall.


Introduction
Mixed convection boundary layer flow through a porous medium is of primary importance due to its applications in industry such as geothermal operations (for example dynamics of hot springs), terrestrial heat flow through an aquifer, flow of moisture through porous industrial material etc. Engineering applications include porous thermal insulation and heat exchangers with fluidized beds etc.
Reviews on convective heat transfer in porous medium are well documented in the books by Nield and Bejan, [1] and Pop and Ingham [2], Kaviany [3] and in the handbook of porous media edited by Vafai [4].Numerous studies on mixed convection boundary through a porous media using the Darcy model and extended Darcy models have been carried out, by others, by Karabis et al. [5] investigated a mixed forced-free convection boundary layer flow over a vertical wall in a highly porous medium, Chen [6] who looked at a boundary layer flow over a horizontal in a saturated porous medium.Using the model developed by Vafai and Tien [7], Hong et al. [8] investigated the non-Darcian effects on a flow over a vertical plate in porous media.
Most fluids used in engineering have temperature varying viscosity.A study by Hossain and Munir [9] investigated a mixed convection boundary layer with temperature dependent viscosity where it was shown that the viscosity variation parameter has an effect on the wall skin friction coefficient and wall heat transfer coefficient.Jayanthi and Kumari [10] studied the effects of variable viscosity on a non-Darcy flow of a non-Newtonian fluid over a vertical surface in porous medium.They showed that variable viscosity parameter has a significant effect on the fluid flow properties.
The current study focuses on a boundary layer flow over a vertical wall embedded in a highly porous medium and with variable viscosity.The non-similar boundary layer equations are solved using the local non-similarity method due to Sparrow et al. [11] and Sparrow and Yu [12].In Section 2 the mixed convection boundary layer flow equations are formulated and transformed to the non-similar form, in Section 3 the equations are solved and the results presented and discussed in Section 4. Conclusions are drawn in Section 5.

Mixed Convection Boundary Layer Flow in a Porous Medium
A steady, two-dimensional flow of a viscous incompressible fluid over a heated semi-infinite vertical wall embedded in a highly porous medium is studied.The x -axis is taken along the wall and the y -axis as the normal axis to x -axis.The free-stream velocity and temperature are taken as U ∞ and T ∞ respectively and the wall temperature as . The fluid viscosity µ is assumed to decrease exponentially with temperature.Under the boundary layer approximation for a flow through a porous medium and the Boussinesq's approximation, the mixed convection boundary layer equations through a porous media take the form: 3) The appropriate boundary conditions for the flow are: Here g is the gravitational acceleration, α is the coefficient of thermal expansion, κ is the thermal conductivity, p c is the specific heat at constant pressure, 0 ρ is the fluid density at reference temperature 0 T and K is the permeability of the porous medium.The wall temperature w T is assumed to be a constant.Most com- mon fluids have dynamic viscosity that varies exponentially with temperature.Here we use Arrhenius model which takes the exponential form ( ) where 0 µ is reference viscosity at reference temperature 0 T and γ is a viscosity variation number defined as 0 1 ln where 0 w µ µ > with be w µ being the fluid viscosity at the wall temperature w T .For gases 0 γ < and for liquids it is positive.For small values of γ Taylor expansion leads to the linear or the inverse relation of vis- cosity with temperature.

Non-Similar Boundary Layer Equations
The equations are rendered non-dimensional by introducing a boundary layer pseudo-similarity variable ( ) We use the scaling where ( ) is the dimensionless viscosity and Assuming power law variation in the free-stream velocity , we have n β = .The non-similar boun- dary layer equations take the form: where x S is the dimensionless local permeability parameter, Pr is the Prandtl number and Ec is the local Eckert number, defined respectively as ( ) The boundary conditions, in non-dimensional form, are: The shear stress and the heat transfer at the wall can be represented using the local skin friction coefficient and the local Nusselt number, defined by ,0 , ,0

Numerical Solutions of the Non-Similar Boundary Layer Equations
The system of Equations (2.4)-(2.6) is solved using the non-similar method of Sparrow et al. [11].This method has been used extensively in solving non-similar mixed convection boundary layer equations by among other Sparrow and Yu [12], Gorla et al. [13], Mureithi and Mason [14], Bég et al. [15] etc.The method involves approximating the solutions at different levels of truncation.We introduce the notation: The first level of approximation involves solving the local similarity equations, where all terms involving the ξ and its derivatives (that is the non-similar terms, 1 1 , f θ ) in the Equations (2.4)-(2.6)are neglected.
At the second level of truncation, the method involves solving the system (2.4)-(2.6)together with the auxiliary system obtained by differentiating (2.4)-(2.6)partially with respect ξ and neglecting all terms involving second order differentiation with respect to  (that is 2 2 , etc. f θ The auxiliary system takes the form: The system (2.4)-(2.9)solved as if it is an ordinary differential system, with ξ being treated as a parameter.
The results obtained are a better approximation than those obtained at the first level of approximation.An even better approximation can be obtained at the third order level of approximation, where the system (2.4)-(2.9)are solved together with an of auxiliary system obtained by further differentiating Equations (2.7)-(2.9)partially with respect to ξ and all terms containing third order partials with respect to ξ are neglected.

Discussion of Results
We approximate the solution of (2.4)-(2.6)using the second level of approximation.The system of Equations (2.4)-(2.9) is solved using the fourth order Runge-Kutta method combined with a Newton root refinement scheme.Computations were carried out for the case when Pr 0.7 = (air) and viscosity that increases exponentially with increase in temperature ( ) 0 ε < .Figure 1 show the effect varying the non-similar parameter ξ on the scaled stream-wise velocity, f ′ , and viscosity, ( ) µ θ , distribution in the boundary layer.At a fixed stream-wise location, ξ becomes a buoyancy- related parameter.As ξ increases velocity overshoots (exceeding free-stream velocity) are observed within the boundary layer.The maximum values of the super-velocities increase with increase in ξ Results in figure 1also show that increasing ξ has the effect of decreasing the thickness of the viscosity boundary layer.
The effect of varying the viscous dissipation parameter Ec was investigated and the results are presented in Figure 2. The results indicate that as Ec increases, the amplitude of the super-velocities increase.Temperature is also shown to increase with increase in Ec and for high enough values of Ec (in this case 0.6 Ec ≥ ) there is a small region adjacent to the wall, where temperature exceeds its wall value of one.It follows that viscosity function also has positive gradient in that small region.This is an indication that increasing Ec results in increased heat transfer from the wall to the fluid.
Figure 3 shows stream-wise velocity and viscosity distributions within the boundary layer for different values of ε .The results show that at a fixed value of ξ velocity overshoots increase with increase in ε .On the con- trary, viscosity decreases with increase in ε .Heat transfer at the wall (represented by the Nusselt number,

ε = −
cient and a decrease in the Nusselt number.The results for the heat transfer coefficient indicate that at a fixed value of 0 Ec > there is a critical value of ξ beyond which the heat transfer coefficient changes sign from positive (heat flows from fluid to the wall) to negative (heat flows from wall to the fluid in the boundary layer).

Conclusions
A non-Darcy mixed convection boundary layer flow has been investigated.The results have shown that the non-similarity parameter ξ , which at particular stream-wise location represents thermal buoyancy, has the ef- fect of accelerating the scaled free-stream velocity, resulting in super-velocities within the boundary layer.This effect is exacerbated by increasing viscous dissipation parameter and the viscosity variation parameter.The temperature dependent viscosity decreases with increase in buoyancy-related parameter and viscosity variation parameter but increases with increase in the viscous dissipation parameter.
The skin friction coefficient increases with increase in the viscosity variation parameter, the viscous dissipation parameter and the buoyancy-related parameter.For flows with viscous dissipation effects, there is a critical value of the buoyancy related parameter, above which the Nusselt number changes sign from positive (wall extracts heat from fluid) to negative (fluid extracts heat from wall).

x
Re being the local Grash of number and the local Reynolds number, defined as

Figure 5 .Figure 1 .Figure 2 .Figure 3 .Figure 4 .Figure 5 .
Figure5.The results show that the skin friction increases linearly with increase in ξ while the Nusselt number decrease with increase in Ec .On the other hand, increase Ec results in increase in the skin friction coeffi-