Effects of Conjugate Heat Transfer on Steady MHD Mixed Convective Heat Transfer Flow over a Thin Vertical Plate Embedded in a Porous Medium with High Porosity

Abstract

This study investigates mixed convection heat transfer about a thin vertical plate in the presence of magneto and conjugate heat transfer effects in the porous medium with high porosity. The fluid is assumed to be incompressible and dense. The nonlinear coupled parabolic partial differential equations governing the flow are transformed into the nonsimilar boundary layer equations, which are then solved numerically using the Keller box method. The effects of the conjugate heat transfer parameter p, the porous medium parameter k₁, the Forchheimer parameter F*, the mixed convection parameter Ri, the magnetic parameter Mn, and the electric field parameter E₁ on the velocity and temperature profiles as well as on the local skin friction and local heat transfer are presented and analyzed. The validity of the methodology and analysis is checked by comparing the results obtained for some specific cases with those available in the literature.

1. Introduction

Thermal buoyancy-induced flow and convective heat transfer in fluid-saturated porous media have been the subject of numerous publications. This interest in the subject stems from various engineering applications in geothermal reservoirs, petroleum industries, transpiration cooling, storage of radioactive nuclear waste materials, separation processes in chemical industries, building thermal insulation, and solar heating systems [1]. Kaviany [2], Pop and Ingham [3], Ingham and Pop [4, 5], and Nield and Bejan [6] have made comprehensive reviews of the studies of heat transfer in relation to the above applications.

The conjugate heat transfer problem, in which the coupled heat transfer processes between the solid body (conduction mechanisms) and the fluid flow (convection mechanisms) are considered simultaneously, has been investigated by several researchers for with and without porous medium. Miyamoto et al. [7] studied two-dimensional conjugate heat transfer problems of free convection from a vertical flat plate with a uniform temperature or a uniform heat flux at the outsidem surface of the plate. Sparrow and Chyu [8] investigated the conjugate problem for a vertical plate fin with various heat transfer coefficients under forced convection. They assumed that the heat conduction in the fin was to be one-dimensional. Char et al. [9] employed the cubic spline collocation numerical method to analyze the conjugate heat transfer in the laminar boundary layer on a continuous, moving plate. Wang [10] studied the thermo-fluid-dynamic field resulting from the coupling of wall conduction with laminar mixed convection heat transfer of micropolar fluids along a vertical flat plate. Pop et al. [11] presented a detailed numerical study of the conjugate mixed convection flow along a vertical flat plate. Chang [12] presented a numerical analysis of the flow and heat transfer characteristics of mixed convection in a micropolar fluid flowing along a vertical flat plate with conduction effects. He assumed that the heat conduction in the wall was only in the transversal direction. Mamun et al. [13] investigated the heat generation effect on natural convection flow along and conduction inside a vertical flat plate. Hsiao and Hsu [14] studied a conjugate mixed convection heat transfer problem of a second-grade viscoelastic fluid past a horizontal flat-plate fin.

For porous medium, Pop and Na [15] reported a numerical study of the steady conjugate free convection over a vertical slender, hollow circular cylinder with the inner surface at a constant temperature and embedded in a porous medium. Vaszi et al. [16] investigated two-dimensional conjugate free convection from an inclined flat plate in a semi-infinite porous medium under the boundary-layer approximation. Méndez et al. [17] studied steady state heat transfer characteristics of a thin vertical strip with internal heat generation placed in a porous medium. Saeid [18] studied steady conjugate natural convection in two-dimensional vertical porous layer sandwiched between two equal-thickness walls. The Darcy model was used in the mathematical formulation for the porous layer and finite volume method was used to solve the dimensionless governing equations. Vaszi et al. [19] investigated two-dimensional conjugate free convection in a porous medium from a vertical plate fin. Hossain et al. [20] investigated conjugate free convection along a vertical cylindrical fin in a non-Newtonian fluid-saturated porous medium. The boundary layer equations based on the power law model appropriate for the Darcy flows were solved numerically.

