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

Authors: Ahmed M. Salem; Mohamed Abd El-Aziz
Journal: Mathematical Problems in Engineering
 (2013)

---

## 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.


---

## Body

## 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. Thex-axis is directed along the plate and points to the direction of motion. The y-axis is perpendicular to the plate. The fluid is permeated by an imposed uniform magnetic field of strength B0, which acts in the positive y-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
(1)∂u∂x+∂v∂y=0,(2)∂∂y(μun)+ρbk*∂∂y(u2)=k*(ρgβ∂T∂y-σB02∂u∂y),(3)u∂T∂x+v∂T∂y=α∂2T∂y2+μk*ρcp(un-1+ρbk*μu)u2.
The associated boundary conditions are
(4)v=0,T=Tw+N∂T∂yaty=0,u⟶0,T⟶T∞asy⟶∞.
Here x,y are the Cartesian coordinates, u and v the velocity components, n the power-law index, µ the viscosity, T the temperature, α the effective thermal diffusivity of the saturated porous medium, β the volumetric thermal expansion, σ the electric conductivity, b an empirical constant associated with the Forchheimer inertia term, ρ the density, cp the specific heat at constant pressure, g the acceleration due to gravity, k* is the modified permeability for the flow of non-Newtonian power law fluid, and N is the thermal slip factor. Christopher and Middleman [13] and Dharmadhikari and Kale [14] proposed the following relationships for the permeability k* as a function of the power-law index n:
(5)k*=(625)(nε3n+1)n(εd3(1-ε))n+1
(see Christopher and Middleman [13]),
(6)k*=6ε(6n+110n+3)n(7516)3(10n-3)/(10n+11)(ε2d8(1-ε))n+1
(see Dharmadhikari and Kale [14]).In the above equation,ε is the porosity and d 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])
(7)μ=μ∞(11+γ(T-T∞))
or
(8)μ=1a(T-T∞),
where
(9)a=γμ∞,Tr=T∞-1γ.
In the above relations, both a and Tr 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, a>0 for liquids and a<0 for gases.To proceed, we define the dimensionless variables(10)η=yxRax,ψ=αRaxf(η),θ(η)=T-TrTw-T∞+θr,
where Rax=(x/α)(k*ρgβ(Tw-T∞)/μ∞)1/n is the local rayleigh number, ψ is the physical stream function defined so that u=∂ψ/∂y, v=-∂ψ/∂x, and the continuity equation is satisfied automatically, θr=(Tr-T∞)/(Tw-T∞)is the viscosity/temperature parameter, f is a dimensionless stream function, η is a similarity space variable, and θ is the dimensionless temperature. The velocity components are readily obtained as
(11)u=αxRaxf′,v=-α2xRax(f-ηf′).
The mathematical problem defined in (2) and (3) is then transformed into the coupled nonlinear ordinary differential equations and their associated boundary conditions:
(12)θr(θr-θ)2(f′)nθ′+nθr(θr-θ)(f′)n-1f′′+2Gr*f′f′′=θ′-Mf′′,(13)θ′′+Ecf′[θθr-θ(f′)n+Gr*(f′′)2]=-12fθ′,(14)f(0)=0,θ(0)=1+Λθ′(0),(15)f′(∞)=0,θ(∞)=0,
where prime denotes ordinary differentiation with respect to η, Rax is the local rayleigh number, M=(σB02k*/μ∞Raxn-1)(x/α)n-1 is the magnetic field parameter, Gr*=f0Rad2-n is the non-Darcian parameter withf0=(α/d)2-n(ρbk*/μ∞) representing the structural and thermophysical properties of the porous medium and Rad the pore diameter-dependent Rayleigh number, which describes the relative intensity of the buoyancy force and the dispersion parameter, Ec=xgβ/cp is the Eckert number, and Λ=N0(a/α)1/2 is the temperature jump parameter.The physical quantity of interest is the local Nusselt numberNux which is defined by
(16)Nux=-xTw-T∞(∂T∂y)y=0=-Rax1/2θ′(0).


