稀薄流过渡区气固两相喷流的建模与数值模拟

从描述粒子运动的细观层次出发,运用DSMC方法,构造了稀薄流过渡区气固相互作用模型和颗粒碰撞模型,数值模拟了有气固两相横向喷流干扰的二维有限平板高超声速流场,分析了流场结构和物面气动系数分布等特性.结果表明,对于不同的来流Knudsen数和颗粒质量载荷等参数,喷流干扰流场存在差异,同时反映了基于DSMC方法的气固相互作用模型是研究稀薄及过渡条件下的气固两相流动的有效方法.     

热辐射计内稀薄气体流动的蒙特卡罗模拟

本文应用直接仿真蒙特卡罗(DSMC)方法,研究了Crookes热辐射计内由于叶片两表面的温度不同而引起的稀薄气体的流动,给出了热辐射计叶片高温面温度900K,低温表面温度为300 K,四种不同克努森数(Kn)下流速矢量场、温度场和叶片两面压力差的仿真结果.并与热辐射计在不同Kn数下的转速实验测量结果进行了比较,两者结果基本一致.     

稀薄气体圆柱绕流的并行化DSMC方法研究

本文针对ＤＳＭＣ方法的内在并行特征构造了两种并行算法,并在ＰＶＭ环境下实现了并行化计算。计算结果表明,实现ＤＳＭＣ方法的并行化能够显著缩短计算时间,提高计算效率,不同并行算法能够达到的加速比存在差异。     

高超声速稀薄流的气粒多相流动DSMC算法建模研究

基于直接模拟Monte Carlo(DSMC)方法,构造适用于DSMC算法的固态和液态颗粒碰撞、聚合和分离模型,发展稀薄条件下双向耦合作用的气粒多相流的DSMC算法,在此基础上初步实现高超声速稀薄流环境中的气粒多相喷流流场数值模拟.算例结果表明该方法能为稀薄过渡区气粒多相流动提供一种新的应用研究手段.     

固体火箭发动机内气粒两相流动的SPH-FVM耦合方法数值模拟

为解决现有数值模型在求解固体火箭发动机内颗粒的燃烧流动等问题时所遇到的难题,基于拟流体模型,采用光滑粒子流体动力学方法 (SPH)求解离散颗粒相,追踪颗粒运动轨迹,采用有限体积方法 (FVM)对气体连续相进行描述,捕捉流场特性,两相间通过曳力、压力梯度、热传导、体积分数等参量进行耦合,建立了两种不同坐标系下方法间的耦合框架。耦合方法模拟了喷气推进实验室(JPL)喷管中气粒两相流动过程,得到了气相和颗粒相的参数分布规律,与相关实验结果及离散颗粒模型(DPM)数值模拟结果对比,吻合较好。针对Φ315mm实验发动机工作状态和结构特点,进行了发动机燃烧室内过载条件下的气粒两相流动数值模拟,分析了纵、横加速度载荷对发动机燃烧室内粒子场和聚集带的影响,与已有的数值模拟结果基本吻合。结果表明,新方法求解发动机气粒两相流问题准确可靠,适用于发动机内更为复杂的两相流问题数值模拟。     

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

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

稀薄过渡流区横向喷流干扰效应数值模拟研究

为了准确预测稀薄过渡流区横向喷流与稀薄大气的干扰流动特征,本文在建立直角与表面非结构网格混合结构的DSMC数值算法以及碰撞网格自适应算法基础上,提出基于MPI的静态随机负载平衡技术,构建了DSMC并行计算代码。计算的不同压力比条件下三维平板模型横向喷流与稀薄大气干扰的分离长度与低密度风洞试验有较好的一致性,验证了本文算法的可靠性。开展了细长钝双锥外形高超声速稀薄来流与超声速/高超声速横向喷流干扰效应的研究,计算分析了不同飞行高度、不同飞行速度、不同飞行攻角、不同喷流推力下复杂流场结构和对气动力特性的影响规律。考察了RCS喷管出口参数不同(均匀/非均匀)对喷口附近分离涡和分离长度的影响。     

可压稀相两相流场的数值模拟

针对可压两相流动，在阻力系数中考虑了可压缩性及颗粒聚积所造成的影响，用TVD格式对沙层，平板及JPL喷管内的气固稀相两相流进行了数值模拟，所获结果和已有的实验数据及有关文献吻合很好，在JPL喷管中，随着质量分数的增加，由于气体和颗粒间的相互作用逐渐增强，气相的温度增大，马赫数减小，但压力变化很小，而随着颗粒直径的增加，气相轴线上的温度减小。     

叶片前缘两相流冲击冷却的耦合数值模拟

