基于宏观方程数值本构关系的气体动理论加速收敛方法
Method of accelerating convergence for gas kinetic algorithm based on digital constitutive relation of macroscopic equations作者机构：中国空气动力研究与发展中心超高速所绵阳621000 英国思克莱德大学机械与航空工程系詹姆斯·维尔流体实验室格拉斯哥G11XJ 国家计算流体力学实验室北京100191
出 版 物：《物理学报》 (Acta Physica Sinica)
年 卷 期：2020年第69卷第20期
页 面：281-292页
核心收录：
学科分类：07[理学] 070104[理学-应用数学] 0701[理学-数学]
基 金：国家自然科学基金青年科学基金(批准号:11902339)资助的课题
摘 要：在跨流域复杂流动问题的模拟中,基于求解速度分布函数演化方程的气体动理论方法的效率问题一直受到工程应用领域关注.研究提升气体动理论方法在定常流动模拟中的计算效率具有重要意义.为了提升定常流动计算收敛速度,本文提出了一种耦合宏观方程数值本构关系的气体动理论加速收敛方法.通过求解Boltzmann模型方程,将应力、热流高阶项的数值解与宏观方程耦合,实现了宏观方程的封闭;另一方面,宏观方程的计算结果被用来更新Boltzmann模型方程的当地平衡态速度分布函数中的宏观物理量,以此构造求解Boltzmann模型方程的全隐式数值格式.通过跨流域方腔流动、超声速圆柱绕流及双圆柱干扰绕流案例的数值模拟,对方法进行了广泛考核.计算结果与常规气体动理论统一算法、直接模拟蒙特卡罗法符合良好,证明该方法很好地描述了稀薄流动中的非线性本构关系,以及激波、强壁面剪切、流动分离等强非平衡特征.进一步,对于低努森数Kn的流动,方法能显著加速收敛过程,提升计算效率;随着努森数Kn增加,气体对流输运效应减弱,方法的加速收敛效果降低.与此同时,如何减少内迭代耗时,进一步提升效率有待更多研究.
维普期刊数据库
同方期刊数据库