Effect of Constant Magnetic Field on Convective Heat Transfer through Anisotropic River Beds

Affiliation(s)

Laboratoire d’Energétique et de Mécanique Appliquées (LEMA), Ecole Polytechnique d’Abomey-Calavi, Cotonou, Bénin.

Laboratoire d’Energétique et de Mécanique Appliquées (LEMA), Ecole Polytechnique d’Abomey-Calavi, Cotonou, Bénin.

ABSTRACT

An analytical investigation is conducted to study the effect of magnetic field on convection heat transfer through packed porous beds which consists of a horizontal fluid layer (river bed) and a porous zone with anisotropic permeability and underlined by a surface heated by a constant temperature* T1*. The free surface of the fluid layer overlying the horizontal porous layer receives solar rays to length of day and is then considered heated isothermally at temperature *T2* such as *T1* < *T2*. Flow in porous medium is assumed to be governed by the generalized Brinkman-extended Darcy law and in the fluid layer by the Navier-Stokes model. The Beavers-Joseph condition is applied at the interface between the two layers. The influence of Hartmann number and hydrodynamic anisotropy on the convective phenomenon is investigated analytically. It is found that the magnetic field, the anisotropic permeability and the thickness of the porous lining, ε, have a strong influence of the geothermal convective flow and the heat transfer rate.

An analytical investigation is conducted to study the effect of magnetic field on convection heat transfer through packed porous beds which consists of a horizontal fluid layer (river bed) and a porous zone with anisotropic permeability and underlined by a surface heated by a constant temperature

1. Introduction

The first study concerning the effect of a magnetic field on the natural convection heat transfer in a rectangular porous cavity seems to be due to [1] . The vertical walls of the enclosure were maintained in isothermal and isoelectrical conditions. It was demonstrated that, for large Hartmann numbers, the electromagnetic magnetic force retards considerably the convective heat transfer. The influence of the Hartmann number on the Nusselt number was investigated numerically. The stability of a conducting fluid saturating a porous layer, in the presence of a uniform magnetic field, was investigated analytically by [2] . On the basis of the linear stability theory, the critical Rayleigh numbers for the onset of motion were obtained for various types of thermal and hydrodynamic boundary conditions. The case of a shallow cavity heated isothermally from the sides was considered by [3] . The effect of a magnetic field on the convective heat transfer was investigated analytically using matched asymptotic expansions. The results indicate that the retarding effect of the electromagnetic body Lorentz force decreases the strength of convection in the enclosure.

Several studies have been made on the effect of the magnetic field on the thermal convection through porous cavity. Our present research concerns the study of the thermal convection in a horizontal fluid-superposed porous layer and few studies have been made in this sense (see our literature magazine [4] ). Like that, we consider the convective heat transfer through a parallel-plate horizontal system consisting of a homogeneous porous bed underlying a single-component fluid layer whose upper surface is free and isothermally heated. **A** magnetic field, ** B** , is applied perpendicularly to the long sides. The lower plate bounding the porous layer is impermeable and is maintained at a constant temperature. The porous medium is homogeneous and anisotropic in permeability whose principal axes are arbitrary oriented, as it is seen in nature and for many realistic applications. Beavers and Joseph boundary condition model is applied at the permeable bounding interface between the two layers [5] . On the basis of the generalized Brinkman-extended Darcy model, of Navier-Stokes equations and of energy equation which takes into account the viscous dissipation, the effects of magnetic field, of anisotropic parameters of the porous matrix and of the influence of the depth ratio on velocity and temperature fields and heat transfer rate are investigated in detail.

2. Mathematical Formulation and Resolution

The physical model illustrating the problem under different considerations is shown in Figure 1. The system of height, h, consists of a horizontal parallel-plate porous channel of thickness, h_{m}, underlying a fluid layer whose upper plate (free surface) is permeable and exposed to a constant temperature T_{2}. The lower impermeable plate lining the non-erodible porous layer is maintained at a constant temperature T_{1} < T_{2}.

The axial and transverse coordinates are respectively x' and y', the latter being measured vertically upwards from the lower impermeable wall. The porous medium is anisotropic in flow permeability, the permeabilities along the two principal axes of the porous matrix are denoted by K_{1} and K_{2}. The anisotropy of the porous layer is characterized by the permeability ratio
${K}^{*}={K}_{1}/{K}_{2}$ and the orientation angle φ, defined as the angle between the horizontal direction and the principal axis with the permeability K_{2}.

Thus, the flow regime is divided into two zones:

Zone 1 (fluid layer) from the free surface of the fluid to the surface of the

Figure 1. Physical model and coordinate system.

porous layer (called the nominal surface) and

Zone 2 (porous layer) from the nominal surface to the impermeable lower plate.

The porous bed is saturated with an incompressible viscous fluid that is in local thermodynamic equilibrium with the solid matrix.

The equations governing the conservation of mass, momentum, energy and electric charge transfer [4] [6] [7] can be written in each Zone as follows

Zone 1 (fluid layer):

Equation governing the conservation of mass

$\nabla \cdot {{V}^{\prime}}_{f}=0,$ (1)

Equation governing the conservation of momentum (Navier-Stokes model with the presence of magnetic and gravitational fields).

${\rho}_{0}\left[\frac{\partial {{V}^{\prime}}_{f}}{\partial {t}^{\prime}}+\left({{V}^{\prime}}_{f}\cdot \nabla \right){{V}^{\prime}}_{f}\right]=-\nabla {p}_{f}+{\rho}_{f}g+{\mu}_{eff}{\nabla}^{2}\left({{V}^{\prime}}_{f}\right)+{j}_{f}\wedge B,$ (2)

Equation governing the conservation of energy

$\frac{\partial {{T}^{\prime}}_{f}}{\partial {t}^{\prime}}+\left({{V}^{\prime}}_{f}\cdot \nabla \right){{T}^{\prime}}_{f}={\alpha}_{f}{\nabla}^{2}\left({{T}^{\prime}}_{f}\right)+\frac{\mu}{{\left(\rho {C}_{p}\right)}_{f}}{\Phi}_{f},$ (3)

