Numerical investigation of the impact of asymmetric fuel injection on shock train characteristics

Numerical simulations are carried out to investigate the impact of asymmetric fuel injection on shock train characteristics using the commercial-code FLUENT. The asymmetry of fuel injection is examined by changing the fuel flow rates of the upper and lower wall fuel injectors. The numerical approach solves the two-dimensional Reynolds-averaged Navier–Stokes (RANS) equations, supplemented with a k-ω model of turbulence. As a result, different ways of fuel injections will always lead to shock train transitions, with the variations of shock train structure, strength and leading edge position. For symmetric fuel injection, the flowfield of the isolator is quite asymmetric with the boundary layer of the upper wall side developing much stronger than that of the lower wall, which is due to the heterogeneity of the incoming flow. Regarding to asymmetric fuel injection with more of lower wall side, though the pressures in the combustor are nearly the same, the first shock of the shock train converts between ‘Distinct symmetric X type shock’ and ‘Obscure and weaker asymmetric shock’ and the shock train leading edge moves upstream with the increase of the asymmetry level. With regard to asymmetric fuel injection with more of upper wall side, ‘incomplete asymmetric X type shock’ occurs and the shock train structures keep nearly the same with low level of fuel injection asymmetry. Unexpected results like unstart will happen when increasing the level of fuel injection asymmetry. And the isolator will come back to normal state by decreasing the differential of upper and lower wall sides fuel injections.     

三维侧压高超声速进气道不启动流场试验与数值模拟研究

对某构型三维侧压高超声速进气道开展了Ma4的自由射流试验和数值仿真,研究 了低马赫数下不启动流场的流动机理。观测到了具有“对涡”结构的底板油流图案,并得到 了分离区的范围和内部流动特征。分析得出,“对涡”结构油流图案和流场中部分离区的形 成是唇口激波和内收缩形成的逆压梯度作用于侧板激波形成的流场中部低能流区的结果。根 据试验和数值模拟结果给出了流场结构示意图,并为下一步进气道构形设计和性能改善工作 提出了建议。     

超声速进气道及冲压发动机动态特性分析

基于小偏差线性化思想,利用超声速进气道动力学模型计算得到,进气道激波位置和波后压力的响应幅值随频率增大整体趋于减小,但在各阶纵向谐振频率上存在谐振峰。并进一步考虑了燃烧室加质燃烧,分析了冲压发动机气路动态特性,推导出适用于冲压发动机的集中燃烧模型,研究表明在燃油喷注流量的扰动下,冲压发动机幅频响应谐振峰显著。     

冲压发动机燃烧引起激波运动过程一维分析

分析了冲压发动机喷油燃烧引起内流道内正激波运动的机理,采用一维激波捕捉方法,建立了燃油喷入对正激波运动位置影响的一维仿真模型。通过仿真发现:喷入燃油并逐步增大燃油-空气当量比时,正激波逐步向上游运动;燃油-空气当量比越大,正激波越接近进气道喉道;当燃油-空气当量比增大到一定程度时,正激波距离进气道喉道最近,但并未越过喉道;进一步增大燃油-空气当量比,正激波开始向下游回退进一步分析发现:冲压发动机流道及燃烧组织匹配设计直接影响到正激波在流道内的运动位置,需要在设计中格外重视。燃油-空气当量比与激波位置的关系分析可为冲压发动机设计提供一定的理论参考。     

超声速内流场横向喷流的流动与混合特性实验研究

为研究超声速内流场中横向喷流的流动与混合特性,将丙酮蒸汽加入喷流介质,用平面激光诱导荧光(PLIF)技术对流场中流向中心截面和横截面上的丙酮进行成像,研究了喷流的运动轨迹、流场结构、混合方式,以及参数对喷流流动与混合的影响。结果表明:喷流柱的波动失稳及喷流剪切层中生成的大尺度结构有助于增强喷流与主流在近场的混合;提高出口马赫数会导致剪切层失稳以及出现大尺度结构的位置移向下游,不利于改善近场的混合;增大喷口直径能增加喷流在展向的扩展,升高喷流总压能增加喷流在展向和横向的扩展,并使出现大尺度结构的位置靠近上游;在喷注流量相同条件下,采用小喷注面积高总压喷注更利于增强混合。     

加扰流环双组元LAE的流场数值仿真

