燃气源复杂通道三维冷流数值模拟

采用两方程RNGκ-ε湍流模型，通过非结构网格上的SIMPLE算法，对某远程火箭弹控制系统燃气源复杂通道内流场进行了冷流数值模拟，获得了通道内各气流参量的三维分布，了解了燃气源内各处特别是控制气流流经的迷宫形多层导热套环结构内的气流流动特性，计算结果与实验结果基本吻合。     

塞式喷管冷流实验和数值模拟研究

开展了塞式喷管冷流实验和数值模拟研究。实验模型采用直排的点膨胀构型，实验测量了两种塞体长度的塞式喷管在不同压比条件下的表面压强分布，拍摄了典型工况下流场的纹影照片。数值模拟采用二阶NND格式求解二维N-S方程，计算得到的壁面压强分布与实验结果吻合较好。通过实验和数值模拟验证了塞式喷管高度补偿特性，得到的结论为塞式喷管设计和性能预报提供了依据。     

灭火弹气体发生器内流场数值模拟

为了提高灭火弹灭火性能,优化灭火弹结构设计,利用Fluent对灭火弹气体发生器内流场进行仿真计算,通过改变气体发生器开孔情况及发生剂参数调整弹内流场分布。计算结果显示:当首孔距离等于15倍气体发生器排气孔孔径时,灭火弹壳体近壁面存在5个大小相近且数值较高的压强峰,压力分布最利于灭火弹壳体预制孔槽破裂;灭火弹内气体发生器上排气孔对称性越强,发生器内压力分布越均匀。经过试验验证,仿真所得压力-时间曲线与试验测得数据吻合良好,说明该模型适用于灭火弹内流场仿真。      

非均匀超声来流矩形隔离段内流场实验

针对超燃冲压发动机隔离段的非均匀进口条件设计了隔离段实验风洞，通过测量隔离段壁面压力和拍摄流场纹影照片研究了矩形隔离段内激波/紊流附面层相干流场。研究发现，隔离段进口的非均匀流使隔离段流场压升特征与附面层发展规律与均匀进口的隔离段流动有显著差异。用截面当量直径取代Waltrup公式中的圆管直径可以取得较好的吻合效果。在进口马赫数小于2时，升高同样的压力，非均匀进口隔离段产生的激波串长度比Waltrup公式预测的长度要长。纹影仪观察发现隔离段内激波存在严重的振荡现象。     

低速高温燃气流热模拟试验方法和设备

对比分析了两种气流状态参数和两种加热情况下典型前缘部件表面热流密度的相似性,论证了利用亚声速高温燃气流加热方式进行近地空间高超声速飞行工况气动热模拟试验的可行性.针对高超声速飞行器典型钝头锥结构提出“小喷口低速高温燃气流+石英灯”组合热试验方案.通过采用新型高效双腔蒸发管型燃气发生器、新型带保温夹层和耐高温陶瓷内衬的水冷不锈钢高温管道结构,同时引入电加热器预热及燃烧室两路供油方案,使所建低速高温燃气流热试验设备产生燃气流温度达到2100K,φ250mm喷口处平均径向温度分布梯度约3K/mm,具有线性温度控制功能且稳态控制温差约46K,满足24km、马赫数为6典型高超声速飞行器工况驻点区域高温/大热流密度气动热试验要求.      

小型化催化重整煤油燃气发生器的数值模拟

催化重整燃气发生器是在线产生氢气的有效途径,有助于碳氢燃料超燃冲压发动机的着火和燃烧.通过对新设计的小型化燃气发生器乳化煤油的燃烧及催化重整进行数值模拟,在计算结果和试验数据吻合基础上,分析了出口燃气的总体状态和影响催化重整各参数的变化规律,确定了研究涉及到的参数范围内的最佳工况点,并揭示了内部的燃烧状态和一些难测量成分含量变化趋势,从而进一步指导设计和试验.      

提升机构微型燃气发生器的探索研究

介绍了一种宇航用的利用产生燃气做功并通过提升机构完成提升动作的微型燃气发生器的设计。建立了描述燃气发生器内部流场的数学物理方程。通过数值方法模拟了燃气发生器内部的燃烧和流动以及提升机构的运动过程。由模拟结果及理论分析来指导药柱的设计, 使得燃气发生器能满足高性能和高可靠性的要求     

离心式喷嘴内流场特性的数值模拟

基于单流体模型及气液两相混合分数的概念，建立了离心式喷嘴内流场的数学模型，并进行了数值计算，计算结果给出了离心式喷嘴内部液膜与气涡共存的流场结构及气液交界面的几何形状，计算得到的壁面压力分布、出口处液膜内轴向速度分布以及雾角大小都与实验结果吻合的很好，与常用的跟踪气液界面的ALE方法相比，本文建立的数值模拟方法具有简单，计算量小和易于实现的特点。     

高空模拟试验舱中燃气流动过程的分析与计算

以一维等熵流动理论为基础,在充分考虑激波效应的情况下,分析了高空模拟试验舱中发动机燃气可能出现的两种流动过程,并计算了燃气流动过程的压力变化曲线,其结果为高空模拟试验舱的结构设计和真空机组的选配提供了参考。     

变比热容对射流推力矢量喷管内流场影响的数值模拟

利用时间推进的有限体积法求解二维雷诺平均Navier-Stokes(N-S)方程,分别在变比热容和定比热容的情况下研究了射流推力矢量喷管的内流性能,并对两种情况下计算结果进行了比较.结果表明,与定比热容情况相比,工质为变比热容的燃气时,喷管内部由二次流产生的斜激波位置向喷管下游移动,回流区减小,喷管的推力矢量角和推力系数减小,温度越高,两种方法的计算结果相差越大.因此,在采用计算流体力学(computational fluidic dynamic,简称CFD)方法研究射流推力矢量喷管内流场性能时应该考虑变比热的影响.      