Equation governing the conservation of electric

$\nabla \cdot {j}_{f}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{j}_{f}=\sigma \left(-\nabla {\varnothing}_{f}+{{V}^{\prime}}_{f}\wedge B\right).$ (4)

Zone 2 (porous layer):

Equation governing the conservation of mass

$\nabla \cdot {{V}^{\prime}}_{m}=0,$ (5)

Equation governing the conservation of momentum (Brinkman-extended Darcy law with the presence of magnetic and gravitational fields).

${{V}^{\prime}}_{m}=\frac{\stackrel{\xaf}{\stackrel{\xaf}{K}}}{\mu}\left(-\nabla {{p}^{\prime}}_{m}+{\rho}_{m}g+{\mu}_{eff}{\nabla}^{2}\left({{V}^{\prime}}_{m}\right)+{j}_{m}\wedge B\right),$ (6)

Equation governing the conservation of energy

${\left(\rho {C}_{p}\right)}_{m}\frac{\partial {{T}^{\prime}}_{m}}{\partial {t}^{\prime}}+{\left(\rho {C}_{p}\right)}_{f}\nabla \left({{V}^{\prime}}_{m}\cdot {{T}^{\prime}}_{m}\right)=k{\nabla}^{2}\left({{T}^{\prime}}_{m}\right)+\mu {\Phi}_{m},$ (7)

Equation governing the conservation of electric

$\nabla \cdot {j}_{m}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{j}_{m}=\sigma \left(-\nabla {\varnothing}_{m}+{{V}^{\prime}}_{m}\wedge B\right).$ (8)

In these equations, $V$ , denotes the velocity vector, p' the pressure and T' the temperature of the fluid, $g$ the gravitational acceleration, t' the time, and

${\left(\rho {C}_{p}\right)}_{f}$ , ${\left(\rho {C}_{p}\right)}_{m}$ the heat capacity of the fluid and the saturated porous

medium, respectively. The subscript “f” denotes the fluid layer, “m” the porous medium. Moreover, μ the dynamic viscosity, ${\rho}^{\prime}={\rho}_{0}\left[1-\beta \left({T}^{\prime}-{T}_{0}\right)\right]$ the density, β the thermal-expansion coefficient, ${T}_{0}$ the constant reference Kelvin temperature, ${\rho}_{0}$ the density of the fluid at ${T}_{0}$ , ${C}_{p}$ the specific heat of the

fluid, k the thermal conductivity and $\alpha =k/{\left(\rho {C}_{p}\right)}_{f}$ the thermal diffusivity. In

Equation (3), $j$ is the electric current density, σ the electrical conductivity of the fluid, ∅ the electric potential and $-\nabla \varnothing $ the associated electric field. As discussed by [8] , for a two dimensional situation Equation (4) and (8) for the electric potential reduces to ${\nabla}^{2}\varnothing =0$ . The unique solution is $\nabla \varnothing =0$ . It follows that the electric field vanishes everywhere. The energy Equations (3) and (7) take into account the viscous dissipations ${\Phi}_{f}$ and ${\Phi}_{m}$ for the fluid layer and the porous medium respectively. It is important to mention here that the viscous dissipation ${\Phi}_{m}$ for the porous layer is the sum of viscous dissipation of the saturated fluid and the Darcy dissipation through the term ${{V}^{\prime}}_{m}^{2}/\stackrel{\xaf}{\stackrel{\xaf}{K}}$ . In Equation (6), ${\mu}_{eff}$ denotes the apparent dynamic viscosity for Brinkman model and the symmetrical second-order permeability tensor $\stackrel{\xaf}{\stackrel{\xaf}{K}}$ is defined as

$\stackrel{\xaf}{\stackrel{\xaf}{K}}=\left[\begin{array}{cc}{K}_{1}{\mathrm{sin}}^{2}\phi +{K}_{2}{\mathrm{cos}}^{2}\phi & \left({K}_{2}-{K}_{1}\right)\mathrm{sin}\phi \mathrm{cos}\phi \\ \left({K}_{2}-{K}_{1}\right)\mathrm{sin}\phi \mathrm{cos}\phi & {K}_{2}{\mathrm{sin}}^{2}\phi +{K}_{1}{\mathrm{cos}}^{2}\phi \end{array}\right]$ (9)