采用交错网格系统SIMPLE算法和二维两相流场燃烧模型，对有不同形状(矩形、圆弧三角形及其混合结构)和位置扰流环的某双组元液体远地点发动机(LAE)流场进行了数值仿真计算。结果表明，扰流环对提高燃烧效率作用明显。环的位置应适中(不能过于靠近头部和喷管出口处)，其高度越大，燃烧效率越高。另外，矩形环的压力损失大于三角形环。     

空气涡轮火箭发动机内外涵气流掺混研究

通过无化学反应、均匀进气条件下肼单组元空气涡轮火箭发动机混流燃烧室内流场的数值计算,得到了流向涡与正交涡系产生、衰减演变过程及其对内外涵气流掺混效率的影响规律。结果表明,大尺度阵列二次环流诱导形成的流向涡对内外涵气流掺混起主导作用,大波瓣穿透率的斜切波瓣混流器的综合性能较优。结合热试车结果,分析了包括波瓣混流器在内的两类掺混方案的强化掺混效率。分析表明,非均匀进气条件对小尺寸空气涡轮火箭发动机掺混燃烧效率影响很大。     

进口水平投影可控的流线追踪内收缩进气道设计

为了满足两侧进气布局飞行器的乘波前体与进气道一体化设计要求,提出了一种进口水平投影可控的流线追踪内收缩进气道设计方法。基于马赫数分布可控的轴对称基准流场,在指定进口水平投影为椭圆的条件下,采用该方法设计了内收缩进气道并在设计点(Ma=5.4)和接力点(Ma=4.0)对其进行数值研究。结果表明,设计点时进气道都能保持基准流场的波系结构和沿程压力分布,无粘时可以全捕获自由来流,喉道性能与基准流场几乎相等。有粘条件下,设计点和接力点时进气道具有较高的压缩效率和良好的流量捕获能力,接力点的流量系数高达0.85。该设计方法为内收缩进气道与乘波前体的一体化设计提供了新途径。     

Numerical investigation on the shock wave transition in a three-dimensional scramjet isolator

The scramjet isolator, which is used to prevent the hypersonic inlet from disturbances that arise from the pressure rise in the scramjet combustor due to the intense turbulent combustion, is one of the most critical components in hypersonic airbreathing propulsion systems. Any engineering error that is possible in the design and manufacturing procedure of the experimental model, and the intense heat release in the scramjet combustor, may cause the performance of the isolator to decrease, leading to its lack of capability in supporting the back pressure. The coupled implicit Reynolds Averaged Navier–Stokes (RANS) equations and the two-equation standard k?ε turbulent model have been employed to numerically simulate the flow fields in a three-dimensional scramjet isolator. The effects of the divergent angle and the back pressure on the shock wave transition and the location of the leading edge of the shock wave train have been estimated and discussed. The obtained results show that the present numerical results are in very good agreement with the available experimental shadow-pictures, and the numerical method is more suitable for capturing the shock wave train and predicting the location of the leading edge of the shock wave train in the scramjet isolator than the present two-dimensional numerical methods. This is due to the small width-to-height ratio of the isolator and the intense three-dimensional flow structures. On increasing the divergent angle of the scramjet isolator, the static pressure along the central symmetrical line of the isolator decreases sharply. This is due to the strong expansion wave generated at the entrance of the isolator, and when the divergent angle of the isolator is sufficiently large, namely 1.5°, a zone of negative pressure is formed just ahead of the leading edge of the shock wave train. At the same time, the shock wave train varies from being oblique to being normal, and then back to oblique. With an increase in the prescribed back pressure at the exit of the scramjet isolator, the leading edge of the shock wave train moves forward towards the entrance of the isolator, and when the back pressure is sufficiently large, unstart conditions in the hypersonic inlet can take place if the shock train reaches the inlet.     

斜激波串与上游激波干扰特性试验研究

为研究斜激波串在背压条件下前移与上游激波相互干扰的流场结构和运动规律，在来流为马赫数 2.7 的直管道内设计一种等宽度斜楔，采用动态压力测量、高速纹影和粒子图像测速(PIV)技术等手段进行了试验。研究结果表明：内置斜楔在管道内产生入射激波、分离激波、膨胀波、再附激波和激波诱导分离等复杂上游激波流场，在分离区附近形成有顺压梯度和逆压梯度的区域。当增大下游压比时，斜激波串逐渐向上游激波流场移动；经过斜楔产生的分离区时，斜激波串的移动速度急剧提升，同时出现非对称分离偏转方向的切换。对比了三种长度尺寸的等楔角斜楔所产生的上游激波流场的差异性，发现在相同的斜楔前缘起始点和楔角时，随着斜楔长度的增加，上游激波流场中激波诱导的分离尺度逐渐变大。     

