旋流燃烧器复合小火焰模型的大涡模拟

对低旋燃烧器内的流动和燃烧进行了大涡模拟,其中化学反应分别采用传统的扩散小火焰模型和所提出的复合小火焰模型描述。复合小火焰模型借助于燃烧区索引的概念区分当地的燃烧模式,若当地的燃烧为扩散燃烧,则调用扩散小火焰库,否则调用预混小火焰库。数值结果与实验数据的对比表明,采用大涡模拟方法结合小火焰模型能够很好地模拟实验室尺度的低旋燃烧器,且采用复合小火焰模型能够得到更加符合实际的结果。     

LES of supersonic combustion in a scramjet engine model

In this study, Large Eddy Simulation (LES) has been used to examine supersonic flow and combustion in a model scramjet combustor. The LES model is based on an unstructured finite volume discretization, using total variational diminishing flux reconstruction, of the filtered continuity, momentum, enthalpy, and passive/reactive scalar equations, used to describe the combustion process. The configuration used is similar to the laboratory scramjet at the Institute for Chemical Propulsion of the German Aerospace Center (DLR) and consists of a one-sided divergent channel with a wedge-shaped flameholder at the base of which hydrogen is injected. Here, we investigate supersonic flow with hydrogen injection and supersonic flow with hydrogen injection and combustion. For the purpose of validation, the LES results are compared with experimental data for velocity and temperature at different cross-sections. In addition, qualitative comparisons are also made between predicted and measured shadowgraph images. The LES computations are capable of predicting both the non-reacting and reacting flowfields reasonably well—in particular we notice that the LES model identifies and differentiates between peculiarities of the flowfields found in the experiments.     

小火焰模型在贫燃预混火焰中的研究

由层流小火焰库引入详细化学反应机理,通过简化的PDF方法计算组分浓度、平均温度和密度等变量,以钝体火焰稳定燃烧室和某燃气轮机上的燃烧室为例,模拟甲烷/空气贫燃条件下预混燃烧的平均火焰位置和火焰厚度,计算结果与实验结果吻合良好,这表明此方法能够较好计算出平均湍流火焰的主要特征。     

超燃冲压发动机推进性能理论分析

超燃冲压发动机的正推力问题和超声速燃烧的稳定性问题是制约超燃冲压发动机发展的两个关键气动物理问题.虽然经过50多年的研究,但是目前国内外对这两个关键问题的机理还没有研究清楚.文章首次将CJ爆轰理论应用于超燃冲压发动机推进性能分析,给出了这两个关键气动问题的理论分析结果.分析结果表明,燃烧室入口空气静温对发动机的推进性能产生重要影响.当爆轰波的爆速大于隔离段内空气来流的速度时,会向隔离段上游传播,导致发动机不起动.飞行Mach数Ma=6~8是超燃发动机的临界不稳定范围,飞行Mach数Ma>9,超声速燃烧将变得稳定.      

CFD analysis of the HyShot II scramjet combustor

The development of novel air-breathing engines such as supersonic combustion ramjets (scramjets) depends on the understanding of supersonic mixing, self-ignition and combustion. These aerothermochemical processes occur together in a scramjet engine and are notoriously difficult to understand. In the present study, we aim at analyzing the HyShot II scramjet combustor mounted in the High Enthalpy Shock Tunnel Göttingen (HEG) by using Reynolds Averaged Navier Stokes (RANS) and Large Eddy Simulation (LES) models with detailed and reduced chemistry. To account for the complicated flow in the HEG facility a zonal approach is adopted in which RANS is used to simulate the flow in the HEG nozzle and test-section, providing the necessary inflow boundary conditions for more detailed RANS and LES of the reacting flow in the HyShot combustor. Comparison of predicted wall pressures and heat fluxes with experimental data show good agreement, and in particular does the LES agree well with the experimental data. The LES results are used to elucidate the flow, mixing, self-ignition and subsequent combustion processes in the combustor. The combustor flow can be separated into the mixing zone, in which turbulent mixing from the jet-in-cross flow injectors dominates, the self-ignition zone, in which self-ignition rapidly takes place, and the turbulent combustion zone, located towards the end of the combustor, in which most of the heat release and volumetric expansion takes place. Self-ignition occurs at some distance downstream of the injectors, resulting in a distinct pressure rise further downstream due to the volumetric expansion as observed in the experiments. The jet penetration is about 30% of the combustor height and the combustion efficiency is found to be around 83%.     

