基于MRT－LBM的方腔流动局部加密算法-JournalofNortheastern.PDF

第３９卷第４期 东 北 大 学 学 报 ( 自 然 科 学 版 ) Ｖｏｌ.３９ꎬＮｏ.４ ２０１８年 ４ 月 Ｊｏｕｒｎａｌ ｏｆ Ｎｏｒｔｈｅａｓｔｅｒｎ Ｕｎｉｖｅｒｓｉｔｙ(Ｎａｔｕｒａｌ Ｓｃｉｅｎｃｅ) Ａｐｒ. ２０ １８ 　 ｄｏｉ:１０.１２０６８/ｊ.ｉｓｓｎ.１００５－３０２６.２０１８.０４.０１７ 基于ＭＲＴ－ＬＢＭ的方腔流动局部加密算法 陈颂英ꎬ汪　 超ꎬ曲延鹏ꎬ王润堃 (山东大学 高效洁净机械制造教育部重点实验室ꎬ山东 济南　 ２５００６１) 摘　 　 　 要: 使用多参数弛豫模型的格子玻尔兹曼方法及Ｄ２Ｑ９模型ꎬ对方腔流顶盖下方左右奇异角落处应 用局部网格加密模拟方腔流动特性.在粗细网格界面上ꎬ建立了分布函数的转换公式ꎬ采用３次样条插值计算 加密界面上插值点的参量ꎬ根据Ｃｈａｐｍａｎ－Ｅｎｓｋｏｇ分析ꎬ给出了应力计算公式.数值计算结果表明:对雷诺数 １０００ 的方腔流动ꎬ沿腔中心线的速度分布与经典文献结果对比效果良好ꎬ压力、涡量轮廓图的噪声明显降低ꎬ 应力振荡明显减少ꎬ证实了所提局部加密方法的有效性. 关　 键　 词: 格子Ｂｏｌｔｚｍａｎｎꎻ多参数弛豫模型ꎻ局部网格加密ꎻ方腔ꎻＤ２Ｑ９模型 中图分类号: Ｏ ３５　 　 　 文献标志码: Ａ　 　 　 文章编号: １００５－３０２６(２０１８)０４－０５３８－０５ Ｌｏｃａｌ Ｇｒｉｄ Ｒｅｆｉｎｅｍｅｎｔ Ｉｍｐｌｅｍｅｎｔａｔｉｏｎ ｏｆ Ｌｉｄ Ｄｒｉｖｅｎ Ｃａｖｉｔｙ Ｆｌｏｗ Ｂａｓｅｄ ｏｎ ＭＲＴ￣ＬＢＭ ＣＨＥＮ Ｓｏｎｇ￣ｙｉｎｇꎬＷＡＮＧ ＣｈａｏꎬＱＵ Ｙａｎ￣ｐｅｎｇꎬＷＡＮＧＲｕｎ￣ｋｕｎ (Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｈｉｇｈ￣Ｅｆｆｉｃｉｅｎｃｙ ａｎｄ Ｃｌｅａｎ Ｍｅｃｈａｎｉｃａｌ Ｍａｎｕｆａｃｔｕｒｅꎬ Ｓｈａｎｄｏｎｇ ＵｎｉｖｅｒｓｉｔｙꎬＪｉｎａｎ ２５００６１ꎬ Ｃｈｉｎａ. Ｃｏｒｒｅｓｐｏｎｄｉｎｇ ａｕｔｈｏｒ:ＣＨＥＮ Ｓｏｎｇ￣ｙｉｎｇꎬＥ￣ｍａｉｌ:ｃｈｅｎｓｙ６６＠ｓｄｕ.ｅｄｕ.ｃｎ) Ａｂｓｔｒａｃｔ: Ｅｍｐｌｏｙｉｎｇｔｈｅ Ｄ２Ｑ９ ｍｏｄｅｌꎬｇｒｉｄｓ ａｔ ｔｈｅ ｄｏｕｂｌｅ ｓｉｎｇｕｌａｒ ｃｏｒｎｅｒｓ ｗｅｒｅ ｌｏｃａｌｌｙ ｒｅｆｉｎｅｄ ｔｏ ｉｎｖｅｓｔｉｇａｔｅ ｔｈｅ ｆｌｏｗ ｐｅｃｕｌｉａｒｉｔｙ ｏｆ ｔｈｅ ｌｉｄ ｄｒｉｖｅｎ ｃａｖｉｔｙ ｂａｓｅｄ ｏｎ ｔｈｅ ｍｕｌｔｉ￣ｒｅｌａｘａｔｉｏｎ ｔｉｍｅ ｏｆ ｌａｔｔｉｃｅ Ｂｏｌｔｚｍａｎｎ ｍｅｔｈｏｄ (ＭＲＴ￣ＬＢＭ). Ｔｈｅ ｔｒａｎｓｆｅｒｒｉｎｇ ｆｏｒｍｕｌａｅ ｏｆ ｔｈｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ｆｕｎｃｔｉｏｎｓ ｗｅｒｅ ｅｓｔａｂｌｉｓｈｅｄ ｂｅｔｗｅｅｎ ｔｈｅ ｃｏａｒｓｅ ａｎｄ ｆｉｎｅ ｇｒｉｄｓ. Ｐａｒａｍｅｔｅｒｓ ｗｅｒｅ ｉｎｔｅｒｐｏｌａｔｅｄ ｗｉｔｈ ｃｕｂｉｃ ｓｐｌｉｎｅ ｉｍｐｌｅｍｅｎｔａｔｉｏｎ ｏｎ ｔｈｅ ｒｅｆｉｎｅｄ ｌａｙｅｒｓꎬ ａｎｄ ｔｈｅ ｓｔｒｅｓｓ ｅｘｐｒｅｓｓｉｏｎｓ ｗｅｒｅ ｇｉｖｅｎ ｏｎ ｔｈｅ ｂａｓｉｓ ｏｆ ｔｈｅ Ｃｈａｐｍａｎ￣Ｅｎｓｋｏｇ. Ｆｏｒ ｔｈｅ ｃａｖｉｔｙ ｆｌｏｗｗｉｔｈ ｔｈｅ Ｒｅｙｎｏｌｄｓｎｕｍｂｅｒ １０００ꎬｔｈｅ ｓｉｍｕｌａｔｉｏｎｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ ｔｈｅｖｅｌｏｃｉｔｙｄｉｓｔｒｉｂｕｔｉｏｎａｌｏｎｇｔｈｅｃａｖｉｔｙｍｉｄｐｌａｎｅａｇｒｅｅｓｗｉｔｈｔｈａｔｏｆｔｈｅｃｌａｓｓｉｃａｌｂｅｎｃｈｍａｒｋｖｅｒｙ ｗｅｌｌꎬｔｈｅｎｏｉｓｅｓｏｎｔｈｅｐｒｅｓｓｕｒｅａｎｄｖｏｒｔｅｘｃｏｎｔｏｕｒｓａｒｅｓｉｇｎｉｆｉｃａｎｔｌｙｒｅｄｕｃｅｄꎬａｎｄｔｈｅｓｔｒｅｓｓｆｌｕｃｔｕａｔｉｏｎ ｏｂｖｉｏｕｓｌｙ ｄｅｃｒｅａｓｅｓ. Ｔｈｅ ｃｏｍｐｕｔａｔｉｏｎａｌ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄｖａｌｉｄａｔｉｏｎ ｏｆ ｔｈｅｐｒｅｓｅｎｔ ｍｅｔｈｏｄ. Ｋｅｙ ｗｏｒｄｓ: ｌａｔｔｉｃｅ Ｂｏｌｔｚｍａｎｎꎻ