为了更好地强化叶片前缘的换热,对叶片前缘进行两相流强化冲击冷却的耦合数值模拟研究。利用水滴/蒸汽冲击曲面冷却系统进行数值模拟程序的校核。在此基础上,利用欧拉-拉格朗日颗粒追踪方法对叶片前缘的两相流冲击冷却进行耦合数值研究。研究了水滴加湿量、水滴直径和水滴与壁面的作用边界条件对换热性能的影响。结果表明当加湿量从3%增加到8%时,最大的平均传热增强系数从1.671增加到4.913;水滴直径从5μm增加到20μm时,最大的平均传热增强系数从2.128降低到1.164;反弹和反弹+破碎边界条件预测的传热效果好于粘附和粘附+破碎边界条件。     

大瑞利数下二维Benard对流演化过程的DSMC仿真研究

利用DSMC（Direct Simulation Monte Carlo）方法研究了大Rayleigh（瑞利）数条件下稀薄气体的二维Benard对流的特性。考察了流动的非定常演化过程及若干流动特征量随时间的变化特征，得到了与实验较为一致的结论。计算结果表明，当流动控制参数不断增大时，二维Berard对流从热传导状态过渡到热对流状态表现出明显的非定常演化特征。最后当瑞利数增大到一定程度后，注动显示出混沌的特性。     

电离对高超声速稀薄流飞行器气动热影响

将电离反应模型扩展到（direct simulation Monte Carlo,DSMC）方法中,研究了电离反应效应对高超声速稀薄流飞行器气动热的影响特性.针对稀薄流场中电子出现带来的实际困难,引入“捆绑法”思想处理电子在流场中的运动,并给出了电离反应模型及电离反应处理方法.在以RAM-C Ⅱ飞行器外形为例对增加了电离反应的DSMC代码进行验证的基础上,以“星尘号”探测器外形为研究对象,针对不同飞行高度下5组元混合气体模型（无电离）和11组元混合气体模型（含电离）的化学非平衡流动开展了数值模拟,细致分析和对比了电离反应效应对探测器气动热的影响规律.研究结果表明:采用的电离反应处理方案能够模拟带电离反应的高超声速化学非平衡稀薄流动.在飞行高度为60km时电离反应对探测器气动热的影响最为强烈,使探测器的驻点热流密度降低了5.12%,电离反应对探测器气动热的影响随气体稀薄程度增加而减弱.      

稀薄气体条件下Rayleigh-Benard系统稳定性的DSMC仿真研究

本文利用DSMC（Direct Simulation Monte Carlo）方法研究了稀薄气体条件下Rayleigh-Benard系统的稳定性问题。考察了在经典的Boussinesq假设不再成立时，Rayleigh-Benard系统稳定性问题的特征。研究指出在气体的稀薄条件下，初场密度梯度的方向和Rayleigh数是影响Rayleigh-Benard系统不稳定性的主要因素。     

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

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

过渡流区N-S/DSMC耦合计算研究

在已有CFD和DSMC方法程序研究基础上,采用MPC耦合处理原理,发展适于流场分区信息交换的亚松弛技术,建立了N-S/DSMC耦合算法.通过对近连续过渡流区不同Kn数条件下钝锥绕流计算与低密度风洞试验结果比较分析,一方面证明了本文所建立的N-S/DSMC耦合计算模型和方法在低Kn数过渡流区的有效性,另一方面剖析了N-S/DSMC耦合计算误差,可望发展一个工程适用的近连续过渡流区高超声速飞行器气动特性预测分析工具.     

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

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

Maxwell气固相互作用模型对稀薄高超声速凹腔绕流流动特征和热环境的影响

针对高空航天飞机等再入飞行器表面缺陷或防热瓦缝隙导致的局部压力过高和气动加热问题,采用直接模拟Monte Carlo （DSMC）方法研究了飞行高度为80 km的稀薄流区高超声速凹腔绕流问题,考虑气固相互作用（GSI）模型对凹腔流场特征和表面压力、热流的影响。结果表明：稀薄流条件（80 km）下,GSI为完全漫反射时,在凹腔前缘分离的剪切层再次附着在后缘,在凹腔内部形成一个充满腔体的单涡结构;随着GSI从完全漫反射向镜面反射变化,气体与凹腔表面之间的切向动量交换减弱,即黏性剪切作用减弱,外部气流被卷入凹腔的程度减弱,导致涡结构不断减小直至消失,凹腔底部逐渐出现所谓的"死水区"。与完全漫反射相比,镜面反射或近镜面反射会导致凹腔上游侧面的峰值压力和峰值热流以及下游侧面的峰值压力剧烈增大,在飞行器设计中,应特别留意上述表面的压力载荷和热载荷。