Transport Phenomena in Porous Media and Fractal GeometryView this Special Issue
Research Article | Open Access
Numerical Simulation of Driven Convective Heat Transfer Based Lattice Boltzmann Method in a Porous Cavity
A lattice Boltzmann model of the uniform velocity, driven convective thermal conductivity in a porous cavity is studied. The Darcy, Richardson, and Reynolds numbers are shown to have a significant influence on the heat transfer behavior and the horizontal velocity of the flow field, while the porosity has little influence on either. The model is validated by the average Nusselt number at different Reynolds numbers, and the numerical results are in good agreement with available published data.
Driven convective heat transfer in porous cavities has become a topic of interest in textile fiber and engineering fields [1, 2] through efforts to improve chemical and petroleum engineering, cooling of electronic devices , food processing, and functional clothing design.
Over the past several decades, convection heat transfer in porous media has been studied both experimentally and theoretically, but there have been relatively few numerical investigations. A coupled lattice Bhatnagar-Gross-Krook (CLBGK) model using the Boussinesq equations to simulate fluid flow in porous media was proposed by Guo et al.  and by other researchers. Guo and Zhao used the representative elementary volume (REV) scale  to study driven convection at a uniform velocity. The REV scale is much larger than the pore scale but much smaller than the domain scale. A number of other models have also been proposed, such as the Darcy model, the Brinkman-extended Darcy model, and the Forchheimer-extended Darcy model . These differ from conventional numerical methods such as the finite-difference, finite-volume, and finite-element methods. Recently, a generalized lattice Boltzmann equation (GLBE) was proposed for isothermal incompressible flows in porous media .
The aim of the present paper is to study the effects of Darcy number (Da), porosity (), Reynolds number (Re), and Richardson number (Ri) on the driven convection using GLBE-based simulations. We adopt a model in which the sidewalls of the porous cavity are at different temperatures. To the best of our knowledge, previous simulations of driven flow have been limited to isothermal problems solved using lattice Boltzmann methods (LBMs); there have been no reports of studies using an LBE model. The flow and heat transfer in porous media are coupled by rather complicated mechanisms by extending the GLBE  proposed to thermal flows. We use the GLBM to simulate the velocity field and use a newly defined temperature distribution function LBM to simulate the temperature field.
We discuss the physical problem and mathematical formulation of porous flow in Section 2, and, in Section 3, the proposed LBM formulation including the proposed treatments of the velocity and temperature boundary conditions is presented. The numerical results and a conclusion are given in Sections 4 and 5, respectively.
2. Mathematical Formulation
Over the past decades, the lattice Boltzmann method (LBM) has been developed considerably for simulating fluid flows and heat transfer phenomena in porous media [5, 8–16]. The physical problem is illustrated schematically in Figure 1.
A square cavity is filled with an isotropic, homogeneous porous medium that is heated by the two vertical sidewalls. The fluid is assumed to flow into the porous cavity from the left-hand side. The (dimensionless) temperatures of the hot and cold sidewalls are and , respectively. The medium is also assumed to have a constant porosity. The governing equation of the generalized model for the incompressible fluid flow and heat transfer in such a porous media is given aswhere and are the volume-averaged velocity and pressure, respectively. is the porosity of the medium. Consider represents the ratio between the heat capacities of the solid and fluid phases with and being the density and capacity of the solid (fluid) phase, respectively. In this work, the Reynolds number () is assumed to be equal to the effective Reynolds number . The Prandtl number is . The last term () on the right-hand side of (1b) represents the total body force due to the presence of a porous medium and other external force fields and is expressed aswhere is Darcy number defined as , is the Grashof number, is the unit vector in the -direction, and is a geometric function defined as 
It can be seen that as , that is, in the absence of a porous medium, the governing equations, given in (1a)–(1c), reduce to the standard equations for free flows. In addition, some nondimensional parameters based on the governing equations can also be defined. The Richardson number is the dimensionless number that expresses the ratio of the buoyancy term, and we can define it in terms of the Grashof and Reynolds numbers as
The Prandtl number can be defined aswhere and are the relaxation times and the viscosity ratio and the ratio between the heat capacities of the solid and fluid phases are set to 1.
3. Thermal Lattice Boltzmann Model for Porous Flows
The evolution of the single-particle density distribution of the velocity field and the temperature distribution of the temperature field can be defined in the LBM as [5, 18]where is the time increment. Here, and are the equilibrium distribution functions (EDF) of and , respectively, and are given bywhere is the weight coefficient and is the speed of sound. In the present work, we adopt the D2Q9 model, and the discrete velocities are given by
Here, is the lattice spacing, the weight coefficients are , (, and (), and the speed of sound is defined as .
The forcing term in (6) is expressed as  which gives the total force due to the presence of the porous medium and the other external force fields. Thus, the fluid density and velocity can be expressed aswhere , , and are given by
By applying the Chapman-Enskog expansion , the macroscopic equations, (6) and (7), become equivalent to those in . For the problem discussed here, the velocity at the inlet is uniform; that is, , and the velocity on the (adiabatic) horizontal sidewalls is equal to zero.
4. Numerical Results and Discussion
To evaluate the accuracy of the present numerical method, we assume the following parameter values: (); (Da = 106); Pr = 0.71; ; ; ; Re = 100, 400, 1000; and Ri = 1.0 × 10−2, 6.25 × 10−4, and 1.0 × 10−4.
The average Nusselt numbers at the hot wall for the present LBM are in good agreement with those reported in previous studies. A comparison is given in Table 1.
To investigate the effect of the porosity on the isothermals and average Nusselt numbers, calculations were carried out under conditions of Ri = 0.0001, Da = 0.01, and Re = 100 by varying the porosity from 0.2 to 0.9. The average Nusselt number at the isothermal walls is defined as 
This number is an important dimensionless parameter for describing the convective heat transport. Figure 2 shows that varying has little impact on the temperature field. There is a steep temperature gradient in the horizontal direction at the right wall due to the driven convection and a weak uniform temperature gradient in the left-hand side of the cavity. As shown in Table 2, the average Nusselt numbers at the cold wall increase slightly with increasing porosity but those of the hot wall decrease, which is consistent with the behavior of the temperature field.
The effect of the Darcy number on the temperature field is shown in Figure 3, where values of Re = 100, Ri = 10−4, and were assumed, and we see that the density of the boundary layer near the hot wall becomes thinner as increases. The data in Table 3 reveal that the average Nusselt numbers at the hot wall increase and those on the cold wall decrease as increases from 10−4 to 10−1. These results indicate that driven convection across the cavity is the dominant mechanism and that conduction heat transfer is strong at high Da values.
A similar phenomenon is also observed as the Reynolds number increases (Da = 10−2, Ri = 10−4, and ), as shown in Figure 4. As Re increases, the density of the isothermals near the hot wall reduces, and bending is observed due to the convection effect near the horizontal sidewalls. As can be seen in Table 4, the average Nusselt numbers near the hot wall increase rapidly with an increasing Re value, while there is only a small increase in the numbers near the cold wall.
The Richardson number represents the importance of natural convection relative to the forced convection, and its effect on the temperature field is shown in Figure 5. We see that there is little change in the field as Ri increases, but a strong temperature gradient appears suddenly in the upper right-hand region of the cell when . This is consistent with the average Nusselt numbers listed in Table 5.
Some additional insight can be gained by examining the horizontal velocity in the vertical midplane of the cavity (Figure 6). We see that there are two peaks in the velocity near the top and bottom walls that decrease as , , or decreases, as shown in Figures 6(a), 6(b), and 6(c), respectively. This decrease in the velocity peaks is slightly weaker in the case of and Re. Figure 6(d) shows the significant effect of Ri on the velocity.
We have proposed an LBM for solving the driven convective heat transfer problem in porous media with a temperature field. The temperature distribution was obtained by solving an additional LBE, and the simulation results showed that variations in and Ri had little impact on the temperature field and velocity profile, whereas significant changes were seen when . The influences of the Reynolds and Darcy numbers on the temperature field and the average Nusselt numbers near the hot wall were significant. Compared with traditional numerical methods, the LBM is more suitable for solving practical application problems in engineering.
Conflict of Interests
The authors declare that they have no financial or personal relationships with other people or organizations that can inappropriately influence their work; there is no professional or other personal interest of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in, or the review of, this paper.
This work was supported by the National Nature Science Foundation of China under Grants nos. U1262109, 51476145, and 51476146 which are gratefully acknowledged.
- Z. H. Chai, Z. L. Guo, and B. C. Shi, “Lattice Boltzmann simulation of mixed convection in a driven cavity packed with porous medium,” in Computational Scienc—ICCS 2007, vol. 4448 of Lecture Notes in Computer Science, pp. 802–809, Springer, 2007.
- P. N. Shankar and M. D. Deshpande, “Fluid mechanics in the driven cavity,” Annual Review of Fluid Mechanics, vol. 32, pp. 93–136, 2000.
- H. F. Oztop, “Combined convection heat transfer in a porous lid-driven enclosure due to heater with finite length,” International Communications in Heat and Mass Transfer, vol. 33, no. 6, pp. 772–779, 2006.
- Z. Guo, B. Shi, and C. Zheng, “A coupled lattice BGK model for the Boussinesq equations,” International Journal for Numerical Methods in Fluids, vol. 39, no. 4, pp. 325–342, 2002.
- Z. L. Guo and T. S. Zhao, “Lattice Boltzmann model for incompressible flows through porous media,” Physical Review E, vol. 66, no. 3, Article ID 036304, 2002.
- K. M. Khanafer and A. J. Chamkha, “Mixed convection flow in a lid-driven enclosure filled with a fluid-saturated porous medium,” International Journal of Heat and Mass Transfer, vol. 42, no. 13, pp. 2465–2481, 1999.
- J. Fan and H. W. K. Tsang, “Effect of clothing thermal properties on the thermal comfort sensation during active sports,” Textile Research Journal, vol. 78, no. 2, pp. 111–118, 2008.
- L. Durlofsky and J. F. Brady, “Analysis of the Brinkman equation as a model for flow in porous media,” Physics of Fluids, vol. 30, no. 11, p. 3329, 1987.
- V. Prasad and F. A. Kulacki, “Natural-convection in porous-media bounded by short concentric vertical cylinder,” ASME Journal of Heat Transfer, vol. 107, no. 1, pp. 147–154, 1985.
- R. Rajamani, C. Srinivas, P. Nithiarasu, and K. N. Seetharamu, “Convective heat transfer in axisymmetric porous bodies,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 5, no. 9, pp. 829–837, 1995.
- P. Nithiarasu, K. N. Seetharamu, and T. Sundararajan, “Natural convective heat transfer in a fluid saturated variable porosity medium,” International Journal of Heat and Mass Transfer, vol. 40, no. 16, pp. 3955–3967, 1997.
- P. Nithiarasu and K. Ravindran, “A new semi-implicit time stepping procedure for buoyancy driven flow in a fluid saturated porous medium,” Computer Methods in Applied Mechanics and Engineering, vol. 165, no. 1-4, pp. 147–154, 1998.
- C. E. Hickox and D. K. Gartling, “A numerical study of natural convection in a vertical, annular, porous layer,” International Journal of Heat and Mass Transfer, vol. 28, no. 3, pp. 720–723, 1985.
- N. Tatsuo, T. Toru, S. Mitsuhiro, K. Yuji, and O. Hiroyuki, “Numerical analysis of natural convection in a rectangular enclosure horizontally divided into fluid and porous regions,” International Journal of Heat and Mass Transfer, vol. 29, no. 6, pp. 889–898, 1986.
- S. Chen and G. D. Doolen, “Lattice Boltzmann method for fluid flows,” Annual Review of Fluid Mechanics, vol. 30, pp. 329–364, 1998.
- S. Succi, The Lattice Boltzmann Equation for Fluid Dynamics and Beyond, Numerical Mathematics and Scientific Computation, Oxford University Press, New York, NY, USA, 2001.
- S. Ergun, “Fluid flow through packed columns,” Chemical Engineering Progress, vol. 48, pp. 89–94, 1952.
- Z. Guo and T. S. Zhao, “A lattice boltzmann model for convection heat transfer in porous media,” Numerical Heat Transfer, Part B: Fundamentals, vol. 47, no. 2, pp. 157–177, 2005.
Copyright © 2015 You-Sheng Xu et al. 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.