基于CARS技术的超燃冲压发动机点火过程温度测量

相干anti-Stokes Raman散射（coherent anti-Stokes Raman scattering,CARS）技术作为一种非接触测量手段,已广泛应用于多种发动机模型燃烧室温度测量及地面试验.然而,目前的工作主要集中在稳态燃烧场温度的测量,缺乏用高分辨率的单脉冲来测量瞬变的燃烧火焰温度及组分浓度的研究.基于CARS理论,结合多参数拟合算法,开发了基于MATLAB的CARS光谱计算和拟合程序CARSCF;利用McKenna平面火焰炉在不同工况下进行了温度测量,并与DLR测量结果进行对比,结果显示开发的CARSCF具有较高的测量重复性和准确性;最后将CARS技术应用于测量超燃冲压发动机点火过程中的温度测量,获取了点火过程中的温度.结果显示,在来流Mach数为3的条件下,H₂/air点火过程中温度呈现急剧上升然后缓慢下降,而CARS信号则呈现急剧上升然后急剧下降随后又缓慢上升的趋势,并且在点火过程中最高温度为1 511 K.      

Numerical investigation of low-frequency instability and frequency shifting in a scramjet combustor

The combustion instability in a laboratory-scale direct-connect hydrogen-fueled scramjet combustor is investigated numerically. The numerical simulation has been carried out using a delayed detached eddy simulation (DDES) with a detailed reaction mechanism. The computational framework has high fidelity by applying multi-dimensional high order accurate schemes for handling convective and viscous fluxes. The field data were accumulated up to 100 milliseconds on each case to capture sufficiently the repetitive behavior of low-frequency instability of order of 100 Hz. The numerical results exhibit the formation/dissipation of pressure and shock wave induced by continuous heat release in the combustor. This motion of pressure/shock wave, so-called upstream-traveling shock wave, presents repeated dynamics between isolator and combustor with a period of several milliseconds. With this periodic hydrodynamic characteristic, the upstream-traveling shock wave interacts with the boundary layer and injected fuel stream affecting fuel/air mixing and burning, and finally inducing the combustion instability in a scramjet combustor. Frequency analysis derived major instability frequencies of 190 Hz and 450 Hz in the isolator and combustor for low and high equivalence ratios, respectively. Current numerical results present the underlying flow physics on the shifting of the instability frequency by changing the equivalence ratio observed by the previous experimental studies. The fact that an instability frequency exists homogeneously from isolator to combustor informs that the combustion instability of scramjet engine is the fully coupled flow/combustion dynamics throughout the engine on a macroscopic scale.     

A spray flamelet/progress variable approach combined with a transported joint PDF model for turbulent spray flames

