Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 761023, 10 pages

http://dx.doi.org/10.1155/2013/761023

## Hydromagnetic Non-Darcian Free-Convective Flow of a Non-Newtonian Fluid with Temperature Jump

^{1}Mathematics Department, Faculty of Science and Arts, King Khalid University, Belqarn, P.O. Box 61985, Saudi Arabia^{2}Mathematics Department, Faculty of Science, Helwan University, P.O. Box 11795, Cairo, Egypt^{3}Mathematics Department, Faculty of Science, King Khalid University, Abha, P.O. Box 9004, Saudi Arabia

Received 16 February 2013; Accepted 30 March 2013

Academic Editor: Mohamed Seddeek

Copyright © 2013 Ahmed M. Salem and Mohamed Abd El-Aziz. 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

In the present study, the effect of viscous dissipation on magnetohydrodynamic (MHD) non-Darcian free-convection flow of a non-Newtonian power-law fluid past a vertical flat plate in a saturated porous medium with variable viscosity and temperature jump is considered. The fluid is permeated by a transverse magnetic field imposed perpendicularly to the plate on the assumption of a small magnetic Reynolds number. The fluid viscosity is assumed to vary as a reciprocal of linear function of temperature. The governing boundary layer equations and boundary conditions are cast into a dimensionless form and simplified by using a similarity transformation into a system of nonlinear ordinary differential equations and solved numerically. The effects of the governing parameters on the flow fields and heat transfer are shown in graphs and tabular form.

#### 1. Introduction

The analysis of convective heat transfer over a vertical plate in a non-Newtonian fluid saturated porous medium is important for the thermal design of industrial equipment dealing with molten plastics, polymeric liquids, oil extraction, ground water pollution, thermal insulation, solid matrix heat exchangers, the storage of nuclear waste, packed-bed catalytic reactors, foodstuffs and slurries, and so forth. Thus, considerable efforts have been directed toward major aspects of this coupled, nonlinear boundary-layer problem. Most of the earlier studies were based on Darcy’s law which states that the volume-averaged velocity is proportional to the pressure gradient. The Darcy model is shown to be valid under the conditions of low velocities and small porosity, Cheng and Hsu [1]. In many practical situations, the porous medium has higher flow rates and reveals nonuniform porosity distribution in the near wall region, making the Darcy’s law inapplicable. To model the real physical situation better, it is therefore necessary to include the aforementioned non-Darcian terms in the analysis of convective transport in a porous medium. The inertia effect is expected to be important at a higher flow rate, and it can be accounted for through the addition of a velocity-squared term in the momentum equation, which is known as the Forchheimer extension. The boundary effect may become significant when heat transfer is considered in a region very close to a solid boundary. The Brinkman’s extension, which incorporates a viscous shear stress term into the momentum equation, together with the no-slip boundary condition, is usually used to shed light on the importance of boundary effects. A comment was made by Fand and Brucker [2] that the effect of viscous dissipation might be significant in the case of natural convection in porous medium in connection with their experimental correlation for heat transfer in external flows. The validity of the comment was tested for the Darcy model by Fand et al. [3], both experimentally and analytically while estimating the heat transfer coefficient from a horizontal cylinder embedded in a saturated porous medium. Their mathematical analysis is confined to studying the dissipation effect using a steady, 1D energy equation; the basis of the equation is from the analogy given by Bejan [4] for the inclusion of viscous dissipation effects. The influence of viscous dissipation can be seen from the analogy given by Tucker and Dessenberger [5] to model the heat transfer and fluid flow through porous media in order to study the resin transfer molding (R.T.M.) for producing fiber-reinforced polymeric parts in final shape. Ibrahim et al. [6] studied the problem of natural convective heat transfer for a non-Newtonian fluid from an impermeable vertical plate embedded in a fluid-saturated porous medium. El-Amin [7] investigated the influences of viscous dissipation on buoyancy-induced flow over (a horizontal or a vertical flat plate) embedded in a non-Newtonian fluid-saturated porous medium. Mansour and El-Shaer [8] studied the effect of magnetic field on non-Darcy axisymmetric free convection in a power-law fluid-saturated porous medium with variable permeability. Mahdy [9] studied the flow, heat, and mass transfer characteristics of non-Darcian mixed convection of a non-Newtonian fluid from a vertical isothermal plate embedded in a homogeneous porous medium with the effect of Soret and Dufour. It is known that the fluid properties change with temperature especially the fluid viscosity. So in order to predict accurately the flow behavior, it is necessary to take into account this variation of viscosity since the recent results on the flow of non-Newtonian fluid have shown that when this effect is included, the flow characteristics may be substantially changed compared to the constant viscosity case. The effect of variable viscosity on non-Darcy free or mixed convective heat transfer along a vertical surface embedded in a porous medium saturated with a non-Newtonian fluid is analysed by Jayanthi and Kumari [10]. Salem [11] studied the problem of flow and heat transfer of an electrically conducting viscoelastic fluid over a continuously stretching sheet in the presence of a uniform magnetic field.