Assuming that when the flow is fully developed in the system, the axial (x'-direction) velocity depends on the transverse coordinate y' (i.e., ${{u}^{\prime}}_{f}={{u}^{\prime}}_{f}\left({y}^{\prime}\right)$ for the fluid layer and ${{u}^{\prime}}_{m}={{u}^{\prime}}_{m}\left({y}^{\prime}\right)$ for the porous layer), and then from the continuity equation, the transverse velocity component must be zero (i.e., ${{v}^{\prime}}_{f}=0$ and ${{v}^{\prime}}_{m}=0$ ). The temperature is assumed to be a function of y'. No assumptions are made with regard to the pressure variation (which, in fact, is found to be a function of x' and y' in Zone 2 [9] ). So, governing Equations (1)-(8) may be reduced as

Zone 1 (fluid layer):

$\frac{\text{d}{{u}^{\prime}}_{f}}{\text{d}{x}^{\prime}}=0$ (10)

$-\frac{1}{{\rho}_{0}}\frac{\partial {{p}^{\prime}}_{f}}{\partial {x}^{\prime}}+\nu \frac{{\text{d}}^{2}{{u}^{\prime}}_{f}}{\text{d}{{y}^{\prime}}^{2}}-\frac{\sigma {B}^{2}{{u}^{\prime}}_{f}}{{\rho}_{0}}=0,$ (11)

$-\frac{1}{{\rho}_{0}}\frac{\partial {{p}^{\prime}}_{f}}{\partial {y}^{\prime}}=0,$ (12)

${\alpha}_{f}\frac{{\text{d}}^{2}{{T}^{\prime}}_{f}}{\text{d}{{y}^{\prime}}^{2}}+\frac{\mu}{{\left(\rho {C}_{p}\right)}_{f}}{\left(\frac{\text{d}{{u}^{\prime}}_{f}}{\text{d}{y}^{\prime}}\right)}^{2}=0.$ (13)

Zone 2 (porous layer):

$\frac{\text{d}{{u}^{\prime}}_{m}}{\text{d}{x}^{\prime}}=0$ (14)

$b{{u}^{\prime}}_{m}=-\frac{{K}_{1}}{\mu}\frac{\partial {{p}^{\prime}}_{m}}{\partial {x}^{\prime}}+{K}_{1}\lambda \frac{{\text{d}}^{2}{{u}^{\prime}}_{m}}{\text{d}{{y}^{\prime}}^{2}}-\frac{{K}_{1}\sigma {B}^{2}{{u}^{\prime}}_{m}}{\mu},$ (15)

$c{{u}^{\prime}}_{m}=\frac{{K}_{1}}{\mu}\frac{\partial {{p}^{\prime}}_{m}}{\partial {y}^{\prime}},$ (16)

$k\frac{{\text{d}}^{2}{{T}^{\prime}}_{m}}{\text{d}{{y}^{\prime}}^{2}}=-\mu \left[{\left(\frac{\text{d}{{u}^{\prime}}_{m}}{\text{d}{y}^{\prime}}\right)}^{2}+\frac{b{\left({{u}^{\prime}}_{m}\right)}^{2}}{{K}_{1}}\right].$ (17)

where

$\{\begin{array}{l}b={\mathrm{sin}}^{2}\phi +{K}^{*}{\mathrm{cos}}^{2}\phi ,\\ c=\left(1-{K}^{*}\right)\mathrm{sin}\phi \mathrm{cos}\phi .\end{array}$ (18)

And $\lambda ={\mu}_{eff}/\mu $ the relative viscosity for which the value in the present study is taken, as a first approximation, equal to unity (i.e., ${\mu}_{eff}\approx \mu $ ).

The appropriate boundary conditions prevailing on the lower impermeable boundary and the upper free surface and at the interface of the two layers (y = h_{m}) of the channel are [4]

${y}^{\prime}=0:{{u}^{\prime}}_{m}=0,\text{\hspace{0.17em}}{{T}^{\prime}}_{m}={T}_{1},$ (19)

${y}^{\prime}=h:\frac{\text{d}{{u}^{\prime}}_{f}}{\text{d}{y}^{\prime}}=0,\text{\hspace{0.17em}}{{T}^{\prime}}_{f}={T}_{2},$ (20)

${y}^{\prime}={h}_{m}:{{u}^{\prime}}_{f}={{u}^{\prime}}_{m}={u}_{B}/\frac{\text{d}{{u}^{\prime}}_{f}}{\text{d}{y}^{\prime}}=\frac{{\beta}_{1}}{\sqrt{{K}_{1}}}\left({u}_{B}-{{{u}^{\prime}}_{m}|}_{y={h}_{m}-\sqrt{{K}_{1}}}\right),$ (21)

${y}^{\prime}={h}_{m}:{{T}^{\prime}}_{f}={{T}^{\prime}}_{m}={T}_{B}/\frac{\text{d}{{T}^{\prime}}_{f}}{\text{d}{y}^{\prime}}=\frac{{\beta}_{2}}{\sqrt{{K}_{1}}}\left({T}_{B}-{{{T}^{\prime}}_{m}|}_{y={h}_{m}-\sqrt{{K}_{1}}}\right).$ (22)

where ${u}_{B}$ is the slip velocity at the nominal surface (interface) which changes to the constant Darcy velocity through the porous, the existence of the temperature slip layer whose thickness given by Beavers and Joseph is of order $\sqrt{{K}_{1}}$ . Indeed, this thickness subsequently has been shown to be equal to $\sqrt{{K}_{1}}$ [10] [11] . According to these authors, the parameter ${\beta}_{1}$ (and ${\beta}_{2}$ ) denotes a constant depending on the material property of the porous medium, which have can be determined only experimentally. The slip temperature ${T}_{B}$ at the nominal surface changes to the constant ambient temperature ${T}_{0}$ through the porous layer whose thickness considered to be the same as that of the velocity slip layer namely $\sqrt{{K}_{1}}$ . So, the existence of the temperature slip layer depends on that of the velocity slip layer and in view of this, ${T}_{0}$ is identified with ${{T}^{\prime}}_{m}$ at a distance $\sqrt{{K}_{1}}$ below the nominal.

When the following variables are introduced

$\begin{array}{l}x={x}^{\prime}/\left(hRe\right),y={y}^{\prime}/h,\\ \left({u}_{f},{u}_{m}\right)=\left({{u}^{\prime}}_{f},{{u}^{\prime}}_{m}\right)/\stackrel{\xaf}{u},\text{\hspace{0.17em}}\left({T}_{f},{T}_{m}\right)=\left[\left({{T}^{\prime}}_{f},{{T}^{\prime}}_{m}\right)-{T}_{1}\right]/\Delta T,\\ \Delta T={T}_{2}-{T}_{1},\left({P}_{f},{P}_{m}\right)=\left({{p}^{\prime}}_{f},{{p}^{\prime}}_{m}\right)/\left(\rho {\stackrel{\xaf}{u}}^{2}\right).\end{array}\}$ (23)

the governing Equations (10)-(13) corresponding to the fluid layer may be written in non-dimensional form as

$\frac{{\text{d}}^{2}{u}_{f}}{\text{d}{y}^{2}}-H{a}^{2}{u}_{f}=\frac{\text{d}{P}_{f}}{\text{d}x},$ (24)

$\frac{{\text{d}}^{2}{T}_{f}}{\text{d}{y}^{2}}=-PrEc{\left(\frac{\text{d}{u}_{f}}{dy}\right)}^{2}.$ (25)

and the governing Equations (14)-(17) corresponding to the anisotropic porous layer may be written in non-dimensional form as

$\frac{{\text{d}}^{2}{u}_{m}}{\text{d}{y}^{2}}-\frac{b+H{a}_{m}^{2}}{Da}{u}_{m}=\frac{\partial {P}_{m}}{\partial x},$ (26)

$\frac{\partial {P}_{m}}{\partial y}=\frac{c}{ReDa}{u}_{m},$ (27)

$\frac{{\text{d}}^{2}{T}_{m}}{\text{d}{y}^{2}}=-PrEc\left[{\left(\frac{\text{d}{u}_{m}}{\text{d}y}\right)}^{2}+{\gamma}^{2}{\left({u}_{m}\right)}^{2}\right].$ (28)

where

${\gamma}^{2}=\frac{b}{Da}.$ (29)

In the above equations, $Pr=\mu {C}_{p}/k$ is the Prandtl number, $Ec={\stackrel{\xaf}{u}}^{2}/\left({C}_{p}\Delta T\right)$ the Eckert number, $Re=\stackrel{\xaf}{u}h/\nu $ the Reynolds number, $Da={K}_{1}/{h}^{2}$ the Darcy number, $Ha=B\cdot h\sqrt{\sigma /\mu}$ the Hartmann number for the fluid layer and $H{a}_{m}=B\sqrt{{K}_{1}\sigma /\mu}$ the Hartmann number for the anisotropic porous layer.

The boundary conditions associated with the non-dimensional Equations (19)-(22) are [4]

$y=0:{u}_{m}=0,{T}_{m}=0,$ (30)

$y=1:\frac{\text{d}{u}_{f}}{\text{d}y}=0,{T}_{f}=1,$ (31)

$y=\epsilon :{u}_{f}={u}_{m}={u}_{B}/\frac{\text{d}{u}_{f}}{\text{d}y}=\frac{{\beta}_{1}}{\sqrt{Da}}\left({u}_{B}-{{u}_{m}|}_{y=\epsilon -\sqrt{Da}}\right),$ (32)

$y=\epsilon :{T}_{f}={T}_{m}={T}_{B}/\frac{\text{d}{T}_{f}}{\text{d}y}=\frac{{\beta}_{2}}{\sqrt{{K}_{1}}}\left({T}_{B}-{{T}_{m}|}_{y=\epsilon -\sqrt{Da}}\right).$ (33)

where $\epsilon ={h}_{m}/h$ is the thickness ratio defining as the ratio of saturated porous layer thickness to packed porous beds thickness (which consists of a horizontal fluid layer (river bed) and a porous zone).

Using the hydrodynamic conditions (30) and (32), the resolution of Equation (24) yields the velocity distribution in Zone 1 (for the fluid layer) expressed as

${u}_{f}\left(y\right)={A}_{1}{\text{e}}^{Ha\cdot y}+{A}_{2}{\text{e}}^{-Ha\cdot y}-\frac{Q}{H{a}^{2}}$ (34)

where $Q=\text{d}{P}_{f}/\text{d}x$ .

Taking into account Equation (34) and considering the boundary conditions for ${T}_{f}$ , Equation (25) can be integrated to give the following fully developed temperature profile in fluid Zone:

${T}_{f}\left(y\right)=-PrEc\left[\frac{{A}_{1}^{2}}{4}{\text{e}}^{2\cdot Ha\cdot y}+\frac{{A}_{2}^{2}}{4}{\text{e}}^{-2\cdot Ha\cdot y}-{A}_{1}{A}_{2}H{a}^{2}y\right]+{A}_{3}\cdot y+{A}_{4},$ (35)

where

$\begin{array}{l}{A}_{1}=\frac{{U}_{B}+Q/H{a}^{2}}{\left[{\text{e}}^{Ha\cdot \epsilon}+{\text{e}}^{Ha\left(2-\epsilon \right)}\right]},\\ {A}_{2}={A}_{1}{\text{e}}^{2\cdot Ha},\\ {A}_{3}=\frac{{T}_{B}-1}{\epsilon -1}+\frac{PrEc}{4\left(\epsilon -1\right)}\left[\left({\text{e}}^{2\cdot Ha\cdot \epsilon}+{\text{e}}^{2\cdot Ha\left(2-\epsilon \right)}\right)+2\left[2H{a}^{2}\left(1-{\epsilon}^{2}\right)-1\right]{\text{e}}^{2\cdot Ha}\right]{A}_{1}^{2},\\ {A}_{4}=\frac{\epsilon -{T}_{B}}{\epsilon -1}+\frac{PrEc}{4\left(\epsilon -1\right)}\left[2\cdot \epsilon \cdot {\text{e}}^{2\cdot Ha}\left(1-2H{a}^{2}\left(1-\epsilon \right)\right)-\left({\text{e}}^{2\cdot Ha\cdot \epsilon}+{\text{e}}^{2\cdot Ha\left(2-\epsilon \right)}\right)\right]{A}_{1}^{2}.\end{array}\}$ (36)

