格子Bonzmann方法用于多孔介质和自由流体开口腔体内自然对流的数值模拟研究.pdfVIP

中国地球物理2011 ·6l· (2)地壳流体与地震预报、成矿成藏及碳封存 格子Boltzmann方法用于多孔介质 与自由流体开口腔体内自然对流的数值模拟研究 LatticeBoltzmannMethodfor NaturalConvection Modeling inPorousand FluidDomains Homogenous 戴传山刘学章’ DaiChuanshanLiu Xuezhang 天津大学机械工程学院地热研究培训中心 天津300072 1．引言 多孔介质空间内的对流与传热是自然界中普遍存在的一种复杂物理现象，也在工程实际中有着广泛的 强有力的优势。本文用格子Boltzmann方法数值模拟了二维开口方腔内垂直填充部分多孔介质的自然对流 换热问题，着重研究了有多孔介质与自由流体空间同时存在情况下，界面的热质传递问题，以及冷壁面自 由空间和多孔介质横向宽度比刀的关系。 2．模型和边界条件 该二维方腔的三边均为无滑移边界，左侧壁面的温度设为0，右侧为开口，流入的流体温度定常为1，流 出的流体边界设为绝热，上下壁面均为绝热。方腔内右部分区域为多孔介质，其他为自由空间，自由流体和多 面上的应力跳跃的边界条件【2】。应力跳跃的条件意味着在界面上平均速度是连续的，平均速度的斜率会有变化。 3．数值模拟结果 (1)经过对比发现，所得的结果与前人的结论基本吻合13】，为了得到不依赖以网格的解，还应使用适 当的网格数。本文在X和Y方向均划分了100个网格。随着的竖直方向】，的增加，Nu数一般先由小变大， 再由大变小，在方腔中间的斜率最大，即Nu数的变化最快。 向增加斜率会变大。Da数的变化也影响的Nu数变化的关键因素。其他因素一定时，随着Da的增加，Nu数 的增加斜率同样会变大，在y方向上的变化几乎呈直线状。方腔左侧的Nu数在随着不同孔隙率有不同程度的 变化。而随着孔隙率的增加，Nu数增加的斜率也会变大，但是这种影响程度要比l匕数和Da数来的小。 几乎未流动，左侧没有形成有效的传热。随着行逐渐变大，Nu数逐渐变大，在竖直，，方向的斜率也有所增加。 (4)对数值模拟得出的结果进行整理，得到了该模型下方腔左侧Nu的准则关联式，如下： Nu．0．076J配o·¨Dao·14rlo·18 103s 10-4sDa‘10—2，疗=0．25—4。 其中， Ra墨3×105， 参考文献 D in 【1】NieldA，Bcj鲫A．Convectionsporommedia,2ndedition[M]．New York：Springer-Verlag,1999． S．Momentumtransferatthe bL-蚋qeena end I： boundary a fluid 【2】Ochoa-TapiaJA，Whiteker poro璐mediumhomogeneous Transfer theoretical JournalofHeatandMass development[J]．International 1995；38：2635—2646． Nasr M．