This article illustrates the effect of wall conduction on non-Darcy MHD mixed convection flow over a thin vertical plate embedded in a porous medium with high porosity. The boundary layer equations governing the flow are reduced to local nonsimilarity equations which are solved using the implicit finite difference method (Keller box). Variation in the fluid-solid interfacial temperature distribution, the local skin friction and local heat transfer parameters as well as the velocity and temperature profiles are presented to highlight the influence of the wall conduction, porous medium parameter, buoyancy, electric field, and magnetic effects parameters.

2. Analysis

This study considers two-dimensional, steady, laminar, non-Darcy, incompressible electrically conducting fluid flow over a thin vertical plate of length 𝐿 and finite thickness 𝑏 (𝐿≫𝑏) embedded in a porous medium. The physical model and coordinate system are shown in Figure 1. Far above/below the surfaces of the thin vertical plate, the velocity and the temperature of the free stream are 𝑢_∞ and 𝑇_∞, respectively. The temperature of the inside surface of the plate is maintained at a constant temperature of 𝑇₀, where 𝑇₀ > 𝑇_∞. The flow region is exposed under uniform transverse magnetic fields ^→𝐵=(0,𝐵₀,0) and uniform electric field ^→𝐸=(0,0,−𝐸₀) [21]. Magnetic and electric fields are known from Maxwell’s equation that ∇⋅^→𝐵=0 and ∇×^→𝐸=0. When magnetic field is not so strong then electric field and magnetic field obey Ohm’s law ^→𝐽=𝜎(^→𝐸+^→𝑞×^→𝐵), where ^→𝐽 is the Joule current, 𝜎 is the magnetic permeability and ^→𝑞 is the fluid velocity. It is assumed that magnetic Reynolds number of the fluid is small so that induced magnetic field and Hall effect may be neglected and takes into account of magnetic field effect as well as electric field in momentum boundary layer equation [21]. The porous medium is assumed to be transparent and in thermal equilibrium with the fluid. Both the fluid and the porous medium are opaque for self-emitted thermal radiation. The properties of the fluid and the porous media, such as viscosity, thermal conductivity, specific heat, and permeability, are assumed to be constant. The porous medium is considered to be homogeneous and isotropic (i.e., uniform with a constant porosity and permeability).

Figure 1 

The schematic of the problem.

The governing equations for this investigation are based on the usual boundary-layer equations modified to include the porous medium effects, the magnetic and electric field effects, and the thermal buoyancy effects. These equations (with the Boussinesq and non-Darcy approximations) can be written as (see Chamkha et al. [1] and Chamkha [22]) 𝜕𝑢+𝜕𝑥𝜕𝜐1𝜕𝑦=0,𝜀²𝑢𝜕𝑢𝜕𝑥+𝜐𝜕𝑢=𝑣𝜕𝑦𝜀𝜕²𝑢𝜕𝑦²+𝜎𝜌𝐸₀𝐵₀−𝐵²₀𝑢−𝑢_∞−𝑣𝐾𝑢−𝑢_∞−𝐹𝑢²−𝑢²_∞+𝑔𝛽𝑇−𝑇_∞𝑢𝜕𝑇𝜕𝑥+𝜐𝜕𝑇=𝑘𝜕𝑦𝜌𝑐^𝑒_𝑝𝜕²𝑇𝜕𝑦². (2.1)(2.2)(2.3)

