超燃冲压发动机燃烧室流动问题的数值验证

基于三维可压缩RANS方程,结合有限速率化学反应模型,发展了一套模拟多组元化 学 反应与湍流流动的计算程序。通过对超燃冲压发动机燃烧室内的几种典型流动问题,如横 向喷流问题、后向台阶流动问题和氢／空气化学反应问题进行数值模拟,验证计算程序的可 靠性,将数值模拟结果与实验结果进行了对比,结果吻合较好。     

Numerical investigation on hydrogen combustion in a scramjet with 3D sidewall compression inlet

Numerical simulations were employed to analyze the flowfield of a scramjet with three-dimensional (3D) sidewall compression inlet, and the effect of inlet distortion on the mixing and combustion process was examined. The numerical approach solved the compressible Reynolds Averaged Navier–Stokes (RANS) equations supplemented with a finite rate chemical reacting model for the combustion of hydrogen fuel and air. Turbulence closure was achieved using Menter shear-stress transport (SST) model. To verify the accuracy of the simulation, the computed wall pressure was compared with the experimental data of the direct-connect combustor test. The metrics employed in the simulations included qualitative assessments related to flow structure as well as quantitative values of fuel mixing efficiency, combustion efficiency and static pressure distribution. Intake sidewalls were found to strongly affect the inlet flow structure, which became more complex in the nonuniform flowfield. The shock train system affected the combustion region located upstream of the injection and led to pairs of asymmetric separation bubbles. Nevertheless, the shock train system dissipated due to the reactions, the combustion patterns of each fuel jets in downstream region were nearly identical, and the degree of improvement of mixing and combustion efficiency near the downstream injectors was less than that near the upstream injectors.     

一种新颖的基于总能量守恒的化学平衡流算法

在推进剂燃烧的建模上,传统的热力计算方法一般基于总焓守恒求解定压绝热燃烧温度和平衡组分,不能考虑壁面传热;在燃气流动的建模上,通常采用的冻结流模型认为本地的组分及热物理性质与燃烧室瞬时一致,忽略了这些参数因来流气体与本网格滞留气体掺混带来的随时间的缓变效应。提出了一种新颖的可以考虑壁面传热的基于总能量守恒的化学平衡流计算方法,运用Fortran2008语言,采用面向对象编程方法建立了化学平衡流燃气发生器管道的模块化仿真模型,并将该模型应用到一个包含42个组件的涡轮试验台气路系统的建模与仿真中。与早期模型仿真结果及试验数据的对比发现,新模型的仿真结果有一定改进,更加接近试验数据。     

Igniter jet dynamics in solid fuel ramjets

The dynamics of a two dimensional plane jet injected at the base of a step, parallel to the wall, in backward facing step flow geometry is numerically studied. The objective of this work is to gain insight into the dynamics of the igniter flow field in solid fuel ramjet motors. Solid fuel ramjets operate by ingestion of air and subsequent combustion with a solid fuel grain such as polyethylene. The system of governing equations is solved with a finite volume approach using a structured grid in which the AUSM⁺ scheme is used to calculate the convective fluxes. The Spalart and Allmaras turbulence model is used in these simulations. Experimental data have been used to validate the flow solver and turbulence model simulation results. The comparison of the numerical results and experimental data has validated the use of the adopted turbulence model for the study of this type of problem. A special attention is paid to the igniter jet exit location. It is shown that the wall jet igniter, issuing from the base of the step, drastically changes the structure of recirculating region of backward facing step flow and produces large and damaging shear stress on the fuel surface. Moving the igniter jet exit location to the top of the backward facing step changes the flow field favorably, by reducing the fuel surface shear stress by an order of magnitude and maintaining the recirculating region behind the step, which can provide proper residence time for the fuel–air mixture chemical reactions.     

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

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

固液混合火箭发动机燃烧室和喷管流动数值模拟

固液混合火箭发动机是采用液体作为氧化剂,固体作为燃料的一种典型的混合火箭发动机.固液混合火箭发动机中的燃烧和流动问题是固液混合火箭发动机设计中的关键问题,对固液混合火箭发动机的燃烧室和喷管进行一体化计算很有必要.利用二维轴对称N-S方程和组分方程对选用液氧/端羟基聚丁二烯推进剂的固液混合火箭发动机的燃烧室和喷管进行了一体化计算.计算采用LU时间隐式格式、MUSCL空间离散和Van Leer矢通量分裂方法,采用有限速率化学反应模型,对化学源相进行了点隐式处理.计算中分别采用了一步化学反应模型和两步化学反应模型方案,计算了多个氧化剂流速和燃烧室压强下的燃烧室和喷管流场分布,对化学模型进行了选择,为固液混合火箭发动机的设计提供了依据.     

