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

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

二维过渡区含化学反应的喷流干扰流场DSMC/EPSM仿真

构造可仿真含化学反应的稀薄／连续不同流动区域的DSMC／EPSM混合算法，应用该算法，数值模拟过渡区有横向喷流干扰的二维有限平板高超声速化学反应流场，分析流场结构和物面气动系数分布等特性，充分反映稀薄／连续混合区域高超声速流动中的化学反应效应，同时通过与DSMC算法比较，论证DSMC／EPSM混合算法的可行性与计算效率相对提高的优越性。     

未来航空空气动力学的发展(上)

介绍了与航空空气动力学相关的技术,包括湍流理论、涡结构、转捩和分离机制及主动涡控制技术,大型高速民航机和军用运输机减阻技术,可重复使用的高超声速空天飞行器以及与未来微型航空器相关的昆虫飞行动力学.     

基于简化Boltzmann方程的超声速流稀薄效应分析

从近连续流到近自由分子流,计算分析超声速流场的稀薄效应。模型Boltzmann方程先后经简化分布函数和离散速度坐标变换后,采用一个隐式通量修正二阶迎风TVD格式差分求解。壁面取漫反射气体分子模型。在Knudsen数为0.001,0.01,0.1,1.0,10等5种情形下,数值模拟二维圆柱氩气体绕流,观察到了不同的波系及尾迹结构。总阻力系数计算值同实验数据基本吻合,碰撞项Shakov动力学模型反映的流场细节略好于修正BGK模型。稀薄效应趋于弱化激波等超声速流动结构,物体影响域增大。     

稀薄效应对可倾瓦动压气体轴承性能的影响

以微型三可倾瓦动压气体轴承为研究对象,采用连续模型、一阶滑移模型和WU新滑移模型速度边界条件,建立考虑稀薄效应的滑移修正雷诺方程。结合牛顿迭代法和有限差分法求解修正雷诺方程,计算考虑稀薄效应的三可倾瓦动压气体轴承各瓦块压力分布和承载力。计算结果表明:三可倾瓦动压气体轴承偏心率越大,转速越大,压缩系数越大,轴承承载力越大;考虑稀薄效应后,计算出的可倾瓦动压气体轴承承载力明显下降,且随着克努森数的增大,WU新滑移模型计算得到的承载力明显低于一阶滑移模型。     

返回器稀薄区气动特性工程计算方法的应用研究

中国航天返回器再入过程采用非弹道式再入轨道,其高空稀薄区滞留时间显著增长,稀薄效应对再入飞行的影响显著增强,因此快速而准确地预测返回器稀薄区气动特性变得非常重要。文章分析归纳了多种稀薄区气动特性工程计算方法即桥函数方法的发展和应用,并以 Stardust 返回舱为对象,对比分析了三种不同桥函数的预测精度,给出了在航天返回器气动预测中更为合适的工程方法。结果显示：不同桥函数预测结果差别很大;在不同攻角下,与数值模拟的对比分析表明,局部桥函数方法气动特性预测结果基本优于其它桥函数,尤其是在力矩系数预测上。因此稀薄区气动特性的预测采用局部桥函数较为合适。     

一类二维非结构网格DSMC方法的实现策略及其应用

研究了一类二维非结构网格DSMC方法的实现策略。在数据结构方面,设计了局部化的数据组织方式,节约了内存与计算时间。发展了一种跟踪模拟分子迁移的算法,该算法仅需少量的逻辑运算与代数运算,不仅可以快速跟踪模拟分子在网格之间的迁移,而且可以准确判别分子与物面是否相互作用,搜索过程中的附带信息给出了分子与物面碰撞的精确时间与位置,避免了重新计算。引入碰撞距离的思想,既减少存储又保证正确的模拟结果。在程序编制过程中,我们充分展现了Fortran90高级语言的主要特性,引入动态数组、指针、链表以及派生类型数据,编制了计算程序。最后对过渡流域高超声绕流进行了数值试验。     

旋转不变的数据驱动稀薄流非线性本构计算方法

DNCR方法基于Navier-Stokes方程和UGKS方法的数值模拟结果作为流场样本训练数据,借助机器学习构建热流/应力与流场特征参数的高维复杂非线性回归关系模型,最终通过耦合数据驱动的非线性本构关系求解宏观量守恒方程得到待预测稀薄非平衡流数值解。但现有DNCR方法的特征参数(速度、压力、密度等)和标记值(热流、应力张量等)不具备旋转不变特征,所得训练模型无法适用于坐标系旋转或平移后的计算网格。针对这一缺陷,本文构建全新且具有旋转不变性的特征参数与标记值,并结合典型算例预测结果与特征参数权重反馈优化所有训练集特征参数;同时瞄准回归模型预示范围与泛化性能提升,针对极端随机树开展选参与调参研究,发展了一种基于旋转不变量的改进DNCR方法。针对不同来流条件、不同几何外形条件下的二维高超声速圆柱绕流与顶盖方腔驱动流,评估了改进DNCR方法对比原始方法的计算精度提升效果。计算结果表明：使用旋转不变量能够显著提升训练模型对坐标系旋转、外形变化的适应能力,使DNCR方法具备更好的泛化性能。     

模拟空间环境下小发动机羽流压力场

为研究空间发动机的羽流特性及其对卫星的污染,在模拟空间环境下,进行了电加热CO2气体模拟发动机羽流流场的试验测量和理论计算。试验台主要包括空间环境模拟系统,电热气体模拟发动机,稳压气源,测量系统及其温控。理论计算采用计算流体动力学(CFD)和直接模拟Monte Carlo(DSMC)相结合的方法,完成了羽流场的数值模拟。测试结果表明,试验系统满足在10-3Pa量级的真空压力和93±5 K的背景温度下稳定运行和精确测量。不同坐标位置处试验数据和数值结果的比较表明,试验测量和理论计算吻合较好。