Journal of Energy

Volume 2013 (2013), Article ID 327179, 7 pages

http://dx.doi.org/10.1155/2013/327179

## Numerical Investigation of Forced Convection in a Channel with Solid Block inside a Square Porous Block

Islamic Azad University Science and Research Ayatollah Amoli Branch, Amol, Iran

Received 13 February 2013; Accepted 7 April 2013

Academic Editor: S. A. Kalogirou

Copyright © 2013 Neda Janzadeh and Mojtaba Aghajani Delavar. 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

The fluid flow and heat transfer in a porous medium have received considerable attention due to its importance in many engineering applications. In this study numerical investigation of fluid flow and heat transfer over a hot solid block inside a square porous block located in a channel was carried out. The lattice Boltzmann method with nine velocities, D2Q9, was used for numerical simulations. Brinkman-Forchheimer model was successfully used to simulate fluid flow in porous media. The effects of parameters such as porosity, Reynolds number on flow pattern, and heat transfer were studied. The different effects of mentioned parameters were discussed in the paper.

#### 1. Introduction

The fluid flow and transport phenomena in a porous medium have been used in many fields of science and engineering, such as geothermal energy extraction, solar collector, solar absorbers food processing, fuel cells, petroleum processing, catalytic and chemical particle beds, transpiration cooling, electronic cooling, drying processes, porous bearing, nuclear reactors, and many others.

*Thermal Insulation*. Fluid flow and convective heat transfer in conduits fully and partially filled with porous medium have been investigated analytically, experimentally, and numerically by many researchers. Poulikakos and Kazmierczak [1] theoretically studied the convection heat transfer in the developed region of two parallel plates and in a circular pipe with porous media attached to the wall. Fu et al. [2] numerically simulated the convection heat transfer through a single porous block mounted on the heated wall, and they studied the effect of different parameters including the porous block size and length and porosity of porous block and Reynolds number on convection heat transfer enhancement. Li et al. [3] numerically investigated the effect of staggered porous blocks mounted on the wall in a laminar flow channel on the velocity field and local heat transfer. Huang et al. [4] numerically studied steady-state heat transfer to a fluid passing through multiple heated block mounted on the wall of a channel with porous covers.

Guo et al. [5] investigated numerically the heat transfer of pulsating flow through a pipe with a porous medium attached to the pipe wall. In their work they used conventional control volume method. Hamdan et al. [6] used an implicit finite-difference method to simulate forced convection in a porous substrate inserted in the core of a parallel plate channel. Jiang et al. [7] investigate the forced convection heat transfer in plate channels saturated with air and water which contain sintered bronze porous media.

The lattice Boltzmann method (LBM) is a powerful numerical technique based on kinetic theory for simulation of fluid flows and modeling the physics in fluids [8–10]. The problem of convective heat in porous medium was investigated numerically by Guo and Zhao [11]. In their work, they used the lattice Boltzmann method to simulate the temperature field. Seta et al. [12] applied the lattice Boltzmann method to simulation of natural convection in porous media. Shokouhmand et al. [13] simulated the laminar flow and convective heat transfer in a porous medium inside the core of two parallel plates of a conduit using lattice Boltzmann method.

Hot solid block inside a square porous block may be applied in some engineering application such as porous media as an active layer in reacting chemical flows, fuel cell, heat exchangers, and electronic devices cooling.

#### 2. Numerical Model

The general form of lattice Boltzmann equation with nine velocities, D2Q9, with external force can be written as [9]
where is the discrete lattice velocity in direction is the external force, denotes the lattice time step, is the lattice relaxation time, denotes the equilibrium distribution function, is weighting factor, and is the lattice fluid density. To consider both the flow and the temperature fields, the thermal LBM utilizes two distribution functions, and , for flow and temperature fields, respectively. The distribution function is as same as discussed previously; the distribution function is as follows [9]:
The flow properties are defined as (*i* denotes the component of the Cartesian coordinates)
The Brinkman-Forchheimer equation was used for flow in porous regions that is written as [9, 12]
where is the porosity, is the permeability, is the effective viscosity, is the kinematic viscosity, and is the acceleration due to gravity. The last term in the right hand in the parenthesis is the total body force, , which was written by using the widely used Ergun’s relation [16]. For porous medium the corresponding distribution functions are as the same as (1). But the equilibrium distribution functions and the best choice for the forcing term are
The forcing term defines the fluid velocity as
According the previous equations is related to , so (6) is nonlinear for the velocity. A temporal velocity is used to solve this nonlinear problem [13] as follows:
The effective thermal conductivity, of the porous medium should be recognized for proper investigation of conjugate convection and conduction heat transfer in porous zone, which was calculated by [16]