A spray flamelet/progress variable approach is developed for use in spray combustion with partly pre-vaporised liquid fuel, where a laminar spray flamelet library accounts for evaporation within the laminar flame structures. For this purpose, the standard spray flamelet formulation for pure evaporating liquid fuel and oxidiser is extended by a chemical reaction progress variable in both the turbulent spray flame model and the laminar spray flame structures, in order to account for the effect of pre-vaporised liquid fuel for instance through use of a pilot flame. This new approach is combined with a transported joint probability density function (PDF) method for the simulation of a turbulent piloted ethanol/air spray flame, and the extension requires the formulation of a joint three-variate PDF depending on the gas phase mixture fraction, the chemical reaction progress variable, and gas enthalpy. The molecular mixing is modelled with the extended interaction-by-exchange-with-the-mean (IEM) model, where source terms account for spray evaporation and heat exchange due to evaporation as well as the chemical reaction rate for the chemical reaction progress variable. This is the first formulation using a spray flamelet model considering both evaporation and partly pre-vaporised liquid fuel within the laminar spray flamelets. Results with this new formulation show good agreement with the experimental data provided by A.R. Masri, Sydney, Australia. The analysis of the Lagrangian statistics of the gas temperature and the OH mass fraction indicates that partially premixed combustion prevails near the nozzle exit of the spray, whereas further downstream, the non-premixed flame is promoted towards the inner rich-side of the spray jet since the pilot flame heats up the premixed inner spray zone. In summary, the simulation with the new formulation considering the reaction progress variable shows good performance, greatly improving the standard formulation, and it provides new insight into the local structure of this complex spray flame.     

Dual/scram-mode combustion limits of ethylene and surrogate endothermically-cracked hydrocarbon fuels at Mach 8 equivalent high-enthalpy conditions

This paper examines the scram/dual-mode combustion limits of hydrocabon fuels within a Mach 8, scramjet combustor. Flight-equivalent flows were delivered to the axisymmetric, cavity combustor via a reflected shock tunnel. Two scramjet fuels were examined: ethylene and a surrogate mixture representing endothermically cracked n-dodecane. Combustion modes were examined via static pressure sensors and through both chemiluminescence imaging, and planar laser induced fluorescence (PLIF) of the OH combustion radical in the combustor exhaust plume. Ethylene-fuelled experiments developed scram-mode combustion under reduced fuelling conditions, experiencing shock wave dominated flowfields. OH PLIF diagnostics indicated such combustion modes developed a ring-like structure of combustion products, primarily axisymmetrically adjacent to the combustor wall. Increased fuelling anchored combustion downstream of the fuel injector, while further increases instigated dual-mode combustion. In this mode, subsonic combustion regions combine with the supersonic coreflow to permit the transfer of information upstream with substantially increased pressure encountered. Optical diagnostics indicate broadly asymmetric, unsteady combustion features. The surrogate mixture representing endothermically cracked n-dodecane experienced rapid onset from no-combustion (optically confirmed) to fully developed dual-mode combustion at critical fuelling rates. OH PLIF signals and chemiluminescence of this fuel were weaker than comparable ethylene cases, indicating potential differences in combustion pathways.     

烧蚀支板对超燃冲压发动机燃烧室性能的影响

为了提高超燃冲压发动机燃烧室的性能,本文提出了燃料喷注支板与烧蚀支板组合的燃烧室新方案,并研究了新方案对超燃冲压发动机燃烧室性能的影响。相比于单燃料喷注支板方式而言,加入烧蚀支板后,虽然燃烧室内的总压恢复系数有所下降,但燃烧室内燃料与空气的混合效率、燃烧效率均有显著提高,燃烧效率的提高弥补了燃烧室内总压损失所带来的机械能损失,使得燃料喷注支板和烧蚀支板组合方式下的燃烧室比冲高于单燃料喷注支板时的比冲。     

JF12激波风洞高Mach数超燃冲压发动机实验研究

针对高Mach数（Ma ≥ 7）超燃冲压发动机高气动阻力下的燃烧组织问题,提出一种双突扩燃烧室结构方案.使用数值模拟方法考察了射流与双突扩燃烧室组合方式的混合燃烧特性.设计了双突扩超燃冲压发动机模型,在力学研究所JF12长试验时间激波风洞内,开展了Ma=7.0和Ma=9.5的氢燃料点火和燃烧试验对比.在风洞有效试验时间100 ms内,实现了Ma=7.0和Ma=9.5超燃冲压发动机的成功点火与稳定燃烧.在Ma=7.0情况下,进气道采用三维压缩,燃烧室入口设计Mach数Ma_c=2.5,壁面压力分布实验结果显示燃烧放热靠近燃烧室扩张段上游;在Ma=9.5情况下,进气道采用二维压缩,燃烧室入口设计Mach数Ma_c=3.5,由于燃烧室流动速度特别高,燃烧放热靠近燃烧室扩张段下游.      