where ${u}_{B}$ is the velocity profile at the interface that must be determined by making use of the condition (32) in which we have to know first the velocity ${u}_{m}$ for the porous layer. It is seen that ${T}_{f}$ depends on ${T}_{B}$ that has to be calculated from the boundary condition (33) for which the temperature profile ${T}_{m}$ (in Zone 2) should be expressed.

Now, the velocity and temperature profiles ${u}_{m}$ and ${T}_{m}$ for the porous layer will be determined.

As $H{a}_{m}^{2}=Da\cdot H{a}^{2}$ , Equation(26) may be written as

$\frac{{\text{d}}^{2}{u}_{m}}{\text{d}{y}^{2}}-{\xi}^{2}{u}_{m}=\frac{\partial {P}_{m}}{\partial x},$ (37)

where ${\xi}^{2}={\gamma}^{2}+H{a}^{2}$ .

Eliminating the pressure from Equation (27) and (37) in the usual way, one may have

$\frac{{\text{d}}^{3}{u}_{m}}{\text{d}{y}^{3}}-{\xi}^{2}\frac{\text{d}{u}_{m}}{dy}=0.$ (38)

It is clear that the bulk mean velocity ${\stackrel{\xaf}{u}}^{\prime}$ (in Zone 2) defined as

${\stackrel{\xaf}{u}}^{\prime}=\left(1/\left[{\rho}^{\prime}\left({h}_{m}\times 1\right)\right]\right){\displaystyle {\int}_{0}^{{h}_{m}}{\rho}^{\prime}{{u}^{\prime}}_{m}\text{d}{y}^{\prime}}$ is calculated in dimensionless terms by