#### 3. Boundary Conditions

From the streaming process the distribution functions out of the domain are known. The unknown distribution functions are those toward the domain. In Figure 1 the unknown distribution function which needs to be determined, are shown as dotted lines. Regarding the boundary conditions of the flow field, the solid walls are assumed to be no slip, and thus the bounce-back scheme is applied. This scheme specifies the outgoing directions of the distribution functions as the reverse of the incoming directions at the boundary sites. For example, for flow field in the north boundary the following conditions is used: In this channel the inlet velocity at west boundary is known . In LBM method the inward distribution functions at the boundaries have to be specified. So the , and are unknowns for west boundary condition. The unknown distribution functions are calculated as In this problem the outlet velocity is unknown. The east boundary condition in Figure 1 represents the outlet condition. Then , and need to be calculated at the east boundary as Furthermore, for the temperature field, the local temperature is defined by (3). For isothermal boundaries such as a bottom hot wall the unknown distribution functions are evaluated as where is the lattice on the boundary. For adiabatic boundary condition such as a bottom wall the unknown distribution functions are evaluated as

#### 4. Computational Domain

Computational domain consists of a hot solid block inside a square porous block located in a channel (Figure 2). The channel height, length, and other simulation parameters are illustrated in the Figure 2 and Table 1.

#### 5. Validation and Grid Independent Check

In this study, flow pattern and thermal field were simulated in a channel with a hot solid block inside a porous block by using lattice Boltzmann method. Figure 3 compares well the velocity profile for flow in a clear channel captured by LBM and results presented by Nield [14]. In this figure the result of simulation is compared well with Seta et al. [15]. Table 2 shows the computed Nusselt number by LBM and Kays and Crawford [17], and good agreement is obtained.

#### 6. Results and Discussion

In the present work the thermal lattice Boltzmann model with nine velocities was used to solve the force convection heat transfer in a channel containing hot solid block inside a square porous block. The effects of porosity and Reynolds number, on the flow field and convective heat transfer, were investigated.

Figure 4 shows the velocity and temperature contours for different values of porosity at . With increasing the porosity, velocity rises slightly in the top and bottom of porous block, because fluid tends to flow in clear passages with lower pressure drop. At higher values of porosity, fluid changes easier its path to clear passages, so velocity increases in these passages.

According to (8) with increasing the porosity, the effective thermal conductivity in the porous block reduces. At higher value of porosity, the heat transfer between fluid and solid block decreases. So the fluid temperature will increase.

The velocity and temperature contours for different Reynolds number have been drown in Figure 5.

It can be seen that as Reynolds number increases the flow velocity increases. With increasing the velocity, the convective coefficient increased according to (), and the magnitude of heat transfer increases. But According to () with increasing the velocity, the temperature gradient of the fluid reduces. So the outlet fluid temperature will decrease.

The temperature and velocity profiles for various Reynolds numbers are represented in Figure 6. As expected the flow velocity increases with increasing the Reynolds number. The higher values of Reynolds number cause the temperature of zones to be decreased due to the decrease in heat transfer process. In Figure 7 the average fluid temperature for various Reynolds numbers at different porosity is drown that shows the effect of Reynolds number on temperature more than porosity. As mentioned before at higher porosity obtains a lower average temperature for fluid.

#### 7. Conclusion

Fluid flow and heat transfer in a channel with solid block inside a square porous block were carried out in this paper. The effect of parameters including porosity and Reynolds number on the velocity and thermal field was investigated numerically by using lattice Boltzmann method. Decreasing the porosity leads to a higher temperature of fluid due to increasing which causes the heat transfer enhancement in porous block. With increasing the Reynolds number the fluid velocity increased, and the average fluid temperature will be reduced.

#### References

