^{*}

^{*}

In this paper, a detailed numerical investigation of two-dimensional laminar forced convection in a porous channel with inlet and outlet slot is presented. A uniform heat flux is applied on one wall of channel and another wall is isolated. The flow in the porous medium is modeled using the Brinkman-Forchheimer-extended Darcy model in which the inertia and boundary effects are taken into consideration and thermal dispersion effects are not included in the energy equation. Parametric studies are conducted to evaluate the effects of particle diameter, Reynolds number on the heat transfer and friction factor. Nusselt number and friction factor are developed for efficient design of a porous heat exchanger based on the present configuration. In order to solve this problem Fluent software was used.

Recent advances in high-technology fields such as electronic system have lead to dramatic increases in heat flux levels, for thermal management of such systems, more efficient heat transfer enhancement techniques are being explored. There has been a growing interest in heat transfer enhancement using a porous medium and several studies have reported that the use of porous medium for heat transfer yields higher heat transfer performance than existing techniques such an pinned-fin array or twisted tape inserts. Hunt and Tien [

In the present study, a porous parallel-plate channel containing a fluid-saturated sintered porous metal matrix that a uniform heat flux is applied on one of the walls and another wall of channel is isolated. To perform a detailed study of the hydrodynamic and heat transfer characteristics within porous channel, the numerical model for a steady-state, incompressible, and two-dimensional flow in metal porous channel is developed. The flow in the porous medium is modeled using the brinkman-forchheimerextended Darcy model that takes into account the boundary and inertia effects, thermal dispersion are neglected in the energy equation. Channeling effects near the walls are neglected in the present model and porosity is assumed to be constant. The numerical model is validated by comparing the present numerical results with the experimental data. Geometrical configuration of channel is shown in

Cf = inertia coefficient

h = heat transfer coefficient

K = permeability of porous medium

Dp = particle diameter

Nu = nusselt number

k = thermal conductivity

Dh = hydraulic diameter

F = friction factor

T = temperature

P = pressure

q” = heat flux

u = fluid velocity along x-direction

v = transverse fluid velocity along the duct height

Greek symbols

µ = dynamics viscosity of the fluid

= fluid density

= porosity

= performance parameter

Subscripts

f = fluid

s = solid

w = wall

o = inlet

e = effective

A schematic configuration of the physical model of porous channel is illustrated in

porous medium is composed of sintered metal particles such as aluminum spherical beads. Steady state, incompressible, and two-dimensional flow in a channel filled with a homogeneous and isotropic sintered porous metal and heated from the bottom wall is considered. Channeling effects near the walls are neglected in the present model and porosity is assumed to be constant. The thermo physical properties of fluid (water) and of the porous medium are assumed to be constant. Based on the physical model discussed above, the governing equation can be written as:

Continuity equation: considering a laminar, incompressible Newtonian fluid flow, the conservation equation of mass is derived as:

Momentum equation: the microscopic momentum equation are given by the navier-stokes equation:

Similarly, integrating equation (2) with respect to a REV (representation volume) and the no-slip condition on the solid surface gives:

For steady-state situations, dividing by on both sides of Equation (3) leads to:

(5)

That in the momentum equation Darcy term, brinkman term and forchheimer term is calculated in respect with

and

The permeability K and the inertia coefficient C_{f} are evaluated with using following empirical correlations:

And

Because of local equilibrium between fluid and solid phases in porous media, we have:

Ke is the effective thermal conductivity.

The boundary conditions for the configuration considered are as follows: no slip condition occurs on all the solid surface (top and bottom walls), the velocity v and temperature T of flow are uniform over the cross section of the inlet port. The dimensionless averaged pressure difference between the inlet and outlet ports is calculated in terms of the friction factor which is defined as:

The local Nusselt number for the heated boundary is defined as:

The governing equation (1-5) with the associated boundary condition was solved numerically utilizing a control volume method. The momentum equation for velocity and pressure drop was solved first and the energy equation for temperature and heat transfer rate was then solved based on the corrected velocity results. After a staggered grid was arranged for the discretization process a second—order upwind was employed to model the interaction between the convection and the diffusion terms in governing equations. An iterative procedure based on the Simple algorithm was employed to deal with the linkage between pressure and velocity in the momentum equations, and temperature in the energy equation. In order to solve this problem with all considered boundary conditions was used. But because doesn’t include equations related to terms of brinkman and forchheimer ,we add up these equations using udf (user defined function) to fluent but before solving this problem using fluent, geometrical configuration are meshed with gambit software. Gambit is the preprocessing software for fluent used for creating and meshing 2-d model. The units are not important while creating the grid in gambit, but the simulation in fluent always requires the input parameters to be in S.I units. For 2-d geometers, the z-coordinate is assigned a default value of zero. Both structured and unstructured meshes are available in gambit. in the present study, structured meshing was employed. Boundary conditions types for inlet, outlet and channel walls were specified and the mesh-file is saved in a binary format for fluent to recognize. In

To investigate the effects of particle diameter on heat transfer and pressure drop, the particle diameter Dp is selected. It is shown in Figures 3 and 4 that in the relatively small particle diameter region (Darcy regime), the length-averaged Nusselt number (Nu) and the pressure drop (f_{Re}_{ }) decrease sharply with increasing . For a large particle diameter (non-Darcy regime), the nusselt number and the pressure drop both decrease slowly as they approach their corresponding asymptotic values for the nonporous case. As would be expected, a higher heat

transfer rate is achieved at a larger penalty in pressure drop. Although both the Nusselt number and pressure drop increase as the particle diameter decreases, the increasing rates of both parameters are different. For example, as dp decrease from 0.2 to 0.1, the pressure drop increase by 52% but the Nusselt number increase only by 30%. Thus a relatively larger penalty for pressure drop will occur when a smaller particle diameter is selected to enhance heat transfer.

The combined effects of the Reynolds number and the particle diameter on heat transfer are illustrated in

channel is also increased as discussed earlier.

In the design of heat exchangers in which a porous medium is not present, a dimensionless parameter is usually used to evaluate the performance based on both the heat transfer enhancement and the associated pressure drop. This parameter is defined as:

In the present study, a similar parameter is used to evaluate the performance of the porous channel heat exchanger .in

as Re increases from 180 to 600, performance parameter decrease by 42% for Dp = 0.169, and by only 28% for Dp = 0.072

A numerical study of laminar forced convection in a two-dimensional porous channel with inlet and outlet ports was conducted. a smaller particle diameter can be used to achieve higher heat transfer enhancement , but a larger particle diameter leads to a more efficient performance based on heat transfer enhancement.