三维热化学非平衡流动非结构网格DSMC方法及其应用

研究了过渡流域三维热化学非平衡流动DSMC方法实现的过程。以四面体非结构网格为基本单元，提出一种新型的高效搜索算法，该算法不仅可以快速跟踪模拟分子在网格之间的迁移，而且可以准确判别分子与物面是否相互作用，避免了原有算法中分子表面反射的非确定论判据。设计了适合三维DSMC方法的动态局部时间步长技术，节约了计算时间。将碰撞距离的思想引人到非结构网格上来，有效地消除了网格尺度小于三分之一平均自由程的限制。利用Fortran90的动态分配内存技术编制了适用于任意外形的通用计算程序。最后对全尺寸航天飞机高超声速绕流进行了数值试验，计算结果验证了该算法的可行性及高效性。     

高超声速钝体尾流的积分电子密度与尺寸度律

本文提出了一个高超声速非烧蚀钝体化学非平衡尾流流场的简化模型和简化的空气化学模型。应用该模型计算尾流的积分电子密度。计算结果与实验数据符合得很好。在此基础上,以P_∞D_N=375mmHg·mm为例,初步论证了尾流双尺度律的适用程度与范围,为弹道靶模型实验提供正确的导向。     

高温圆管内多孔介质耦合换热的辐射作用

通过数值模拟圆管内填充多孔介质的辐射对流耦合换热,研究高温固体骨架辐射效应对温度分布及换热的影响.基于局部非热平衡模型分别建立流体和固相能量方程,采用蒙特卡罗法求解固体骨架的辐射换热;对不考虑热辐射引起的温度场偏差和管壁发射率以及多孔结构参数的影响进行讨论.结果表明:多孔固体骨架的辐射效应对入口段温度场的影响明显,不考虑辐射将导致较大偏差,壁温为1500K时最大偏差为16%.管壁发射率对温度场影响较小,壁温为1500K时影响小于3%.孔隙率或孔径增大,壁面辐射热流密度比例增加,辐射效应体现明显.      

湿空气凝结对ONERA M6机翼气动特性影响的数值研究

大气飞行条件下机翼附近的跨声速湿空气流动中空气中含有的水蒸汽可能越过饱和线而发生非平衡凝结。水蒸汽凝结潜热加热气流会改变跨声速流动的特性,从而对机翼气动特性造成显著影响。本文建立了湿空气非平衡凝结流动的数学物理模型,对ONERA M6机翼在跨声速条件下的湿空气非平衡凝结流动进行了分析。结果表明,与干空气流动相比,在攻角为3.06°和空气相对湿度为50%时,ONERA M6机翼表面压力系数有显著变化。造成机翼气动特性显著变化的原因在于:湿空气中水蒸汽凝结放热对跨声速气流加热,导致机翼表面附近的流速、压力与流场结构发生了显著变化。     

ESA再入探测器的某些空气热动力问题

给出了两种不同的再入飞行器:充气的再入和下滑演示器(IRDT)和Huygens.两者的任务都涉及在靠近轨道的条件下再入并均已于2005年完成.概述了IRDT的使命和具体设计问题,综述了IRDT最近飞行的初步结果和ESTEC估计Huygens受到的辐射通量采用的方法.     

超声速非平衡电离磁流体动力技术实验系统

介绍了超声速非平衡电离磁流体动力技术实验系统的设计思想和组成,设计、制作了马赫数为3.5的吸气式双喉道风洞,采用以陶瓷板为阻挡介质的电容耦合射频阻挡放电,实现了超声速流场中大体积、连续、稳定、均匀等离子体的产生.主要结论有:实验测试风洞稳定工作时间为18s,稳定运行时实验段静压为650Pa;电容耦合射频放电典型工作状态下,通过伏安诊断法测出超声速条件下等离子体的电导率为1.27×10-3S/m.      

基于湍流输运特性对S-A模型在压气机角区流动模拟中的改进研究

在非设计工况下,压气机叶栅角区的三维分离加剧,流场内的湍流处于较强的非平衡状态,因此必须对于非平衡状态的湍流的输运过程给出准确的预估.基于不同工况下压气机内角区分离流动的湍流输运过程的分析,提出了一种修正Spalart-Allmaras(S-A)模型中模拟湍流黏性生成和耗散关系的改进方法,经过数值研究发现改进的S-A模型对于静子角区分离的模拟更接近物理实际,从而提高了对非设计工况的预估能力.      

高超声速非平衡流场多个振动温度模型的数值研究

利用数值求解热化学非平衡Navier-Stokes方程的方法,建立了考虑一个和多个振动温度的高超声速热化学非平衡流场的CFD计算程序,并对半球模型和球锥模型进行了数值计算,研究结果表明:(1)各个分子组分的振动温度分布于考虑一个振动温度模型的振动温度周围,即一个振动温度模型的振动温度是多个振动温度模型的各个分子组分振动温度的一种平均值;(2)在模型身部,振动温度的峰值高于平动温度的峰值;在模型近尾出现振动温度冻结并高于平动温度。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。