A computational study of the HyShot II combustor performance

Experimental and flight data for hypersonic air-breathing engines are both difficult and extremely expensive to obtain, motivating the use of computational models to enhance the understanding of the complex physics involved. Here, a comprehensive numerical study has been carried out for the HyShot II scramjet combustor. This study makes use of Reynolds Average Navier Stokes (RANS) based models and Large Eddy Simulation (LES) based models with semi-detailed reaction kinetics. In this investigation we focus on the underlying flow-mixing-combustion physics at different operating conditions tested in the High Enthalpy Shock Tunnel Göttingen (HEG). To account for the complex flow in the HEG facility a zonal approach is employed in which RANS is used to simulate the flow in the HEG nozzle and test-section, providing the necessary inflow boundary conditions for the combustor RANS and LES, being the focus of this analysis. Specifically, we here combine results from RANS and LES computations with data from the HEG experiments and the target HyShot II flight-tests at two different flight-altitudes (28 and 33 km). The LES model is observed to capture the experimental wall-pressure and heat-flux data very well for both the 33 and 28 km altitude cases, whereas the RANS model is only able to predict the wall-pressure and heat flux data for the 28 km altitude case. Based on the LES results, the flow at both altitudes is found to be unsteady, but with unsteady transitional flow features dominating the 33 km case. Moreover, these results show that the equivalence ratio is of key importance to the resulting flow, mixing and combustion physics, with richer mixtures being prone to transitional flow features. The LES results are also used to describe the flow physics in detail for both altitudes, and the key processes responsible for the transition between the two combustion modes observed.     

Injection and mixing in a scramjet combustor: DES and RANS studies

The low momentum flux ratio jet in the HyShot II scramjet combustor is studied by DES (Detached Eddy Simulation) and RANS (Reynolds-Averaged Navier–Stokes) methods. The flow structure near the injector, shock pattern in the symmetry plane as well as the instantaneous coherent structures are presented and explained. Further insight into the flow physics is obtained by visualizing instantaneous coherent structures. The formation of Ω-shaped vortices, which was previously observed in experiments but never well-studied numerically, is discussed in detail. A new schematic of flow physics is proposed to enhance the understanding of the low momentum flux ratio jet. Compared to the DES result, the RANS method is unable to capture the dynamics of turbulent structures. The DES method provides much detailed information about mixing patterns and a more reliable mixing efficiency than the RANS result. The RANS method over-predicts the eddy-viscosity during turbulence modeling and suppresses unsteady turbulent fluctuations by time averaging, which results in a 25% over-estimation of the mixing efficiency.     

Vortex dynamics in different combustion regions of a cavity-based scramjet

A hybrid RANS/LES study of a cavity-based scramjet was performed and reasonable agreements were found between simulation results and experimental measurements. In the current case, the flame was stabilized by the subsonic cavity shear layer and propagated downstream into the supersonic flow. The vortex dynamic in the flow, mixing, and combustion regions was comparatively investigated. The averaged vorticity in the combustion regions was lower by 55% compared to the mixing region, primarily due to dilatation as a result of the heat release. Furthermore, the combustion zone was decomposed into four regions based on premixed/diffusion flame and subsonic/supersonic combustion. Then the vorticity and its transport in the four regions were compared. The averaged vorticity in the premixed combustion regions was only slightly larger than that in the diffusion combustion regions. However, the averaged heat release rate was nearly 3 times larger in the premixed regions, leading to higher contributions of dilatation and baroclinic torque in the premixed regions, with an overall weak positive impact on the vorticity generation. In the subsonic combustion regions, the vorticity was three times larger than that in the supersonic combustion regions, despite similar heat release rates on average. It could be explained by the relatively large magnitude of dilatation and baroclinic torque in the supersonic flow. Vortex stretching and dilatation were comparable in the supersonic flame but the former became two times larger than the latter in the subsonic flame. Moreover, the baroclinic torque had larger contributions than diffusion in the supersonic flame whereas the opposite trend was found in the subsonic flame. The results highlight that the subsonic combustion regions in the cavity shear layer and near the walls significantly contribute to the vortex dynamics and mixing process, in addition to flame stabilization.     

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.     

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.     

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,由于燃烧室流动速度特别高,燃烧放热靠近燃烧室扩张段下游.      