The above equations are called Brinkman-Forchheimer-extended-Darcy equations [23]. Here 𝑢 and 𝑣 are the velocity components in the 𝑥 and 𝑦 directions, respectively; 𝑣 is the kinematic viscosity; 𝑔 is the acceleration due to gravity; 𝑘_𝑒 is the effective thermal conductivity of porous medium; 𝜌 is the density of the fluid; 𝛽 is the coefficient of thermal expansion; 𝑇 is the temperature of the fluid inside the thermal boundary layer; 𝐾 is the permeability of the porous medium; 𝐹 is the empirical constant (Forchheimer number) in the second-order resistance and setting 𝐹=0 in (2.2), the equation is then reduced to the Darcy’s law. The third and fourth terms on the right hand side of (2.2) stand for the first-order (Darcy) resistance and second-order porous inertia resistance, respectively [21]. The thermal dispersion effect is minimal when the thermal diffusivity (𝑘_𝑒/𝜌𝑐_𝑝) of the porous matrix is of the same order of magnitude as that of the working fluid [1].

The appropriate boundary conditions for the velocity and temperature of this problem are 𝑥=0,𝑦>0,𝑇=𝑇_∞,𝑢=𝑢_∞,𝑥>0,𝑦=0,𝑇=𝑇_𝑤(𝑥),𝑢=0,𝜐=0,𝑦⟶∞,𝑇⟶𝑇_∞,𝑢⟶𝑢_∞, (2.4) where subscripts 𝑤 and ∞ refer to the wall and the boundary layer edge, respectively. In addition, 𝑇_𝑤(𝑥) is the surface temperature of the plate, which is not known a priori.

One objective of the current study is to predict the surface temperature of the plate 𝑇_𝑤(𝑥). Therefore, an additional governing equation is required for the plate based on the simplification that the plate steadily transfers its heat to the surrounding fluid. Since the thickness of the plate, 𝑏, is small compared with its length, 𝐿, the axial conduction term in the heat conduction equation of the thin plate can be omitted [24, 25]. The governing equation for the temperature distribution within the plate is given by 𝑑²𝑇𝑑𝑦²_𝑠=0;0≤𝑥≤𝐿;−𝑏<𝑦≤0. (2.5)

The boundary conditions for the wall of plate are given At𝑦=−𝑏,𝑇_𝑠=𝑇₀,Attheinterface𝑦=0,𝑇_𝑠=𝑇(𝑥,0)∶−𝑘_𝑠𝑑𝑇𝑑𝑦_𝑦=0,𝑠=−𝑘_𝑓𝜕𝑇(𝑥,0)𝜕𝑦_𝑦=0,𝑓, (2.6a)(2.6b)where 𝑘_𝑠 and 𝑘_𝑓 are the thermal conductivity of the plate and the fluid, respectively. The boundary conditions given in (2.6b) state the physical requirements that the temperature and heat flux of the plate must be continuous across the solid-fluid interface. From (2.5) and (2.6a), (2.6b), the temperature distribution 𝑇_𝑤 at the interface is shown to be 𝑇_𝑤(𝑥)()𝑘=𝑇𝑥,0=𝑏^𝑓𝑘_𝑠𝜕𝑇(𝑥,0)𝜕𝑦+𝑇₀. (2.7)

To seek a solution, the following dimensionless variables are introduced: 𝜉(𝑥)=𝑥𝐿,𝜓𝑥,𝑦=𝑣𝑢_∞𝑥^1/2𝑓,𝑢𝜉,𝜂𝜂=𝑦^∞𝑣𝑥^1/2,𝜃=𝑇−𝑇_∞𝑇₀−𝑇_∞, (2.8)

where 𝜓(𝑥,𝑦) is the free stream function that satisfies (2.1) with 𝑢=𝜕𝜓/𝜕𝑦 and 𝜐=−𝜕𝜓/𝜕𝑥.

In terms of these new variables, the velocity components can be expressed as 𝑢=𝑢_∞𝑓′,𝜐=−𝑣𝑢_∞𝑥12𝑥12𝑓+𝜉𝜕𝑓−𝜂𝜕𝜉2𝑓^′. (2.9)(2.10)

