卫星姿控发动机混合物羽流场分区耦合计算研究

研究求解喷管内流场N-S方程数值计算方法,发展基于N-S方程物面边界滑移流理论计算技术。提出求解羽流核心区轴对称DSMC模拟方法与远场三维DSMC仿真方案,发展多组元混合物羽流DSMC仿真方法。研究求解卫星姿控发动机内外近场、远场、倒流区和物面相互作用影响区多流域流场分区耦合计算技术,建立了一套用于求解混合物燃气羽流及对太阳电池帆板与卫星体表面撞击污染影响数值模拟方法。通过对分别安装于某在轨卫星不同位置两个典型姿控发动机燃气五组元混合物羽流计算研究及相关结果对比分析,证实本文数值方法可靠性。     

姿控发动机高空羽流流场干扰效应的DSMC方法研究

在研究权函数设置及热力学非平衡模型的基础上，将DSMC方法应用于求解三维问题，数值模拟了轴线平行的两姿控发动机喷流干扰形成的高空羽流流场。结果表明，过渡区流场中存在显著的干扰效应，采用将单个流场物理量简单相加的方法是不合适的。     

Hall推力器羽流数值模拟

以稳态等离子体推力器(SPT)羽流流场为研究对象,基于直接蒙特卡罗模拟与粒子单元混合方法建立数值模型,模拟等离子体稳态流动,研究流场特性。将计算结果与实验数据相比较,并与同类推力器羽流参数进行比较,以保证仿真可靠性。研究表明,所建立的仿真模型能够较好地预测SPT羽流电子数密度和离子电流密度等参数。     

高空羽流场的DSMC计算和实验研究

以小推力发动机的高空羽流场为研究对象，实现了高空小喷管氮气流DSMC方法数值模拟，验证了DSMC方法在高空膨胀流中的可用性，并在现有设备基础上，进行一组高空羽流试验，给出试验与计算的对比结果。结果表明：喷管出口附近能得到满意结果，由于实现设备能力所限，远场试验没有达到预期效果。     

姿控发动机内流场及高空羽流流场的DSMC一体化数值模拟

通过引入“平衡抽样方法”有效解决了DSMC求解近连续流区所遇到的碰撞顶计算效率问题;对于过渡流及自由分子流区,构造了适合于DSM仿真的分子碰撞传能模型并配置了一种光滑过渡以减少误差的权函数。在此基础上首次实现了对姿控发动机喷管内流场及包括倒流区在内的高空弱流流场的DSMC方法一体化数值模拟。     

高超声速平头圆柱绕流化学反应流场的DSMC/EPSM混合算法研究

本文发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的EPSM耦合模型，并通过混合参数进行流区的自动识别，将EPSM方法与蒙特卡罗直接模拟方法(DSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的轴对称情况下高超声速平头圆柱绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

模拟真空羽流场的特征线法

介绍了特征线法（MOC法）对真空羽流场进行数值模拟的过程，用MOC法进行了喷管出口流动不均匀情况真空羽流场算例的求解，所得结果与理论分析相符合，与DSMC法计算结果吻合，提高了计算效率，节省了计算时间。结果表明：用MOC法可代替DSMC法对连续介质区进行模拟，并可在用N-S和DSMC耦合数值模拟中确定耦合边界。     

基于DSMC方法的真空羽流不确定量化研究-空间推进专刊*

真实流动环境下的真空羽流必然存在着多种不确定性,那么确定条件输入的数值模拟会存在偏差,因此需进一步研究不确定性对羽流流动特征的影响规律。本文采用直接模拟蒙特卡洛（DSMC）方法对不确定性输入的羽流流场进行模拟,采用稀疏的概率配置点方法对来流、壁面及模型参数等输入不确定量进行描述,对不确定性的传播和输出目标的平均值、方差及不确定度进行计算。研究表明,流场不确定性沿流线传播至流速最大处之后迅速增强,并在声速线前出现骤减的现象,传播至声速线之后,挡板壁面输入不确定性的影响凸显。最突出的是,压力不确定度在挡板驻点位置达到全场最大值,约为输入不确定量（3.54%）的2.1倍。此外,温度跳跃不确定度受到壁面温度不确定性输入的限制而近似保持为一个恒定值,约为输入不确定度的0.8倍。进而,壁面热流不确定度（5.54%）比壁面正应力不确定度（6.25%）略小,切应力不确定度最小（5.07%）。Sobol’全局敏度分析表明,喉道速度和喉道压力的输入不确定性对气动力/热不确定度的贡献是最大的,且远远超过了壁面温度和模拟分子直径不确定性输入的贡献。     

基于PIC/MCC/DSMC方法霍尔推力器热分析

以提高霍尔推力器性能和使用寿命为目的对霍尔推力器进行了热分析。建立了等离子体在通道陶瓷壁面和阳极的能量沉积计算模型,并将模型耦合到PIC/MCC/DSMC流场计算程序中,计算加速通道壁面上的能量沉积分布。将计算得到的壁面能量沉积作为霍尔推力器温度计算的热流边界条件,考虑结构间的热传导与热辐射,计算温度分布。为了考察壁面的热流分布方式不同对温度场的影响,在总热流相同的情况下,将通道热流假定为线性分布与平均分布,计算霍尔推力器的温度分布。结果表明,壁面能量沉积占总功率的20.4%,陶瓷壁面能量沉积沿轴向位置先增大后减小,最大值在电离区。阳极的能量沉积,沿径向在中间位置达到最大,两端靠近壁面处较小。三种热流边界条件下的温度结果表明,三者最高温度都为700K左右,但高温位置不同,PIC热流边界条件下,最高温在电离区,而线性与平均热流边界条件下,高温区分别在出口区与近阳极区。通过与实验结果比较表明,PIC计算热流边界条件下高温区温度与测量者吻合更好,误差小于0.8%。     

壁面辐射平衡DSMC方法及其在双锥构型中应用

临近空间高超声速飞行器在高空长时间飞行,受气流加热影响,飞行器表面温度显著升高,依赖地面试验和传统DSMC仿真预测的热流值明显高于飞行观测值,导致飞行器防热系统的保守设计.本文发展了一种基于壁面辐射平衡的DSMC边界模型,通过热流值反算辐射平衡壁面温度,并以此温度作为下一个时间步DSMC计算的边界条件,迭代更新至给出壁面温度的收敛值.基于该温度边界条件,开发了适用于轴对称构型的DSMC求解器,并以钝锥构型对计算模型和求解器进行了验证.重点针对激波风洞试验条件下的双锥构型,开展数值模拟研究,结果表明:该构型恒温冷壁条件得到的壁面压力分布和热流与风洞试验结果吻合,两种温度条件下的压力峰值差异约为15.4％,但是整体气动力特性差异仅约为0.33％;相对于冷壁,辐射平衡计算得到的前缘处热流峰值降低约50％,再附点处的热流峰值降低约三分之二;两种条件相结合,可以给出壁面热流的预测范围.