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.
1. 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.
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.
Schematic of the physical model and coordinate for the problem with boundary conditions.
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 Th=0.5 and Tc=-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(1a)∇·u=0,(1b)∂u∂t+u·∇uε=-∇εp+1Ree∇2u+F,(1c)σ∂T∂t+u·∇T=1PrRe∇2T,where u and p are the volume-averaged velocity and pressure, respectively. ε is the porosity of the medium. Consider σ=ε+(1-ε)ρscps/ρfcpf represents the ratio between the heat capacities of the solid and fluid phases with ρs(ρf) and cps(cpf) being the density and capacity of the solid (fluid) phase, respectively. In this work, the Reynolds number (Re=LU/υ) is assumed to be equal to the effective Reynolds number Ree. The Prandtl number is Pr. The last term (F) 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 as(2)F=-εDaReu-εFεDaPu∣u+εGrRee2kT,where Da is Darcy number defined as Da=K/L2, Gr is the Grashof number, k is the unit vector in the y-direction, and Fε is a geometric function defined as [17](3)Fε=1.75150ε3.
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(4)Ri=GrRee2.
The Prandtl number can be defined as(5)Pr=Jeτ-0.5στ1-0.5,where τ and τ1 are the relaxation times and the viscosity ratio Je 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 fi(x,t) of the velocity field and the temperature distribution gi(x,t) of the temperature field can be defined in the LBM as [5, 18](6)fix+eiδt,t+δt-fix,t=-1τfix,t-fieqx,t+δtFi,(7)gix+eiδt,t+δt-gix,t=-1τ′gix,t-gieqx,t,where δ is the time increment. Here, fieq(x,t) and gieq(x,t) are the equilibrium distribution functions (EDF) of fi(x,t) and gi(x,t), respectively, and are given by(8)fieqx,t=ωiρ1+ei·ucs2+e·u22s4-u22s2,gieqx,t=ωiTσ+e·ucs2,where ωi is the weight coefficient and cs is the speed of sound. In the present work, we adopt the D2Q9 model, and the discrete velocities ei are given by(9)ei=0,fori=0;cosi-1π2,sini-1π2,fori=1:4;2cosi-5π2+π4,sini-5π2+π4,fori=5:8.
Here, c=δx/δt is the lattice spacing, the weight coefficients are ω0=4/9, ωi=1/9 (i=1–4), and ωi=1/36 (i=5–8), and the speed of sound is defined as cs=c/3.
The forcing term Fi in (6) is expressed as [5] (10)Fi=ωiρ1-12τ·ei-uε·Fcs2+u·ei×F·eiεcs4,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(11)ρ=∑fi,u=vb0+b02+b1v,where b0, b1, and v are given by(12)b0=121+δt2εDaRe,b1=12εFεDa,v=∑fieiρ+δt2GrReekT.
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, u=0.1, 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: (ε=0.999); (Da = 106); Pr = 0.71; τ=0.833; σ=1; Nx×Ny=257×257; 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.
Comparison of the average Nusselt number at the hot wall in this work and that reported in previous studies for Da = 106, Ri = 10−2, ε = 0.999, and τ = 0.833.
Re
Present
Reference [1]
Reference [6]
Reference [7]
100
2.04
1.97
1.94
2.01
400
4.05
4.03
3.84
3.91
1000
6.56
6.56
6.33
6.33
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] (13)Nu¯=∫01∂Tx,0∂ydx.
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.
Average Nusselt numbers at the hot and cold walls for different porosities.
ε
0.2
0.4
0.6
0.9
Cold wall
0.091
0.111
0.134
0.167
Hot wall
57.673
56.980
56.052
54.662
The isothermals for porosity: ε = 0.2, 0.4, 0.6, 0.9.
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.
Average Nusselt numbers at the hot and cold wall for different Darcy numbers.
Da
10-4
10-3
10-2
10-1
Cold wall
1.238
0.393
0.122
0.028
Hot wall
36.120
44.833
56.524
60.358
Isothermals for the Darcy numbers: Da = 10−4, 10−3, 10−2, and 10−1.
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.
Average Nusselt numbers at the hot and cold wall for different Reynolds numbers.
Re
10
100
200
400
Cold wall
0.091
0.122
0.139
0.146
Hot wall
6.254
56.524
101.764
171.784
Isothermals for the Reynolds numbers: Re = 10, 100, 200, and 400.
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.
Average Nusselt numbers at the hot and cold wall for different Richardson numbers.
Ri
10-4
10-2
10-1
1.0
Cold wall
0.122
0.122
0.122
0.134
Hot wall
56.524
56.524
56.524
53.188
Isothermals for the Richardson numbers: Ri = 10−4, 10−2, 10−1, and 1.0.
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. Figure 6(d) shows the significant effect of Ri on the velocity.
Horizontal velocity profiles in the vertical midplane for different (a) Darcy numbers, (b) porosities, (c) Reynolds numbers, and (d) Richardson numbers.
5. 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.
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.
Acknowledgment
This work was supported by the National Nature Science Foundation of China under Grants nos. U1262109, 51476145, and 51476146 which are gratefully acknowledged.
ChaiZ. H.GuoZ. L.ShiB. C.Lattice Boltzmann simulation of mixed convection in a driven cavity packed with porous medium20074448Springer802809Lecture Notes in Computer ScienceShankarP. N.DeshpandeM. D.Fluid mechanics in the driven cavity20003293136OztopH. F.Combined convection heat transfer in a porous lid-driven enclosure due to heater with finite length200633677277910.1016/j.icheatmasstransfer.2006.02.0032-s2.0-33747025304GuoZ.ShiB.ZhengC.A coupled lattice BGK model for the Boussinesq equations200239432534210.1002/fld.337MR19047862-s2.0-0037054101GuoZ. L.ZhaoT. S.Lattice Boltzmann model for incompressible flows through porous media200266303630410.1103/physreve.66.0363042-s2.0-34548358628KhanaferK. M.ChamkhaA. J.Mixed convection flow in a lid-driven enclosure filled with a fluid-saturated porous medium199942132465248110.1016/s0017-9310(98)00227-02-s2.0-0345535644FanJ.TsangH. W. K.Effect of clothing thermal properties on the thermal comfort sensation during active sports200878211111810.1177/07316844080800462-s2.0-67651095761DurlofskyL.BradyJ. F.Analysis of the Brinkman equation as a model for flow in porous media19873011332910.1063/1.8664652-s2.0-0023140832PrasadV.KulackiF. A.Natural-convection in porous-media bounded by short concentric vertical cylinder1985107114715410.1115/1.32473712-s2.0-0022012033RajamaniR.SrinivasC.NithiarasuP.SeetharamuK. N.Convective heat transfer in axisymmetric porous bodies19955982983710.1108/EUM0000000004126ZBL0860.760862-s2.0-0029395262NithiarasuP.SeetharamuK. N.SundararajanT.Natural convective heat transfer in a fluid saturated variable porosity medium1997401639553967ZBL0925.7666010.1016/S0017-9310(97)00008-22-s2.0-0031250142NithiarasuP.RavindranK.A new semi-implicit time stepping procedure for buoyancy driven flow in a fluid saturated porous medium19981651-4147154MR166353610.1016/s0045-7825(98)00036-x2-s2.0-0032476850HickoxC. E.GartlingD. K.A numerical study of natural convection in a vertical, annular, porous layer198528372072310.1016/0017-9310(85)90196-62-s2.0-0022028986TatsuoN.ToruT.MitsuhiroS.YujiK.HiroyukiO.Numerical analysis of natural convection in a rectangular enclosure horizontally divided into fluid and porous regions198629688989810.1016/0017-9310(86)90184-52-s2.0-0022738546ChenS.DoolenG. D.Lattice Boltzmann method for fluid flows199830329364SucciS.2001New York, NY, USAOxford University PressNumerical Mathematics and Scientific ComputationMR1857912ErgunS.Fluid flow through packed columns1952488994GuoZ.ZhaoT. S.A lattice boltzmann model for convection heat transfer in porous media200547215717710.1080/104077905908834052-s2.0-13844272699