The transformed momentum and energy equations together with the boundary conditions, (2.2), (2.3), and (2.5), can be written as 1𝜀𝑓^′′′+12𝜀²𝑓𝑓^′′+Mn𝜉𝐸₁−𝑓^′−1−𝑘₁𝜉𝑓^′−1−𝐹^∗𝜉𝑓^′2𝜉−1+Ri𝜉𝜃=𝜀²𝑓^′𝜕𝑓^′𝜕𝜉−𝑓^′′𝜕𝑓1𝜕𝜉𝜃Pr^′′+12𝑓𝜃^′=𝜉𝑓^′𝜕𝜃𝜕𝜉−𝜃^′𝜕𝑓,𝜕𝜉 (2.11) with the boundary conditions: 𝑓(𝜉,0)+2𝜉𝜕𝑓𝜕𝜉=0,𝑓^′(𝜉,0)=0,𝜃(𝜉,0)−1=𝑝𝜉^−1/2𝜃^′(),𝑓𝜉,0_′(𝜉,∞)=1,𝜃(𝜉,∞)=0, (2.12) where 𝑝=(𝑘_𝑓/𝑘_𝑠)(𝑏/𝐿)Re_𝐿^1/2 is the conjugate heat transfer parameter. It should be noticed that for the limiting case of 𝑝=0, the thermal boundary condition in (2.12) on the wall becomes isothermal. Hence, the magnitude of p determines the importance of the wall heat conduction effect [25].

The corresponding dimensionless groups that appeared in the governing equations defined as: Mn=HaRe,Ha=𝜎𝐵²₀𝐿²𝜇,𝐸₁=𝐸,𝐸Re𝐸=⁰𝐵₀𝑣/𝐿,𝑘₁=𝑣𝐿𝐾𝑢_∞,𝐹^∗=𝐹𝐿,Ri=Gr_𝐿Re_𝐿,Gr_𝐿=𝑔𝛽𝑇_𝑤−𝑇_∞𝐾𝐿𝑣²,Re_𝐿=𝑢_∞𝐿𝑣,Pr=𝜇𝑐_𝑝𝑘=𝑣𝛼, (2.13) where Mn is the magnetic parameter; Ha is the Hartman number; 𝐸₁ is the electric field parameter; 𝑘₁ is the porous medium parameter; 𝐹^∗ is the Forchheimer parameter expressing the relative importance of the inertia effect. Ri is the Richardson number, which measures the relative importance of free to forced convection. Ri = 0 corresponds to the case of purely forced convection condition. Ri → ∞ corresponds to the case of purely free convection condition. It is noted that Ri is not the function of 𝑥. Gr_𝐿 is the average Grashof number; Re_𝐿 is the average Reynolds number; Pr is the Prandtl number.

From the definition of the dimensionless wall temperature, it can be shown that 𝜃_𝑤=𝑇_𝑤−𝑇_∞𝑇₀−𝑇_∞. (2.14)

3. Numerical Solution

The system of transformed equations under the boundary conditions, Equations (2.11), (2.12), has been solved numerically using the Keller box scheme, which is proved to be an efficient and accurate finite-difference scheme [30]. Readers are referred to Cebeci and Bradshaw [30] for the details of the numerical methods. This is a very popular implicit scheme, which demonstrates the ability to solve systems of differential equations of any order as well as featuring second-order accuracy (which can be realized with arbitrary non-uniform spacing), allowing very rapid 𝑥 or 𝜉 variations [31].

A set of nonlinear finite-difference algebraic equations derived are then solved by using the Newton quazi-linearization method. The same methodology as that followed by Hossain et al. [20] is followed. Therefore, for the finite-difference forms of the equations, the reader is referred to Takhar and Beg [31] for the brevity of the paper.

In the calculations, a uniform grid of the step size 0.01 in the η-direction and a nonuniform grid in the ξ-direction with a starting step size 0.001 and an increase of 0.1 times the previous step size were found to be satisfactory in obtaining sufficient accuracy within a tolerance better than 10⁻⁶ in nearly all cases. The value of η_∞ = 16 is shown to satisfy the velocity to reach the relevant stream velocity.