The present work is devoted to study the problem of free-convective heat transfer from a vertical flat plate embedded in a non-Darcian porous medium saturated with a power law non-Newtonian fluid. The fluid is subjected to a transverse uniform magnetic field. It is assumed that the magnetic Reynolds number is small so that the induced magnetic field is neglected. Viscous dissipation due to mechanical work is included in the temperature field. The fluid viscosity is assumed to vary as a reciprocal of linear function of temperature. Temperature jump is also considered which gives interesting features regarding such flow.

#### 2. Analysis

Let us consider the flow of steady, laminar, incompressible, two-dimensional, non-Darcy, non-Newtonian power-law fluid flow with heat transfer along a vertical flat plate in a saturated porous medium. It is assumed that the fluid and the solid matrix are everywhere in local thermal equilibrium. The -axis is directed along the plate and points to the direction of motion. The -axis is perpendicular to the plate. The fluid is permeated by an imposed uniform magnetic field of strength , which acts in the positive -axis. The magnetic Reynolds number is taken to be small enough so that the induced magnetic field can be neglected. With these assumptions and the application of the Boussinesq and boundary layer approximations, the governing equations, namely, the equation of continuity, the non-Darcy flow model (i.e., the model given by Shenoy [12]), and the energy equation, are The associated boundary conditions are Here are the Cartesian coordinates, and the velocity components, the power-law index, µ the viscosity, the temperature, the effective thermal diffusivity of the saturated porous medium, the volumetric thermal expansion, the electric conductivity, an empirical constant associated with the Forchheimer inertia term, the density, the specific heat at constant pressure, the acceleration due to gravity, is the modified permeability for the flow of non-Newtonian power law fluid, and is the thermal slip factor. Christopher and Middleman [13] and Dharmadhikari and Kale [14] proposed the following relationships for the permeability as a function of the power-law index : (see Christopher and Middleman [13]), (see Dharmadhikari and Kale [14]).

In the above equation, is the porosity and is the pore diameter.

The fluid properties are assumed to be isotropic and constant, except for the fluid viscosity which is assumed to be reciprocal of a linear function of temperature and can be written as (Kumari [15], Lai and Kulacki [16]) or where In the above relations, both and are constants, and their values depend on the reference state and the thermal property of the fluid . The viscosity of a liquid usually decreases with increasing temperature and it increases with gases in general, for liquids and for gases.

To proceed, we define the dimensionless variables where is the local rayleigh number, is the physical stream function defined so that , , and the continuity equation is satisfied automatically, is the viscosity/temperature parameter, is a dimensionless stream function, is a similarity space variable, and is the dimensionless temperature. The velocity components are readily obtained as The mathematical problem defined in (2) and (3) is then transformed into the coupled nonlinear ordinary differential equations and their associated boundary conditions: where prime denotes ordinary differentiation with respect to , is the local rayleigh number, is the magnetic field parameter, is the non-Darcian parameter with representing the structural and thermophysical properties of the porous medium and the pore diameter-dependent Rayleigh number, which describes the relative intensity of the buoyancy force and the dispersion parameter, is the Eckert number, and is the temperature jump parameter.

The physical quantity of interest is the local Nusselt number which is defined by

#### 3. Numerical Solution

The coupled non-linear differential equations (12) and (13) subject to the boundary conditions (14) and (15) constitute a two-point boundary value problem. In order to solve these equations numerically, we follow the most efficient numerical shooting technique with fifth-order Runge-Kutta-Fehlberg integration scheme. In this method, it is most important to choose the appropriate finite values of . To select , we begin with some initial guess value and solve the problem with some particular set of parameters to obtain and . The solution process is repeated with another large value of until two successive values of and differ only after desired digit signifying the limit of the boundary along . The last value of is chosen as appropriate value of the limit for that particular set of parameters. Equations (12) and (13) were first formulated as a set of four first-order simultaneous equations of four unknowns following the method of superposition [17]. Thus, we set , , , and ; (12) and (13) then reduced into a system of ordinary differential equations; that is,
where and are determined such that they satisfy and . The shooting method is used to guess and until the boundary conditions and are satisfied. Then the resulting differential equations can be integrated by fifth-order Runge-Kutta-Fehlberg integration scheme. The above procedure is repeated until we get the results up to the desired degree of accuracy, 10^{−6}.

