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 k1, the Forchheimer parameter F*,
the mixed convection parameter Ri, the magnetic parameter Mn, and the electric
field parameter E1 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 L and finite thickness b (L≫b) 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 u∞ and T∞, respectively. The temperature of the inside surface of the plate is maintained at a constant temperature of T0, where T0 > T∞. The flow region is exposed under uniform transverse magnetic fields B→=(0,B0,0) and uniform electric field E→=(0,0,-E0) [21]. Magnetic and electric fields are known from Maxwell’s equation that ∇·B→=0 and ∇×E→=0. When magnetic field is not so strong then electric field and magnetic field obey Ohm’s law J→=σ(E→+q→×B→), where J→ is the Joule current, σ is the magnetic permeability and q→ 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).
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])
(2.1)∂u∂x+∂υ∂y=0,(2.2)1ε2[u∂u∂x+υ∂u∂y]=vε∂2u∂y2+σρ[E0B0-B02(u-u∞)]-vK(u-u∞)-F(u2-u∞2)+gβ(T-T∞)(2.3)u∂T∂x+υ∂T∂y=keρcp(∂2T∂y2).
The above equations are called Brinkman-Forchheimer-extended-Darcy equations [23]. Here u and v are the velocity components in the x and y directions, respectively; v is the kinematic viscosity; g is the acceleration due to gravity; ke is the effective thermal conductivity of porous medium; ρ is the density of the fluid; β is the coefficient of thermal expansion; T is the temperature of the fluid inside the thermal boundary layer; K is the permeability of the porous medium; F is the empirical constant (Forchheimer number) in the second-order resistance and setting F=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 (ke/ρcp) 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
(2.4)x=0,y>0,T=T∞,u=u∞,x>0,y=0,T=Tw(x),u=0,υ=0,y⟶∞,T⟶T∞,u⟶u∞,
where subscripts w and ∞ refer to the wall and the boundary layer edge, respectively. In addition, Tw(x) 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 Tw(x). 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, b, is small compared with its length, L, 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
(2.5)d2Tdy2|s=0;0≤x≤L;-b<y≤0.
The boundary conditions for the wall of plate are given (2.6a)Aty=-b,Ts=T0,(2.6b)Attheinterface(y=0),Ts=T(x,0):-ksdTdy|y=0,s=-kf∂T(x,0)∂y|y=0,f,where ks and kf 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 Tw at the interface is shown to be
(2.7)Tw(x)=T(x,0)=bkfks∂T(x,0)∂y+T0.
To seek a solution, the following dimensionless variables are introduced:
(2.8)ξ(x)=xL,ψ(x,y)=(vu∞x)1/2f(ξ,η),η=y(u∞vx)1/2,θ=T-T∞T0-T∞,
where ψ(x,y) is the free stream function that satisfies (2.1) with u=∂ψ/∂y and υ=-∂ψ/∂x.
In terms of these new variables, the velocity components can be expressed as
(2.9)u=u∞f′,(2.10)υ=-(vu∞x)12x{12f+ξ∂f∂ξ-η2f′}.
The transformed momentum and energy equations together with the boundary conditions, (2.2), (2.3), and (2.5), can be written as
(2.11)1εf′′′+12ε2ff′′+Mnξ[E1-(f′-1)]-k1ξ(f′-1)-F*ξ[(f′)2-1]+Riξθ=ξε2(f′∂f′∂ξ-f′′∂f∂ξ)1Prθ′′+12fθ′=ξ(f′∂θ∂ξ-θ′∂f∂ξ),
with the boundary conditions:
(2.12)f(ξ,0)+2ξ∂f∂ξ=0,f′(ξ,0)=0,θ(ξ,0)-1=pξ-1/2θ′(ξ,0),f′(ξ,∞)=1,θ(ξ,∞)=0,
where p=(kf/ks)(b/L)ReL1/2 is the conjugate heat transfer parameter. It should be noticed that for the limiting case of p=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:
(2.13)Mn=HaRe,Ha=σB02L2μ,E1=ERe,E=E0B0v/L,k1=vLKu∞,F*=FL,Ri=GrLReL,GrL=gβ(Tw-T∞)KLv2,ReL=u∞Lv,Pr=μcpk=vα,
where Mn is the magnetic parameter; Ha is the Hartman number; E1 is the electric field parameter; k1 is the porous medium parameter; F* 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 x. GrL is the average Grashof number; ReL 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)θw=Tw-T∞T0-T∞.
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 x 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−6 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.
Comparison of the values −θ′(0, 0) for various values Pr at Mn=0.0, E1=0.0, ε=1.0, F*=0.0, Ri=0.0, k1=1.0, and p=0.0.
Pr
Kuznetsov and Nield [26]
Present study
0.1
0.158
0.148
1
0.332
0.332
5
0.570
0.577
10
0.730
0.728
20
0.910
0.918
30
1.050
1.055
40
1.150
1.155
50
1.245
1.244
60
1.320
1.321
70
1.390
1.390
80
1.450
1.452
90
1.510
1.514
100
1.570
1.573
Comparison of the values −θ′(ξ, 0) for various values ξ with Pr=0.7, Mn=1.0, Ri=1.0, ε=1.0, F*=0.0, k1=0.0, and p=0.0.
ξ
Chamkha et al. [27]
Aydin and Kaya [28]
Present study
0.0
0.293
0.293
0.293
0.01
0.322
0.330
0.330
0.02
0.325
0.350
0.350
0.03
0.363
0.367
0.367
0.04
0.380
0.379
0.379
0.05
0.390
0.390
0.390
0.06
0.400
0.400
0.400
0.07
0.409
0.408
0.408
0.08
0.410
0.415
0.415
Comparison of the values −θ′(ξ, 0) for various values ξ and ε with Pr=1,F*=0.5,Mn=0.0,E1=0.0,k1=1.0,Ri=0.0, and p=0.0.
ε=0.5
ε=0.75
ε=1.0
ξ
Aydin and Kaya [29]
Present study
Aydin and Kaya [29]
Present study
Aydin and Kaya [29]
Present study
0.0
0.367
0.367
0.346
0.346
0.332
0.332
0.01
0.368
0.368
0.350
0.350
0.338
0.338
0.02
0.369
0.369
0.353
0.353
0.343
0.343
0.03
0.371
0.371
0.355
0.355
0.348
0.348
0.04
0.372
0.372
0.358
0.358
0.352
0.352
0.05
0.374
0.374
0.361
0.361
0.356
0.356
0.06
0.375
0.375
0.364
0.364
0.360
0.360
0.07
0.376
0.376
0.366
0.366
0.364
0.364
0.08
0.377
0.377
0.369
0.369
0.368
0.368
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 p=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 E1 = 0.1, 0.5, 1.0, and 2.0; porous medium parameter k1 = 0.1, 0.5, 1.0 and 2; inertia parameter F* = 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 p, Ri, Mn, E1, k1, and F* 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, E1 = 0.1, k1 = 0.5, F* = 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 ks or higher convective cooling effect due to greater kf increases the value of p 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.
Dimensionless velocity (a) and temperature (b) profiles for different p while ε=0.98, Pr=1.0, Ri=1, Mn=1, E1=0.1, k1=0.5, F*=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 p with ξ are shown in Figures 3(a), 3(b), and 3(c), respectively where Ri = 1.0, Mn = 0.5, E1 = 0.1, k1 = 0.5, F* = 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 p (Figure 3(a)). Comparing with isothermal wall (p=0), an increase in the conjugate heat transfer parameter, p, causes a reduction in the interfacial temperature. This is because an increased value of p corresponds to a lower wall conductance ks and promotes a greater surface temperature variations (Figure 3(a)). The increases value of p 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].
Effect of conjugate heat transfer parameter p 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, E1=0.1, k1=0.5, and F*=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 (p=0) and nonisothermal plate (p>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.
Dimensionless velocity (a) and temperature (b) profiles for different Ri while ε =0.98, Pr=1.0, Mn=1, E1=0.1, k1=0.5, F*=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 p=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 (p=0) and a nonisothermal plate (p>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].
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, E1=0.1, k1=0.5, and F*=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 (p=0) and nonisothermal plate (p>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 p decreases velocity and temperature gradients at the wall.
Dimensionless velocity (a) and temperature (b) profiles for different Mn while ε=0.98, Pr=1.0, Ri=5, E1=0.1, k1=0.5, F*=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 (p=0) and nonisothermal plate (p=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.
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, E1=0.1, k1=0.5, and F*=0.5.
The dimensionless velocity (a) and temperature (b) profiles inside the boundary layer for different values of the electric field parameter E1 for the cases of isothermal plate (p=0) and nonisothermal plate (p>0.0) are illustrated in Figure 8. Increasing of the electric field parameter E1 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 E1 decreases the interfacial temperature (Figure 9(a)).
Dimensionless velocity (a) and temperature (b) profiles for different E1 while ε=0.98, Pr=1.0, Ri=5, Mn=0.5, k1=0.5, F*=0.5, and ξ=0.5.
Effect of E1 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, k1=0.5, and F*=0.5.
The increasing of the porous medium parameter k1 decreases momentum and thermal boundary layer thickness and increases velocity and temperature profiles for both isothermal (p=0.0) and nonisothermal (p=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 k1 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.
Dimensionless velocity (a) and temperature (b) profiles for different k1 while ε=0.98, Pr=1.0, Ri=5, Mn=0.5, E1=0.1, F*=0.5 and ξ=0.5.
Effect of k1 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, E1=0.1, and F*=0.5.
Figure 12 shows the effect of drag (inertia) coefficient of porous medium (or Forchheimer parameter) F* on the dimensionless velocity and temperature profiles at ξ=0.4. From this figure it is observed that the effect of drag coefficient F* 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.
Dimensionless velocity (a) and temperature (b) profiles for different F* while ε=0.98, Pr=1.0, Ri=5, Mn=0.5, E1=0.1, k1=0.5, and ξ=0.4.
Effect of F* 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, E1=0.1, k1=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:
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.
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, E1, and Ri decreases the velocity and temperature gradients at wall for nonisothermal cylinder (i.e., p>0).
An increase in the porous medium parameter k1 and Forchheimer parameter F* increases the local skin friction and local heat transfer parameters and decreases the dimensionless interfacial temperature distributions.
NomenclatureB0:
Magnetic field strength
cp:
Specific heat of the convective fluid
E1:
Electric field parameter
f:
Dimensionless stream function
F*:
Forchheimer parameter
Gr:
Grashof number
Ha:
Hartman number
ke:
Effective thermal conductivity of porous medium
K:
Porous medium permeability, m2
k1:
Porous medium parameter
L:
Plate length
Mn:
Magnetic parameter, M=Ha/Re
Pr:
Prandtl number
Re:
Reynolds number
T:
Temperature
u,υ:
Velocities in x and y directions, respectively
x,y:
Coordinates in horizontal and vertical directions, respectively.
Greek Symbolsσ:
Magnetic permeability
β:
Coefficient of thermal expansion
η:
Pseudosimilarity variable, yRex1/2/x
ξ:
Nonsimilarity variable, x/L
ρ:
Fluid density
μ:
Dynamic viscosity
ν:
Kinematic viscosity
θ:
Dimensionless temperature profile in(2.8).
Subscriptss:
Solid
w:
Wall
∞:
Free stream.
ChamkhaA. J.IssaC.KhanaferK.Natural convection from an inclined plate embedded in a variable porosity porous medium due to solar radiation200241173812-s2.0-003614856210.1016/S1290-0729(01)01305-9KavianyM.1995New York, NY, USASpringerPop andI.InghamD. B.2001Amsterdam, The NetherlandPergamonInghamD. B.Pop andI.1998New York, NY, USAPergamonInghamD. B.Pop andI.2002New York, NY, USAPergamonNieldD. A.BejanA.1999New York, NY, USASpringer1656781ZBL0947.76078MiyamotoM.SumikawaJ.AkiyoshiT.NakamuraT.Effects of axial heat conduction in a vertical flat plate on free convection heat transfer19802311154515532-s2.0-001908069910.1016/0017-9310(80)90158-1SparrowE. M.ChyuM. K.Conjugated forced convection-conduction analysis of heat transfer in a plate fin1982104204206CharM. I.ChenC. K.CleaverJ. W.Conjugate forced convection heat transfer from a continuous, moving flat sheet19901132572612-s2.0-0025493804WangT. Y.The coupling of conduction with mixed convection of micropolar fluids past a vertical flat plate1998258107510842-s2.0-003221524610.1016/S0735-1933(98)00098-0PopI.InghamD. B.YuanY.Mixed convective conjugate heat transfer from a vertical flat plate19967652812892-s2.0-21344471657ChangC. L.Numerical simulation of micropolar fluid flow along a flat plate with wall conduction and buoyancy effects2006396, article 019113211402-s2.0-3375116419810.1088/0022-3727/39/6/019MamunA. A.ChowdhuryZ. R.AzimM. A.MalequeM. A.Conjugate heat transfer for a vertical flat plate with heat generation effect2008132213223HsiaoK. L.HsuC. H.Conjugate heat transfer of mixed convection for viscoelastic fluid past a horizontal flat-plate fin200929128362-s2.0-5374909818310.1016/j.applthermaleng.2008.01.027PopI.NaT. Y.Conjugate free convection over a vertical slender hollow cylinder embedded in a porous medium20003653753792-s2.0-003427144610.1007/s002310000124VasziA. Z.InghamD. B.LesnicD.MunslowD.PopI.Conjugate free convection from a slightly inclined plate embedded in a porous medium20018174654792-s2.0-23044528531MéndezF.TreviñoC.PopI.LiñánA.Conjugate free convection along a thin vertical plate with internal nonuniform heat generation in a porous medium2002387-86316382-s2.0-003669330410.1007/s002310100244SaeidN. H.Conjugate natural convection in a vertical porous layer sandwiched by finite thickness walls20073422102162-s2.0-3384662653510.1016/j.icheatmasstransfer.2006.11.003VasziA. Z.ElliottL.InghamD. B.PopI.Conjugate free convection from a vertical plate fin with a rounded tip embedded in a porous medium20044712-13278527942-s2.0-184248236510.1016/j.ijheatmasstransfer.2004.01.001HossainM. A.NakayamaA.PopI.Conjugate free convection of non-Newtonian fluids about a vertical cylindrical fin in porous media19953031491532-s2.0-002921550110.1007/BF01476523PalD.MondalH.Hydromagnetic non-Darcy flow and heat transfer over a stretching sheet in the presence of thermal radiation and Ohmic dissipation2010155119712092-s2.0-7044952118510.1016/j.cnsns.2009.05.051ChamkhaA. J.Hydromagnetic natural convection from an isothermal inclined surface adjacent to a thermally stratified porous medium19973510-119759862-s2.0-0031200244LauriatG.GhafirR.VafaiK.HadimH.Forced convective heat transfer in porous media2000New York, NY, USAMarcel Dekker201204PozziA.LupoM.The coupling of conduction with laminar natural convection along a flat plate1988319180718142-s2.0-0024072374CharM.ChangC. L.Effect of wall conduction on natural convection flow of micropolar fluids along a flat plate19974015364136522-s2.0-0041630679KuznetsovA. V.NieldD. A.Boundary layer treatment of forced convection over a wedge with an attached porous substrate2006976836942-s2.0-3384715779610.1615/JPorMedia.v9.i7.70ChamkhaA. J.MujtabaM.QuadriA.IssaC.Thermal radiation effects on MHD forced convection flow adjacent to a non-isothermal wedge in the presence of a heat source or sink20033943053122-s2.0-0038617483AydinO.KayaA.Radiation effect on MHD mixed convection flow about a permeable vertical plate20084522392462-s2.0-5414909197710.1007/s00231-008-0428-yAydinO.KayaA.Non-Darcian forced convection flow of viscous dissipating fluid over a flat plate embedded in a porous medium20087321731862-s2.0-4254916747710.1007/s11242-007-9166-8CebeciT.BradshawP.1977Washington DC, USAHemisphereTakharH. S.BegO. A.Effects of transverse magnetic field, prandtl number and reynolds number on non-darcy mixed convective flow of an incompressible viscous fluid past a porous vertical flat plate in a saturated porous medium1997211871002-s2.0-0030990865AydinO.KayaA.Mixed convection of a viscous dissipating fluid about a vertical flat plate20073158438532-s2.0-3384559051310.1016/j.apm.2005.12.015ChangC. L.Buoyancy and wall conduction effects on forced convection of micropolar fluid flow along a vertical slender hollow circular cylinder20064925-26493249422-s2.0-3375034061310.1016/j.ijheatmasstransfer.2006.05.037ChangC. L.Numerical simulation for natural convection of micropolar fluids flow along slender hollow circular cylinder with wall conduction effect20081336246362-s2.0-3484886187410.1016/j.cnsns.2006.05.012KumariM.NathG.Radiation effect on mixed convection from a horizontal surface in a porous medium20043144834912-s2.0-234245331510.1016/j.mechrescom.2003.11.006