In order to verify the accuracy of the present method, the present results were compared with those of Kuznetsov and Nield [26] (Table 1), Chamkha et al. [27], and Aydin and Kaya [28] (Table 2) and Aydin and Kaya [29] (Table 3). The comparison is found to be in good agreement, as shown in Tables 1, 2, and 3, respectively.

Table 1 

Comparison of the values −θ′(0, 0) for various values Pr at Mn=0.0, 𝐸₁=0.0, 𝜀=1.0, 𝐹^∗=0.0, Ri=0.0, 𝑘₁=1.0, and 𝑝=0.0.

Table 2 

Comparison of the values −θ′(ξ, 0) for various values 𝜉 with Pr=0.7, Mn=1.0, Ri=1.0, 𝜀=1.0, 𝐹^∗=0.0, 𝑘₁=0.0, and 𝑝=0.0.

Table 3 

Comparison of the values −θ′(ξ, 0) for various values 𝜉 and 𝜀 with Pr=1,𝐹^∗=0.5,Mn=0.0,𝐸₁=0.0,𝑘₁=1.0,Ri=0.0, and 𝑝=0.0.

4. Results and Discussion

The aim of this study was to investigate the flow and heat transfer characteristics for the non-Darcy mhd mixed convection flow over a thin vertical plate with wall conduction effect. The following ranges of the main parameters are considered: conjugate heat transfer parameter 𝑝=0.0, 0.1, 0.2, and 0.35; mixed convection parameter Ri = 0.0, 1.0, 5.0, and 10.0; magnetic interaction parameter Mn = 0.0, 0.5, 1.0, and 2.0; electric field parameter 𝐸₁ = 0.1, 0.5, 1.0, and 2.0; porous medium parameter 𝑘₁ = 0.1, 0.5, 1.0 and 2; inertia parameter 𝐹^∗ = 0.1, 0.5, 1.0, and 2.0; Pr = 1.0 and ε = 0.98. Hence, the numerical computations were performed.

The combined effects of 𝑝, Ri, Mn, 𝐸₁, 𝑘₁, and 𝐹^∗ on the momentum and heat transfer are analyzed. The Richardson number, Ri represents a measure of the effect of the buoyancy in comparison with that of the inertia of the external forced or free stream flow on the heat and fluid flow. Outside the mixed convection region, either the pure forced convection or the free convection analysis can be used to describe accurately the flow or the temperature field. Forced convection is the dominant mode of transport when Ri→0, whereas free convection is the dominant mode when Ri→∞. Buoyancy forces can enhance the surface heat transfer rate when they assist the forced convection [32].

The effect of conjugate heat transfer parameter p on the velocity (a) and temperature (b) profiles within the boundary layer with Ri = 1.0, Mn = 0.5, 𝐸₁ = 0.1, 𝑘₁ = 0.5, 𝐹^∗ = 0.5, and Pr = 1.0 is shown in Figures 2(a) and 2(b), respectively. The increasing of the conjugate heat transfer parameter decreases velocity and temperature gradients at the wall. A lower wall conductance 𝑘_𝑠 or higher convective cooling effect due to greater 𝑘_𝑓 increases the value of 𝑝 as well as causes greater temperature difference between the two surfaces of the plate. The temperature at the solid-fluid interface is reduced since the temperature at the outside surface of the plate is kept constant.

(a)

(b)

Figure 2 

Dimensionless velocity (a) and temperature (b) profiles for different 𝑝 while 𝜀=0.98, Pr=1.0, Ri=1, Mn=1, 𝐸1=0.1, 𝑘1=0.5, 𝐹∗=0.5 and 𝜉=0.5.

