Flashback-induced flame shape transition in a two-stage LPP aeronautical combustor

Large-Eddy Simulations were performed to study the flashback-induced flame shape transition of a lean premixed M flame in a staged liquid-fuelled aeronautical lean-burner, as observed experimentally. The BIMER combustor is a Lean Premixed Prevapourised (LPP) burner composed of two stages, each with its own injector and swirler: the main outer stage, called multipoint, uses jet-in-crossflow injection to achieve the LPP regime, while the central stage, called pilot, uses a pressure swirl injector to create a hollow cone spray to stabilise the flame. During LPP operation, this M flame presents a strong acoustic activity, promoting a periodic flashback of its leading edge. When, aiming to stabilise the flame, the pilot injection is increased and the multipoint injection decreased, the oscillating leading edge (due to the longitudinal acoustic perturbations) attaches to the pilot spray, changing the flame into a Tulip shape. Two phenomena were identified as being the most relevant causes of this flame shape transition. First, the leading edge position and the thermoacoustic instability amplitude are directly linked to the combustion chamber final temperature. The higher the temperature in the chamber, the more upstream the leading edge stabilises, and the higher the acoustic oscillation amplitude, both increasing the risk of a successful flashback. Second, the injection regime with high pilot injection allows the leading edge to attach to the pilot spray, as the flame only transitions when the pilot spray is sufficiently high. The higher the pilot fuel flow, the higher the amount of fuel sprayed in the critical region where the flame might attach for a transition to the Tulip shape. Therefore, as the change in injection regime is the main mechanism lean staged burners use to reduce emissions while increasing operability, this works shows that an M flame is unsuitable to such burners with similar aerodynamic topology and properties.     

LES-CMC simulations of a turbulent bluff-body flame

The large Eddy simulations (LES)-conditional moment closure (CMC) method with detailed chemistry is applied to a bluff-body stabilized flame. Computations of the velocity and mixture fraction fields show good agreement with the experiments. Temperature and major species are well-predicted throughout the flame with the exception of the flow regions in the outer shear layer close to the nozzle where the pure mixing between hot recirculating products and fresh oxidizer cannot be captured. LES-CMC generally improves on results obtained with RANS-CMC and on LES that uses one representative flamelet to model the dependence of reactive species on mixture fraction. Simulated CO mass fractions are generally in good agreement with the experimental data although a 10% overprediction can be found at downstream positions. NO predictions show a distinct improvement over the flamelet approach, however, simulations overpredict NO mass fractions at all downstream locations due to an overprediction of temperature close to the nozzle. The potential of LES-CMC to predict unsteady finite rate effects is demonstrated by the prediction of endothermic—or “flame cooling”—regions close to the neck of the recirculation zone that favours ethylene production via the methane fuel decomposition channel.     

Examination of mode shapes in an unstable model combustor

The coupling between the fluid dynamics, heat addition, and the acoustics of a combustor system determine whether it is prone toward combustion instability. This paper presents results from a benchmark study of the eigenmodes in an unstable experimental combustor. The axisymmetric combustor configuration is representative of a number of practical systems and comprises an injector tube, geometric expansion into a combustion chamber, and a short converging nozzle. Instability limit cycle amplitudes ranged from 5% to nearly 50% of the mean 2.2 MPa pressure. Multiple harmonics were measured for the highly unstable cases. The model combustor was designed to provide a fairly comprehensive set of tested effects: sonic vs subsonic inlets; oxidizer tube lengths that were either quarter-wave, half-wave, or off-resonant acoustic equivalents to the combustion chamber; a significant injector mean flow with Ma∼0.4; and a varied combustion chamber length. The measured mode shape data were analyzed and reduced to provide comparison with results from a linearized one-dimensional Euler model, which included the effects of real boundary conditions, entropy generation, area change, and heat and mass addition, but did not include a model for unsteady heat addition. For low-amplitude instabilities, the measured resonance frequencies agreed with those calculated by the model for the injector tube-combustion chamber system. Resonance frequencies for the high-amplitude oscillation cases corresponded to the first longitudinal frequency of the combustion chamber and its integer multiples. Good quantitative agreement was obtained between computed and measured phase difference profiles, and mode envelopes agreed qualitatively. These results provide a basis for subsequent combustion response studies on the effects of unsteady heat addition.     

