Numerical Simulation of Driven Convective Heat Transfer Based Lattice Boltzmann Method in a Porous Cavity

A lattice Boltzmannmodel of the uniform velocity, driven convective thermal conductivity in a porous cavity is studied.TheDarcy, 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.


Introduction
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 [3], 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. [4] and by other researchers. Guo and Zhao used the representative elementary volume (REV) scale [5] 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 [5]. 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 [5].
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 [5] 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. and heat transfer phenomena in porous media [5,[8][9][10][11][12][13][14][15][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 ℎ = 0.5 and = −0.5, 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 as where and are the volume-averaged velocity and pressure, respectively. is the porosity of the medium. Consider = + (1 − ) / 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 (Re = / ) is assumed to be equal to the effective Reynolds number Re . The Prandtl number is Pr. The last term ( ) on the righthand side of (1b) represents the total body force due to the presence of a porous medium and other external force fields and is expressed as where Da is Darcy number defined as Da = / 2 , Gr is the Grashof number, is the unit vector in the -direction, and is a geometric function defined as [17] It can be seen that as = 1, 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 as where and 1 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.

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, eq ( , ) and eq ( , ) are the equilibrium distribution functions (EDF) of ( , ) and ( , ), respectively, and are given by where is the weight coefficient and is the speed of sound.
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 as where 0 , 1 , and V are given by By applying the Chapman-Enskog expansion [5], the macroscopic equations, (6) and (7), become equivalent to those in [5]. For the problem discussed here, the velocity at the inlet is uniform; that is, = 0.1, and the velocity on the (adiabatic) horizontal sidewalls is equal to zero.
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 [7] Nu = ∫ 1 0 ( , 0) . 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 = 0.5 were assumed, and we see that the density of the boundary layer near the hot wall becomes thinner as Da 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 Da 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 = 0.5), 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 Ri ≤ 0.1. 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 Da, , or Re 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.

Conclusions
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 Ri > 0.1. 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.