The variation of the interfacial temperature, the local skin friction, and the local heat transfer parameters for different values of 𝑝 with ξ are shown in Figures 3(a), 3(b), and 3(c), respectively where Ri = 1.0, Mn = 0.5, 𝐸₁ = 0.1, 𝑘₁ = 0.5, 𝐹^∗ = 0.5, and Pr = 1.0. It can be seen that the temperature of the fluid on the wall increases with 𝜉 for a given value of 𝑝 (Figure 3(a)). Comparing with isothermal wall (𝑝=0), an increase in the conjugate heat transfer parameter, 𝑝, causes a reduction in the interfacial temperature. This is because an increased value of 𝑝 corresponds to a lower wall conductance 𝑘_𝑠 and promotes a greater surface temperature variations (Figure 3(a)). The increases value of 𝑝 increases the momentum and thermal boundary layer thickness and therefore the local skin friction and the local heat transfer parameters decrease as shown in Figures 3(b) and 3(c). Also, increasing ξ increases the interfacial temperature, the local skin friction, and the local heat transfer parameters (Figures 3(a), 3(b), and 3(c)). Similar results were found in the literture [33, 34].

(a)

(b)

(c)

Figure 3 

Effect of conjugate heat transfer parameter 𝑝 on the dimensionless interfacial temperature (a), local skin friction (b), and local heat transfer (c) parameters against the streamwise distance 𝜉 while 𝜀=0.98, Pr=1.0, Ri=1, Mn=1, 𝐸1=0.1, 𝑘1=0.5, and 𝐹∗=0.5.

Figure 4 shows the dimensionless velocity (a) and temperature (b) profiles inside the boundary layer for different values of the buoyancy parameter Ri for the cases of isothermal plate (𝑝=0) and nonisothermal plate (𝑝>0). The increasing of Ri increases velocity and temperature gradients at the wall. And also, increasing the conjugate heat transfer parameter decreases the velocity and temperature profiles in the boundary layer.

(a)

(b)

Figure 4 

Dimensionless velocity (a) and temperature (b) profiles for different Ri while ε =0.98, Pr=1.0, Mn=1, 𝐸1=0.1, 𝑘1=0.5, 𝐹∗=0.5, and 𝜉=0.5.

Figure 5 shows the variation of the dimensionless interfacial temperature distributions (a), the local skin friction (b), and the local heat transfer (c) parameters in the boundary layer. It can be seen that as the values of ξ increases, the interfacial temperature rises. Compared with the limiting case of Ri = 0.0, an increase in the value of Ri gives rise to a reduced interfacial temperature since a greater value of Ri indicates a greater buoyancy effect, which increases the convection cooling effect and hence reduces the wall temperature (Figure 5(a)). Figure 5(b) illustrates the effect of the buoyancy force on the local skin friction factor for 𝑝=0 (solid lines) and 0.2 (dashed lines). It is observed that the local skin friction parameter increases with the buoyancy effect. The reason for this is that an increase in the buoyancy effect in mixed convection flow leads to an acceleration of the fluid flow, which increases the local skin friction parameter. Additionally, the higher the value of the buoyancy effect, the more the sensitivity of the wall conduction effects influences the skin friction factor. In Figure 5(c), the effect of the buoyancy force on the local heat transfer parameter is illustrated. It is noted that as the value of the Ri increases, the local heat transfer parameter also increases, both for the case of an isothermal plate (𝑝=0) and a nonisothermal plate (𝑝>0). This is because an increased buoyancy effect generates a greater buoyancy force, which increases the fluid velocity and hence raises the local heat transfer parameter [33].

(a)

(b)

(c)

Figure 5 

Effect of Ri on the dimensionless interfacial temperature (a), local skin friction (b), and local heat transfer (c) parameters against the streamwise distance 𝜉 while 𝜀=0.98, Pr=1.0, Mn=1, 𝐸1=0.1, 𝑘1=0.5, and 𝐹∗=0.5.

Figure 6 shows the dimensionless velocity (a) and temperature (b) profiles inside the boundary layer for different values of the magnetic parameter Mn for the cases of isothermal plate (𝑝=0) and nonisothermal plate (𝑝>0.0). The increasing of the magnetic parameter Mn increases velocity and temperature gradients at the wall due to magnetic field effects on external flow field. As mentioned above, increasing the conjugate heat transfer parameter 𝑝 decreases velocity and temperature gradients at the wall.