$\stackrel{\xaf}{u}={\displaystyle {\int}_{0}^{\epsilon}{u}_{m}\text{d}y}$ (39)

which equals to ε.

Making use of Equation (39) and the hydrodynamic conditions (31) and (32), the resolution of Equation (38) yields the velocity distribution (for the porous layer) expressed as

${u}_{m}\left(y\right)=\frac{1}{\xi}\left({B}_{1}{\text{e}}^{\xi y}-{B}_{2}{\text{e}}^{-\xi y}\right)+{B}_{3}$ (40)

where

$\begin{array}{l}{B}_{1}=\frac{{\xi}^{2}\epsilon \left(1-{\text{e}}^{-\xi \epsilon}\right)-{u}_{B}\left({\text{e}}^{-\xi \epsilon}+\xi \epsilon -1\right)\xi}{4\cdot \mathrm{cosh}\left(\xi \epsilon \right)-2\cdot \xi \cdot \epsilon \cdot \mathrm{sinh}\left(\xi \epsilon \right)-4}\\ {B}_{2}=\frac{{B}_{1}\left(1-{\text{e}}^{\xi \epsilon}\right)+\xi {u}_{B}}{1-{\text{e}}^{-\xi \epsilon}}\\ {B}_{3}=\frac{{B}_{2}-{B}_{1}}{\xi}\end{array}\}$ (41)

By substituting Equation (40) into Equation (28) and using the thermal boundary conditions (32) and (33), the temperature profile for the porous medium is obtained as

