贫油预混预蒸发燃烧室回火现象的大涡模拟

为研究贫油预混预蒸发(LPP)燃烧室的回火现象及形成机理,基于大涡模拟(LES)及部分预混燃烧模型,对Schneider试验进行数值计算,分别将冷态及燃烧温度场与试验结果进行对比,验证了模型的准确性。将该模型应用于LPP燃烧室,并利用详细的化学机理对航空煤油火焰传播速度进行修正,成功重现了回火现象的发生并向主燃级上游传递的动态过程,通过流场分析发现局部形成的不规则涡结构在回火现象中扮演着重要角色。进一步从涡量方程的角度对涡结构进行分析,发现燃烧引起的斜压及体积膨胀是诱发回火现象的主要原因。     

基于混合网格的三维Navier-Stokes方程并行算法

提出了一种基于混合网格的三维Navier-Stokes方程的并行计算方法。Navier-Stokes的求解采用了基于面的有限体积方法,该方法适用于任何网格类型。采用一方程Spalart-Allmaras模型来计算紊流黏性。并行计算采用区域分裂的方法,利用METIS网格分区系统实现了各节点的加载平衡。节点间的数据交换通过调用MPI库函数来实现,采用非阻断通讯的方式来减少数据交换时间。充分利用FORTRAN90的动态存储特性来减少对内存的需求。最后,通过对绕DLR-F6外形(翼身组合体+挂架+发动机短舱)黏性流动的数值模拟,验证了该并行程序的准确性,高性能并行计算以及处理复杂几何外形的能力。     

非均匀网格湍流大涡模拟高精度有限体积解法

为准确预测不可压复杂湍流,提出了一种可用于大涡模拟均匀或非均匀网格上的高精度有限体积法。该方法空间离散采用有限体四阶紧致格式,时间推进采用四阶Runge-Kutta法,压力-速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动证实了该方法具有近四阶的空间精度;并在此基础上,采用动态Smagor-insky亚格子应力模式,成功地实现了充分发展槽道湍流和后台阶湍流流动的大涡模拟计算,所得结果与直接数值模拟结果吻合良好,且采用非均匀网格可在比均匀网格数少的离散系统上得到同样满意的结果。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

三维非结构贴体直角网格的N-S方程数值模拟

基于八叉树数据结构,采用直接分层光顺的方法,解决了非结构直角网格在物面不正交贴体的难点,发展了任意复杂外形的三维非结构贴体直角网格,并且可以根据流场特征进行自适应处理。在此基础上采用最小二乘法对N-S方程进行了MUSCL有限体积法的非结构网格的粘性求解,提出了带约束边界条件的最小二乘法的处理办法。计算表明,三维非结构贴体直角网格的粘性求解,可以适应各种高超声速的复杂外形,计算结果与试验吻合较好。     

机载中空纤维膜富氮性能实验

通过建立机载中空纤维膜环境模拟性能实验平台,对与国内某技术中心合作研制的机载中空纤维膜进行了系统的环境模拟性能实验及分析.实验结果显示:①渗余富氮空气氧体积分数随输入空气压力的上升而下降,随输入流量的上升而上升,而与输入空气温度基本无关;②渗余富氮空气氧体积分数随富氮空气流量的上升而上升,且输入空气压力低时明显;③海拔高度对机载中空纤维膜富氮性能基本不影响,仅启动时渗余富氮空气氧体积分数略有下降.通过采用多元回归方法,得到有一定普遍意义的机载膜富氮性能公式.      

体积分数法测量几种转静系中的最小封严流量

为防止涡轮盘腔转静系统高温“燃气入侵”现象发生,需要确定最小封严流量.在不同的转速和不同出气间隙下,采用体积分数法分别对6种典型封严结构进行实验研究,以确定它们在不同工况下的最小封严流量并得出其变化规律:在固定的出气间隙值下,最小封严冷气流量随转盘转速的增加而增大,且几乎呈线性变化;对于固定的转速,随着出气间隙的增加,最小封严冷气流量也随之增加.拟合出不同封严结构的级间封严准则函数关系式,并探寻出封严效果最优的封严结构.      

姿控发动机羽流的可见光/近红外特性

为了研究发动机羽流辐射特性对光学探测成像的影响,基于不同直径粒子散射理论、辐射传递方程和分子吸收线数据库,应用有限体积法,建立了火箭发动机羽流在不同波段的辐射特性计算程序.针对实际工作情况进行建模,开展某姿控发动机工作在环月球轨道时的可见及近红外波段羽流的辐射特性计算.在0.4~0.9μm波长范围内,针对羽流气体组分、波长、观察天顶角以及太阳辐射对羽流光谱的影响作了数值计算及分析.研究结果表明:应用有限体积法开发的程序较好地模拟火箭发动机羽流的辐射特性,并具有广泛适用性;在可见光和近红外波段,太阳辐射作用对火箭发动机喷流辐射特性影响最大.      

压气机进气畸变数值模拟技术研究

发展了一种针对进气畸变条件下的风扇/压气机进行性能预估和稳定性分析的计算方法。首先研究了将叶片作用力简化为体积力源项的建模方法,在此基础上开发出一套基于体积力的三维进气畸变数值模拟程序,使用该程序对NASA Rotor 35在均匀进气、进口存在稳态总压畸变及同时存在总压和总温畸变的流场进行了模拟分析。结果表明,该程序获得的压气机特性及参数分布与雷诺平均Navier-Stokes(RANS)计算吻合得很好,同时正确地模拟出了压气机转子与上游畸变来流的耦合作用及其对压气机性能和稳定工作裕度的影响。     

热压罐工艺仿真技术

随着FEM和CFD仿真手段的发展,利用仿真手段替代部分试模,预报试模的结果已成为可能。通过仿真手段可以模拟热压罐工艺过程中罐内的流场情况、温度场分布、预浸料的固化过程,以及最终工件的变形和残余应力等。而在进行了大量的虚拟仿真试验之后,则可利用神经网络建立热压罐工艺的知识库和专家系统,从而指导工装工件摆放、工装设计以及诸多工艺参数的优化,从根本上改变热压罐工艺方案的设计方式。     

HLLC方法在三维非结构网格上的应用

以三维非结构网格的显式有限体积法为基础，采用格心方法以及基于近似黎曼解的二阶精度Godonov扩展方法求解三维Euler方程，使用HLLC(Harten，Lax，van Leer,Contact)近似黎曼解的方法计算网格单元边界处的守恒量通量。为了验证方法的可行性，利用三维解算器对NACA0012翼型绕流进行了数值计算，结果令人满意。在此基础上对某型弹的流场进行了数值模拟，取得了比较好的结果。