(a)

(b)

Figure 6 

Dimensionless velocity (a) and temperature (b) profiles for different Mn while 𝜀=0.98, Pr=1.0, Ri=5, 𝐸1=0.1, 𝑘1=0.5, 𝐹∗=0.5 and 𝜉=0.5.

The variation of the dimensionless interfacial temperature distributions, the local skin friction, and the local heat transfer parameters in the boundary layer for different magnetic parameter are shown in Figures 7(a), 7(b), and 7(c), respectively. Increasing the magnetic parameter Mn decreases the interfacial temperature (Figure 7(a)). The local skin friction and the local heat transfer parameters with different values of Mn for isothermal plate (𝑝=0) and nonisothermal plate (𝑝=0.2) are illustrated in Figures 7(b) and 7(c). The magnetic force aiding the flow and increases the local skin friction (i.e., shear stress) and the local heat transfer (i.e., heat transfer rate) parameters at the wall.

(a)

(b)

(c)

Figure 7 

Effect of Mn on the dimensionless interfacial temperature (a), local skin friction (b), and local heat transfer (c) parameters against the streamwise distance 𝜉 while 𝜀=0.98, Pr=1.0, Ri=5, 𝐸1=0.1, 𝑘1=0.5, and 𝐹∗=0.5.

The dimensionless velocity (a) and temperature (b) profiles inside the boundary layer for different values of the electric field parameter 𝐸₁ for the cases of isothermal plate (𝑝=0) and nonisothermal plate (𝑝>0.0) are illustrated in Figure 8. Increasing of the electric field parameter 𝐸₁ decreases the momentum and thermal boundary layers thickness (Figures 8(a) and 8(b)) because Lorentz force arising due to electric field acts as an accelerating force in reducing the frictional resistance [21]. Therefore, increasing the velocity and temperature gradients at wall increases the local skin friction and local heat transfer parameters as shown in Figures 9(b) and 9(c). Also increasing the electric field parameter 𝐸₁ decreases the interfacial temperature (Figure 9(a)).

(a)

(b)

Figure 8 

Dimensionless velocity (a) and temperature (b) profiles for different 𝐸1 while 𝜀=0.98, Pr=1.0, Ri=5, Mn=0.5, 𝑘1=0.5, 𝐹∗=0.5, and 𝜉=0.5.

(a)

(b)

(c)

Figure 9 

Effect of 𝐸1 on the dimensionless interfacial temperature (a), local skin friction (b), and local heat transfer (c) parameters against the streamwise distance 𝜉 while 𝜀=0.98, Pr=1.0, Ri=5, Mn=0.5, 𝑘1=0.5, and 𝐹∗=0.5.

The increasing of the porous medium parameter 𝑘₁ decreases momentum and thermal boundary layer thickness and increases velocity and temperature profiles for both isothermal (𝑝=0.0) and nonisothermal (𝑝=0.2) plate as shown in Figure 10. Increasing the porous medium parameter decreases the interfacial temperature (Figure 11(a)) and increases local skin friction (Figure 11(b)) and local heat transfer (Figure 11(c)) parameters. The increasing the porous medium parameter 𝑘₁ increases the temperature gradients at wall for 0.1<𝜉<0.55 (as expected from (2.9) and as observed in Aydin and Kaya [29] and Kumari and Nath [35]) whereas decreases the temperature gradient at wall for 0.55<𝜉<1.0.

(a)

(b)

Figure 10 

Dimensionless velocity (a) and temperature (b) profiles for different 𝑘1 while 𝜀=0.98, Pr=1.0, Ri=5, Mn=0.5, 𝐸1=0.1, 𝐹∗=0.5 and 𝜉=0.5.

(a)

(b)

(c)

Figure 11 

