基于宏观方程数值本构关系的气体动理论加速收敛方法

摘要

在跨流域复杂流动问题的模拟中,基于求解速度分布函数演化方程的气体动理论方法的效率问题一直受到工程应用领域关注.研究提升气体动理论方法在定常流动模拟中的计算效率具有重要意义.为了提升定常流动计算收敛速度,本文提出了一种耦合宏观方程数值本构关系的气体动理论加速收敛方法.通过求解Boltzmann模型方程,将应力、热流高阶项的数值解与宏观方程耦合,实现了宏观方程的封闭;另一方面,宏观方程的计算结果被用来更新Boltzmann模型方程的当地平衡态速度分布函数中的宏观物理量,以此构造求解Boltzmann模型方程的全隐式数值格式.通过跨流域方腔流动、超声速圆柱绕流及双圆柱干扰绕流案例的数值模拟,对方法进行了广泛考核.计算结果与常规气体动理论统一算法、直接模拟蒙特卡罗法符合良好,证明该方法很好地描述了稀薄流动中的非线性本构关系,以及激波、强壁面剪切、流动分离等强非平衡特征.进一步,对于低努森数Kn的流动,方法能显著加速收敛过程,提升计算效率;随着努森数Kn增加,气体对流输运效应减弱,方法的加速收敛效果降低.与此同时,如何减少内迭代耗时,进一步提升效率有待更多研究.