多松弛格子Boltzmann方法的曲边界数值模拟.pdfVIP

多松弛格子Boltzmann方法的曲边界数值模拟 李凯1，钟诚文1，2，卓丛山1，郑佩l，陈效鹏3 (1．西北工业大学翼型、叶栅空气动力学国防科技重点实验室，陕西西安710072) (2．西北工业大学高性能计算研究与发展中心，陕西西安710072) (3．西北工业大学力学与土木建筑学院，陕西西安710072) 格，通过引入曲边插值边界处理方法YMS格式，分别对Re=40，Ma=0．1下的 和对比研究，两者的结果与文献中的数值结果均吻合得很好。证明了这种方法 适用于低雷诺数粘性流动数的值模拟，为以后复杂流动的模拟提供了一种有效 的途径。 关键词：格子Boltzmann，多松弛模型，YMS格式，笛卡尔网格，曲边界 1引言 Boltzmann 经过20多年的发展，格子Boltzmann方法(LatticeMethod，LBM)已经成 为流体数值计算的一种重要方法。作为一种介观方法，LBM可以比较直观、方便地处 理流体内部以及流体与周围环境的相互作用，而且凭借其计算简单，容易实现以及天然 的并行性等优势在流体力学的诸多领域都得到了成功的应用。传统的单松弛模型简单实 与单松弛模型相比具有很大的优势，因此MRT模型被越来越多地应用于高雷诺数流动 模拟。另外，在LBM中边界处理非常重要，直接影响到计算格式的精度和稳定性。但 是复杂流场的边界往往是不规则的，对于基于笛卡尔网格的标准LBM，物理边界点不 提出一种基于虚平衡态分布的插值方法，而后Meil41等对该方法作了改进，由于该方法 中插值的松弛因子与LBM中的无量纲松弛时间有关，因而稳定性受到一定的影响。Guot5】 等人结合宏观量插值和非平衡态外推的思想提出应用非平衡态外推方法处理曲边界。 Bouzidil6】等人利用反弹格式的处理方法，按照反弹方向直接对分布函数作插值。 易实现，而且具有二阶数值精度，也可用于运动边界条件的处理。由于该方法中需要参 照q使用不同的计算格式，从而导致分布函数的跳变，产生非物理的震荡。对此Yuf81 16 型，引入YMS格式来进行边界处理，分别对圆柱和翼型绕流进行了数值模拟和对比研 究。 2数值方法 2．1LBM模型 首先在本文中采用了D2Q9格子模型，其离散速度矢量定义为： c·， 平衡分布函数定义为： 烈咖wp[·等蜡丢] ㈤ 其中Wo=4／9，嵋=…=w4=1／9，w5=…=％=1／36，巳=1／4j为声速。 LBM的MRT模型中分别定义了由9个矩和分布函数构成的矩空间和分布函数空间 【9】： l聊)=(p，邮，六…q矗，q，，P。，P掣)2 (3) If)--(yo(x,f)，石(x，f)，…，五(x，f))2(4) 在矩空间中p是密度，P是能量，s是能量的平方，正和矗分别是动量的x和Y方向的 分量，qx和g。分别是能量通量的x和Y方向的分量，‰和％分别是粘性应力张量的对 角和非对角分量。矩空间和分布函数空间之间的关系为： I所)=M·I厂) (5) 其中 17 l l 4 2 ●2 4 1 1 O ●o乏o ●o乏o ●o乏o 1 o，2●o M=0 2 o 1 o (6) 0 ●o之●之o O o1