基于等效理论的RQL燃烧室掺混特性研究

采用RNG k-ε湍流模型对RQL（rich—burn／quick—quench／lean—burn）燃烧室内气流的掺混过程进行了数值模拟,运用等效思想,通过对混合流场内部温度、近壁区温度以及出口温度的分析,讨论不同掺混角度和掺混射流速度对RQL燃烧室混合特性的影响,进而了解RQL燃烧室工作时热流流场的结构状况。研究结果表明,不同的掺混气流入射角度和速度,对RQL燃烧室内气流掺混的高温区位置、壁面温度以及出口温度分布的影响明显。     

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.     

耦合凝相-气相动力学机理的二硝酰胺铵燃烧模型

为了研究二硝酰胺铵ADN(Ammonium Dinitramide)固体推进剂燃烧物理化学过程并预测其燃 烧特性,建立一个耦合凝相-气相动力学机理的ADN燃烧模型。该模型基于凝相与气相的总连 续方程、组元连续方程、能量守恒方程及有限速率化学动力学方程而建立,并引入多组元系 统状态方程封闭方程组。模型中包含34种组元,1个固相（凝相）ADN分解总化学反应和165 个气相细节（基元）化学反应,并使用以温度函数表示的物性参数进行计算。应用气相燃烧 模型 对0.6 MPa下ADN燃烧火焰温度、组元摩尔浓度分布进行预测;应用耦合凝相-气相的燃烧模 型对0.2 MPa～36 MPa压强区域内柱状端燃ADN推进剂燃速、燃烧表面温度进行预测,计算结 果与文献报道试验数据较吻合。说明该燃烧模型能够较准确描述ADN气相燃烧波结构和ADN固 体推进剂燃速特性。
     

Supercomputer modeling of hydrogen combustion in rocket engines

Hydrogen being an ecological fuel is very attractive now for rocket engines designers. However, peculiarities of hydrogen combustion kinetics, the presence of zones of inverse dependence of reaction rate on pressure, etc. prevents from using hydrogen engines in all stages not being supported by other types of engines, which often brings the ecological gains back to zero from using hydrogen. Computer aided design of new effective and clean hydrogen engines needs mathematical tools for supercomputer modeling of hydrogen–oxygen components mixing and combustion in rocket engines.The paper presents the results of developing verification and validation of mathematical model making it possible to simulate unsteady processes of ignition and combustion in rocket engines.     

Influences of H2O mass fraction and chemical kinetics mechanism on the turbulent diffusion combustion of H2–O2 in supersonic flows

Hydrogen is one of the most promising fuels for the airbreathing hypersonic propulsion system, and it attracts an increasing attention of the researchers worldwide. In this study, a typical hydrogen-fueled supersonic combustor was investigated numerically, and the predicted results were compared with the available experimental data in the open literature. Two different chemical reaction mechanisms were employed to evaluate their effects on the combustion of H₂–O₂, namely the two-step and the seven-step mechanisms, and the vitiation effect was analyzed by varying the H₂O mass fraction. The obtained results show that the predicted mole fraction profiles for different components show very good agreement with the available experimental data under the supersonic mixing and combustion conditions, and the chemical reaction mechanism has only a slight impact on the overall performance of the turbulent diffusion combustion. The simple mechanism of H₂–O₂ can be employed to evaluate the performance of the combustor in order to reduce the computational cost. The H₂O flow vitiation makes a great difference to the combustion of H₂–O₂, and there is an optimal H₂O mass fraction existing to enhance the intensity of the turbulent combustion. In the range considered in this paper, its optimal value is 0.15. The initiated location of the reaction appears far away from the bottom wall with the increase of the H₂O mass fraction, and the H₂O flow vitiation quickens the transition from subsonic to supersonic mode at the exit of the combustor.     

水反应金属燃料发动机三维两相燃烧数值模拟

用三维湍流N-S方程及颗粒轨道模型描述水反应金属燃料发动机内部喷雾两相湍流燃烧过程。通过耦合求解气液两相流模型方程,得到发动机燃烧流场。通过模拟Mg与水的反应,分析比较了一次和二次进水方案的不同流场特性。研究结果表明,二次进水方案更有利于火焰稳定和提高燃烧性能。     

富氢/富氧燃气温度对同轴直流气-气喷嘴性能的影响

通过求解使用k-ε湍流模型的Navier-Stokes方程组对采用同轴直流气-气单喷嘴燃烧室的燃烧流场进行数值模拟,对比分析了富氢/富氧燃气推进剂与常温氢气/氧气推进剂条件下的燃烧流场、燃烧室室壁和喷注面板处的燃气温度,研究了富氢/富氧燃气温度变化对燃烧流场和燃烧室热载的影响。数值结果表明:富氢/富氧燃气气-气喷嘴的燃烧性能较好,但热载较高;富氢/富氧燃气温度一定范围内提高对燃烧性能影响不明显,而热载增加。     