Three dimensional investigation of the shock train structure in a convergent–divergent nozzle

Three-dimensional computational fluid dynamics analyses have been employed to study the compressible and turbulent flow of the shock train in a convergent–divergent nozzle. The primary goal is to determine the behavior, location, and number of shocks. In this context, full multi-grid initialization, Reynolds stress turbulence model (RSM), and the grid adaption techniques in the Fluent software are utilized under the 3D investigation. The results showed that RSM solution matches with the experimental data suitably. The effects of applying heat generation sources and changing inlet flow total temperature have been investigated. Our simulations showed that changes in the heat generation rate and total temperature of the intake flow influence on the starting point of shock, shock strength, minimum pressure, as well as the maximum flow Mach number.     

下颌式进气道固冲发动机二次燃烧性能研究

针对采用下颌式进气道的固体火箭冲压发动机,建立了二次燃烧性能计算模型,对掺混燃烧性能进行了仿真研究。研究表明,采用掺混装置可大幅提升下颌式进气道的固冲发动机补燃室一次燃气和空气的掺混均匀度,并通过数值仿真对掺混装置进行了优化。结合数值仿真优化结果,通过地面直连试验,验证了不采用与采用掺混装置的补燃室二次燃烧性能。试验结果表明,合理设计掺混装置,可显著提高补燃室二次燃烧性能,特征速度燃烧效率均在93%以上;空燃比在6~20之间的发动机高空比冲提升了55%以上,空燃比在20~30之间的发动机高空比冲提升了75%以上。     

大流量、高压航天反力式涡轮最佳通道气动设计

基于闭式补燃循环液体火箭发动机流量大的特点，欲设计效率水平高的航天反力式涡轮。因涡轮进出口压力均很高，而膨胀比小、载荷系数大，为保证较高的涡轮效率水平，对涡轮气动设计方法进行了优化。在涡轮进口总温、总压、转速和功率一定条件下，以AMDC／KQ涡轮叶栅损失模型为基础，依据涡轮中径的一维气动计算，对涡轮子午通道、叶栅通道及叶栅造型几组参数组合分别进行了气动设计的优化，研究了涡轮中径、叶高、叶栅稠度、导动叶喉宽匹配及动叶进口构造角对涡轮效率的影响，实现了涡轮效率水平最高。     

Numerical validation and parametric investigation on the cold flow field of a typical cavity-based scramjet combustor

The three-dimensional coupled implicit Reynolds Averaged Navier–Stokes (RANS) equations and the two equation standard k–ε turbulence model has been employed to numerically simulate the cold flow field in a typical cavity-based scramjet combustor. The numerical results show reasonable agreement with the schlieren photograph and the pressure distribution available in the open literature. The pressure distribution after the first pressure rise is under-predicted. There are five shock waves existing in the cold flow field of the referenced combustor. The first and second pressure rises on the upper wall of the combustor are predicted accurately with the medium grid. The other three shock waves occur in the core flow of the combustor. The location of the pressure rise due to these three shock waves could not be predicted accurately due to the presence of recirculation zone downstream of the small step. Further, the effect of length-to-depth ratio of the cavity and the back pressure on the wave structure in the combustor has been investigated. The obtained results show that there is an optimal length-to-depth ratio for the cavity to restrict the movement of the shock wave train in the flow field of the scramjet combustor. The low velocity region in the cavity affects the downstream flow field for low back pressure. The intensity of the shock wave generated at the exit of the isolator depends on the back pressure at the exit of the combustor and this in turn affects the pressure distribution on the upper wall of the combustor.     

燃烧室结构对固液火箭发动机燃烧与流动的影响研究

建立了85%H2O2-PE固液火箭发动机氧化剂H2O2催化分解、PE燃料热解以及热解气体与氧化剂分解气体扩散燃烧的综合模型,计算了固液火箭发动机燃烧室轴对称二维内流场,对不同结构燃烧室内流场的计算结果进行了对比,研究了补燃室和氧化剂入口突扩结构对发动机燃烧性能的影响.结果表明,增加氧化剂入口突扩段有利于发动机稳定工作和充分燃烧,增加补燃室长度可以提高发动机平均燃烧温度,使燃烧更加充分.     

