格子Boltzmann亚网格模型及其流动模拟研究.pdfVIP

中国工程热物理学会 热机气动热力学 第十届年会 编号:012059 格子Boltzmann亚网格模型及其流动模拟 ‘ 晏一光 ‘郑楚光 郭照立 (华中科技大学煤嫩烧国家重点实脸室 武汉 430074) (联系电话:0278403 E-mail:yig_yanB263.net;co) 摘 要 格子Boltzman唯B)方法是一种模拟轴运现象的并行算法，它的发展为解决流体力学问刃提供了一 种新的途径.本文结合大祸模拟 (LES)和亚网格摸型探讨了一种格子Boltzmann亚网格的湍流模型 并利用这种模型对后台阶流动进行了高Reynolds致情况下的棋拟.并将此摸型所得的结果与实验和其 他方法的棋拟结果进行定性较。结果表明此模型能较好的模拟高Reynolds数湍流流动. ，关镇词韶子Boltzmann方念大涡棋拟，亚网格模型，后台阶流动 ) Abstract LatticeBoltzmannmethod,anewmeansforsolvingfluiddynamicsproblems,isaparallelalgorithmfor simulationoftransportphenomena.Combininglargeeddysimulationangsubgridmodel,wediscussalattice Boltzmannsubgidmodelforsimulatingturbulence.Asforbackwardfacingstepflowwecomparethe simulationresultsformthismodelandtheresultsfromothermodelsandexperimentalresultsqualitively?he comparisonindicatesthat也islatticeBoltzmannsubgridmodelcansimulateturbulengceandsomecharacters offlowverywell KeywordsLatticeBoltzmannmethod,Largeeddysimulation,Subgridmodel,Backwardfacingstepflow 1.简介 格子Boltztnann方法是一种从LGA方法发展而来的局部介观动力学方法，它可以 在微观上模拟宏观输运现象，如粘性流动、反应扩散问题等。可以选取合适的平衡分布 函数和松弛因子。并通过Chapman-Enskog展开和多尺度技术在二阶精度上恢复到宏观 的Navier-Stok“方程.自从1986年Frisch,Hasslacher,Pomeau(FHP)Ill提出第一个二 维模型并再现了不可压的N-S方程后，LB方法对流动的模拟不断取得进展，LB方法对 各种流动现象的模拟结果与试验以及其他一些方法的模拟结果均相符很好，而且，由于 LB方法独有的特点。它在模拟某些复杂流动现象上有无可比拟的优势，如多孔介质流 动。 本文在二维九点LB模型的基础上，通过考察大祸模拟 (LES)和亚网格模型中的 一些物理概念和处理方法，在LB方法中引入了湍流粘性的概念，并由此导出了模型中 总的松弛时间，即松弛因子，从而得到了一种模拟高Reynolds数流动的一种湍流模型. 我们利用此模型模拟了高Reynolds数下的二维后台阶流动，并与三维的大涡模拟和实验 结果进行了比较，结果相符较好，说明此模型能较好的模拟高Reynolds数流动。下面要 。本文得到国家重点基础研究专项经费(1999022207)资助 用到二维九点LB模型，但由于此模型己被广泛地探讨过[:l此处不再累述。 2.格子Boltamann亚网格模型 在前人的研究工作中，LB方法作为一种直接的模拟方法.其最小的捕获尺度为一 个网格单位，而其最大的捕获尺度取决于其特征尺度，这经常受到计算能力的限制，因 而到目前为止，LB方法仅限于解决低Re数流动问题。如果进一步减小流体枯性而提高 Re数，模拟中将会产生计算不稳定性，其根源在于流场中存在较大的局部速度梯度。 这意味着现有的LB模型中的粘性存在