Hydrodynamic Anisotropy Effects on Radiation-Mixed Convection Interaction in a Vertical Porous Channel

The effects of hydrodynamic anisotropy on the mixed-convection in a vertical porous channel heated on its plates with a thermal radiation are investigated analytically for fully developed flow regime. The porous medium is anisotropic in permeability whose principal axes are oriented in a direction that is oblique to the gravity. The generalized Brinkman-extended Darcy model which allows the no-slip boundary-condition on solid wall is used in the formulation of the problem. The flow reversal, the thermal radiation influence for natural, and forced convection are considered in the limiting cases for low and high porosity media. It was found that the anisotropic permeability ratio, the orientation angle of the principal axes of permeability and the radiation parameter affected significantly the flow regime and the heat transfer.


Introduction
Convection heat transfer in porous media is a fundamental importance in such technologies as geothermal exploitation, oil recovery, radioactive waste management, insulation of building and cold storage, drying processes, transpiration cooling, powder metallurgy, agricultural engineering, solidification and binary alloys, etc. It is also important to geophysics and environmental sciences. Much of this activity, both numerical and experimental, has been summarized by Nield and Bejan (1999) [1]. Moreover, a relevant recent work on the topic is done by Chen et al. (2000) [2].
Thermal radiation always exists and can strongly interact with convection in many situations of engineering interest. The convection heat transfer in a porous channel (or in an enclosed space) in the presence of thermal radiation continues to receive considerable attention because of its importance in many practical applications such as furnaces, combustion chambers, cooling towers, rocket engines and solar collectors. During the last decade, many experimental and numerical investigations on the phenomenon of the interaction of natural or mixed convection with thermal radiation in vertical porous channels or enclosures have been presented. Mahmud and Fraser (2003) [3] have examined analytically the effects of radiation heat transfer on magnetohydrodynamic mixed convection through a vertical channel packed with fluid saturated porous substances. Basing their analysis on the first and second laws of thermodynamics, these authors have given special attention to entropy generation characteristics and their dependency on the various governing parameters of the problem. Jbara et al. (2011) [4] have conducted a numerical study to investigate the effects of the radiation heat transfer on the transient natural convection in a vertical channel opened at both ends and filled with a fluid saturated porous medium. Taking into account the Rosseland approximation in the energy equation, they have showed that approximation is only valid for large Planck number values for temperature ratios closed to unity and/or for single scattering albedo near or equals to one. An unsteady natural convection magnetohydrodynamic flow of a rotating viscous electrically fluid in a vertical channel with radiation effects has been investigated by Chauban and Rastogi (2010) [5] who obtained analytical expressions for velocity and temperature fields and discussed the effects of radiation on the phenomenon. A transient study of coupled natural convection and radiation in a porous vertical channel using the finitevolume method has been conducted by Slimi et al. (2004) [6]. The results indicated that the controlling parameters of the problem have significant effects on the flow and the thermal fields and on the transient process of heating and cooling of the medium. They have demonstrated that the volumetric flow rate and the convected heat flux at the channel's exit increase when the Planck number is decreased and/or the optical thickness and the wall emissivity are increased.
Most of the existing theories and experimental investigations on the topic, are concerned isotropic porous media. However, in several applications the porous materials are anisotropic. Such porous media are in fact encountered in numerous systems in industry and nature. As examples, we can cite fibrous materials, biological materials, geological formations, and oil extraction. The inclusion of more physical realism in the matrix properties of the medium is important for the accurate modeling of the anisotropic media. Anisotropy is generally a consequence of a preferential orientation and asymmetric geometry of the grain or fibbers which constitute the porous medium. Despite its vast range of applications, convection in such anisotropic porous media has received relatively little attention.
Thermal convection in a porous medium with anisotropic permeability was first considered by Castinel and Combarnous (1974) [7] who conducted an experimental and theoretical investigation in a horizontal layer heated from below. Their results were extended by Epherre (1975) [8] and Kvernvold and Tyvand (1979) [9] who considered a more general type of anisotropy. Convection attributable to side-heating in a vertical layer has been considered by several authors. For example, the case of square cavity with homogeneous media that is both thermally and hydrodynamically anisotropic has been investigated analytically by Kimura et al. (1993) [10] and numerically by Ni and Beckermann (1991) [11]. The effect of both anisotropy of permeability and thermal conductivity on the overall Nusselt number was found to be equally significant. A few studies have also been concerned with the case when the principal axes of anisotropy of the porous structure are inclined with respect to the gravity force. For this situation, the onset of motion in a porous layer heated from below was predicted by Tyvand and Storesletten (1991) [12] and Zhang et al. (1993) [13]. It was demonstrated that the influence of the anisotropy orientation considerably modifies the stability limit. Convection heat transfer in a vertical cavity heated from the side with various thermal boundary conditions has been investigated by Zhang (1993) [14], Degan et al. (1995) [15] and Vasseur (1997, 2002) [16] [17]. It was demonstrated that the heat transfer was maximum when the principal axis with higher permeability is parallel to the vertical direction and minimum when it is perpendicular.
The contemporary trend in the field of heat transfer and thermal design is to apply a second law of thermodynamics analysis, and its design-related concept of entropy generation minimization (see Bejan (1996) [18]). Entropy generation is associated with thermodynamic irreversibility which is present in all heat transfer processes. Because of the abundance of publications on entropy, no attempt is made here to review the literature.
The present work deals with an analytical study of coupled fluid flow and heat transfer by mixed convection and radiation in a vertical channel opened at both ends and filled with a fluid-saturated porous medium. The bounding walls of the channel are isothermal and gray. The effects of hydrodynamic anisotropy of the porous medium will be investigated, since the physical problem is of significant importance to many engineering-related applications.

