Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 261871, 19 pages

doi:10.1155/2012/261871

Research Article

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

Ahmet Kaya

Department of Mechanical Engineering, Kahramanmaras Sutcu Imam University, 46100 Kahramanmaras, Turkey

Received 8 May 2012; Accepted 4 June 2012

Academic Editor: Serge Prudhomme

Copyright © 2012 Ahmet Kaya. 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.

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 and uniform electric field [21]. Magnetic and electric fields are known from Maxwell’s equation that and . 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])

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

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

The boundary conditions for the wall of plate are given

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

To seek a solution, the following dimensionless variables are introduced:

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

The transformed momentum and energy equations together with the boundary conditions, (2.2), (2.3), and (2.5), can be written as

with the boundary conditions:

where is the conjugate heat transfer parameter. It should be noticed that for the limiting case of , 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:

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 . is the average Grashof number; is the average Reynolds number; Pr is the Prandtl number.

From the definition of the dimensionless wall temperature, it can be shown that

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

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

Table 2: Comparison of the values −θ′(ξ, 0) for various values with , , , , , , and .

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

Table 3: Comparison of the values −θ′(ξ, 0) for various values and with , and .

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 , 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 , whereas free convection is the dominant mode when . 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 , , , , , , and .

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

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 () and nonisothermal plate (). 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 while ε =0.98, , , , , , and .

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 (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 () and a nonisothermal plate (). 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 , , , , , and .

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 () and nonisothermal plate (). 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 , , , , , and .

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 () and nonisothermal plate () 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 , , , , , and .

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 () and nonisothermal plate () 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 while , , , , , , and .

(a)

(b)

(c)

Figure 9: Effect of on the dimensionless interfacial temperature (a), local skin friction (b), and local heat transfer (c) parameters against the streamwise distance while , , , , , and .

The increasing of the porous medium parameter decreases momentum and thermal boundary layer thickness and increases velocity and temperature profiles for both isothermal () and nonisothermal () 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 (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 .

(a)

(b)

Figure 10: Dimensionless velocity (a) and temperature (b) profiles for different while , , , , , and .

(a)

(b)

(c)

Figure 11: Effect of on the dimensionless interfacial temperature (a), local skin friction (b), and local heat transfer (c) parameters against the streamwise distance while , , , , , and .

Figure 12 shows the effect of drag (inertia) coefficient of porous medium (or Forchheimer parameter) on the dimensionless velocity and temperature profiles at . 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 , , , , , , and .

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

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

: | Magnetic field strength
: | Specific heat of the convective fluid
: | Electric field parameter
: | Dimensionless stream function
: | Forchheimer parameter
: | Grashof number
: | Hartman number
: | Effective thermal conductivity of porous medium
: | Porous medium permeability,
: | Porous medium parameter
: | Plate length
: | Magnetic parameter,
: | Prandtl number
: | 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,  
: | Nonsimilarity variable,  
: | Fluid density
: | Dynamic viscosity
: | Kinematic viscosity
: | Dimensionless temperature profile in(2.8).

Subscripts

: | Solid
: | Wall
: | Free stream.

References

A. J. Chamkha, C. Issa, and K. Khanafer, “Natural convection from an inclined plate embedded in a variable porosity porous medium due to solar radiation,” International Journal of Thermal Sciences, vol. 41, no. 1, pp. 73–81, 2002.

M. Kaviany, Principles of Heat Transfer in Porous Media, Springer, New York, NY, USA, 1995.

I. Pop and and D. B. Ingham, Convective Heat Transfer, Pergamon, Amsterdam, The Netherland, 2001.

D. B. Ingham and I. Pop and, Transport Phenomena in Porous Media. I, Pergamon, New York, NY, USA, 1998.

D. B. Ingham and I. Pop and, Transport Phenomena in Porous Media. II, Pergamon, New York, NY, USA, 2002.

D. A. Nield and A. Bejan, Convection in Porous Media, Springer, New York, NY, USA, 1999.

M. Miyamoto, J. Sumikawa, T. Akiyoshi, and T. Nakamura, “Effects of axial heat conduction in a vertical flat plate on free convection heat transfer,” International Journal of Heat and Mass Transfer, vol. 23, no. 11, pp. 1545–1553, 1980.

E. M. Sparrow and M. K. Chyu, “Conjugated forced convection-conduction analysis of heat transfer in a plate fin,” ASME Journal of Heat Transfer, vol. 104, pp. 204–206, 1982.

M. I. Char, C. K. Chen, and J. W. Cleaver, “Conjugate forced convection heat transfer from a continuous, moving flat sheet,” International Journal of Heat and Fluid Flow, vol. 11, no. 3, pp. 257–261, 1990.

T. Y. Wang, “The coupling of conduction with mixed convection of micropolar fluids past a vertical flat plate,” International Communications in Heat and Mass Transfer, vol. 25, no. 8, pp. 1075–1084, 1998.

I. Pop, D. B. Ingham, and Y. Yuan, “Mixed convective conjugate heat transfer from a vertical flat plate,” Zeitschrift fur Angewandte Mathematik und Mechanik, vol. 76, no. 5, pp. 281–289, 1996.

C. L. Chang, “Numerical simulation of micropolar fluid flow along a flat plate with wall conduction and buoyancy effects,” Journal of Physics D, vol. 39, no. 6, article 019, pp. 1132–1140, 2006.

A. A. Mamun, Z. R. Chowdhury, M. A. Azim, and M. A. Maleque, “Conjugate heat transfer for a vertical flat plate with heat generation effect,” Nonlinear Analysis. Modelling and Controlling, vol. 13, no. 2, pp. 213–223, 2008.

K. L. Hsiao and C. H. Hsu, “Conjugate heat transfer of mixed convection for viscoelastic fluid past a horizontal flat-plate fin,” Applied Thermal Engineering, vol. 29, no. 1, pp. 28–36, 2009.

I. Pop and T. Y. Na, “Conjugate free convection over a vertical slender hollow cylinder embedded in a porous medium,” Heat and Mass Transfer, vol. 36, no. 5, pp. 375–379, 2000.

A. Z. Vaszi, D. B. Ingham, D. Lesnic, D. Munslow, and I. Pop, “Conjugate free convection from a slightly inclined plate embedded in a porous medium,” Zeitschrift fur Angewandte Mathematik und Mechanik, vol. 81, no. 7, pp. 465–479, 2001.

F. Méndez, C. Treviño, I. Pop, and A. Liñán, “Conjugate free convection along a thin vertical plate with internal nonuniform heat generation in a porous medium,” Heat and Mass Transfer, vol. 38, no. 7-8, pp. 631–638, 2002.

N. H. Saeid, “Conjugate natural convection in a vertical porous layer sandwiched by finite thickness walls,” International Communications in Heat and Mass Transfer, vol. 34, no. 2, pp. 210–216, 2007.

A. Z. Vaszi, L. Elliott, D. B. Ingham, and I. Pop, “Conjugate free convection from a vertical plate fin with a rounded tip embedded in a porous medium,” International Journal of Heat and Mass Transfer, vol. 47, no. 12-13, pp. 2785–2794, 2004.

M. A. Hossain, A. Nakayama, and I. Pop, “Conjugate free convection of non-Newtonian fluids about a vertical cylindrical fin in porous media,” Heat and Mass Transfer, vol. 30, no. 3, pp. 149–153, 1995.

D. Pal and H. Mondal, “Hydromagnetic non-Darcy flow and heat transfer over a stretching sheet in the presence of thermal radiation and Ohmic dissipation,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 5, pp. 1197–1209, 2010.

A. J. Chamkha, “Hydromagnetic natural convection from an isothermal inclined surface adjacent to a thermally stratified porous medium,” International Journal of Engineering Science, vol. 35, no. 10-11, pp. 975–986, 1997.

G. Lauriat and R. Ghafir, “Forced convective heat transfer in porous media,” in Handbook of Porous Media, K. Vafai and H. Hadim, Eds., pp. 201–204, Marcel Dekker, New York, NY, USA, 2000.

A. Pozzi and M. Lupo, “The coupling of conduction with laminar natural convection along a flat plate,” International Journal of Heat and Mass Transfer, vol. 31, no. 9, pp. 1807–1814, 1988.

M. Char and C. L. Chang, “Effect of wall conduction on natural convection flow of micropolar fluids along a flat plate,” International Journal of Heat and Mass Transfer, vol. 40, no. 15, pp. 3641–3652, 1997.

A. V. Kuznetsov and D. A. Nield, “Boundary layer treatment of forced convection over a wedge with an attached porous substrate,” Journal of Porous Media, vol. 9, no. 7, pp. 683–694, 2006.

A. J. Chamkha, M. Mujtaba, A. Quadri, and C. Issa, “Thermal radiation effects on MHD forced convection flow adjacent to a non-isothermal wedge in the presence of a heat source or sink,” Heat and Mass Transfer, vol. 39, no. 4, pp. 305–312, 2003.

O. Aydin and A. Kaya, “Radiation effect on MHD mixed convection flow about a permeable vertical plate,” Heat and Mass Transfer, vol. 45, no. 2, pp. 239–246, 2008.

O. Aydin and A. Kaya, “Non-Darcian forced convection flow of viscous dissipating fluid over a flat plate embedded in a porous medium,” Transport in Porous Media, vol. 73, no. 2, pp. 173–186, 2008.

T. Cebeci and P. Bradshaw, Momentum Transfer in Boundary Layers, Hemisphere, Washington DC, USA, 1977.

H. S. Takhar and O. A. Beg, “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 medium,” International Journal of Energy Research, vol. 21, no. 1, pp. 87–100, 1997.

O. Aydin and A. Kaya, “Mixed convection of a viscous dissipating fluid about a vertical flat plate,” Applied Mathematical Modelling, vol. 31, no. 5, pp. 843–853, 2007.

C. L. Chang, “Buoyancy and wall conduction effects on forced convection of micropolar fluid flow along a vertical slender hollow circular cylinder,” International Journal of Heat and Mass Transfer, vol. 49, no. 25-26, pp. 4932–4942, 2006.

C. L. Chang, “Numerical simulation for natural convection of micropolar fluids flow along slender hollow circular cylinder with wall conduction effect,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 3, pp. 624–636, 2008.

M. Kumari and G. Nath, “Radiation effect on mixed convection from a horizontal surface in a porous medium,” Mechanics Research Communications, vol. 31, no. 4, pp. 483–491, 2004.