进口条件对超燃冲压发动机氢气燃烧流场特性的影响分析

对不同进口条件下的超燃冲压发动机燃烧室内氢气喷流超声速燃烧流动特性进行了数值模拟与分析.宽范围超燃冲压发动机是吸气式高超声速飞行器推进系统设计中的热点问题之一,受实验设备硬件条件及实验技术限制,数值模拟技术仍然是超燃冲压发动机燃烧室内燃气燃烧特性及流场特性的主要研究手段.采用基于混合网格技术的多组元N-S方程有限体积方法求解器,在不同进口Mach数及压强条件下,对带楔板/凹腔结构的燃烧室模型氢气喷流燃烧流场进行了数值模拟,对比分析了氢气喷流穿透深度、喷口前后回流区结构、掺混效率及燃烧效率等流场结构与典型流场参数的变化特性及影响规律.研究成果可为宽范围超燃冲压发动机喷流燃烧流动特性分析提供参考.      

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.     

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

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

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.     

Combustion characteristics of a dual-mode scramjet combustor with cavity flameholder

Combustion characteristics of a laboratory dual-mode ramjet/scramjet combustor were studied experimentally. The combustor consists of a sonic fuel jet injected into a supersonic crossflow upstream of a wall cavity pilot flame. These fundamental components are contained in many dual-mode combustor designs. Experiments were performed with an isolator entrance Mach number of 2.2. Air stagnation temperatures were varied from 1040 to 1490 K, which correspond to flight Mach numbers of 4.3–5.4. Both pure hydrogen and a mixture of hydrogen and ethylene fuels were used. High speed imaging of the flame luminosity was performed along with measurements of the isolator and combustor wall pressures. For ramjet mode operation, two distinct combustion stabilization locations were found for fuel injection a sufficient distance upstream of the cavity. At low T₀, the combustion was anchored at the leading edge of the cavity by heat release in the cavity shear layer. At high T₀, the combustion was stabilized a short distance downstream of the fuel injection jet in the jet-wake. For an intermediate range of T₀, the reaction zone oscillated between the jet-wake and cavity stabilization locations. Wall pressure measurements showed that cavity stabilized combustion was the steadiest, followed by jet-wake stabilized, and the oscillatory case. For fuel injection close to the cavity, a hybrid stabilization mode was found in which the reaction zone locations for the two stabilization modes overlapped. For this hybrid stabilization, cavity fueling rate was an important factor in the steadiness of the flow field. Scramjet mode combustion was found to only exist in the cavity stabilized location for the conditions studied.     

Large-eddy simulation of equilibrium plasma-assisted combustion in supersonic flow

A large-eddy simulation (LES) model with a new localized dynamic subgrid closure for the magnetohydrodynamics (MHD) equations is used to investigate plasma-assisted combustion in supersonic flow. A 16-species and 74-reactions kinetics model is used to simulate hydrogen-air combustion and high-temperature air dissociation. The numerical model is validated with experimental data for non-reacting and reacting supersonic flow over a rearward-facing step. The creation of a plasma source near the step corner is shown to have a strong localized effect with the high temperature region resulting in an increase of the radical species concentration in the mixing region. This has the potential for enhancing combustion. In addition, downstream fuel–air mixing is improved, primarily by the creation of a strong baroclinic torque effect in the near field of the plasma source. Furthermore, by adding an uniform external magnetic field, the Lorentz force effect helps to further enhance mixing by lifting the shear layer and increasing fuel penetration by approximately 20%.     

Planar laser-induced fluorescence imaging of OH in a supersonic combustor fueled with ethylene and methane