The impact of dilatation,scrambling, and pressure transport in turbulent premixed flames

Premixed turbulent flames feature strong interactions between chemical reactions and turbulence that affect scalar and turbulence statistics. The focus of the present work is on clarifying the impact of pressure dilatation/flamelet scrambling effects with a comprehensive second-moment closure used for evaluation purposes. Model extensions that take into account flamelet orientation and molecular diffusion are derived. Isothermal pressure transport is included with an additional variable density contribution derived for the flamelet regime of combustion. Full closure is assessed by comparisons with Direct Numerical Simulations (DNSs) of statistically ‘steady’ fully developed premixed turbulent planar flames at different expansion ratios. Subsequently, the prediction of lean premixed turbulent methane–air flames featuring fractal grid generated turbulence in an opposed jet geometry is considered. The overall agreement shows that ‘dilatation’ effects contribute to counter-gradient transport and can also increase the turbulent kinetic energy significantly. Levels of anisotropy are broadly consistent with the DNS data and key aspects of opposed jet flames are well predicted. However, it is also shown that complications arise due to interactions between the imposed pressure gradient and combustion and that redistribution is affected along with the scalar flux at the leading edge. The latter is strongly affected by the reaction rate closure and, potentially, by pressure transport. Overall, the derived models offer significant improvements and can readily be applied to the modelling of premixed turbulent flames at practical rates of heat release.     

Modelling compression ignition engines by incorporation of the flamelet generated manifolds combustion closure

Tabulated chemistry models allow to include detailed chemistry effects at low cost in numerical simulations of reactive flows. Characteristics of the reactive fluid flows are described by a reduced set of parameters that are representative of the flame structure at small scales so-called flamelets. For a specific turbulent combustion configuration, flamelet combustion closure, with proper formulation of the flame structure can be applied. In this study, flamelet generated manifolds (FGM) combustion closure with progress variable approach were incorporated with OpenFOAM® source code to model combustion within compression ignition engines. For IC engine applications, multi-dimensional flamelet look-up tables for counter flow diffusive flame configuration were generated. Source terms of non-premixed combustion configuration in flamelet domain were tabulated based on pressure, temperature of unburned mixture, mixture fraction, and progress variable. A new frozen flamelet method was introduced to link one dimensional reaction diffusion space to multi-dimensional Computational Fluid Dynamics (CFD) physical space to fulfill correct modelling of thermal state of the engine at expansion stroke when charge composition was changed after combustion and reaction rates were subsided. Predictability of the developed numerical framework were evaluated for Sandia Spray A (constant volume vessel), Spray B (light duty optical Diesel engine), and a heavy duty Diesel engine experiments under Reynolds averaged Navier Stokes turbulence formulation. Results showed that application of multi-dimensional FGM combustion closure can comprehensively predict key parameters such as: ignition delay, in-cylinder pressure, apparent heat release rate, flame lift-off , and flame structure in Diesel engines.     

Prediction of soot and thermal radiation in a model gas turbine combustor burning kerosene fuel spray at different swirl levels