- D. Poulikakos and M. Kazmierczak, “Forced convection in a duct partially filled with a porous material,”
*Journal of Heat Transfer*, vol. 109, no. 3, pp. 653–662, 1987. View at Google Scholar · View at Scopus - W. S. Fu, H. C. Huang, and W. Y. Liou, “Thermal enhancement in laminar channel flow with a porous block,”
*International Journal of Heat and Mass Transfer*, vol. 39, no. 10, pp. 2165–2175, 1996. View at Publisher · View at Google Scholar · View at Scopus - H. Y. Li, K. C. Leong, L. W. Jin, and J. C. Chai, “Analysis of fluid flow and heat transfer in a channel with staggered porous blocks,”
*International Journal of Thermal Sciences*, vol. 49, no. 6, pp. 950–962, 2010. View at Publisher · View at Google Scholar · View at Scopus - P. C. Huang, C. F. Yang, J. J. Hwang, and M. T. Chiu, “Enhancement of forced-convection cooling of multiple heated blocks in a channel using porous covers,”
*International Journal of Heat and Mass Transfer*, vol. 48, no. 3-4, pp. 647–664, 2005. View at Publisher · View at Google Scholar · View at Scopus - Z. Guo, S. Y. Kim, and H. J. Sung, “Pulsating flow and heat transfer in a pipe partially filled with a porous medium,”
*International Journal of Heat and Mass Transfer*, vol. 40, no. 17, pp. 4209–4218, 1997. View at Google Scholar · View at Scopus - M. O. Hamdan, M. A. Al-Nimr, and M. K. Alkam, “Enhancing forced convection by inserting porous substrate in the core of a parallel-plate channel,”
*International Journal of Numerical Methods for Heat and Fluid Flow*, vol. 10, no. 5, pp. 502–517, 2000. View at Google Scholar · View at Scopus - P. X. Jiang, M. Li, T. J. Lu, L. Yu, and Z. P. Ren, “Experimental research on convection heat transfer in sintered porous plate channels,”
*International Journal of Heat and Mass Transfer*, vol. 47, no. 10-11, pp. 2085–2096, 2004. View at Publisher · View at Google Scholar · View at Scopus - S. Succi,
*The Lattice Boltzmann Equation for Fluid Dynamics and Beyond*, Clarendon Press, Oxford, UK, 2001. - A. A. Mohammad,
*Lattice Boltzmann Method, Fundamentals and Engineering Applications with Computer Codes*, Springer, London, UK, 2011. - P. H. Kao, Y. H. Chen, and R. J. Yang, “Simulations of the macroscopic and mesoscopic natural convection flows within rectangular cavities,”
*International Journal of Heat and Mass Transfer*, vol. 51, no. 15-16, pp. 3776–3793, 2008. View at Publisher · View at Google Scholar · View at Scopus - Z. Guo and T. S. Zhao, “A lattice Boltzmann model for convection heat transfer in porous media,”
*Numerical Heat Transfer, Part B*, vol. 47, no. 2, pp. 157–177, 2005. View at Publisher · View at Google Scholar · View at Scopus - T. Seta, E. Takegoshi, and K. Okui, “Lattice Boltzmann simulation of natural convection in porous media,”
*Mathematics and Computers in Simulation*, vol. 72, no. 2-6, pp. 195–200, 2006. View at Publisher · View at Google Scholar · View at Scopus - H. Shokouhmand, F. Jam, and M. R. Salimpour, “Simulation of laminar flow and convective heat transfer in conduits filled with porous media using lattice Boltzmann Method,”
*International Communications in Heat and Mass Transfer*, vol. 36, no. 4, pp. 378–384, 2009. View at Publisher · View at Google Scholar · View at Scopus - D. A. Nield, “Forced convection in a parallel plate channel with asymmetric heating,”
*International Journal of Heat and Mass Transfer*, vol. 47, no. 25, pp. 5609–5612, 2004. View at Publisher · View at Google Scholar · View at Scopus - T. Seta, E. Takegoshi, and K. Okui, “Thermal lattice Boltzmann model for incompressible flows through porous media,”
*Journal of Thermal Science and Technology*, vol. 1, no. 2, pp. 90–100, 2006. View at Publisher · View at Google Scholar - S. Ergun, “Fluid flow through packed columns,”
*Chemical Engineering Progress*, vol. 48, no. 2, pp. 89–94, 1952. View at Google Scholar - W. M. Kays and M. E. Crawford,
*Solutions Manual, Convective Heat and Mass Transfered*, McGraw-Hill, New York, NY, USA, 3rd edition, 1993.