Numerical investigation of the nonreacting and reacting flow fields in a transverse gaseous injection channel with different species

The mixing and combustion process has an important impact on the engineering realization of the scramjet engine. The nonreacting and reacting flow fields in a transverse injection channel have been investigated numerically, and the predicted results have been compared with the available experimental data in the open literature, the wall pressure distributions, the separation length, as well as the penetration height. Further, the influences of the molecular weight of the fuel and the jet-to-crossflow pressure ratio on the wall pressure distribution have been studied. The obtained results show that the predicted results show reasonable agreement with the experimental data, and the variable trends of the penetration height and the separation distance are almost the same as those obtained in the experiment. The vapor pressure model is suitable to fit the relationship between the penetration height, the separation distance and the jet-to-crossflow pressure ratio. The combustion process mainly occurs upstream of the injection port, and it makes a great difference to the wall pressure distribution upstream of the injection port, especially when the jet-to-crossflow pressure ratio is large enough, namely 17.72 and 25.15 in the range considered in the current study. For hydrogen, the combustion downstream of the injection port occurs more intensively, and this may be induced by its smaller molecular weight.     

Multi-objective design optimization of the transverse gaseous jet in supersonic flows

The mixing process between the injectant and the supersonic crossflow is one of the important issues for the design of the scramjet engine, and the efficiency mixing has a great impact on the improvement of the combustion efficiency. A hovering vortex is formed between the separation region and the barrel shock wave, and this may be induced by the large negative density gradient. The separation region provides a good mixing area for the injectant and the subsonic boundary layer. In the current study, the transverse injection flow field with a freestream Mach number of 3.5 has been optimized by the non-dominated sorting genetic algorithm (NSGA II) coupled with the Kriging surrogate model; and the variance analysis method and the extreme difference analysis method have been employed to evaluate the values of the objective functions. The obtained results show that the jet-to-crossflow pressure ratio is the most important design variable for the transverse injection flow field, and the injectant molecular weight and the slot width should be considered for the mixing process between the injectant and the supersonic crossflow. There exists an optimal penetration height for the mixing efficiency, and its value is about 14.3 mm in the range considered in the current study. The larger penetration height provides a larger total pressure loss, and there must be a tradeoff between these two objection functions. In addition, this study demonstrates that the multi-objective design optimization method with the data mining technique can be used efficiently to explore the relationship between the design variables and the objective functions.     

固体火箭冲压发动机的若干技术问题

简述了固体火箭冲压发动机类型及工作原理,总体评价了固体火箭冲压发动机发展时快时慢的原因,弹-机一体化设计、贫氧推进剂、进气道、转级机构、补燃室等设计中应注意的问题,提出了应加强的研究工作,即开展高能低沉积燃烧产物贫氧推进剂研究;完善多种燃气流量调节装置方案,提高其可靠性;进一步开展一次燃烧和二次掺混燃烧的理论和实验研究,提高燃烧效率;尽快建立自由射流等试验研究手段,开展相关的研究工作;适时开展固体燃料冲压、固体超燃冲压及膏体冲压等发动机的研究,不断拓宽应用领域。     

Ignition delay time in swirling supersonic flow

Supersonic hydrogen–air cylindrical mixing layer is numerically analyzed to investigate the effect of inlet swirl on ignition delay time in scramjets. Combustion is treated using detailed chemical kinetics. One-equation turbulence model of Spalart and Allmaras is chosen to study the problem and advection upstream splitting method is used as computational scheme. The results show that swirling both fuel and oxidizer streams may drastically decrease the ignition distance in supersonic combustion, unlike using the swirl just in fuel stream which has no helpful effect.     

采用新型基准流场的高超声速内收缩进气道性能分析

通过改变中心体形状,设计了新型轴对称基准流场,可显著降低反射激波强度,明显提高压缩效率。基于该基准流场和传统基准流场,分别设计了两个圆形出口内收缩进气道,并对二者的流场及总体性能进行了数值研究。结果表明,新的进气道设计点和接力点肩点附近激波附面层相互作用减弱,流场结构优于传统进气道,压缩效率明显提高,同时进气道起动性能得到改善。