Combustion of kerosene fuel spray has been numerically simulated in a laboratory scale combustor geometry to predict soot and the effects of thermal radiation at different swirl levels of primary air flow. The two-phase motion in the combustor is simulated using an Eulerian–Lagragian formulation considering the stochastic separated flow model. The Favre-averaged governing equations are solved for the gas phase with the turbulent quantities simulated by realisable k–? model. The injection of the fuel is considered through a pressure swirl atomiser and the combustion is simulated by a laminar flamelet model with detailed kinetics of kerosene combustion. Soot formation in the flame is predicted using an empirical model with the model parameters adjusted for kerosene fuel. Contributions of gas phase and soot towards thermal radiation have been considered to predict the incident heat flux on the combustor wall and fuel injector. Swirl in the primary flow significantly influences the flow and flame structures in the combustor. The stronger recirculation at high swirl draws more air into the flame region, reduces the flame length and peak flame temperature and also brings the soot laden zone closer to the inlet plane. As a result, the radiative heat flux on the peripheral wall decreases at high swirl and also shifts closer to the inlet plane. However, increased swirl increases the combustor wall temperature due to radial spreading of the flame. The high incident radiative heat flux and the high surface temperature make the fuel injector a critical item in the combustor. The injector peak temperature increases with the increase in swirl flow mainly because the flame is located closer to the inlet plane. On the other hand, a more uniform temperature distribution in the exhaust gas can be attained at the combustor exit at high swirl condition.     

支板凹腔一体化超燃冲压发动机实验研究

本文针对以凹腔支板一体化燃烧室为基本结构的超燃冲压模型发动机在自由射流风洞中的性能,主要研究了燃料在不同位置喷入时,燃烧室几何结构/气动性能/燃料混合及燃烧特性的相互耦合,以及对发动机推力性能的影响.结果表明支板与凹腔的一体化在合理配置燃料分布情况下可以获得较好的发动机性能.     

Investigation on the flameholding mechanisms in supersonic flows: backward-facing step and cavity flameholder

Abstract  

As effective devices to extend the fuel residence time in supersonic flow and prolong the duration time for hypersonic vehicles cruising in the near-space with power, the backward-facing step and the cavity are widely employed in hypersonic airbreathing propulsive systems as flameholders. The two-dimensional coupled implicit RANS equations, the standard k-ε turbulence model, and the finite-rate/eddy-dissipation reaction model have been used to generate the flow field structures in the scramjet combustors with the backward-facing step and the cavity flameholders. The flameholding mechanism in the combustor has been investigated by comparing the flow field in the corner region of the backward-facing step with that around the cavity flameholder. The obtained results show that the numerical simulation results are in good agreement with the experimental data, and the different grid scales make only a slight difference to the numerical results. The vortices formed in the corner region of the backward-facing step, in the cavity and upstream of the fuel injector make a large difference to the enhancement of the mixing between the fuel and the free airstream, and they can prolong the residence time of the mixture and improve the combustion efficiency in the supersonic flow. The size of the recirculation zone in the scramjet combustor partially depends on the distance between the injection and the leading edge of the cavity. Further, the shock waves in the scramjet combustor with the cavity flameholder are much stronger than those that occur in the scramjet combustor with the backward-facing step, and this causes a large increase in the static pressure along the walls of the combustor.     

湍流扩散燃烧的数值研究-PDF方法和火焰面模型的性能比较

采用标量概率密度函数(PDF)方法、稳态和非稳态火焰面模型三种方法对一个值班湍流CH_4/O_2/N_2射流扩散火焰(Sandia Flame D)进行数值计算,以比较不同燃烧模型的性能。PDF方法通过计算反应标量的PDF输运方程来得到标量分布,而火焰面模型只求解单标量混合物分数的PDF方程,组分和温度分布通过火焰面方程的求解或者火焰面数据库的插值得到。计算结果和实验数据对比表明PDF方法计算结果最好但计算量相当大,稳态火焰面模型则反之。综合而言,非稳态火焰面模型的预测结果相对稳态模型有了非常大的改进,而计算量仍然容易接受,非常适合工程应用。     

涡轮导向器对旋转爆轰波传播特性影响的实验研究