$\begin{array}{c}{T}_{m}\left(y\right)=-PrEc[\frac{1}{4\cdot {\xi}^{2}}\left(1+\frac{{\gamma}^{2}}{{\xi}^{2}}\right)\left({B}_{1}^{2}{\text{e}}^{2\xi y}+{B}_{2}^{2}{\text{e}}^{-2\xi y}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{2{B}_{1}{\gamma}^{2}}{{\xi}^{3}}\left({B}_{1}{\text{e}}^{\xi y}-{B}_{2}{\text{e}}^{-\xi y}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\left(2\left(1-\frac{{\gamma}^{2}}{{\xi}^{2}}\right){B}_{1}{B}_{2}+{\gamma}^{2}{B}_{3}^{2}\right)\frac{{y}^{2}}{2}]+{B}_{4}\cdot y+{B}_{5}.\end{array}$ (42)

where

$\begin{array}{l}{B}_{4}=\frac{{T}_{B}-{B}_{5}}{\epsilon}+\frac{PrEc}{\epsilon}[\frac{1}{4\cdot {\xi}^{2}}\left(1+\frac{{\gamma}^{2}}{{\xi}^{2}}\right)\left({B}_{1}^{2}{\text{e}}^{2\xi \epsilon}+{B}_{2}^{2}{\text{e}}^{-2\xi \epsilon}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{2{B}_{3}{\gamma}^{2}}{{\xi}^{3}}\left({B}_{1}{\text{e}}^{\xi \epsilon}-{B}_{2}{\text{e}}^{-\xi \epsilon}\right)+\left(2\left(1-\frac{{\gamma}^{2}}{{\xi}^{2}}\right){B}_{1}{B}_{2}+{\gamma}^{2}{B}_{3}{}^{2}\right)\frac{{\epsilon}^{2}}{2}],\\ {B}_{5}=PrEc\left[\frac{1}{4\cdot {\xi}^{2}}\left(1+\frac{{\gamma}^{2}}{{\xi}^{2}}\right)\left({B}_{1}^{2}+{B}_{2}^{2}\right)+\frac{2{B}_{3}{\gamma}^{2}}{{\xi}^{3}}\left({B}_{1}-{B}_{2}\right)\right]\end{array}$ (43)

It is noticed that both the velocity profiles ${u}_{f}$ and ${u}_{m}$ depend on ${u}_{B}$ . the same holds good for both the temperature profiles ${T}_{f}$ and ${T}_{m}$ which depend on ${T}_{B}$ . The profiles ${u}_{B}$ and ${T}_{B}$ at the interface are determined by using boundary conditions (32) and (33) respectively. After manipulations, one can obtain

${u}_{B}=\frac{{a}_{5}\cdot Ha\cdot {\beta}_{1}-Q\cdot {a}_{7}\cdot \sqrt{Da}}{Ha\left(Ha\cdot {a}_{7}\sqrt{Da}-{\beta}_{1}\left(1+{a}_{6}\right)\right)}$ (44)

where

$\begin{array}{l}{a}_{5}={a}_{1}\frac{1-{\text{e}}^{\xi \left(\epsilon -\sqrt{Da}\right)}}{\epsilon}-{a}_{3}\frac{1-{\text{e}}^{-\xi \left(\epsilon -\sqrt{Da}\right)}}{\epsilon}\\ {a}_{6}={a}_{2}\frac{1-{\text{e}}^{\xi \left(\epsilon -\sqrt{Da}\right)}}{\epsilon}-{a}_{4}\frac{1-{\text{e}}^{-\xi \left(\epsilon -\sqrt{Da}\right)}}{\epsilon}\\ {a}_{7}=\frac{{\text{e}}^{Ha\cdot \epsilon}-{\text{e}}^{Ha\left(2-\epsilon \right)}}{{\text{e}}^{Ha\cdot \epsilon}+{\text{e}}^{Ha\left(2-\epsilon \right)}}\\ {a}_{1}=\frac{{\xi}^{2}\epsilon \left(1-{\text{e}}^{-\xi \epsilon}\right)}{{a}_{0}}\\ {a}_{0}=4\cdot \mathrm{cosh}\left(\xi \epsilon \right)-2\cdot \xi \cdot \epsilon \cdot \mathrm{sinh}\left(\xi \epsilon \right)-4\\ {a}_{2}=\frac{1-\xi \epsilon -{\text{e}}^{-\xi \epsilon}}{{a}_{0}}\\ {a}_{3}=\frac{{a}_{1}\left(1-{\text{e}}^{\xi \epsilon}\right)}{1-{\text{e}}^{-\xi \epsilon}}\\ {a}_{4}=\frac{\xi +{a}_{2}\left(1-{\text{e}}^{\xi \epsilon}\right)}{1-{\text{e}}^{-\xi \epsilon}}\end{array}\}$ (45)

Similarly, after development and reduction of algebraic terms, the characteristic temperature distribution ${T}_{B}$ at the interface is calculated by

${T}_{B}=\frac{\left[{b}_{0}\sqrt{Da}+\left({b}_{1}+{b}_{2}\right){\beta}_{2}\right]\epsilon \left(1-\epsilon \right)}{{\beta}_{2}\cdot \epsilon \left(1-\epsilon \right)+\left(\sqrt{Da}-\epsilon \right){\beta}_{2}\left(1-\epsilon \right)+\sqrt{Da}}$ (46)

where

$\begin{array}{l}{b}_{0}=\frac{PrEc}{4\left(\epsilon -1\right)}\left[\left(1-2\cdot Ha\left(\epsilon -1\right)\right){\text{e}}^{2\cdot Ha\cdot \epsilon}+\left(1+2\cdot Ha\left(\epsilon -1\right)\right){\text{e}}^{2\cdot Ha\left(2-\epsilon \right)}+{b}_{3}\right]{A}_{1}^{2}+\frac{1}{\epsilon -1}\\ {b}_{1}=-PrEc\left[\frac{1}{4\cdot {\xi}^{2}}\left(1+\frac{{\gamma}^{2}}{{\xi}^{2}}\right)\left({B}_{1}^{2}{\text{e}}^{2\xi \left(\epsilon -\sqrt{Da}\right)}+{B}_{2}^{2}{\text{e}}^{-2\xi \left(\epsilon -\sqrt{Da}\right)}\right)+{b}_{3}\right]\\ {b}_{2}=\frac{\epsilon -\sqrt{Da}}{\epsilon}\left[PrEc\left(\frac{1}{4\cdot {\xi}^{2}}\left(1+\frac{{\gamma}^{2}}{{\xi}^{2}}\right)\left({B}_{1}^{2}{\text{e}}^{2\xi \epsilon}+{B}_{2}^{2}{\text{e}}^{-2\xi \epsilon}\right)+{b}_{5}\right)-{B}_{5}\right]\\ {b}_{3}=\left[8\cdot \epsilon H{a}^{2}\left(\epsilon -1\right)+2\cdot \left(2\cdot H{a}^{2}\left[1-{\epsilon}^{2}\right]-1\right)\right]{\text{e}}^{2\cdot Ha}\\ {b}_{4}=\frac{2{B}_{3}{\gamma}^{2}}{{\xi}^{3}}\left({B}_{1}{\text{e}}^{\xi \left(\epsilon -\sqrt{Da}\right)}-{B}_{2}{\text{e}}^{-\xi \left(\epsilon -\sqrt{Da}\right)}\right)+\left(2\left(1-\frac{{\gamma}^{2}}{{\xi}^{2}}\right){B}_{1}{B}_{2}+{\gamma}^{2}{B}_{3}^{2}\right)\frac{{\left(\epsilon -\sqrt{Da}\right)}^{2}}{2}\\ {b}_{5}=\frac{2{B}_{3}{\gamma}^{2}}{{\xi}^{3}}\left({B}_{1}{\text{e}}^{\xi \epsilon}-{B}_{2}{\text{e}}^{-\xi \epsilon}\right)+\left(2\left(1-\frac{{\gamma}^{2}}{{\xi}^{2}}\right){B}_{1}{B}_{2}+{\gamma}^{2}{B}_{3}^{2}\right)\frac{{\epsilon}^{2}}{2}\end{array}\}$ (47)

The heat transfer rates through the upper free surface and the lower wall bounding the porous layer underlining the fluid zone are expressed by the Nusselt numbers given in dimensionless terms respectively by

$N{u}_{u}={\frac{\text{d}{T}_{f}}{\text{d}y}|}_{y=1}=2\cdot PrEc\cdot H{a}^{2}{A}_{1}^{2}{\text{e}}^{2\cdot Ha}+{A}_{3}$ (48)

and

$N{u}_{\mathcal{l}}={\frac{\text{d}{T}_{m}}{\text{d}y}|}_{y=0}=-PrEc\left(1-\frac{{\gamma}^{2}}{{\xi}^{2}}\right)\frac{{B}_{1}-{B}_{2}}{2\cdot \xi}+\frac{2{B}_{3}{\gamma}^{2}}{{\xi}^{3}}\left({B}_{1}+{B}_{2}\right)+{B}_{4}$ (49)

Consequently, one can deduce that the mean heat transfer rate $N{u}_{moy}$ through the superposed layers is calculated by the following expression

3. Results and Discussion

3.1. Horizontal Velocity Distribution and Temperature Distribution for the Entire System Width

Figure 2 illustrates the analytical results obtained for the Horizontal velocity distribution for the entire system width when Da = 0.05, Q = −2, β_{1} = 0.1, ε = 0.2, K^{*} = 0.1, φ = 0˚ and various values of Ha. From this figure it is clearly seen that the velocity distribution decreases as the Hartmann number increases. For K^{*} = 0.1 (When the permeability in the horizontal direction, K_{2}, is higher than the permeability in the vertical direction, K_{1}.) the convective flow is high when Ha = 0.8 (and for the lowest values of Ha) and is reduced when Ha = 5 (and for the highest values of Ha). It is noticed that the velocity fields have the same behavior as described previously [4] .

Figure 3 illustrates the analytical results obtained for Temperature distribution for the entire system width when K^{*} = 0.1, β_{1} = β_{2} = 0.1, Da = 005, Q = −2, ε = 0.4, Pr = 7, Ec = 0.05, φ = 0˚. It is noticed that the temperature fields have the same behavior as described previously for the velocity distribution, revealing that the effect of varying convection in the two superposed layers depends strongly on the magnetic field, i.e., Ha.

Figure 2. Horizontal velocity distribution for the entire system width.

Figure 3. Horizontal temperature distribution for the entire system width.

3.2. Velocity Profile, u_{B}, at the Interface

Figure 4 show the effects of various values of the Hartmann number, Ha, of varying the anisotropic ratio, K^{*}, and the thickness ratio, ε, on the velocity profile u_{B} at the interface when φ = 10˚, β_{1} = β_{2} = 0.1, Pr = 7, Q = −3, Ec = 5 ´ 10^{−3}, Da = 10^{−3}. Because of the fact that
${u}_{B}$ is not physically defined at ε = 0 (corresponding to the situation when the two layers are assimilated into a pure fluid layer) and at ε = 1 (case for which the system would become a pure porous layer), all curves observed tend asymptotically to the pure fluid case when
$\epsilon \to 0$ . As expected, in this limit, the effects of the anisotropic properties of the almost negligible porous layer are irrelevant. For intermediate values of the thickness ratio, it is seen that, as the porous layer becomes more and more important, the anisotropic effects of the porous medium become more significant (please refer [4] ). It is also noticed that each curve for a given value of K^{*} is distinguishable from others at a special point
${\epsilon}_{c}$ that depend on the Hartmann number and below which the effects of anisotropy and magnetic field is irrelevant. For example,
${\epsilon}_{c}=0.0744$ when Ha = 1,
${\epsilon}_{c}=0.0731$ when Ha = 0.5 and
${\epsilon}_{c}=0.0722$ when Ha = 0.005. For each value of K^{*}, Figure 4 indicates that the velocity profile u_{B} at the interface decreases as the Hartmann number increases.

3.3. The Mean Heat Transfer Rate, Nu_{moy}

The variations of Nusselt number with the Hartmann number and the anisotropic ratio is shown in Figure 5 as a function of the thickness ratio, ε, when Da = 75 × 10^{−3}, φ = 0˚, Pr = 7.0, Ec = 0.05, β_{1} = β_{2} = 0.001, and Q = 50. When the porous lining is negligible with respect to the fluid layer (i.e., when the thickness ratio is small), the heat transfer given by Equation (50), does not vary with the anisotropic ratio, K^{*}. Then the Nusselt number decreases with the increase of the thickness ratio and the effect of anisotropy is predominant. Figure 5 indicates

Figure 4. Velocity profile, u_{B}, at the interface.

Figure 5. Effect of the Hartmann number on Nusselt number, Nu_{moy}.

that the heat transfer through the entire system is enhanced when K^{*} = 2.5 with respect to the isotropic case for which K^{*} = 1.0. That is damaging to the aquatic species life. We also note that the increase of the Hartmann number corresponds to the decrease of the heat transfer when the porous layer becomes more and more important. What encourages the aquatic species life.

In Figure 6, we have plotted the average Nusselt number as a function of the thickness ratio, ε, with the variations of Hartmann number and the anisotropic orientation φ for Da = 75 × 10^{−3}, K* = 0.1, Pr = 7.0, Ec = 0.05, β_{1} = β_{2} = 0.001, and Q = 50. It is noticed that the Nusselt number have the same behavior as described previously, revealing that the heat transfer through the entire system for

Figure 6. Effect of the Hartmann number and φ on Nusselt number, Nu_{moy}.

a given value of Hartmann number is found to be maximum when φ = 90˚ and minimum when φ = 0˚.

The average Nusselt number is presented in Figure 7 as a function of the Darcy number, of the Hartmann number and of varying the anisotropic ratio when φ = 0˚, Pr = 7, Ec = 0.05, β_{1} = β_{2} = 10^{−6}, ε = 0.25 and Q = 5. From this figure it is clearly seen that the average Nusselt number decreases as the Hartmann number and the Darcy number increase, the highest decrease of the average Nusselt number occurs for the region of the low values of the Darcy number (i.e., in the porous zone where the effect of anisotropy is predominant). For high values of the Darcy number (case for which the system would become a pure fluid layer), no significant decrease occurs on the average Nusselt number which is not affected by the anisotropy ratio.

In Figure 8, the average Nusselt number is plotted as a function of the Darcy number for K^{*} = 0.5, Pr = 7, Ec = 0.05, β_{1} = β_{2} = 10^{−6}, ε = 0.25 and Q = 5 and various values of the Hartmann number and of φ. It is noticed that the average Nusselt number have the same behavior as described previously in Figure 7 and is found to be maximum when φ = 90˚ and minimum when φ = 0˚.

4. Conclusion

In this paper we have examined the effect of magnetic field on convection heat transfer through packed porous beds which consists of a horizontal fluid layer (river bed) and a porous zone with anisotropic permeability. The results obtained show that, the application of a magnetic field on the natural convection in the y-direction has some important changes on the structure of the convective flow and on convection heat transfer. It generates a reduction in the heat transfer and reduces considerably the velocity and temperature distribution. We have also examined the effect of anisotropic permeability ratio and anisotropic angle

Figure 7. Effect of the Hartmann number and Da on Nusselt number, Nu_{moy}.

Figure 8. Effect of the Hartmann number and φ on Nusselt number, Nu_{moy}.

on the velocity profile, u_{B}, at the interface and on the average Nusselt number. The heat transfer, for a given value of the Hartmann number and for φ = 0˚, is found to be maximum (minimum) when the permeability in the y-direction (x-direction) is bigger than the permeability in the x-direction (y-direction) and for a given value of the anisotropic permeability ratio (K^{*} = 0.5), the average Nusselt number is found to be maximum when φ = 90˚ and minimum when φ = 0˚.

Cite this paper

Yovogan, J. , Degan, G. and Fagbemi, L. (2018) Effect of Constant Magnetic Field on Convective Heat Transfer through Anisotropic River Beds.*Journal of Crystallization Process and Technology*, **8**, 57-71. doi: 10.4236/jcpt.2018.82004.

Yovogan, J. , Degan, G. and Fagbemi, L. (2018) Effect of Constant Magnetic Field on Convective Heat Transfer through Anisotropic River Beds.

References

[1] Ni, J., Beckerman, C. and Smith, T.F. (1993) Effect of an Electromagnetic Field on Natural Convection in Porous Medium. Fundam Heat Transf Electromagn Electrost Acoust Field ASME HTD, 248, 23-33.

[2] Alchaar, S., Vasseur, P. and Bilgen, E. (1995) Effects of a Magnetic Field on the Onset of Convection in a Porous Medium. Heat and Mass Transfer, 30, 259-267.

https://doi.org/10.1007/BF01602772

[3] Bian, W., Vasseur, P. and Bilgen, E. (1996) Effect of an External Magnetic Field on Buoyancy-Driven Flow in a Shallow Porous Cavity. Numerical Heat Transfer, Part A, 29, 625-638.

https://doi.org/10.1080/10407789608913810

[4] Yovogan, J. and Degan, G. (2013) Effect of Anisotropic Permeability on Convective Heat Transfer through a Porous River Bed Underlying a Fluid Layer. Journal of Engineering Mathematics, 81, 127-140.

https://doi.org/10.1007/s10665-012-9605-6

[5] Beavers, G.S. and Joseph, D.D. (1974) Boundary Conditions at a Naturally Permeable Wall. Journal of Heat Transfer, 96, 324-330.

[6] Farid, B., Benissad, S., Talbi, K. and Bouttout, A. (2007) Effet d’un champ magnétique externe sur la convection naturelle dans une enceinte rectangulaire tridimensionnelle. Jean-Jacques BEZIAN. JITH, Albi, France. ENSTIMAC, 5 p,.

[7] Banjer, H.M. and Abdullah, A.A. (2010) Convection in Superposed Fluid and Porous Layers in the Presence of a Vertical Magnetic Field. WSEAS Transactions on Fluid Mechanics, 5, 175-185.

[8] Garandet, J.P., Alboussiere, T. and Moreau, R. (1992) Buoyancy-Driven Convection in a Rectangular Enclosure with a Transverse Magnetic Field. International Journal of Heat and Mass Transfer, 35, 741-748.

https://doi.org/10.1016/0017-9310(92)90242-K

[9] Degan, G., Zohoun, S. and Vasseur, P. (2002) Forced Convection in Horizontal Porous Channels with Hydrodynamic Anisotropy. International Journal of Heat and Mass Transfer, 45, 3181-3188.

https://doi.org/10.1016/S0017-9310(02)00032-7

[10] Rudraiah, N. and Veerabhadraiah, R. (1977) Temperature Distribution in Couette Flow past a Permeable Bed. Proceedings of the Indian Academy of Science, Section A, 86, 537.

[11] Rudraiah, N. and Veerabhadraiah, R. (1978) Effect of Buoyancy on the Free Surface Flow past a Permeable Bed. Warme- und Stoffübertragung, 11, 265.

https://doi.org/10.1007/BF02587790

[1] Ni, J., Beckerman, C. and Smith, T.F. (1993) Effect of an Electromagnetic Field on Natural Convection in Porous Medium. Fundam Heat Transf Electromagn Electrost Acoust Field ASME HTD, 248, 23-33.

[2] Alchaar, S., Vasseur, P. and Bilgen, E. (1995) Effects of a Magnetic Field on the Onset of Convection in a Porous Medium. Heat and Mass Transfer, 30, 259-267.

https://doi.org/10.1007/BF01602772

[3] Bian, W., Vasseur, P. and Bilgen, E. (1996) Effect of an External Magnetic Field on Buoyancy-Driven Flow in a Shallow Porous Cavity. Numerical Heat Transfer, Part A, 29, 625-638.

https://doi.org/10.1080/10407789608913810

[4] Yovogan, J. and Degan, G. (2013) Effect of Anisotropic Permeability on Convective Heat Transfer through a Porous River Bed Underlying a Fluid Layer. Journal of Engineering Mathematics, 81, 127-140.

https://doi.org/10.1007/s10665-012-9605-6

[5] Beavers, G.S. and Joseph, D.D. (1974) Boundary Conditions at a Naturally Permeable Wall. Journal of Heat Transfer, 96, 324-330.

[6] Farid, B., Benissad, S., Talbi, K. and Bouttout, A. (2007) Effet d’un champ magnétique externe sur la convection naturelle dans une enceinte rectangulaire tridimensionnelle. Jean-Jacques BEZIAN. JITH, Albi, France. ENSTIMAC, 5 p,

[7] Banjer, H.M. and Abdullah, A.A. (2010) Convection in Superposed Fluid and Porous Layers in the Presence of a Vertical Magnetic Field. WSEAS Transactions on Fluid Mechanics, 5, 175-185.

[8] Garandet, J.P., Alboussiere, T. and Moreau, R. (1992) Buoyancy-Driven Convection in a Rectangular Enclosure with a Transverse Magnetic Field. International Journal of Heat and Mass Transfer, 35, 741-748.

https://doi.org/10.1016/0017-9310(92)90242-K

[9] Degan, G., Zohoun, S. and Vasseur, P. (2002) Forced Convection in Horizontal Porous Channels with Hydrodynamic Anisotropy. International Journal of Heat and Mass Transfer, 45, 3181-3188.

https://doi.org/10.1016/S0017-9310(02)00032-7

[10] Rudraiah, N. and Veerabhadraiah, R. (1977) Temperature Distribution in Couette Flow past a Permeable Bed. Proceedings of the Indian Academy of Science, Section A, 86, 537.

[11] Rudraiah, N. and Veerabhadraiah, R. (1978) Effect of Buoyancy on the Free Surface Flow past a Permeable Bed. Warme- und Stoffübertragung, 11, 265.

https://doi.org/10.1007/BF02587790