Mathematical Formulation
The problem, under different considerations, concerns an optically thin and electrically conducting fluid flowing through a vertical channel opened at both ends and filled with a porous medium as shown by the physical model considered in Figure 1. The bounding walls of the channel are isothermal and kept at the same or different temperatures. Moreover, these bounding walls are submitted to the influence of radiative intensities 1 2 I − and 2 1 I − exerting reciprocally actions on the right and the left walls, respectively. The depth of the channel along the z′axis is assumed sufficiently long compared to other directions, x′ being the direction of the flow and y′direction perpendicular to that latter. A fully developed, laminar and induced flow is maintained through this porous channel upward, so that the natural convection aids the forced flow. As considered in nature in fact, the porous medium is anisotropic in flow permeability, the permeabilities along the two principal axes of the porous matrix are denoted by 1 K and 2 K . The anisotropy in flow permeability of the porous medium is then characterized by the permeability ratio * 1 2 K K K = and the orientation ϕ , defined as the angle between the horizontal direction and the principal axis with permeability 2 K . The principal directions of the thermal conductivity are assumed to coincide to those of the coordinate axes. The saturating fluid is viscous, incompressible and assumed to be in local thermodynamic equilibrium with the solid matrix everywhere and its thermophysical properties are assumed constant, except for the density in the buoyancy term in the momentum equation.
Under the above approximations, the equations governing the conservation of mass, momentum (generalized Brinkman-extended Darcy's law) and energy can be written as follows (see [3] [17]): where ′ V is the superficial flow velocity, p′ is the pressure, µ the dynamic viscosity, eff µ apparent dynamic viscosity for Brinkman's model, ( ) p f c ρ the heat capacity of the fluid, ρ the density, α the thermal diffusivity, T ′ the temperature and R q ′ represents the radiative heat flux. The symmetrical second-order permeability tensor K is defined as:  cos sin sin cos sin cos sin cos The energy Equation (3) takes into account the radiative heat flux term R q ∇ ⋅ ′ . It is a non-linear term which will be evaluate to have approximate analytical solution for the present problem. Due to the fact that the radiative heat flux is normal to the bounding walls, 0 and only the term 0 will be considered in the energy equation, and for which an approximate mathematical model will be estimated. According to Arpaci et al. [19], one can write as follows: where ς , Ω and b E are the monochromatic absorption coefficient, solid angle, and blackbody emissive power, respectively. The second term at the right hand side of Equation (5) represents the effect of boundaries. Assuming hemispherical isotropy, and replacing actual intensity with two-stream intensities (outward and inward) as shown in Figure 1, the second term of the right hand side of Equation (5) can be written as Ω and b E are the monochromatic absorption coefficient, solid angle, and blackbody emissive power, respectively. The second term at the right hand side of Equation (5) represents the effect of boundaries. Assuming hemispherical isotropy, and replacing actual intensity with two-stream intensities (outward and inward) as shown in Figure 1, the second term of the right hand side of Equation (5) can be written as ( ) Considering that the absorption is negligible for a thin gas, 1 2 I − and 2 1 I − remain uniform across the thickness of the gas between the plates (see Figure 1). Consequently, 1 2 I Substituting results obtained for these intensities into Equation (6) and noting that 1 ρ ε + = (the transmissivity being negligible in this study) and π b I E = , one can have the following approximation E are the monochromatic absorption coefficients on the left and the right walls, respectively.
Taking into account that the emissivities of the left and the right walls are identical, (i.e., 1 2 ε ε = ), the substitution of Equation (9) into Equation (5) yields ( ) In the above equation, ( ) can be replaced by the mean emissive power 0 b E , such that the two terms, b E and 0 b E can be further expressed as 4 T σ ′ and 4 0 T σ ′ , where σ is the Stefan-Boltzmann constant. These terms in 4 T ′ appearing in the energy equation have to be linearized in order to obtain analytical solution. Expressing T as a Taylor series about 0 T ′ with the assumption of small Truncating the above series after the second term and using the definitions of b E and 0 b E , Equations (10) and (11) Introducing the Boussinesq approximation and assuming that when the flow is fully developed in the channel, the axial (x'-direction) velocity depends only of the transverse coordinate y' (i.e., ( ) u u y ′ ′ ′ = ), and then from the continuity equation, the transverse velocity component must be zero (i.e., 0 v′ = ). Moreover, accounting for uniform wall temperature, one gets 0 T y ′ ′ ∂ ∂ =, and the governing equations (for steady convective flow) may be written in reduced form as The hydrodynamic and thermal boundary conditions for the vertical channel are Taking dRe (and d), 0 U ′ , 0 U ρ ′ and ΔT as respective dimensional scales for length x′ (and y′ ), velocities u′ and v′ , pressure and temperature, the governing Equations (14) and (17) can be written in nondimensional form as − . The boundary conditions, Equations (19) and (20) become (21) and (24) and the boundary conditions Equations (26) and (27), it is seen that the present problem is governed by eight dimensionless parameters, namely Ra , Da , Re , Pe , λ , Pl , τ and θ . However, in the present study, the value of eff µ in Brinkman's extension is taken, as a first approximation, equal to µ (i.e., 1 λ = ).

Fully Developed Flow Solution
Using the boundary conditions for the temperature, Equation (24) can be integrated to give the following fully developed temperature profile By substituting Equation (28) into Equation (22), and using the boundary conditions, Equations (26) and (27), the velocity profile is obtained as follows In the above expression of the distribution of the velocity, the parameter γ still remains to be evaluated. Substituting of Equation (29) into the continuity Equation (21) then gives directly the pressure variation as It is noticed that when 0 This result is in agreement with that which has been found by Degan and Vasseur [17]. The wall friction is defined by the following expression where the plus and minus signs correspond to the left and the right walls. Hence, on the left wall the friction is expressed as follows while, on the right wall, for the friction one can have such that the average friction defined by Concerning theory of fully developed confined convection including flow reversal in vertical channels, when buoyancy effects are increased (i.e., when the heat flux is increased), the fluid will accelerate near the walls. Then, mass conservation requires that the fluid decelerates in the center of the channel. Consequently, if buoyancy effects are strong enough, a minimum will form in the velocity distribution at the channel centerline. For even stronger buoyancy effects, a flow reversal will form at the channel centerline. In the case of aiding mixed convection through the porous channel, when a reverse flow occurs, the relatively lower velocity negative flow passes along side the cold wall hence carries a lower level of thermal energy. Since a net dimensionless mass flow is fixed, an equal quantity of fluid is added to the fluid flowing in the positive (upward) direction in this fluid flows adjacent to the hot wall thereby carrying a larger amount of energy.
Following Aung and Worku [20] and many authors, it is possible to deduce an expression or a criterion by which to predict whether flow reversal occurs. It is seen that, when 1 θ < , the occurrence of flow reversal is given by the condition Applying Equation (29), the above condition translates into G. Degan et al.
Accordingly, one can deduce the flow reversal function Two cases are of interest, one with 1 ξ  and the other with ): This case corresponds to a weaker anisotropic porous medium situation for which the resistance resulting from the boundary effects is predominant with respect to that due to the solid matrix, as Da → ∞ when 0 ξ → . This situation approaches to the fluid medium case in which the anisotropic effects of the porous medium are irrelevant. As a result, the pressure gradient, from Equation (31), is given by Consequently, as ξ goes to zero, the limit of the pressure gradient is written as Then, the velocity profile, Equation (29) and its limit, as ξ goes to zero, become as follows The average friction and its limit, as ξ goes to zero, are Also, the flow reversal function, Equation (41) and its limit can be written as Consequently, as d R goes to zero, applying Equation (49), the function of the flow reversal function's limit and its limit ( ) Similar results has been obtained by Degan and Vasseur [17]. The same remark is also made when comparing the result above with that one obtained by Aung and Worku [20] for the situation corresponding to a pure fluid medium in absence of radiation condition imposed on the walls.

Limiting Case of Forced Convection Solution
The limiting case of forced convection solution is obtained by setting ( ) 0 Ra Pe = in Equation (29). Then the velocity distribution is given by Here, the pressure variation is determined by the following expression The results presented above will be specified for two cases of interest ( 1 ξ  and 1 ξ  ).
The pressure gradient is written • 1 ξ  : In this situation, one can obtain the following results for the velocity profile and for the pressure gradient

Limiting Case of Natural Convection Solution
This important limiting case will be studied by setting 0 p x γ =∂ ∂ = . For natural convection, the entrance velocity cannot be specified. In this case, new dimensionless variables should be used as and, as 1 λ = in this study, Equation (22) can be rewritten as follow The solution for the velocity distribution may be written in terms of variables utilized herein, giving In the above velocity profile, the corresponding expressions to B and C are those indicated in Equation (30) G. Degan et al.

32
When 0 d R → , Equation (67) gives similar expression to that one obtained by Degan and Vasseur [17] corresponding to the case studied here when the radiation effect on the side walls is not taken into account.
The total heat absorbed by the fluid in traversing the channel is Writting Q′ in its dimensionless form and using Equations (28) and (65), the total heat absorbed is obtained as where B and C are the corresponding expressions indicated in Equation (30). An average Nusselt number may be defined as is the Darcy-modified Rayleigh number and Q the total heat absorbed, calculated in Equation (69).
The flow reversal condition for the limiting case of natural convection is ( ) such that the flow reversal criterion becomes Taking 0 d R → , one obtains from Equation (72), similar flow reversal occurrence criterion to that one obtained by Degan and Vasseur [17].
We notice here two cases of interest, the first one with 1 ξ  and the second with 1 ξ  .
• 1 ξ  :It is found in this case that, the velocity and the flow volume rate and the average Nusselt number are expressed as follows Similar result has been found by Degan and Vasseur [17] for analogous situation when the thermal radiation effect is not taking into account.

Results and Discussion
The effects of varying d R , the radiation parameter, d R and ϕ , the anisotropic parameters on the velocity profile for the case of mixed convection using Equation (29) are presented in Figures 2-4 for 5 Ra Pe = , 0.5 θ = and 0.01 Da = . In Figure 2, the radiation parameter d R affects the flow by promoting the convective circulation within the channel near the hot wall. The inverse behavior is obtained when 0 d R → for which a reverse flow is observed as the fluid near the cold wall sinks and the velocity decreases more and more. Such a conclusion has also been made by Degan and Vasseur [17] when studying mixed-convection through a vertical anisotropic porous channel with oblique principal axes without radiation effect. Aung and Worku [20] explained that, the reverse motion which occurs in the flow field is due to the fact that the ratio of wall temperature differences θ is less than unity.
From Figure 3 and , it is observed that near the hot wall, the velocity increased as the anisotropic parameters of the porous medium * K and ϕ are made higher. Moreover, the curves show that the peak value of velocity increases rapidly near the hot wall as ϕ or * K increases. Contrary to the following behavior, increasing values of these parameters contribute to slow the convective circulation along the warm wall as it is shown in Figure 3 when 0 ϕ →  for * 1 K < ( * 0.1 K = ), all curves sink in the neighborhood of the cold wall but the reversal flow is absent. We concluded that small values of anisotropic parameters of the porous medium contribute to stabilize the fluid motion. Figure 5 shows the dimensionless temperature profile for mixed convection when 0.5 θ = for various values of the radiation parameter d R . Curves plotted in Figure 5 reveal that the temperature decreases with an increase of d R . When 0 d R = , the temperature distribution is a function only of the ratio of wall temperature differences θ as predited by Aung and Worku [21]. Moreover, for all values of 0 d R > , the temperature is lower than that of the absence of thermal radiation. One can deduce that, the thermal radiation contributes to reduce the temperature  in the flow field when 1 θ < . Increasing d R means that the release of the heat energy from the flow region and so the fluid temperature decreases.
The effects of the anisotropic parameters of the porous matrix and the thermal radiation parameter on the gradient of the pressure are presented in Figure 6 and Ra Pe = and 0.2 θ = . From Figure 6 and Figure 7, it is noticed that, when Da is small enough tends asymptotically toward a constant value that depend on * K and d R . The limit 0 Da → corresponds to a pure Darcy medium situation for which ( ) given by Equation (62). In pure Darcy medium, the effects of both radiation and anisotropic parameters are observed to be significant. The increase of both of them diminishes the gradient of the pressure. This affects consequently the fluid motion and contributes to reduce the reversal flow as shown in Figure 2. We concluded that radiation effect contributes to stabilize the fluid motion. For example, when 3 10 Da ≈ for * 0.1 K = and * 2 K = respectively. As expected from equation 43, when the permeability of the porous medium Da is increased, the boundary frictional resistance becomes gradually more important and adds to the bulk frictional drag induced by the solid matrix to slow the convection motion. As consequence, the effects of varying the values of the anisotropic parameters of the porous medium and the thermal radiation parameter become less and less important and the present solution approaches that for a pure viscous fluid indicated in Equation (42).
In Figure 8   varying anisotropic parameters of the porous matrix are observed significant, since ( ) fRe − decreases as ϕ is made larger. When the Darcy number is made large enough, the results indicate that the curves, for a given value of ϕ , tend asymptotically toward the pure fluid situation (Equation (46)). The Darcy number required to reach this limit increases as the value of ϕ is made higher. For example, this happens at The effects of the radiation parameter R d and the permeability ratio K * on the Nusselt number is illustrated in Figure 9 for Ra = 100, Da = 10 −2 , 45 ϕ =  , 0.9 θ = , 0 0.8 T = and 0.8 It is observed that, for a given value of R d , the convective heat transfer increases as R d is made smaller. For a fixed value of R d , Nu tends asymptotically toward a constant value as K * is made small enough *  heat transfer for this limiting case, it is observed that As it can be observed in Figure 9, the convection becomes less and less important as * K is made larger. Thus, for each value of d R considered in Figure 9, the heat transfer rate approaches the pure condition regime, 1 Nu → as * K is made large enough. The value of * K necessary to reach the pure condition regime depends upon d R . For instance, for 0.5 d R = , pure condition is reached when * 5.74 K ≈ while for 2 d R = , * 3.78 K ≈ . The influence of the anisotropic orientation ϕ on Nusselt number Nu , is presented in Figure 10 where the bounding walls are submitted to various values of the radiation parameter R d and K * when 100 Ra = , permeability ratio * K is made smaller than one (i.e., * 1 K < ). Moreover, for a given * K , the higher heat transfer rate is obtained when d R is smaller. For the range 0 90 ϕ < <   and * 1 K < , Figure 10 indicates that Nu is minimum at 0 ϕ =  for which permeability in vertical direction is minimum, but is maximum at 90 ϕ =  .
The inverse is observed for * 1 K > where the convective heat transfer is maximum at 0 ϕ =  and minimum at 90 ϕ =  . For a given d R , when ( ) * * 1 1 K K > < , Nu is maximum (minimum) at 0 ϕ =  and minimum (maximum) at 90 ϕ =  . The maximum (minimum) increased (decreases) when d R is small (great) enough. Therefore, a maximum (minimum) convective heat transfer is reached when orientation of the principal axis with higher permeability of anisotropic porous medium is parallel (perpendicular) to the gravity. A Similar result has been obtained by Degan and Vasseur [17] while studying mixed convection through a vertical anisotropic porous channel with oblique principal axes.

Conclusions
A study has been made of mixed convection through a parallel-plate vertical porous channel submitted to a thermal radiative flux on its wall. The porous medium is assumed to be hydrodynamically anisotropic with its principal axes oriented in a direction that is oblique to the gravity.
Analytical expressions valid for fully developed flow and based on the generalized Brinkman-extended Darcy are obtained. The main conclusions of the present analysis are: • Both thermal radiation and anisotropic parameters have a strong influence on the fluid motion and the heat transfer through the parallel-plate vertical porous channel • In the pure Darcy medium 0 Da → , the velocity, the average wall friction, the pressure gradient and the heat transfer rates are similar to those given by a pure Darcy analysis. In this limit, the orientation angle of the principal axes of permeability and the permeability ratio and the thermal radiation parameter affect considerably the convection heat transfer. • In the fluid medium ( 1 Da > ), the resistance resulting from the boundary effects is predominant and adds to the bulk frictional drag induced by the solid porous matrix to slow down the convection motion. Thus, the effects of anisotropic parameters (i.e., * K and ϕ ) and thermal radiation parameter d R on the convective heat transfer are progressively irrelevant, and the solutions obtained become independent of them.
• The effects of increasing values of the anisotropic parameters and the thermal radiation parameter tend to decrease the temperature and the gradient of the pressure. Moreover, the decrease of values of thermal radiation reduces the fluid velocity and makes the reversal flow occurs. • For a given value of the thermal radiation parameter, a maximum (minimum) heat transfer rate through the parallel-plate vertical porous channel is obtained when the porous matrix is oriented in such a way that the principal axis with higher permeability is parallel (perpendicular) to the gravity.