羽流中固体颗粒在真空环境下分布的数值仿真

通过对气?固两相间动量和能量相互作用解耦处理,建立了一种适于模拟真空环境气固两相混合物羽流的DSMC双向耦合算法与固体颗粒空间输运变化特性的TPMC计算技术。仿真了固体火箭发动机两相羽流流场和固体颗粒在远离发动机喷口数十公里空间扩散运动分布特性,通过将计算结果与典型文献结果及理论分析比较确认,证实本文方法的准确可靠性。结果表明,固体颗粒对气相的扩散有一定的阻滞作用;仅在离发动机喷口一定距离,气相对固体颗粒有较大影响,致颗粒温度下降、速度增加;颗粒温度随发动机喷口距离增大而减小,一直要在远离喷口上百公里颗粒温度才随轴向位置趋于平衡,且不同尺寸的颗粒温度差别较大,对指导外层空间高真空环境气固两相羽流传输影响工程研制具有重要意义。     

一种稀薄两相流动的数值模拟方法

从描述粒子运动的微观层次出发,采用双向耦合技术,建立了一种适用于稀薄条件下两相流动的DSMC数值模拟方法。对相间相互作用进行解耦处理,实现了气固两相间动量和能量相互作用的模拟。相关算例验证了方法的可靠性。结果表明:在大样本统计平均条件下,该模拟方法能够保证相间动量和能量交换时守恒,固体颗粒尺寸、气体稀薄程度等条件对相间的相互作用有很大的影响。     

基于修正N-S方程本构关系的气体动理论耦合方法

随着稀薄程度的增加,Navier-Stokes方程的线性本构关系难以正确描述稀薄气体输运特性,高阶非线性本构关系往往数学形式极为复杂,对数值求解造成稳定性差等问题。为了发展适宜于近空间飞行器气动特性分析的高超声速稀薄流动模拟方法,本文利用求解Boltzmann模型方程的气体动理论统一算法(Gas Kineitc Unified Algorithm,GKUA)对应力张量、热流等宏观量数值积分求解的优势,提出了一种基于数值修正N-S方程本构关系的气体动理论耦合方法。通过将GKUA获得的应力张量及热流用于修正N-S方程的本构关系,实现了存在局部稀薄效应的流动模拟,并且通过可压缩平板边界层、圆柱绕流问题的数值模拟,验证了方法的有效性。     

分区对接网格在跨流域气体运动论统一算法的应用研究北大核心CSCD

为了模拟跨流域复杂物形绕流或多体流场,需要在位置空间建立能描述整个流场特性的网格系统,相应的数值计算方法也需要有所改进。本文采用分区对接网格处理复杂物形绕流或多体干扰流场,得到能准确描述流场特征的贴体网格系统。在此基础上数值求解考虑转动非平衡影响的Boltzmann-Rykov模型方程,在速度空间应用离散速度坐标法、位置空间采用数值差分NND格式。网格对接面上考虑分布函数的传递,包含网格点信息和网格点上速度的信息。通过计算得到复杂物形的流动特性。相关算例表明网格分区之间的流场参数光滑过渡。与文献及DSMC结果的对比证实本文基于分区对接网格的气体运动论统一算法是可靠的。