## 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 f′(0) and θ′(0). The solution process is repeated with another large value of η∞ until two successive values of f′(0) and θ′(0) 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 y1=f, y2=f′, y3=θ, and y4=θ′; (12) and (13) then reduced into a system of ordinary differential equations; that is,
(17)y1′=y2,y1(0)=0,y2′=y4[(θr-y3)2-θry2n][nθr(θr-y3)y2n-1+(θr-y3)2(2Gr*y2+M)],y2(0)=δ1,y3′=y4,y3(0)=1+Λδ2,y4′=-12y1y4-Ec(θr(θr-y3)y2n+1+Gry23),y4(0)=δ2,
where δ1and δ2 are determined such that they satisfy y2(∞)=0 and y4(∞)=0. The shooting method is used to guess δ1and δ2 until the boundary conditions y2(∞)=0 and y4(∞)=0 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 (Gr*=0) and non-Darcy (Gr*=2) cases and different values of power-law index n (0.5≤n≤1.5), magnetic field parameter M (0≤M≤2), viscosity/temperature parameter θr (-∞≤θr≤-1), viscous dissipation parameter (the Eckert number) Ec (0≤Ec≤0.8), and temperature jump parameter Λ (0≤Λ≤0.8).Table1 shows a comparison of the present results for -θ′(0) with those reported by Ibrahim et al. [6] for various values of n and Gr* in the absence of magnetic parameter (M=0), temperature jump (Λ=0), viscous dissipation (Ec=0), and variable viscosity (θr→-∞). It can be seen from this table that the comparison is excellent especially for higher values of the viscosity index (n=1) and smaller values of the non-Darcian parameter (Gr*=0).Table 1
Comparison of the values of-θ′(0)  for  Λ=Ec=M=0,θr→-∞ and various values of n and Gr* with previously published data.
-
θ
′
(
0
)
n
Gr
*
[6]
Present results
0.5
0.0
0.40613
0.399325
0.5
1.0
0.32263
0.319624
1.0
0.0
0.44455
0.444552
1.0
1.0
0.38064
0.379509Typical velocity distributions are presented at selected values ofM and Gr* for a shear thinning fluid of n=0.5 in Figure 1(a) and for a shear thickening fluid of n=1.5 in Figure 1(b) when θr=-1, Λ=0.2, and Ec=0.05. It is seen from these figures that the magnetic field decreases the velocity f′ in both Darcy (Gr*=0) and non-Darcy (Gr*=2) 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 (Gr*=0), the velocity decreases dramatically in the absence of the magnetic field M, but this behavior is more prominent for a shear thinning fluid.(a) Velocity profiles in the Darcy and non-Darcy cases for various values ofM with n=0.5, θr=-1.0, Λ=0.2, and Ec=0.05. (b) Velocity profiles in the Darcy and non-Darcy cases for various values of M with n=1.5, θr=-1.0, Λ=0.2, and Ec=0.05.
(a)
(b)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 n=0.5 and in Figure 2(b) for a shear thickening fluid of n=1.5. 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 M is increased. Further, Figure 2(a) shows that the effect of M 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 M, Figure 2(b) shows that the effect of non-Darcian parameter Gr* 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).(a) Temperature profiles in the Darcy and non-Darcy cases for various values ofM with n=0.5, θr=-1.0, Λ=0.2, and Ec=0.05. (b) Temperature profiles in the Darcy and non-Darcy cases for various values of M with n=1.5, θr=-1.0, Λ=0.2, and Ec=0.05.
(a)
(b)Variations of the dimensionless velocity profilesf′(η) for various values of the Eckert number Ec and non-Darcian parameter Gr* are shown in Figure 3(a) for n=0.5 (shear thinning fluid) and in Figure 3(b) for n=1.5 (shear thickening fluid), while the corresponding temperature θ(η) profiles are presented in Figure 4(a) for n=0.5 and in Figure 4(b) for n=1.5. It is clear from these figures that as compared to the case of no viscous dissipation (Ec=0), it can be seen that the dimensionless velocity and temperature are increased with increasing the viscous dissipation Ec. 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 Ec(Ec=0.8). 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 (Gr*=0) and non-Darcy (Gr*=2) 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 Ec, as expected. Further, viscous dissipation demonstrates a more pronounced influence on the temperature distribution at a larger Gr* for both shear thinning and shear thickening fluids. On the other hand, at a given value of Ec, the impact of the non-Darcian parameter Gr* is to increase the momentum and thermal boundary layer thickness for both shear thinning and shear thickening fluids.(a) Velocity profiles in the Darcy and non-Darcy cases for various values ofEc with n=0.5, θr=-1.0, M=0.2, and Λ=0.2. (b) Velocity profiles in the Darcy and non-Darcy cases for various values of Ec with n=1.5, θr=-1.0, M=0.2, and Λ=0.2.
(a)
(b)(a) Temperature profiles in the Darcy and non-Darcy cases for various values ofEc with n=0.5, θr=-1.0, M=0.2, and Λ=0.2. (b) Temperature profiles in the Darcy and non-Darcy cases for various values of Ec with n=1.5, θr=-1.0, M=0.2, and Λ=0.2.
(a)
(b)The effects of viscosity/temperature parameterθr on the velocity f′(η) and temperature θ for both Darcy (Gr*=0) and non-Darcy (Gr*=2) flow cases are illustrated in Figures 5(a) and 6(a) for a shear thinning (n=0.5) fluid and in Figures 5(b) and 6(b) for a shear thickening (n=1.5) fluid. The results for the case of constant viscosity (i.e., θr→-∞) 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 θr. It is observed from Figures 5(a) and 5(b) that the velocity f′(η) increases at a distance η0 from the plate, where the value of η0 depends on the values of Gr* and n, while the exact opposite trend is noticed for η>η0 as θr increases, and this conduct is veritable for all values of Gr*and n. For example, η0≃1.55 for Gr*=0 and n=0.5, while the corresponding value of η0≃1.8 for Gr*=0 and n=1.5. For a specified value of Gr*, on the contrary, to increase θr 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 θr and Gr* on the velocity and temperature profiles are found to be more intelligible for the case of shear thinning fluid.(a) Velocity profiles in the Darcy and non-Darcy cases for various values ofθr with n=0.5, M=0.2, Λ=0.2, and Ec=0.05. (b). Velocity profiles in the Darcy and non-Darcy cases for various values of θr with n=1.5, M=0.2, Λ=0.2, and Ec=0.05.
(a)
(b)(a) Temperature profiles in the Darcy and non-Darcy cases for various values ofθr with n=0.5, M=0.2, Λ=0.2, and Ec=0.05. (b) Temperature profiles in the Darcy and non-Darcy cases for various values of θr with n=1.5, M=0.2, Λ=0.2, and Ec=0.05.
(a)
(b)Figures7(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 Gr* 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.(a) Velocity profiles in the Darcy and non-Darcy cases for various values ofΛ with n=0.5, θr=-1.0, M=0.2, and Ec=0.05. (b) Velocity profiles in the Darcy and non-Darcy cases for various values of Λ with n=1.5, θr=-1.0, M=0.2, and Ec=0.05.
(a)
(b)(a) Temperature profiles in the Darcy and non-Darcy cases for various values ofΛ with n=0.5, θr=-1.0, M=0.2, and Ec=0.05. (b) Temperature profiles in the Darcy and non-Darcy cases for various values of Λ with n=1.5, θr=-1.0, M=0.2, and Ec=0.05.
(a)
(b)Representative numerical results for the fluid velocity at the platef′(0), the local Nusselt number NuxRax-1/2, and the temperature of the plate θ(0) for various parametric values are presented in Table 2 for a shear thinning fluid with n=0.5 and in Table 3 for a shear thickening fluid with n=1.5. It can be seen from Table 2 that the fluid velocity at the plate f′(0) and the local Nusselt number NuxRax-1/2 increase while the temperature of the plate θ(0) decreases as the viscosity/temperature parameter θr increases. In contrast with θr effect, increases in the magnetic parameter M result in decreases in the fluid velocity at the plate f′(0) and the local Nusselt number NuxRax-1/2, and increases in the temperature of the plate θ(0). Also, from Table 2 we observe that increases in the values of temperature jump parameter Λ cause increases in the values of f′(0), NuxRax-1/2 and θ(0). It can also be observed that as the Eckert number Ec increases, both the fluid velocity and the temperature at the plate f′(0) and θ(0), respectively, increase while the local Nusselt number NuxRax-1/2 decreases. Further, Table 2 shows that the behaviors of f′(0), NuxRax-1/2, and θ(0) due to the effects of θr, Λ, M, and Ec are the same for both Darcy and non-Darcy flow cases, but these behaviors are more pronounced in the case of Darcy flow.Table 2
Variation off′(0),NuxRax-1/2, and θ(0)  with θr,Λ,M, and Ec at n=0.5.
θ
r
Λ
M
Ec
f
′
(
0
)
Nu
x
Ra
x
-
1
/
2
θ
(
0
)
Gr
*
=
0
Gr
*
=
2
Gr
*
=
0
Gr
*
=
2
Gr
*
=
0
Gr
*
=
2
−1.0
0.2
0.2
0.05
1.466160
0.494216
0.391782
0.270713
0.921644
0.945857
−2.0
0.2
0.2
0.05
1.094730
0.450206
0.343294
0.256483
0.931341
0.948703
−∞
0.2
0.2
0.05
0.662177
0.369577
0.280718
0.233396
0.943856
0.953321
−1.0
0.0
0.2
0.05
1.727240
0.518504
0.463863
0.294485
1.000000
1.000000
−1.0
0.5
0.2
0.05
1.209750
0.462126
0.322324
0.241737
0.838838
0.879131
−1.0
0.8
0.2
0.05
1.039080
0.435250
0.276641
0.218961
0.778687
0.824831
−1.0
0.2
0.0
0.05
2.960420
0.533910
0.487462
0.279969
0.902508
0.944006
−1.0
0.2
1.0
0.05
0.561688
0.369883
0.27535
0.238824
0.944930
0.952235
−1.0
0.2
2.0
0.05
0.331220
0.272547
0.224427
0.210505
0.955115
0.957899
−1.0
0.2
2.0
0.0
1.453310
0.492309
0.411001
0.287924
0.917800
0.9424152
−1.0
0.2
2.0
0.5
1.616710
0.513857
0.161406
0.047102
0.967719
0.9905796
−1.0
0.2
2.0
0.8
1.78933
0.546404
−0.094211
−0.318840
1.018842
1.0637680Table 3
Variation off′(0),NuxRax-1/2, and  ;θ(0) with θr,Λ,M, and Ec at n=1.5.
θ
r
Λ
M
Ec
f
′
(
0
)
Nu
x
Ra
x
-
1
/
2
θ
(
0
)
Gr
*
=
0
Gr
*
=
2
Gr
*
=
0
Gr
*
=
2
Gr
*
=
0
Gr
*
=
2
−1.0
0.2
0.2
0.05
1.189271
0.553978
0.424613
0.312988
0.915077
0.937402
−2.0
0.2
0.2
0.05
1.027950
0.532614
0.402245
0.307785
0.919551
0.938443
−∞
0.2
0.2
0.05
0.831943
0.495849
0.373951
0.299693
0.925210
0.940061
−1.0
0.0
0.2
0.05
1.298940
0.576869
0.483828
0.340723
1.000000
1.000000
−1.0
0.5
0.2
0.05
1.067201
0.524817
0.361203
0.27957
0.819399
0.860215
−1.0
0.8
0.2
0.05
0.976295
0.500298
0.315835
0.2531
0.747332
0.797520
−1.0
0.2
0.0
0.05
1.443571
0.591791
0.456514
0.321858
0.908697
0.935628
−1.0
0.2
1.0
0.05
0.657822
0.429802
0.331357
0.279032
0.933729
0.944194
−1.0
0.2
2.0
0.05
0.406623
0.323619
0.267707
0.244694
0.946459
0.951061
−1.0
0.2
2.0
0.0
1.183923
0.552565
0.445510
0.332016
0.910898
0.933597
−1.0
0.2
2.0
0.5
1.255322
0.571986
0.168316
0.0675017
0.966336
0.986500
−1.0
0.2
2.0
0.8
1.341190
0.60079
−0.162245
−0.336549
1.032449
1.067310Values off′(0), NuxRax-1/2, and θ(0) in the case of a shear thickening fluid case of n=1.5 for the same physical parameters presented in Table 2 are listed in Table 3. Effects of all physical parameters on the f′(0), NuxRax-1/2, and θ(0) 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 f′(0) of the shear thinning fluid is higher than that of the shear thickening fluid for all values of θr, Λ, M, and Ec while the opposite result was observed for the non-Darcy case (Gr*=2). For the Darcy flow regime, Tables 2 and 3 show that the local Nusselt number NuxRax-1/2 of the shear thinning fluid is lower than that of the shear thickening fluid for all values of θr, Λ, and Ec while the opposite trend is true for the non-Darcy flow regime. On the contrary, the magnetic field M has an opposite effect to the effects of θr, Λ, and Ec on the NuxRax-1/2 for both Darcy and non-Darcy flow cases. On the other hand, it is observed from these tables that the temperature of the plate θ(0) for a shear thinning fluid is higher than that for a shear thickening fluid for all values of θr, Λ, and M, 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 f′(0) and the local Nusselt number NuxRax-1/2 decrease while the plate temperature θ(0) 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.



---

*Source: 761023-2013-04-24.xml*