Flame dynamics of azimuthal forced spinning and standing modes in an annular combustor

Azimuthal forcing has been applied to flames in a laboratory scale annular combustor in order to accurately control the azimuthal mode of excitation. A new forcing configuration permitted not only the pressure amplitude, but also the spin ratio and mode orientation to be accurately controlled, in order to generate standing modes and for the first time strong spinning modes in both a clockwise (CW) and anti-clockwise (ACW) direction. The phase averaged heat release dynamics of these modes was compared and a number of differences observed depending on the direction of pressure wave propagation, demonstrating characteristic ACW and CW heat release patterns. A new spin compensating averaging method was then introduced to analyse the flame dynamics, and it was shown that through the application of this method the dynamics of standing wave oscillations could be decomposed to recover the characteristic ACW and CW heat release responses. The global heat release response was also assessed during strongly spinning modes, and the magnitude of the response was shown to depend strongly on the direction of propagation, demonstrating the importance of the local swirl direction on the global heat release response, with important implications for the modelling of such flows.     

Assessment of LES of intermittent soot production in an aero-engine model combustor using high-speed measurements

Soot production in turbulent flames is an extremely intermittent phenomenon since it is the result of specific thermochemical conditions occasionally occurring in space and time. In realistic configurations such as the swirling flames used in gas-turbines, the presence of large-scale flow motions can additionally affect soot formation processes, leading to even more pronounced intermittency. Classically, the validation of numerical simulations is performed by comparing time-averaged results with experimental data of the phenomenon under investigation. This comparison can be considered as rigorous only if a statistically converged numerical representation is obtained. In case of sporadic events such as intermittent soot formation in turbulent flames, this means to perform the simulation over thousands of milliseconds of physical time, which can have extremely high CPU demands when performing Large Eddy Simulation (LES). In this work, a possible strategy to overcome this issue is proposed based on the use of high-speed measurements and numerically synthesized signals from LES. To illustrate the approach, numerical and experimental soot light scattering signals are considered here by looking at the model aero-engine combustor developed at DLR for the study of pressurized swirled sooting flames. The light scattering signal is numerically synthesized from an LES. Experimental high-speed measurements are used to statistically account for the high temporal and spatial variability of soot when considering time intervals similar to what is today achievable with LES. The feasibility of this approach is finally demonstrated by comparing numerical results to the ensemble of possible soot production states observed experimentally in the DLR burner allowing to eventually validate the present LES results.     

Detection of early stage damage in carbon fiber reinforced polymers for aeronautical applications using an HTS SQUID magnetometer

We present an experimental characterization of multidirectional fibre composites based on eddy current testing using HTS dc SQUID (Superconducting Quantum Interference Device) magnetometers. The correlation between the mechanical tolerance of CFRPs with different thickness and the phase gradient of the magnetic field generated by damage is shown. The eddy current based SQUID NDE is used to detect damage not visible to the naked eye and to study the effect of an impact far from the point of impact. The aim of this work is to demonstrate the capability of our technique to provide information on early stage damage due to an impact process in an unshielded environment.     

Characteristics of self-excited spinning azimuthal modes in an annular combustor with turbulent premixed bluff-body flames