为了研究涡轮导向器对旋转爆轰波传播特性的影响,以氢气为燃料,空气为氧化剂,在不同当量比下开展了实验研究.基于高频压力传感器及静态压力传感器的信号,详细分析了带涡轮导向器的旋转爆轰燃烧室的工作模式以及涡轮导向器对非均匀不稳定爆轰产物的影响.实验结果表明:在当量比较低时,爆轰燃烧室以快速爆燃模式工作;逐渐增大当量比,爆轰燃烧室开始以不稳定旋转爆轰模式工作;继续增大当量比,爆轰燃烧室以稳定旋转爆轰模式工作,且旋转爆轰波的传播速度和稳定性均随当量比的增大逐渐提高.爆轰波下游的斜激波与涡轮导向器相互作用,涡轮导向器对压力振荡的幅值具有明显的抑制作用,但对压力振荡频率的影响较小.随着当量比的增大,涡轮导向器上下游的静压均同时增大,经过涡轮导向器的作用,涡轮下游静压明显降低.     

Computational investigations of the coupling between transient flame dynamics and thermo-acoustic instability in a self-excited resonance combustor

Combustion instability due to thermo-acoustic interactions is a critical combustion problem that requires a thorough understanding because of its adverse impact on stable and reliable operation of combustors in high-speed propulsion devices like gas turbines and rockets. This work conducts computational investigations of the coupling between the transient flame dynamics such as the ignition delay and local extinction and the thermo-acoustic instability developed in a self-excited resonance combustor to gain deep insights into the mechanisms of thermo-acoustic instability. A 2D modelling framework that employs different flamelet models (the steady flamelet model and the flamelet/progress variable approach) is developed to enable the examination of the effect of the transient flame dynamics caused by the strong coupling of the turbulent mixing and finite-rate chemical kinetics on the occurrence of thermo-acoustic instability. The models are validated by using the available experimental data for the pressure signal. Parametric studies are performed to examine the effect of the occurrence of the transient flame dynamics, the effect of artificial amplification of the Damköhler number, and the effect of neglecting mixture fraction fluctuations on the predictions of the thermo-acoustic instability. The parametric studies reveal that the occurrence of transient flame dynamics has a strong influence on the onset of the thermo-acoustic instability. Further analysis is then conducted to localise the effect of a particular flame dynamic event, the ignition delay, on the thermo-acoustic instability. The reverse effect of the occurrence of the thermo-acoustic instability on the transient flame dynamics in the combustor is also investigated by examining the temporal evolution of the local flame events in conjunction with the pressure wave propagation. The above observed two-way coupling between the transient flame dynamics (the ignition delay) and the thermo-acoustic instability provides a plausible mechanism of the self-excited and sustained thermo-acoustic instability observed in the combustor despite the fact that the results are obtained from 2D simulations. The same analysis is expected to be extensible to fully 3D simulations.     

Effect of cavity location on combustion flow field of integrated hypersonic vehicle in near space

Abstract  

The cavity has been widely employed as the flame holder to prolong the residence time of fuel in supersonic flows since it improves the combustion efficiency in the scramjet combustor, and also imposes additional drag on the engine. In this paper, the two-dimensional coupled implicit Reynolds Average Navier–Stokes equations, the RNG k–ε turbulence model and the finite-rate/eddy-dissipation reaction model have been employed to numerically simulate the combustion flow field of an integrated hypersonic vehicle. The effect of cavity location on the combustion flow field of the vehicle has been investigated, and the fuel, namely hydrogen, was injected upstream of the cavity on the walls of the first stage combustor. The obtained results show that the viscous lift force, drag force and pitching moment of the vehicle are nearly unchanged by varying the cavity location over the location range and designs considered in this article, namely the configurations with single cavity, double cavities in tandem and double cavities in parallel. The variation of the fuel injection strategy affects the separation of the boundary layer, and the viscous effect on the drag force of the vehicle is remarkable, but the viscous effects on the lift force and the pitching moment are both small and they can be neglected in the design process of hypersonic vehicles. In addition to varying the location of the cavities, three fuel injection configurations were considered. It was found that one particular case can restrict the inlet unstart for the scramjet engine.