A supersonic combustor was experimentally investigated using both conventional instrumentation and laser-based diagnostics. Planar laser-induced fluorescence (PLIF) imaging of OH was used in the main section of the combustor to examine flameholding and flame propagation during a series of evaluations at conditions simulating Mach-5.5 flight. Parameters of interest in this study included the angle of the primary fuel injectors, the distribution of fuel throughout the combustor, and the fuel composition. Changes in fuel-injection angle were expected to influence the mixing and combustion processes, and therefore combustor operation. Fuel-distribution variations were expected to modify the flame propagation between flameholding regions. Finally, ethylene and methane were used to examine the suitability of the flameholder designs over a wide range of fuel reactivity. Results suggest that the combustor provides relatively robust flameholding regardless of the fuel used and good flame propagation as long as the fuel distribution provides favorable conditions in the flameholding regions. In addition, the results show that the primary injectors can be useful in controlling certain aspects of combustor operability.     

A multivariate quadrature based moment method for LES based modeling of supersonic combustion

The transported probability density function (PDF) approach is a powerful technique for large eddy simulation (LES) based modeling of scramjet combustors. In this approach, a high-dimensional transport equation for the joint composition-enthalpy PDF needs to be solved. Quadrature based approaches provide deterministic Eulerian methods for solving the joint-PDF transport equation. In this work, it is first demonstrated that the numerical errors associated with LES require special care in the development of PDF solution algorithms. The direct quadrature method of moments (DQMOM) is one quadrature-based approach developed for supersonic combustion modeling. This approach is shown to generate inconsistent evolution of the scalar moments. Further, gradient-based source terms that appear in the DQMOM transport equations are severely underpredicted in LES leading to artificial mixing of fuel and oxidizer. To overcome these numerical issues, a semi-discrete quadrature method of moments (SeQMOM) is formulated. The performance of the new technique is compared with the DQMOM approach in canonical flow configurations as well as a three-dimensional supersonic cavity stabilized flame configuration. The SeQMOM approach is shown to predict subfilter statistics accurately compared to the DQMOM approach.     

Application and theoretical analysis of the flamelet model for supersonic turbulent combustion flows in the scramjet engine

In the framework of Reynolds-averaged Navier–Stokes simulation, supersonic turbulent combustion flows at the German Aerospace Centre (DLR) combustor and Japan Aerospace Exploration Agency (JAXA) integrated scramjet engine are numerically simulated using the flamelet model. Based on the DLR combustor case, theoretical analysis and numerical experiments conclude that: the finite rate model only implicitly considers the large-scale turbulent effect and, due to the lack of the small-scale non-equilibrium effect, it would overshoot the peak temperature compared to the flamelet model in general. Furthermore, high-Mach-number compressibility affects the flamelet model mainly through two ways: the spatial pressure variation and the static enthalpy variation due to the kinetic energy. In the flamelet library, the mass fractions of the intermediate species, e.g. OH, are more sensible to the above two effects than the main species such as H₂O. Additionally, in the combustion flowfield where the pressure is larger than the value adopted in the generation of the flamelet library or the conversion from the static enthalpy to the kinetic energy occurs, the temperature obtained by the flamelet model without taking compressibility effects into account would be undershot, and vice versa. The static enthalpy variation effect has only little influence on the temperature simulation of the flamelet model, while the effect of the spatial pressure variation may cause relatively large errors. From the JAXA case, it is found that the flamelet model cannot in general be used for an integrated scramjet engine. The existence of the inlet together with the transverse injection scheme could cause large spatial variations of pressure, so the pressure value adopted for the generation of a flamelet library should be fine-tuned according to a pre-simulation of pure mixing.     

A quadrature-based LES/transported probability density function approach for modeling supersonic combustion

The joint-scalar probability density function (PDF) approach provides a comprehensive framework for large eddy simulation (LES) based combustion modeling. However, currently available stochastic approaches for solving the high-dimensional PDF transport equation can be error prone and numerically unstable in highly compressible shock-containing flows. In this work, a novel Eulerian approach called the direct quadrature method of moments (DQMOM) is developed for evolving the PDF-based supersonic combustion model. The DQMOM technique uses a set of scalar transport equations with specific source terms to recover the PDF. The new technique is coupled to a compressible LES solver through the energy equation. The DQMOM approach is then used to simulate two practical flow configurations: a supersonic reacting jet and a cavity-stabilized supersonic combustor. Comparisons with experimental data demonstrate the predictive accuracy of the method.     