In this paper we investigate self-excited azimuthal modes in an annular combustor with turbulent premixed bluff-body stabilised flames. Previous studies have shown that both swirl and equivalence ratio influence modal dynamics, i.e. the time-varying nature of the modes. However, self-excited azimuthal modes have not yet been investigated in turbulent flames without bulk swirl, which do not generate any preferential flow in either azimuthal direction, and may therefore lead to different behaviour. Joint probability density functions of the instability amplitudes at various flowrates and equivalence ratios showed a strong bi-modal response favouring both ACW and CW spinning states not previously observed. Operating conditions leading to a bi-modal response provide a unique opportunity to investigate whether the structure of the global fluctuating heat release rate of self-excited spinning modes in both directions exhibit similar dynamics and structure. This was investigated using high-speed OH* chemiluminescence images of the annular combustor and a new rotational averaging method was applied which decomposes the spinning components of the global fluctuating heat release rate. The new rotational averaging, which differs from standard phase-averaging, produces spatial averages in a frame of reference moving with the spinning wave. The results show that the structure of the fluctuating heat release rate for spinning modes is highly asymmetric as characterised by large, crescent shaped regions of high OH* intensity, located on the far side of each flame, relative to the direction of the azimuthally propagating pressure wave. In comparison with interacting swirling flames, these results indicate that the previously observed radial asymmetry of OH* fluctuations may be introduced through advection by local swirl.     

Towards the development of an efficient low-NOx ammonia combustor for a micro gas turbine

Recent studies have demonstrated stable generation of power from pure ammonia combustion in a micro gas turbine (MGT) with a high combustion efficiency, thus overcoming some of the challenges that discouraged such applications of ammonia in the past. However, achievement of low NOx emission from ammonia combustors remains an important challenge. In this study, combustion techniques and combustor design for efficient combustion and low NOx emission from an ammonia MGT swirl combustor are proposed. The effects of fuel injection angle, combustor inlet temperature, equivalence ratio, and ambient pressure on flame stabilization and emissions were investigated in a laboratory high pressure combustion chamber. An FTIR gas analyser was employed in analysing the exhaust gases. Numerical modeling using OpenFOAM was done to better understand the dependence of NO emissions on the equivalence ratio. The result show that inclined fuel injection as opposed to vertical injection along the combustor central axis resulted to improved flame stability, and lower NH₃ and NOx emissions. Numerical and experimental results showed that a control of the equivalence ratio upstream of the combustor is critical for low NOx emission in a rich-lean ammonia combustor. NO emission had a minimum value at an upstream equivalence ratio of 1.10 in the experiments. Furthermore, NO emission was found to decrease with ambient pressure, especially for premixed combustion. For the rich-lean combustion strategy employed in this study, lower NOx emission was recorded in premixed combustion than in non-premixed combustion indicating the importance of mixture uniformity for low NOx emission from ammonia combustion. A prototype liner developed to enhance the control and uniformity of the equivalence ratio upstream of the combustor further improved ammonia combustion. With the proposed liner design, NOx emission of 42?ppmv and ammonia combustion efficiency of 99.5% were achieved at 0.3?MPa for fuel input power of 31.44?kW.     

Leading point behavior during the ignition of an annular combustor with liquid n-heptane injectors

Experimental and numerical investigations of ignition in combustors with multiple burners have recently emerged and have provided new insights on the last phase of ignition in gas turbine-like annular geometries where the flame propagates from burner to burner. Previous comparisons between calculations and experiments of light-round in a laboratory scale annular combustion chamber have demonstrated the ability of large-eddy simulation to predict such processes for perfectly premixed conditions and, more recently, for n-heptane spray injection. The present analysis focuses on two additional operating points with liquid n-heptane sprays and the turbulent flame propagation in the two-phase mixture is examined through the behavior of its leading points. The validation of the light-round process is characterized in terms of ignition delays. The detailed analysis of the propagation through the definition of a leading point enables to highlight some key phenomena responsible for the flame behavior, such as the influence of the liquid droplet spray and its vaporization in the chamber. Calculations indicate that the volumetric expansion due to the chemical reaction at the flame induces a strong azimuthal flow in the fresh stream at a distance of several sectors ahead of the flame, which modifies conditions in this region. This creates heterogeneities in the gas composition and wakes on the downstream side of the swirling jets formed by the injectors, with notable effects on the motion of the leading point and on the absolute flame velocity.     