#### 4. Results and Discussion

The system of transformed governing equations (12) and (13) subject to the boundary conditions (14) and (15) is solved numerically using the shooting method. Solutions are obtained for both Darcy () and non-Darcy () cases and different values of power-law index (), magnetic field parameter (), viscosity/temperature parameter (), viscous dissipation parameter (the Eckert number) (), and temperature jump parameter ().

Table 1 shows a comparison of the present results for with those reported by Ibrahim et al. [6] for various values of and in the absence of magnetic parameter (), temperature jump (), viscous dissipation (), and variable viscosity (). It can be seen from this table that the comparison is excellent especially for higher values of the viscosity index () and smaller values of the non-Darcian parameter ().

Typical velocity distributions are presented at selected values of and for a shear thinning fluid of in Figure 1(a) and for a shear thickening fluid of in Figure 1(b) when , , and . It is seen from these figures that the magnetic field decreases the velocity in both Darcy () and non-Darcy () flow regimes, and this behavior is more noticeable for the Darcy flow case. Also, the momentum boundary layer thickness for shear thickening fluid is larger than that for shear-thinning fluid. Moreover, for the case of Darcy flow (), the velocity decreases dramatically in the absence of the magnetic field , but this behavior is more prominent for a shear thinning fluid.

Representative temperature profiles within the boundary layer for the same values of parameters considered in Figures 1(a) and 1(b) are presented in Figure 2(a) for a shear thinning fluid of and in Figure 2(b) for a shear thickening fluid of . For both Darcy and non-Darcy flow regimes, it is clear from these figures that the temperature increases with an increase in the thermal boundary layer thickness as the magnetic field is increased. Further, Figure 2(a) shows that the effect of in the non-Darcy flow case is seen to be more pronounced for the shear thinning fluid, while the opposite result is observed when the fluid is shear thickening fluid as observed from Figure 2(b). For given , Figure 2(b) shows that the effect of non-Darcian parameter is to increase the temperature for the case of a shear thickening fluid, while the opposite behavior is obtained when the fluid is shear thinning as shown in Figure 2(a).

Variations of the dimensionless velocity profiles for various values of the Eckert number and non-Darcian parameter are shown in Figure 3(a) for (shear thinning fluid) and in Figure 3(b) for (shear thickening fluid), while the corresponding temperature profiles are presented in Figure 4(a) for and in Figure 4(b) for . It is clear from these figures that as compared to the case of no viscous dissipation (), it can be seen that the dimensionless velocity and temperature are increased with increasing the viscous dissipation . This phenomenon is evident from the fact that viscous dissipation produces heat due to drag between the fluid particles, and this extra heat causes an increase in the fluid initial temperature, which in turn gives rise to the velocity of the fluid particles, and hence the fluid velocity increases. It is also observed that a distinctive peak in the temperature profile occurs in the fluid adjacent to the plate for higher values of . The presence of the peak implies that the temperature attains maximum value inside the body of the fluid adjacent to the surface but not at the surface, and therefore heat is expected to transfer from the ambient fluid to the surface. This behavior is true in the Darcy () and non-Darcy () flow cases. Also, the increase in the fluid velocity and temperature due to viscous heating is observed to be more pronounced for a higher value of , as expected. Further, viscous dissipation demonstrates a more pronounced influence on the temperature distribution at a larger for both shear thinning and shear thickening fluids. On the other hand, at a given value of , the impact of the non-Darcian parameter is to increase the momentum and thermal boundary layer thickness for both shear thinning and shear thickening fluids.