Effect of 𝑘1 on the dimensionless interfacial temperature (a), local skin friction (b), and local heat transfer (c) parameters against the streamwise distance 𝜉 while 𝜀=0.98, Pr=1.0, Ri=5, Mn=0.5, 𝐸1=0.1, and 𝐹∗=0.5.

Figure 12 shows the effect of drag (inertia) coefficient of porous medium (or Forchheimer parameter) 𝐹^∗ on the dimensionless velocity and temperature profiles at 𝜉=0.4. From this figure it is observed that the effect of drag coefficient 𝐹^∗ is to decrease the velocity profile, the temperature profile (Figure 12), and the interfacial temperature (Figure 13(a)) and increases local skin friction (Figure 13(b)) and local heat transfer (Figure 13(c)) parameters.

(a)

(b)

Figure 12 

Dimensionless velocity (a) and temperature (b) profiles for different 𝐹∗ while 𝜀=0.98, Pr=1.0, Ri=5, Mn=0.5, 𝐸1=0.1, 𝑘1=0.5, and 𝜉=0.4.

(a)

(b)

(c)

Figure 13 

Effect of 𝐹∗ on the dimensionless interfacial temperature (a), local skin friction (b) and local heat transfer (c) parameters against the streamwise distance 𝜉 while 𝜀=0.98, Pr=1.0, Ri=5, Mn=0.5, 𝐸1=0.1, 𝑘1=0.5.

5. Conclusions

This study has analyzed the effect of wall conduction on the non-Darcy MHD mixed convection flow with high porosity. The nonlinear formulation governing equations and their associated boundary conditions have been obtained and solved using the nonsimilarity transform and the finite difference method (Keller box), respectively. The influences of the conjugate heat transfer parameter, mixed convection parameter, magnetic interaction, electric field parameter, porous medium parameter and inertia parameter on the solid-liquid interfacial temperature distribution, the local skin friction, and the local heat transfer parameters have been systematically examined. From the present numerical investigation, the following conclusions can be drawn:(1)An increase in the conjugate heat transfer parameter decreases the velocity and the temperature gradient and therefore decreases the dimensionless interfacial temperature distribution, the local skin friction, and the local heat transfer parameters.(2)An increase in the magnetic, electric field and buoyancy parameters increases the local skin friction and local heat transfer parameters and decreases the dimensionless interfacial temperature distributions. Increasing Mn, 𝐸₁, and Ri decreases the velocity and temperature gradients at wall for nonisothermal cylinder (i.e., 𝑝>0).(3)An increase in the porous medium parameter 𝑘₁ and Forchheimer parameter 𝐹^∗ increases the local skin friction and local heat transfer parameters and decreases the dimensionless interfacial temperature distributions.

Nomenclature

𝐵0:	Magnetic field strength
𝑐𝑝:	Specific heat of the convective fluid
𝐸1:	Electric field parameter
𝑓:	Dimensionless stream function
𝐹∗:	Forchheimer parameter
Gr:	Grashof number
Ha:	Hartman number
𝑘𝑒:	Effective thermal conductivity of porous medium
𝐾:	Porous medium permeability, m2
𝑘1:	Porous medium parameter
𝐿:	Plate length
Mn:	Magnetic parameter, 𝑀=Ha/Re
Pr:	Prandtl number
Re:	Reynolds number
𝑇:	Temperature
𝑢,𝜐:	Velocities in 𝑥 and 𝑦 directions, respectively
𝑥,𝑦:	Coordinates in horizontal and vertical directions, respectively.

Greek Symbols

𝜎:	Magnetic permeability
𝛽:	Coefficient of thermal expansion
𝜂:	Pseudosimilarity variable,  𝑦Re𝑥1/2/𝑥
𝜉:	Nonsimilarity variable,  𝑥/𝐿
𝜌:	Fluid density
𝜇:	Dynamic viscosity
𝜈:	Kinematic viscosity
𝜃:	Dimensionless temperature profile in(2.8).

Subscripts

𝑠:	Solid
𝑤:	Wall
∞:	Free stream.