Numerical analysis of self-excited tangential combustion instability for an MMH/NTO rocket combustor

This study presented the numerical simulation of the tangential combustion instability in a hypergolic liquid bipropellant rocket thrust chamber, which applied fuel liquid film cooling method and unlike impinging injectors. The liquid spray was modeled using Lagrangian approach, while the gas was regarded as Euler phase. Stress-blended eddy simulation and finite rate/eddy–dissipation model were adopted to simulate the turbulent combustion process. Consistent with the experiment results, this work successfully simulated the transformation of tangential combustion instability from standing mode to spinning mode. The mean pressure, amplitude and frequency of limit cycle oscillation were in good agreement with the experiment. There was a detailed analysis about the flow field, Rayleigh index, and driving mechanism of the combustion instability. It was found that the oscillation began with hot spots of heat release rate due to the interaction between the spray of impinging injectors and cooling fuel jet. More than that, cooling fuel jet also contributed to drive the oscillation. In the standing mode, injectors in the inner and outer rings drive the oscillation together, while the spinning mode is mainly driven by injectors in the outer ring. The pressure wave is subsonic and its Mach number is close to 1. It was shown that the pressure wave contained a complex structure divided into three parts. This led to the in-phase of the pressure along the axial direction and the double-peak characteristic of the downstream pressure signal. Besides, a positive feedback closed-loop system associated with periodic oxidizer/fuel ratio was believed to sustain the combustion instability. The oscillation can be maintained when pressure, heat release and oxidizer/fuel ratio are coupled together. The analysis results indicate that rotating detonation is an implication to tangential combustion instability.     

Scaling of heat transfer in gas-gas injector combustor

The scaling of heat transfer in gas-gas injector combustor is investigated theoretically, numerically and experimentally based on the previous study on the scaling of gas-gas combustion flowfield. The similarity condition of the gas-gas injector combustor heat transfer is obtained by conducting a formulation analysis of the boundary layer Navier-Stokes equations and a dimensional analysis of the corresponding heat transfer phenomenon. Then, a practicable engineering scaling criterion of the gas-gas injector combustor heat transfer is put forward. The criterion implies that when the similarity conditions of inner flowfield are satisfied, the size and the pressure of gas-gas combustion chamber can be changed, while the heat transfer can still be qualitatively similar to the distribution trend and quantitatively correlates well with the size and pressure as q ∝ pc0 .8d t-0.2. Based on the criterion, single-element injector chambers with different geometric sizes and at different chamber pressures ranging from 1 MPa to 20 MPa are numerically simulated. A single-element injector chamber is designed and hot-fire tested at seven chamber pressures from 0.92 MPa to 6.1 MPa. The inner wall heat flux are obtained and analysed. The numerical and experimental results both verified the scaling criterion in gas-gas injector combustion chambers under different chamber pressures and geometries.     

Simulation of spray combustion in a lean-direct injection combustor

Large-eddy simulation (LES) of a liquid-fueled lean-direct injection (LDI) combustor is carried out by resolving the entire inlet flow path through the swirl vanes and the combustor. A localized dynamic subgrid closure is combined with a subgrid mixing and combustion model so that no adjustable parameters are required for both non-reacting and reacting LES. Time-averaged velocity predictions compare well with the measured data. The unsteady flow features that play a major role in spray dispersion, fuel–air mixing and flame stabilization are identified from the simulation data. It is shown that the vortex breakdown bubble (VBB) is smaller with more intense reverse flow when there is heat release. The swirling shear layer plays a major role in spray dispersion and the VBB provides an efficient flameholding mechanism to stabilize the flame.     

Large-eddy simulation of evaporating spray in a coaxial combustor

