He-Luo模型同基于赫米特多项式Shan模型评估.pdf

第13卷第25期2013年9月科学技术与工程 V01．13No．25Sep．2013 Science and ⑥2013Sci．Tech．Engrg． 1671—1815(2013)25-7427—06TechnologyEngineering 航空航天 He．Luo模型同基于赫米特多项式 Shan模型评估 路益嘉 (清华大学航空航天学院工程力学系，北京100084) 缩流动时的优势很明显。 关键词He—Luo模型Shan模型赫米特多项式 中图法分类号V211．3； 文献标志码A 在过去的20年，格子波尔兹曼方法已经成为计 方程，因为在格子玻尔兹曼模拟中密度的空间变化 算流体力学中一个有效数值工具。跟传统的计算 不为零。He和Luo通过忽略可压缩效应项中马赫 流体力学方法相比，格子波尔兹曼方法容易实现， 数的高阶部分来减弱可压缩效应1。 可高度并行化，并能解决复杂边界问题。更重要的 文献证明LBE可以通过对连续波尔兹曼方程 是，由于其动力学特性，格子玻尔兹曼方法可以更 使用小马赫数泰勒展开旧1或赫米特多项式展开一’ 直接地对复杂的流体现象建模。 得到，后者基于流体的动力学表达并具有完整的理 格子波尔兹曼方法发展于研究纳维斯托克斯 论框架。Grad使用赫米特多项式提出了13阶矩来 方程的格子气自动机方法…，作为一个无噪音可替 代LGA的技术，布尔变量被单粒子概率分布函数取 方程在速度空间离散以便把格子波尔兹曼方法同 Grad 代心J。在小马赫数和小努森数时，通过假设分布函 13阶矩系统进行关联。基于这个方法，能够得 数接近于依赖气体密度和速度本地值的局部平衡 到高于纳维斯托克斯方程的波尔兹曼方程高阶近 态来线性化碰撞项目o。若干个小组通过将BGK碰似。但是Shah的模型依然是基于BGK碰撞项的单 撞项H1引入格子波尔兹曼方程得到了标准的格子 松弛方法。 BGK方程‘5’6J。麦克斯韦平衡态分布”1被建议用来 现通过泊肃叶流动和顶盖驱动方腔流两个算 取代Fermi—Dirac分布∽1作为局部平衡态分布函数。例比较和评价了He—Luo方法和Shan方法的性能。 在标准单松弛模型中，可压缩效应可能会产生比较 1模型和方法 严重的数值错误。一般来说，不可压缩可以通过使 密度恒定为一个常数得到，而在波尔兹曼方程中保 1．1格子模型 持一个定值密度是比较困难的。理论上来说格子 格子波尔兹曼方程使用二维九速对称格子 波尔兹曼方程更适合于模拟可压缩纳维斯托克斯 (D2Q9)，在这个模型中，空间被离散为一个方形格 子，有9个离散速度。 2013年4月21日收到，5月15日修改 作者简介：路益嘉，博士后。E—mail：yjia．1u@gmail．corn。 万方数据 7428 科学技术与工程 13卷 (0,0)， 01=0 1．3 Shan单松弛方法 “=1～4 (cos【学k【学])c 对波尔兹曼方程展开得到LBGK方程¨…。 d=5—8 Ⅳ阶Hermite多项式为 (cos[学】,sin[学肛 这里C=统届。，其中艿。为 州扑错V∞l0J 4础