Large-eddy simulation of shock-turbulence interaction in supersonic diffuser flows

LES based on explicit filtering is used to study the shock train phenomenon in turbulent supersonic diffuser flows with circular cross-section and isothermal wall with an incoming pipe flow at friction Reynolds number 245 and centerline Mach number 1.7. Alternate regions of compression and expansion are found in the shock train which is followed by a shock-free ‘mixing’ region as observed in experiments and simulations in the literature. Turbulence amplification and local peaks in pressure-dilatation correlation are observed in the vicinity of the shocks. Low-frequency oscillations of the shock train are also observed.     

LES modelling of turbulent non-premixed jet flames with correlated dynamic adaptive chemistry

Large eddy simulations (LES) for turbulent flames with detailed kinetic mechanisms have received growing interest. However, a direct implementation of detailed kinetic mechanisms in LES modelling of turbulent combustion remains a challenge due to the requirement of huge computational resources. An on-the-fly mechanism reduction method named correlated dynamic adaptive chemistry (CoDAC) is proposed to overcome this issue. A LES was conducted for Sandia Flame-D, with the reaction mechanism of GRI-Mech 3.0 consisting of 53 species and 325 reactions. The reduction threshold used in LES was obtained a priori by using auto-ignition model and partially stirred reactor (PaSR) with pairwise mixing model. LES results with CoDAC are in good agreement with experimental data and those without reduction. The conditional mean of the number of selected species indicates that a large size of locally reduced mechanism is required in the reaction zone where CH₄ is destructed. A computational time analysis shows that the PaSR model predicts better than the auto-ignition model on the wall time reduction with CoDAC in LES.     

Large eddy simulation of piloting effects on turbulent swirling flames

Pilot flames, created by additional injectors of pure fuel, are often used in turbulent burners to enhance flame stabilization and reduce combustion instabilities. The exact mechanisms through which these additional rich zones modify the flame anchoring location and the combustion dynamics are often difficult to identify, especially when they include unsteady hydrodynamic motion. This study presents Large Eddy Simulations (LES) of the reacting flow within a large-scale gas turbine burner for two different cases of piloting, where either 2 or 6% of the total methane used in the burner is injected through additional pilot flame lines. For each case, LES shows how the pilot fuel injection affects both flame stabilization and flame stability. The 6% case leads to a stable flame and limited hydrodynamic perturbations in the initial flame zone. The 2% case is less stable, with a small-lift-off of the flame and a Precessing Vortex Core (PVC) in the cold stabilization zone. This PVC traps some of the lean cold gases issuing from the pilot passage stream, changes the flame stabilization point and induces instability.     

Simulation of diesel spray combustion using LES and a multicomponent vapourisation model

Diesel spray and combustion in a constant-volume engine cylinder was simulated by a large eddy simulation (LES) approach coupling with a multicomponent vapourisation (MCV) modelling. The simulation focused on the inclusion of the interaction between fuel spray and gas-phase turbulence flow at the sub-grid scale. The LES was based on the dynamic structure sub-grid model, and an additional source term was added to the filtered momentum equation to account for the effect of drop motion on the gas-phase turbulence. The multicomponent drop vapourisation modelling was based on the continuous thermodynamics approach using a gamma distribution to describe the complex diesel fuel composition and was capable of predicting a more complex drop vapourisation process. The effect of gas-phase turbulence flow on the fuel drop vapourisation process was evaluated through the solution of the gas-phase moments of the distribution in the present LES framework. A non-evaporative spray in a constant-volume engine cylinder was first simulated to examine the behaviours of LES, in comparison with a Reynolds-averaged Navier–Stokes (RANS) simulation based on the RNG k–? model. More realistic diesel spray structures and improved agreement on liquid penetration length with the corresponding experimental data were predicted by the LES, using a grid resolution close to that of RANS. A more comprehensive simulation of diesel spray and combustion in cylindrical combustor was also performed. Predicted distributions of soot particles were compared to the experimental image, and improved agreement with the experimental data was also observed by using the present LES and MCV models. Consequently, results of the present models proved that improved overall performance of the fuel spray simulation can be achieved by the LES without a significant increase in the computational load compared to the RANS.