Large-eddy simulation of an evaporating isopropyl alcohol spray in a coaxial combustor is performed. The Favre-averaged, variable density, low-Mach number Navier-Stokes equations are solved on unstructured grids with dynamic subgrid scale model to compute the turbulent gas-phase. The original incompressible flow algorithm for LES on unstructured grids by [Mahesh et al., J. Comp. Phys. 197 (2004) 215–240] is extended to include density variations and droplet evaporation. An efficient particle-tracking scheme on unstructured meshes is developed to compute the dispersed phase. Experimentally measured droplet size distribution and size-velocity correlation near the nozzle exit are used as the inlet conditions for the spray. The predictive capability of the LES approach on unstructured grids together with Lagrangian droplet dynamics models to capture the droplet dispersion characteristics, size distributions, and the spray evolution is examined in detail. The mean and turbulent quantities for the gas and particle phases are compared to experimental data to show good agreement. It is shown that for low evaporation rates considered in the present study, a well resolved large-eddy simulation together with simple subgrid models for droplet evaporation and motion provides good agreement of the mean and turbulent quantities for the gas and droplet phases compared to the experimental data. This work represents an important first step to assess the predictive capability of the unstructured grid LES approach applied to spray vaporization. The novelty of the results presented is that they establish a baseline fidelity in the ability to simulate complex flows on unstructured grids at conditions representative of gas-turbine combustors.     

PIV measurement of velocity field in a spray combustor

This paper reports a velocity measurement technique using PIV for application to a luminous flame in a spray combustor. The present system consists of a standard PIV system, a rotary shutter and a band-pass filter, the combination of which removes the influence of the high intensity of the luminous flame. The effectiveness of the rotary shutter is studied by changing the shutter speed from 2 ms to 37 ms. The simultaneous observation of the velocity field and the flame structure was carried out in the combustor model for a boiler. The measured velocity field indicates that the exit velocity from the burner is increased by chemical reactions, but the flow pattern inside the combustor is kept similar to that without combustion.     

Statistical characteristics of pressure oscillations in a premixed combustor

This paper describes an investigation of the statistical characteristics of self-excited and noise-driven pressure oscillations in a premixed combustor. This work was motivated by observations that certain characteristics of these oscillations appear random and cannot be entirely characterized within a deterministic framework (e.g., spontaneous, noise-induced transitions of the combustor from stable to unstable operation or cycle-to-cycle variations in the oscillating pressure). In an effort to elucidate these stochastic elements, we performed an analysis of cycle-to-cycle variations in combustor pressure whose results are described in this paper. Data obtained from our combustor shows that the probability density function of the amplitude of these oscillations transitions from a Rayleigh to a Gaussian-type distribution as the combustor moves from stable to unstable operation. These data also show that the instability phase is nearly uniformly distributed; i.e., there is no phase value with maximum probability of occurrence. We also describe a theoretical analysis of the statistical features of a non-linear combustor model that is forced by random noise. Solutions of this model are presented and shown to be in agreement with measured data. The good agreement between the predictions and measured data suggest that the analysis presented in this paper provides a useful framework for interpreting many other apparently random features of combustor stability characteristics; for example, cyclic variability, “fuzziness” in stability boundaries, or noise-induced transitions.     

高空离轨发动机流场红外辐射特性研究

高超声速飞行器大攻角机动时,其离轨发动机产生的喷流与高速稀薄的大气来流产生强烈干扰,流场情况复杂,流场红外辐射也是天基红外系统探测的标志性事件.本文针对高超声速飞行器发动机喷流与稀薄来流的相互干扰情况,采用数值求解Navier-Stokes方程模拟干扰流场,采用逐线积分法得到气体红外辐射特性,结合反向蒙特卡洛方法计算得...     

化学激光器燃烧室两步反应法的实验研究

在环柱型增益发生器实验系统燃烧室的设计中,将主稀释剂He从两个不同的区域注入燃烧室,其中一部分与D2预混以后通过D2喷注孔注入,另一部分通过主喷管叶片后端的后向喷注孔注入。讨论了不同区域的主稀释剂流量的变化对激光器特性的影响。实验结果表明：后He流量改变时,激光器的功率变化及点火稳定性等明显优于前He流量改变的情况。若要增加主稀释剂比例,应优先考虑调整后He。