The effects of viscosity/temperature parameter on the velocity and temperature for both Darcy () and non-Darcy () flow cases are illustrated in Figures 5(a) and 6(a) for a shear thinning () fluid and in Figures 5(b) and 6(b) for a shear thickening () fluid. The results for the case of constant viscosity (i.e., are also shown in these figures. It is seen from these figures that the velocity and temperature distributions are significantly influenced by the value of . It is observed from Figures 5(a) and 5(b) that the velocity increases at a distance from the plate, where the value of depends on the values of and , while the exact opposite trend is noticed for as increases, and this conduct is veritable for all values of and . For example, for and , while the corresponding value of for and . For a specified value of , on the contrary, to increase is to reduce the temperature for both shear thinning and shear thickening fluids as indicated in Figures 6(a) and 6(b). Comparison of Figures 5 and 6 shows that the effects of and on the velocity and temperature profiles are found to be more intelligible for the case of shear thinning fluid.

Figures 7(a) and 8(a) demonstrate, respectively, the effects of temperature jump parameter on the dimensionless velocity and temperature profiles in the boundary layer for both Darcy and non-Darcy flow regimes for a shear thinning fluid, whereas Figures 7(b) and 8(b) are for a shear thickening fluid. For Darcy and non-Darcy flow cases, the velocity decreases rapidly near the boundary by increasing the temperature jump parameter as observed from Figures 7(a) and 7(b). On the other hand, to increase the magnitude of temperature jump parameter is to reduce the dimensionless temperature profile, and this behavior is conducted for both values of as obvious from Figures 7(b) and 8(b). The temperature jump parameter demonstrates a more pronounced influence on the velocity and temperature distribution of a shear thinning fluid.

Representative numerical results for the fluid velocity at the plate , the local Nusselt number , and the temperature of the plate for various parametric values are presented in Table 2 for a shear thinning fluid with and in Table 3 for a shear thickening fluid with . It can be seen from Table 2 that the fluid velocity at the plate and the local Nusselt number increase while the temperature of the plate decreases as the viscosity/temperature parameter increases. In contrast with effect, increases in the magnetic parameter result in decreases in the fluid velocity at the plate and the local Nusselt number , and increases in the temperature of the plate . Also, from Table 2 we observe that increases in the values of temperature jump parameter cause increases in the values of , and . It can also be observed that as the Eckert number increases, both the fluid velocity and the temperature at the plate and , respectively, increase while the local Nusselt number decreases. Further, Table 2 shows that the behaviors of , , and due to the effects of , , , and are the same for both Darcy and non-Darcy flow cases, but these behaviors are more pronounced in the case of Darcy flow.

Values of , , and in the case of a shear thickening fluid case of for the same physical parameters presented in Table 2 are listed in Table 3. Effects of all physical parameters on the , , and are noticed to be qualitatively similar but with quantitatively increased or decreased magnitude as compared to the case of a shear thinning fluid. Comparison study of Tables 2 and 3 indicates for the Darcy case that of the shear thinning fluid is higher than that of the shear thickening fluid for all values of , , , and while the opposite result was observed for the non-Darcy case (). For the Darcy flow regime, Tables 2 and 3 show that the local Nusselt number of the shear thinning fluid is lower than that of the shear thickening fluid for all values of , , and while the opposite trend is true for the non-Darcy flow regime. On the contrary, the magnetic field has an opposite effect to the effects of , , and on the for both Darcy and non-Darcy flow cases. On the other hand, it is observed from these tables that the temperature of the plate for a shear thinning fluid is higher than that for a shear thickening fluid for all values of , , and , and this is true for both Darcy and non-Darcy flow regimes. Finally, keeping all other parameters fixed, it is shown from Tables 2 and 3 that due to the transition from Darcy regime to non-Darcy regime, the fluid velocity at the plate and the local Nusselt number decrease while the plate temperature increases.

#### 5. Conclusion

The problem of MHD free-convective flow and heat transfer from a vertical flat plate embedded in a non-Darcian porous medium saturated with a power law non-Newtonian fluid is investigated with viscous dissipation variable viscosity. By using similarity transformations, the governing boundary layer equations, for momentum and thermal energy are reduced to coupled non-linear ordinary differential equations which are then solved numerically using a shooting method. The results pertaining to the present study indicate that the flow and the temperature field are significantly influenced by considering the non-Darcy case for the pseudoplastic and dilatant fluids. Increasing the thermal slip parameter causes decreasing the velocity and temperature for Darcy and non-Darcy flow cases. Also increasing the thermal slip parameter causes reducing the heat transfer from the sheet.

#### Acknowledgments

This research is supported by King Khalid University (KKU), Saudi Arabia, under the Project no. KKU-SCI-11-12. This support is highly appreciated and acknowledged. Also the authors are very thankful to the referees for their encouraging comments and constructive suggestions to improve the presentation of this paper.

#### References

- P. Cheng and C. T. Hsu, “Higher-order approximations for Darcian free convective flow about
a semi-infinite vertical flat plate,”
*Journal of Heat Transfer*, vol. 106, no. 1, pp. 143–151, 1984. View at Google Scholar · View at Scopus - R. M. Fand and J. Brucker, “A correlation for heat transfer by natural convection from horizontal cylinders that accounts for viscous dissipation,”
*International Journal of Heat and Mass Transfer*, vol. 26, no. 5, pp. 709–716, 1983. View at Google Scholar · View at Scopus - R. M. Fand, T. E. Steinberger, and P. Cheng, “Natural convection heat transfer from a horizontal cylinder embedded in a porous medium,”
*International Journal of Heat and Mass Transfer*, vol. 29, no. 1, pp. 119–133, 1986. View at Google Scholar · View at Scopus - A. Bejan,
*Convection Heat Transfer*, Wiley, New York, NY, USA, 1984. - C. L. Tucker III and R. B. Dessenberger, “Governing equations for flow and heat transfer in stationary fiber beds,” in
*Flow and Rheology in Polymer Composites Manufacturing*, S. G. Advani, Ed., Elsevier Science, Amsterdam, The Netherlands, 1994. View at Google Scholar - F. S. Ibrahim, M. A. Mansour, and S. M. Abdel-Gaied, “Radiative and thermal dispersion effects on non-Darcy natural convection with lateral mass flux for non-Newtonian fluid from a vertical flat plate in a saturated porous medium,”
*Transport in Porous Media*, vol. 61, no. 1, pp. 45–57, 2005. View at Publisher · View at Google Scholar · View at Scopus - M. F. El-Amin, “Combined effect of magnetic field and viscous dissipation on a power-law fluid over plate with variable surface heat flux embedded in a porous medium,”
*Journal of Magnetism and Magnetic Materials*, vol. 261, no. 1-2, pp. 228–237, 2003. View at Publisher · View at Google Scholar · View at Scopus - M. A. Mansour and N. A. El-Shaer, “Effect of magnetic field on non-Darcy axisymmetric free convection in a power-law fluid saturated porous medium with variable permeability,”
*Journal of Magnetism and Magnetic Materials*, vol. 250, pp. 57–64, 2002. View at Publisher · View at Google Scholar · View at Scopus - A. Mahdy, “Soret and Dufour effect on double diffusion mixed convection from a vertical surface in a porous medium saturated with a non-Newtonian fluid,”
*Journal of Non-Newtonian Fluid Mechanics*, vol. 165, no. 11-12, pp. 568–575, 2010. View at Publisher · View at Google Scholar · View at Scopus - S. Jayanthi and M. Kumari, “Effect of variable viscosity on non-Darcy free or mixed convection flow on a vertical surface in a non-Newtonian fluid saturated porous medium,”
*Applied Mathematics and Computation*, vol. 186, no. 2, pp. 1643–1659, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. M. Salem, “Variable viscosity and thermal conductivity effects on MHD flow and heat transfer in viscoelastic fluid over a stretching sheet,”
*Physics Letters A*, vol. 369, no. 4, pp. 315–322, 2007. View at Publisher · View at Google Scholar · View at Scopus - A. V. Shenoy, “Darcy-Forchheimer natural, forced and mixed convection heat transfer in non-Newtonian power-law fluid-saturated porous media,”
*Transport in Porous Media*, vol. 11, no. 3, pp. 219–241, 1993. View at Publisher · View at Google Scholar · View at Scopus - R. H. Christopher and S. Middleman, “Power-law flow through a packed tube,”
*Industrial and Engineering Chemistry Fundamentals*, vol. 4, no. 4, pp. 422–426, 1965. View at Google Scholar · View at Scopus - R. V. Dharmadhikari and D. D. Kale, “Flow of non-Newtonian fluids through porous media,”
*Chemical Engineering Science*, vol. 40, no. 3, pp. 527–528, 1985. View at Google Scholar · View at Scopus - M. Kumari, “Effect of variable viscosity on non-Darcy free or mixed convection flow on a horizontal surface in a saturated porous medium,”
*International Communications in Heat and Mass Transfer*, vol. 28, no. 5, pp. 723–732, 2001. View at Publisher · View at Google Scholar · View at Scopus - F. C. Lai and F. A. Kulacki, “The effect of variable viscosity on convective heat transfer along a vertical surface in a saturated porous medium,”
*International Journal of Heat and Mass Transfer*, vol. 33, no. 5, pp. 1028–1031, 1990. View at Google Scholar · View at Scopus - T. Y. Na,
*Computational Methods in Engineering Boundary Value Problems*, vol. 145 of*Mathematics in Science and Engineering*, Academic Press, New York, NY, USA, 1979. View at MathSciNet