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.

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 -axis is taken along the wall and the -axis as the normal axis to -axis. The free-stream velocity and temperature are taken as and 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:

The appropriate boundary conditions for the flow are:

Here is the gravitational acceleration, is the coefficient of thermal expansion, is the thermal con- ductivity, is the specific heat at constant pressure, is the fluid density at reference temperature and is the permeability of the porous medium. The wall temperature 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 is reference viscosity at reference temperature and is a viscosity variation number defined as

where with be being the fluid viscosity at the wall temperature. For gases and for liquids it is positive. For small values of Taylor expansion leads to the linear or the inverse relation of viscosity with temperature.

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,) 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 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.

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 (air) and viscosity that increases exponentially with increase in temperature.

Figure 1 show the effect varying the non-similar parameter on the scaled stream-wise velocity, , 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 was investigated and the results are presented in Figure 2. The results indicate that as increases, the amplitude of the super-velocities increase. Temperature is also shown to increase with increase in and for high enough values of (in this case) 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 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, and shear stress at the wall (repre- sented by skin friction coefficient, were investigated as functions of the effect of and various values of Figure 4 shows that the skin friction coefficient and Nusselt number increase with increase in. The skin friction coefficient increases linearly with while the Nusselt number decreases with increase in.

Similar results showing the effect of on both the skin friction and the Nusselt number are presented in Figure 5. The results show that the skin friction increases linearly with increase in while the Nusselt number decrease with increase in. On the other hand, increase results in increase in the skin friction coeffi-

cient and a decrease in the Nusselt number. The results for the heat transfer coefficient indicate that at a fixed value of 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).

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 effect 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).

The author would like to thank the East African Universities Mathematics Programme (EAUMP) for their financial support.