Experimental and numerical study on pressure oscillation in a combustor with rotary valve

A rotary valve coupled in a combustor assembly can generate periodic pressure oscillations inside the combustor and can be used to study the combustion instability in the combustion chamber of a rocket engine. This paper proposes a cold gas flow experimental system based on a rotary valve and a corresponding rotation model. A 3D numerical model is proposed to obtain the transient flow inside the rotary valve, and the dynamic mesh technique and User-Defined Functions (UDFs) are adopted to implement a swing motion instead of a rotary motion. Several cold gas flow experiments are carried out at rotating speeds of 75 225, and 375 rpm to verify the validity of the numerical model. The effects of rotating speed, stroke length ratio and radius of the RED (Rotor Exhaust Duct) on the pressure oscillations are studied using this numerical model. The results show that the maximum and peak-to-peak values of the pressure oscillations gradually decrease with increasing rotating speed. The ratio of the corresponding peak-to-peak value to the maximum pressure (pressure amplitude ratio) is reduced from 1.81 to 0.6%. The stroke length ratio affects the pressure waveform because it leads to a change in the time spent in the non-exhaust stage. When the stroke length ratio is 0, the waveform closely resembles a sine wave. With the increase in the stroke length ratio, the pressure waveform exhibits a more square or a triangular wave shape. Finally, a high-frequency and high-amplitude pressure oscillation can be obtained by appropriately increasing the radius of the RED.     

Experimental and numerical simulation of bird-strike performance of lattice-material-infilled curved plate

The anti-bird-strike performance of a lattice-material-infilled curved plate is investigated herein. Since automatically filling the curved structure by classical lattice material filling methods will cause a large number of manufacturing defects, a space-dependent lattice material filling method for the curved plate is firstly proposed in this paper Next, using a face-centered cubic lattice, a lattice-material-infilled test piece with a hollow ratio of 40.8% is built. The test pieces are manufactured via additive manufacturing using titanium alloy. In bird-strike experimental tests, the test pieces are crashed against gelatin birds at an impact velocity of 200 m/s. Dynamic strain gauges are used to record the crash history and the results are discussed. Furthermore, a numerical analysis to simulate the bird-strike experiment is performed. The results from the experimental tests and numerical simulation agree well. This work shows that the lattice-material-infilled curved plate yields promising bird-strike resistance. Therefore, lattice-infilled materials are feasible for protecting aerospace components against bird-strike as well as for reducing the component weight.     

二维涡轮级的定常/非定常数值模拟

为了解定常与非定常计算差别对工程应用参考价值的影响，本文通过求争二维Navier-Stokes方程，针对一级涡轮叶栅分别进行了定常／非常数值模拟。结果发现，对于叶片表面压力分布，第一排导叶定常计算与非定常计算得到一致的结果，但是在动叶上差别较大，尤其在叶排小间距时。当叶排间距拉大后，定常计算结果趋近于非定常计算结果时间平均值。     

底排装置空腔干扰下的火箭燃气流场数值模拟

从二维轴对称Euler方程出发,应用隐式有限体积TVD格式和分区技术数值模拟来流马赫数Ma=2 0,攻角α=0°时超声速绕流与底排装置空腔内的火箭燃气流相互干扰的复杂流场,计算得到了清晰的流场波系结构,包括底排装置空腔中的激波,火箭燃气流在底排装置空腔外底部与外来绕流作用后形成的相交激波、反射激波等,其结果可为底排火箭复合增程弹丸的结构设计和气动特性研究提供理论依据。     

射流式旋涡发生器对涡轮流动分离控制

通过数值计算,研究了稳态小孔射流式旋涡发生器VGJs(Vortex Generator Jets)对低雷诺数涡轮流动分离的控制。通过对VGJs射流孔径、孔间距的几何参数研究,分析了射流对分离流动的控制机理。在此基础上,提出并详细研究了多种不同孔径射流孔组合排列结构对分离流动的控制效果与作用机理。研究结果表明,射流孔组合排列结构,以大孔射流为主,采用小孔小流量射流对大孔射流控制的分离流场进行"弥补"与调节,可实现总压损失进一步降低22%,射流流量减小25%。     

CRM-WB风洞模型高阶精度数值模拟

通过求解雷诺平均Navier-Stokes （RANS）方程,采用5阶空间离散精度的WCNS和多块对接结构网格技术,开展了CRM翼身组合体风洞试验模型的高阶精度数值模拟,计算构型和计算状态来自AIAA第6届阻力预测研讨会（DPW Ⅵ）。主要目的是采用高阶精度方法,评估静气动弹性变形和模型支撑装置对CRM翼身组合体数值模拟结果的影响。通过与刚性外形的计算结果和NASA NTF风洞试验结果的对比,高阶精度数值模拟结果表明：迎角为3.0°时,静气动弹性变形使得机翼上表面激波位置前移并显著降低了外侧机翼上表面激波位置前的负压;迎角为3.0°时,模型支撑装置使得机翼上表面激波位置进一步前移,并导致翼梢处机翼上表面流动结构发生变化;迎角为4.0°时,计算模型中没有包含模型支撑装置是导致升力系数下降的主要原因;计算模型中包含机翼静气动弹性变形和模型支撑装置的数值模拟结果更加接近试验结果。