粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

固体火箭燃气射流近场形成与发展的数值模拟

固体火箭武器燃气射流的近场区域存在有复杂的波系，膨胀波、压缩波的混合作用是机-弹、弹-舰相容性研究的重要内容，复杂波系压力、温度场的分布是燃气射流动力冲击和热冲击作用的直接因素。因为这一区域流动的复杂性，使得实验研究非常困难。为此应用数值分析的方法，研究了固体火箭高度欠膨胀燃气射流流场起始冲击波、近场波系的形成及发展。计算结果揭示了近场区域的一些特性，与文献提供的数据较为吻合，为工程应用提供了帮助。     

扩张管和燃料部分填充对PDE效能的影响

采用基于有限速率基元化学反应模型的网格自适应有限体积方法,数值模拟以脉冲爆轰发动机为背景的爆轰波通过扩张喷管的流场,研究了喷管形状及喷管充气状况对脉冲爆轰发动机推进性能的影响。结果表明,氢-氧混合物产生爆轰波的等截面爆轰管冲量在某一比例合适的部分填充时达到最大,而对较大扩张角的喷管来说,其冲量一般随喷管内可燃气体填充比例的增加而增加,比冲通常随可燃气体填充比例的减少而增加;合适的喷管和燃料填充系数能使比冲显著提高。     

球面收敛调节片喷管RCS数值模拟研究

结合物理光学迭代法(IPO)和等效棱边电磁流法(EEC),自主开发了电大尺寸腔体电磁散射特性计算程序,并通过对某试验模型雷达散射截面积(RCS)进行计算,完成了程序的可靠性验证。利用该程序对带有电磁信号遮挡板且采用3种不同尾缘修形技术的球面收敛调节片矢量喷管(SCFN)进行了RCS数值模拟研究。结果表明:随着入射方位角的增大,喷管的RCS整体上有减小的趋势;喷管尾缘的修形只有在雷达波以较大方位角入射时对减小RCS效果明显;电磁信号遮挡板的存在很大程度上减弱了喷管的电磁信号特征。     

瓦状特征型面塞式喷管三维计算与实验

为了进一步了解瓦状塞式喷管的性能,采用NND差分格式求解三维N S方程和空气冷流对6单元瓦状特征型面塞式喷管进行了数值模拟和实验研究。研究模型的内喷管面积比为4,总面积比为40,设计压强比为1047。计算得到了流场马赫数和塞锥表面压强分布、喷管推力系数效率,以及不同压强比下中心平面、过渡平面和边缘平面的塞锥表面压强变化规律。计算结果与实验数据吻合得较好,效率数值最大相差1%。实验塞式喷管最大的推力系数效率为0 995,同钟型喷管相比,具有很好的高度补偿能力:从地面到高空,效率在0 93～0 995之间变化。和以前简化型面的4单元瓦状塞式喷管相比,实验和数值模拟均说明塞锥特征型面的优化设计提高了喷管性能,更充分体现了塞式喷管的高度补偿特性,可以成为未来工程应用的选择方案。     

等熵轻质活塞风洞中主要部件尺寸的确定

等熵轻质活塞风洞是一种工作时间很短的风洞,作为气源适合于瞬时热传导测量技术的研究工作。在研究了此风洞工作特性的基础上,提出了确定风洞中主要部件尺寸的方法,为风洞全面设计提供依据。压力比P0/P和trun是确定管道尺寸的主要参数,匹配压力振荡最大振幅的允许值是确定活塞质量的关键。     

双喷管差动式小推力测量装置的设计与分析

为了研究喷管型面对小推力火箭发动机推力性能的影响,设计并建设了双喷管差动式小推力测量装置及其附属设备.在测量装置中,基准喷管和待测喷管同轴反向设置,通过测量结构应变来直接测量两个喷管的推力差.对测量装置进行了简化模型理论分析,导出其灵敏系数理论公式.依据该公式进行关键参数设计,确定合理的灵敏系数,指导试验系统设计,以保证测量精度.验证性试验表明测量装置可以测量牛级推力差,标定曲线线性较好,以基准喷管推力作为比较对象,推力差的相对扩展不确定度为0.5％.推力差与基准喷管推力的比值越小,以基准喷管推力为比较对象的相对测量不确定度越小.测量装置可以用于小推力火箭发动机喷管性能分析及优化的试验研究.     

烧蚀反应动力学参数的确定

本文提出了一种理论模型与实验相结合确定碳基材料烧蚀反应动力学参数(活化能和指前因子)的方法。首先建立了碳基材料热化学烧蚀的湍流边界层模型,运用数值方法对边界层守恒方程组、表面烧蚀和固相热传导进行耦合求解,计算过程中不需要事先假设烧蚀的控制类型。其次,在电弧加热器中模拟火箭发动机喷喉的烧蚀环境进行烧蚀材料与氧化性气体反应的烧蚀实验。最后,运用所建立的模型和在电弧加热器上的实验结果进行预测比较,确定烧蚀反应的动力学参数。作为实例,确定了石墨与CO_2的烧蚀反应的动力学参数。     

Numerical study of a trapezoidal bypass dual throat nozzle

Bypass Dual Throat Nozzle (BDTN) is a novel type of fluidic thrust vectoring nozzle. To improve the infrared stealth performance of BDTN, a nozzle based on BDTN is proposed and numerically simulated. Each cross-section along the x-axis of the novel nozzle becomes a trapezoid, which is named “BDTN-TRA.” The main numerical simulation results show that BDTN-TRA can produce a thrust vectoring angle when the upper or lower bypass valve is open. The angle difference between the two conditions mentioned above is usually approximately 1°–2°. Even if the two bypasses are closed, BDTN-TRA can produce a small thrust vectoring angle at around 3°–5°. When the sidewall angle increases from 60° to 90°, the thrust coefficient and thrust vectoring angle under each work condition usually decrease. A larger aspect ratio indicates better performance. As the aspect ratio increases over 7.2, the performance of BDTN-TRA is quite close to that of BDTN with rectangular cross-sections at the same aspect ratio. These features will benefit the control and trimming for future aircraft design, especially for the flying wing layout aircraft. Last but not least, BDTN-TRA has a more extraordinary mixing performance compared with BDTN. The distributions of static temperature and axial velocity along the x-axis of BDTN-TRA with sidewall angle of 60° decrease faster than those of BDTN. When the total temperature of the inlet equals 1600 K, the static temperature difference between BDTN-TRA with sidewall angles of 60° and 90° is over 360 K at twice the length of the nozzle downstream of the nozzle exit, which is the reflection for excellent infrared stealth for the fighter.     

椭圆截面喷管排气辐射热流数值分析

文中对直升机红外抑制器椭圆截面喷管排气的辐射热流进行数值分析,并预估排气温度分布、速度和燃气成分分布对辐射热流的影响。为了计算辐射热流和考虑尾焰对周围环境的传热,采用形式简单并可用通用方程计算程序计算的六热流辐射模型。由于排气流场内温度和密度变化较大,故紊流模型采用考虑密度变化的修正双方程k-ε模型。本文中所用的比热C_(pm)是随气流温度和成分而变化,温度用迭代法求得。喷管出口截面为椭圆形,排气为三元自由射流,控制方程为三元抛物形偏微分方程,可采用前进积分求解,这样既可满足计算精度要求,又可节省机时。并根据三元喷流周向边界条件周期性变化的特点,采用CTDMA解法。计算结果与试验结果基本一致,说明本计算方法是可行的。