Thrust modulation of solid propellant rockets by a fluidic vortex valve with secondary combustion

An analytical model for calculation of the thrust and pressure modulation of a solid propellant rocket by means of a fluidic vortex valve with secondary combustion has been developed. Thrust control by the vortex valve method was found superior to the axial injection of control flow. Addition of oxygen in the injected flow improves the energetic performance of the system as well as the thrust modulation capability. Experiments have been conducted using a mixture of nitrogen and oxygen as the injection gas. The two main parameters investigated in a series of experiments were the oxygen percentage in the injection gas and the ratio between the mass flow rates of the control gas and the solid propellant combustion products. The results show an increase of thrust by a factor of 2 for a 25% addition in mass flow rate by secondary injection at optimal conditions.     

超燃燃烧室内流场数值模拟研究

对超燃燃烧室的高速流动与燃烧的耦合流场进行了一维和三维数值模拟,揭示了其流场的本质和各个参数的关系。首先阐述了超燃燃烧室一维流动特性,然后分析了一维流动模型,获得了芯流面积的计算公式。将一维流动模型应用到实际超燃燃烧室中,与试验数据对比,误差小于5％,并且选出了最优模型。由于一维模型仅能模拟轴向的一维参数变化规律,不能...     

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

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

Ignition of fuel sprays by shock wave mathematical modeling and numerical simulation

The paper presents the results of developing of physical and mathematical model making it possible to take into account the effect of droplets non-uniformity in space and size distribution on ignition conditions for fuel sprays. The influence of condensed phase volume fraction on ignition and combustion of sprays was studied, physical and mathematical models for multi-phase flows, mixture formation and combustion of liquid fuels based on solving Navier–Stokes equations for gas phase accounting for thermal and mechanical interaction with poly-dispersed droplets array. The problems of particulate phase dynamics are regarded accounting for the interaction with gas phase atomization, evaporation and combustion.It was shown that depending on droplet size distribution and aerosol cloud density different flow scenarios were possible.Several ignition zones could be formed behind incident shock wave depending on mixture properties and initiation parameters. The possibility of numerical simulation permitting variation of definite parameters only made it possible to explain this fact.     

The turbulent combustion of a propane-air mixture

This paper models the combustion of a turbulent homogeneous mixture of propane and air within a duct having a stationary one-dimensional mean flow. The Bray-Moss model is applied to the closure of the chemical production terms, using a probability density function (pdf) of the temperature which is chosen as the characteristic variable. Under the conditions chosen for the study, chemical kinetic factors are important and the conventional assumption, that heat release is controlled by turbulent mixing, is not valid. The chemical model of Edelman and Fortune for the combustion of hydrocarbons is used and simplifying assumptions are made which reduce the systems of unknowns to that of the temperature alone. This leads to the introduction of two chemical production terms which are defined respectively in a “delay zone”, where the heat release is modest, and a “combustion zone”. The required equations for the Favre-averaged temperature, turbulence kinetic energy and the mean square fluctuation of the temperature are solved numerically. In the delay zone, a comparison is made between a second order Borghi type closure and the pdf closure. Good agreement is found in the case of relatively small turbulence intensity. It is shown that the pdf formulation does not require the two zones to be spatially distinct. Differing chemical source terms can be discriminated instantaneously by the reaction progress variable and contributions to the average production terms appropriately apportioned by its pdf. Predictions are made of the profiles of mean temperature and mean square fluctuation under different initial turbulence levels.     

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.     

RBCC燃烧室超声速反应混合层特性的大涡模拟

以飞行马赫数为4.5Ma的RBCC发动机典型工作状态为研究背景,采用大涡模拟研究了支板火箭射流和空气来流形成的超声速反应混合层的掺混燃烧过程,获得了燃烧室内详细的流场结构和流动特征,分析了强射流条件下超声速反应混合层的特性。结果表明由于速度梯度的存在,火箭射流进入燃烧室后与空气来流形成环形剪切层,剪切层内丰富的旋涡结构主导火箭射流和空气来流的掺混燃烧,随着湍流能量的串级输运,化学反应过程中释放的能量将被转化成细观尺度的湍流动能,大尺度旋涡将能量传递给小尺度旋涡并最终耗散,细小尺度的旋涡一方面能够促进燃烧反应物的掺混并强化燃烧过程,另一方面会给化学反应过程带来强烈的脉动,使得局部火焰淬灭,火焰结构表现出明显的非定常性。