In this work, we propose and test a partial premixed fuel injection design of NH 3 H 2 O 2 with double ring-shaped inlets to enhance ammonia combustion in an open-ended combustor by generating and sustaining pulsating combustion oscillations. Emphasis is being placed on determining the minimum amount of hydrogen being blended with ammonia in the presence of such self-excited pulsating oscillations. With the numerical model validated by comparing with experimental and theoretical data, we identify and systemically investigate three key thermodynamic parameters. They are shown to strongly affect the thermal, combustion, and emission performances. These parameters include the following: (1) total fuel mass flow rate m ̇ f; (2) mass fraction of hydrogen ω ̇ H 2; and (3) the temperature TH of a heat exchanger implemented downstream of the combustor. It is interesting to observe that intermittent pulsating oscillations are sustained by such ammonia–hydrogen combustion. Furthermore, comparison is conducted between the present results and those with the classical single ring-shaped fuel inlet under the same flow and operating conditions. It is found that the exothermic heat of the proposed double-ring inlets is increased by 98.7% on average. The frequency of such intermittent oscillations is shown to increase with the decreased NH 3 proportion. When pure hydrogen is supplied and passing through the outer ring inlet, the combustion limit can be greatly expanded, even if the inlet mass fraction of hydrogen is very small. The minimum hydrogen blended with ammonia is shown to be 0.1% to achieve a sustainable combustion and large-amplitude oscillations. The NO emission is found to be decreased, and H 2 O is shown to increase. The present study open ups an approach to enhance ammonia combustion by improving its flammability limit with the minimal hydrogen blended.

The greenhouse effect induced by the combustion of fossil fuel characterized with hydrocarbon components is an urgent environmental concern. Ammonia, NH 3, as a renewable energy with a high proportion of hydrogen content, may be utilized to achieve net-zero coalition by reducing CO 2 emission dramatically or completely. As one of alternative energy sources of the conventional hydrocarbon fuel, ammonia is not only a renewable resource with a high volume and a large energy density but also easy to be produced and processed. The key process for ammonia synthesis is known as the Haber–Bosch process.1,2 It is one of the most promising carbon-free options to supply thermal power for fuel cells, internal combustion engines, ships, next-generation aircrafts, and many other power generation systems.1 However, due to its relatively narrow flammability range, a high ignition energy, a low flame temperature, and a low flame propagation speed,2 it is still quite challenging to be directly utilized as carbon-free fuel in practical combustion applications3,4 on the conventional engines, which are designed to accommodate the classical hydrocarbon fuels. Although the combustion of ammonia does not produce CO 2, the emission of NOx is the main contributor to photochemical pollution. How to reduce NOx production and emission while improving the energy conversion efficiency is another critical problem to be solved.5,6

In order to broaden the flammability range of ammonia fuel7–9 and to improve its flame/combustion speed, many methods/means have been attempted.10–12 Two of them have obvious effects and been widely investigated. One way is to ignite ammonia8 by changing the horizontal flow to the swirling flow11–13 at the combustor inlet. As the vortex decomposition generated, downstream of the swirling combustor at the outlet can generate a recirculating flow region. The heat release (thermal energy) and active radicals can be recirculated to the bottom of the flame, which effectively increases the flame limit. Therefore, a stable NH 3–air flame can be obtained within a wide range of the equivalence ratio.14,15 In order to enable gas turbines being fueled and operated with ammonia, the flame stability and NO x emissions of premix NH 3/air combustion are experimentally studied under ambient temperature and pressure.16,17 Their results show that the flammability limits of such ammonia combustion are close to those of lean and rich combustion limits at different equivalence ratios and inlet velocities.17 Furthermore, the blowout limits of the flame are shown to be related to the swirling number and the geometry shape of the combustor.18 It has been found from the previous analysis of the emission composition that the emissions of NO and NH 3 are in the same quantities under a specific equivalence ratio, when the oxygen is deficient. Furthermore, experimental studies were conducted on the structure and emission of NH 3/air19 swirling flame in single-stage and two-stage combustors. It is found that the combustor wall heat loss inhibited the NO emission. However, it did promote the emission of NH 3 and N 2 O. Furthermore, PLIF (planar laser-induced fluorescence) imaging and FTIR (Fourier transform infrared spectroscopy) spectral analyses showed that the wall heat transfer had a greater impact on the flame than the operating pressure. Although swirling combustion can increase the flammability limit of NH 3,13 it may be challenging to be applied in micro-20 and meso-scale combustors or large-scale practical propulsion systems.21–23 

The other means to significantly improve the flammability range and to reduce the ignition delay speed is mixing ammonia with fuels of high calorific value such as hydrogen,13,23 dimethyl ether (DME),16 methane ( CH 4), or propane ( C 3 H 8). The burning performance of ammonia is found to be improved by co-firing hydrogen with ammonia13 via enhancing the burning velocity and reducing NO x emissions.17 It is also found that the laminar burning velocity (unstretched) is increased nonlinearly with the increased hydrogen fraction. Experimental investigations on the combustion characteristics of NH 3/ H 2/air fuels with different mixing ratios in a high-pressure gas turbine combustor were conducted.14,23 It was found that free radicals are related to the combustion speed. Furthermore, OH is shown to play the most significant role on ammonia decomposition.14 In addition, the heat release rate of NH 3 flame is positively correlated with the concentration of OH and N free radicals.15 The previous studies14,15,17 have shown that if the molar ratio of NH 3 to H 2 is approximately 7:3, then it can be applied as a stable energy source for gas turbine engines. Additionally, experimental studies14 are conducted on the flame propagation limit of NH 3/ CH 4/air in a fan-stirred fixed-volume burner under turbulent conditions. It is found that the flame propagation limit is increased due to the addition of ammonia (gas). In addition, the ammonia flame can still propagate with a high turbulent intensity, even when the equivalence ratio is reduced to 0.9. Finally, the experimental studies on spark ignition engines21 powered by NH 3/ H 2/air have confirmed that ammonia flame stability and combustion heat cycle can be significantly improved, if the molar fraction (volume) of hydrogen being blended is exceeding 20%.

Combustion instability24–27 is an undesirable self-excited thermo-acoustics coupling phenomenon. It is quite often observed to occur in low-emission premixed combustion systems for power generation and propulsion applications.26,27 It is characterized by self-sustained large-amplitude pulsating pressure/velocity oscillations.28–30 For the low-emission premixed combustion systems, self-excited pressure/velocity oscillations may occur, when the unsteady heat release is coupled with one or more acoustic resonant modes31 or periodic flow/mass fluctuations32–34 resulting from flame-holders or the geometric shapes of the combustor.35 Its main characteristics include (1) periodic flame motions and fluctuating thermodynamic properties, and (2) large-amplitude flow and flame reversal and overheating.36–38 These pulsating oscillations can promote heat and mass transfer and so combustion to a certain extent. It may also improve the efficiency of chemical energy conversion. However, at the same time, it will significantly reduce the engines' lifespan and the safe operation range.39–42 To prevent the occurrence of such instability, passive control methods,41 such as perforated liners and an electric heater, may be applied to dampen such pulsating oscillations by changing the acoustic damping or varying the temperature fields of the combustion systems.38,43,44 It has been shown that promising results were achieved, as an electric heater is applied as an “acoustic damper” in an open-ended thermoacoustic combustor.45–47 

To the best knowledge of the present authors, there are few investigations, which focus on enhancing the combustion performance of NH 3 with minimum hydrogen blended in the presence of self-excited pulsating oscillations. In this work, an open–open longitudinal combustor with double ring-shaped inlets is proposed and numerically examined. Emphasis is placed on addressing the key question of how to reduce the blended hydrogen fraction to a minimal value while sustaining a stable ammonia combustion. For this, three  key thermodynamic parameters are identified. Their effects on triggering pulsating combustion oscillations and combustion emissions are evaluated one at a time. The structure of this paper is as follows. In Sec. II, the model and numerical method are described. In additional, the model validation is conducted. In Sec. III, the effect of the total inlet fuel flow rate m ̇ f = m ̇ NH 3 + m ̇ H 2 is examined. In Sec. IV, the effect of the blended hydrogen mass fraction ω ̇ H 2 in the mixture of both ammonia and hydrogen is evaluated to reveal the amplitude and frequency of the pulsating oscillations varied with hydrogen content. In Sec. V, the effect of the downstream heat exchanger surface temperature TH on producing pulsating combustion oscillations is investigated. In Sec. VI, the minimal hydrogen being blended is determined in the presence of such pulsating oscillations and a sustainable combustion. In Sec. VII, the emissions of combustion products such as NO and NO 2 and the reactants are predicted in the presence of such pulsating oscillations, as the minimum hydrogen ω ̇ H 2 min is provided. Finally, the key findings are summarized in Sec. VIII.

An acoustically open–open combustor with a double ring-shaped fuel injector and heat exchangers confined is considered in this work. It is schematically shown in Fig. 1. The length (L) and diameter (D) of the combustor are 1000 and 50 mm, respectively. The geometry and dimensions are chosen by following our previous experimental setup.48 The gases of NH 3 / H 2/air are premixed in a certain proportion and passes into the double-ring inlets. The air (consisting of O 2 and N 2 only) is continuously squeezed into the combustor at the bottom due to the pressure difference. The air inlet width (L1) is 8.75 mm, and the fuel mixture inlet width (L2) is 5 mm. Three heat exchangers with a constant surface temperature (powered by electricity in practice) are fixed 125 mm away from the combustor outlet (Lh) to amplify or attenuate the pulsating oscillations.

FIG. 1.

Schematic of the open–open thermoacoustic combustor with double ring-shaped fuel injectors confined, and the 2D meshes at different locations of the combustor.

FIG. 1.

Schematic of the open–open thermoacoustic combustor with double ring-shaped fuel injectors confined, and the 2D meshes at different locations of the combustor.

Close modal

The 2D (two-dimensional) structure grid is generated by Hypermesh 14.0. In addition, more refined grids are applied near the entrance and the heat exchangers in order to improve the simulation accuracy. In this work, Reynolds-averaged Navier–Stokes (RANS) method and k–ϵ turbulence model are used to solve all coupled governing equations. The third-order MUSCL scheme except for pressure and turbulence dissipation (second-order upwind) is used to solve the pressure–velocity coupling equations. In addition, the eddy dissipation concept (EDC) including detail chemical mechanisms of 19 species and 63-step reactions is selected for the present model. In addition, all cases use the stable solver for calculation first. After obtaining the stable flame, it is changed to the transient solver. The time step is set as Δ t = 1 × 10 4 s. The numerical simulations are demonstrated by ANSYS Fluent R19.3. The mass fraction of oxygen in the mixture is fixed at 11.6% to ensure NH 3 burning at a high concentration.

The flow properties in the combustor involve with a mean (time-average) and a fluctuating part. For this, we study the effects of fluctuations on the mean flow properties by using Reynolds decomposition. The thermodynamic variables consist of a mean and fluctuating components as u = U + u and p = P + p ( t ). Thus, the governing equations for turbulent compressible flow include thermodynamic equations and the unsteady Reynolds-averaged Navier–Stokes ones, which include (1) mass conservation, (2) Reynolds-averaged momentum equations, (3) energy conservation, and (4) species conservation.

The Reynolds stress includes ρ u 2 ¯ , ρ v 2 ¯, ρ w 2 ¯ , ρ u v ¯, etc. Such turbulent stresses are shown to increase with the mean rate of deformation. Boussinesq proposed the equation to describe the deformation relationship. The standard κ ϵ model uses the transport equations for κ and ϵ. More details could be found in Ref. 38.

In order to guarantee that the number of selected grids can lead to sufficiently accurate results, we choose 1.5 × 10 4, 3 × 10 4, and 6 × 10 4 grids for the mesh independence study. The results of the grid sensitivity investigations are shown in Figs. 2(a) and 2(b). It can be observed from the temperature and velocity variation along the axial direction of the combustor that the deviation of the numerical results obtained from the different meshes is shown near the inlet and outlet of the combustor. The largest deviation is observed, when the 1.5 × 10 4 grids are used. However, as the other two numbers of numerical grids are applied, the results are well agreed and consistent.

FIG. 2.

Grid independence study of the predicted (a) temperature and (b) velocity along the centerline of the combustor; (c) validating the model by comparing the numerically predicted longitudinal mode-shapes with the theoretical results obtained from Ref. 34.

FIG. 2.

Grid independence study of the predicted (a) temperature and (b) velocity along the centerline of the combustor; (c) validating the model by comparing the numerically predicted longitudinal mode-shapes with the theoretical results obtained from Ref. 34.

Close modal
In order to further verify the current model, the longitudinal pressure mode-shapes of the combustor are obtained at ω 1 211 Hz and ω 2 422 Hz, which are normalized by its maximum value. This is shown in Fig. 2(c). Comparison is then made between the present results and the classical standing-wave longitudinal mode-shape as theoretically determined34 as
(1)
It can be seen that although there is a slight difference in the axial distribution between the peak pressure and the theoretical value, the overall agreement is acceptable.

Furthermore, comparison is conducted by comparing the numerical predicted frequency of the dominant mode with the experimentally measured ones49 from the open-ended combustor with a single ring-shaped inlet, as the combustor total length is varied. This is summarized in Table I. It can be seen that the difference is less than 10%. More detailed comparison and validations on various turbulence models and combustion models could be found in our previous work.38 It has already been shown38 that the difference is less than 10% in terms of the frequency and amplitude of the dominant mode predicted from the present 2D model with the 3D ones. In general, it is confirmed that the model is applicable to predict the self-excited pulsating combustion oscillations in the open-ended combustors in the presence of the single ring-shaped inlet or two annulus ring-shaped inlets.

TABLE I.

Comparison of the numerical predicted frequency of the dominant mode with the experimentally measured ones.49 

Axial total length (mm) Numerical prediction (Hz) Experimental measurements (Hz)
L = 900  ω d Num / 2 π = 217.4 Hz  ω d Exp / 2 π = 225 Hz 
L = 1200  ω d Num / 2 π = 160.2 Hz  ω d Exp / 2 π = 169 Hz 
Axial total length (mm) Numerical prediction (Hz) Experimental measurements (Hz)
L = 900  ω d Num / 2 π = 217.4 Hz  ω d Exp / 2 π = 225 Hz 
L = 1200  ω d Num / 2 π = 160.2 Hz  ω d Exp / 2 π = 169 Hz 

Figure 3 shows the time evolution of the acoustic pressure fluctuations, as the inlet fuel mass flow rate m ̇ f is set to three different values. It can be seen from the phase diagrams that initial negligible pressure disturbances grow into large-amplitude periodic limit cycles, as indicated by the circle shape. As the mass flow rate m ̇ f is varied, the small-amplitude pressure perturbations take different times to grow up. This indicates that the growth rates are different and related to m ̇ f. More details will be examined in Secs. III and IV.

FIG. 3.

Preliminary results of (a) the pressure fluctuations in time-domain, as the inlet fuel mass flow rate m ̇ f is set to three different values, (b) phase diagrams of the pressure fluctuations at m ̇ f = 0.5 × 10 5 kg/s, (c) phase diagrams of the pressure fluctuations at m ̇ f = 1.0 × 10 5 kg/s, and (d) phase diagrams of the pressure fluctuations at m ̇ f = 2.0 × 10 5 kg/s.

FIG. 3.

Preliminary results of (a) the pressure fluctuations in time-domain, as the inlet fuel mass flow rate m ̇ f is set to three different values, (b) phase diagrams of the pressure fluctuations at m ̇ f = 0.5 × 10 5 kg/s, (c) phase diagrams of the pressure fluctuations at m ̇ f = 1.0 × 10 5 kg/s, and (d) phase diagrams of the pressure fluctuations at m ̇ f = 2.0 × 10 5 kg/s.

Close modal

Figures 4(a) and 4(b) depict the time evolution of the mode-shapes of pressure fluctuations and velocity fluctuations, respectively. It can be seen that longitudinal mode-shapes are present in the combustor, as revealed by both pressure and velocity fluctuations in the normalized axial location x/L. However, the velocity mode-shapes involve with more local and small-amplitude wavy peaks in comparison with that of the pressure mode-shaped. This reveals that nonlinearity does occur to the acoustic velocity fluctuations.

FIG. 4.

Preliminary results of (a) pressure and velocity (b) fluctuations along the axial direction of the combustor at 0.59 t 0.6 s.

FIG. 4.

Preliminary results of (a) pressure and velocity (b) fluctuations along the axial direction of the combustor at 0.59 t 0.6 s.

Close modal

For comparison, the conventional single ring-shaped inlet is chosen as the benchmark case. Figures 5(a)–5(d) illustrates the comparison of the temperature contours of the conventional single ring-shaped and double ring-shaped inlet combustors, as the total inlet fuel mass flow rate m ̇ f is set to the same value. Note that the inlet mass flow ratio for the double ring-shaped inlets is set to 1:1 for simplicity. It can be seen that, under the same m ̇ f, the flame of the double ring-shaped inlets is shorter but “fatter” than that of the single ring-shaped one. However, the downstream temperature is significantly higher. This indicates that the double ring-shaped inlet structure is more conducive to the full combustion of the non-carbon fuel, thus releasing more thermal energy. Same conclusions can be obtained from Figs. 5(c) and 5(d), as m ̇ f is increased to 1.0 × 10 5 kg/s from 0.5 × 10 5 kg/s [see Figs. 5(a) and 5(b)]. The increase in the fuel mass flow rate leads to a longer and larger flame, and so a larger high-temperature region. In addition, a closer observation shows that the flame downstream of the double ring-shaped inlet combustor is detached from the fuel inlet at some distance downstream of the flame-holder.

FIG. 5.

Numerically predicted temperature contour distributions of (a) single ring-shaped inlet and (b) double ring-shaped inlets, as the inlet fuel mass flow rate m ̇ f is set to 0.5 × 10 5 kg/s (a) and (b), 1.0 × 10 5 kg/s (c) and (d); the corresponding mole fraction contours of OH radicals are shown in (e)–(h).

FIG. 5.

Numerically predicted temperature contour distributions of (a) single ring-shaped inlet and (b) double ring-shaped inlets, as the inlet fuel mass flow rate m ̇ f is set to 0.5 × 10 5 kg/s (a) and (b), 1.0 × 10 5 kg/s (c) and (d); the corresponding mole fraction contours of OH radicals are shown in (e)–(h).

Close modal

Figures 5(e)–5(h) shows the corresponding mole fraction distribution contours of OH radical. It can be seen that the mole fraction of OH radical in the double ring-shaped inlet structure is larger, when m ̇ f is set to be the same. This phenomenon becomes more apparent with increased m ̇ f [Figs. 5(e) and 5(g) in comparison with (f) and (h)]. As more OH radicals are formed, the heat of reaction is also increased. This explains why the average temperature in the combustor with the double ring-shaped inlets is always higher than that in the combustor with the single inlet.

The mass fraction of ammonia and hydrogen may greatly affect the combustion heat release, thus changing the frequency and amplitude of self-excited oscillations. Figure 6 shows the pressure fluctuation at the center point of the central axis with different hydrogen mass fraction w ̇ H 2 and the corresponding recurrence plots (RP) [see Figs. 6(b), 6(d), and 6(f)]. It can be seen that the intermittent oscillations became more obvious, as w ̇ H 2 is increased from 15% to 88.4%. In addition, the oscillation period is increased. The minimum amplitude is reduced to 10 Pa. As observed from the recurrence plot [see Figs. 6(b), 6(d), and 6(f)], the combustion oscillations are quickly increased and grow into the limit cycle, even when the hydrogen mass fraction is quite small. With the increase in the hydrogen fraction w ̇ H 2, intermittent pulsating oscillations occur. The minimum amplitude of local peaks is decreased dramatically, as observed from the dark regions variation in the recurrence plots of Figs. 6(b), 6(d), and 6(f), which are corresponding to the limit cycles as shown in Figs. 6(a), 6(c), and 6(e), respectively.

FIG. 6.

Time evolution of the acoustic pressure, as m ̇ f = 2.0 × 10 5 kg/s, and the fuel mass flow rate m ̇ f is set to the mass fraction of H 2, which is set to be (a) w ̇ H 2 = 25 %, (c) w ̇ H 2 = 45 %, and (e) w ̇ H 2 = 88.4 %; (b), (d), and (f) are the corresponding RP (recurrence plot).

FIG. 6.

Time evolution of the acoustic pressure, as m ̇ f = 2.0 × 10 5 kg/s, and the fuel mass flow rate m ̇ f is set to the mass fraction of H 2, which is set to be (a) w ̇ H 2 = 25 %, (c) w ̇ H 2 = 45 %, and (e) w ̇ H 2 = 88.4 %; (b), (d), and (f) are the corresponding RP (recurrence plot).

Close modal

In the absence of limit cycles, the recurrence plots are expected to be involved with nodes distributed randomly. As the periodic oscillations are present, there are a number of tiny “gray” square blocks as observed from Fig. 6(b) along the diagonal (principal) in the RP. It is related to the dominant high-frequency periodic conversion between unsteady heat and sound waves. There are a number of connected horizontal and vertical lines. As multiple frequency tones are generated (i.e., intermittent oscillations), especially with the extreme low-frequency tone, as shown in Figs. 6(d) and 6(f), there are a few large black square blocks along the diagonal direction. The square edges create a network with the nodes entangled with each other, and they intersect each other. The square length is proportional to the period of the extreme low-frequency tone. It is worth noting that this kind of intermittent pulsating oscillations with such large-amplitude periodic variation should be avoided in practical applications. When it occurs, the thermodynamic properties such as density and velocity and pressure are all violently fluctuating, as shown in Fig. 7. These large-amplitude fluctuations are easy to lead to mechanical fatigue problems by inducing deformation or fracture/crack of engine casing material.

FIG. 7.

Phase diagram of the acoustic pressure fluctuation, instantaneous velocity, and instantaneous density, as m ̇ f = 2.0 × 10 5 kg/s and w ̇ H 2 = 45 %; TH of the heat exchanger is set to 300 K.

FIG. 7.

Phase diagram of the acoustic pressure fluctuation, instantaneous velocity, and instantaneous density, as m ̇ f = 2.0 × 10 5 kg/s and w ̇ H 2 = 45 %; TH of the heat exchanger is set to 300 K.

Close modal
Here, we apply the heat exchangers downstream in the combustor, aiming to promote or attenuate the combustion-driven oscillations to affect the less flameable ammonia combustion and NOx emission. Here, the influence of the temperature TH of the heat exchangers on the oscillations amplitude is examined first, as TH is kept at 300 K during t 4.0 s. At t = 4.0 s, TH is linearly raised to 1300 K within 0.1 s. Figures 8(a) and 8(b) illustrate the variations of the pressure amplitude and the acoustical energy inside the burner with the temperature change of the heat exchangers, respectively. It can be seen that with increased TH, the maximum and minimum amplitudes of the pressure intermittent oscillation amplitude are reduced to some extent. The corresponding minimum acoustical energy is observed to reduce to 0, as t > 4 s. Note that the acoustical energy is defined as the sum of the acoustic kinetic and acoustic potential energy [see Eq. (57) in Ref. 45].
(2)
where ρ ¯ is the mean density, u is the velocity fluctuation, c ¯ is the mean speed of sound, and A c = π D 2 / 4 is the cross-sectional area of the cylindrical combustor.
FIG. 8.

Time evolution of the acoustic pressure fluctuations (a) and acoustical energy Ea (b), as m ̇ f = 2.0 × 10 5 kg/s and w ̇ H 2 = 45 %; TH of the heat exchanger is set to 300 K at t 4.0 s and slowly increased to 1300 K over 0.1 s and then remains constant at 4.1 t 8.0 s.

FIG. 8.

Time evolution of the acoustic pressure fluctuations (a) and acoustical energy Ea (b), as m ̇ f = 2.0 × 10 5 kg/s and w ̇ H 2 = 45 %; TH of the heat exchanger is set to 300 K at t 4.0 s and slowly increased to 1300 K over 0.1 s and then remains constant at 4.1 t 8.0 s.

Close modal

Figure 9 compares the frequency spectrum of the heat of reaction released by the premixed chemical combustion and the heat flux of the heat exchangers, before and after changing the temperature TH at t = 4.0 s. It can be seen that their changes are closely related to each other. Moreover, with increased TH, the amplitudes of the peaks of both the heat of reaction and the heat flux are decreased. Furthermore, the variation ranges of the dominant peak frequencies are basically the same.

FIG. 9.

Comparison of the frequency spectrum of the normalized (a) heat of reaction and (b) the heat flux of the heat exchangers, as m ̇ f = 2.0 × 10 5 kg/s and ω ̇ H 2 = 45 %, and the temperature TH of the heat exchanger is varied from 300 K at t 4.0 s to 1300 K at t > 4.0 s.

FIG. 9.

Comparison of the frequency spectrum of the normalized (a) heat of reaction and (b) the heat flux of the heat exchangers, as m ̇ f = 2.0 × 10 5 kg/s and ω ̇ H 2 = 45 %, and the temperature TH of the heat exchanger is varied from 300 K at t 4.0 s to 1300 K at t > 4.0 s.

Close modal

There is a strong interest in determining the minimum hydrogen blended with ammonia to achieve a stable combustion. For this, we have reduced the hydrogen mass flow rate ω ̇ H 2 with respect to the ammonia mass flow rate to 0.1%, as self-excited pulsating oscillations are present. Here, the total flow rate m ̇ f of both ammonia and hydrogen is set to 1.0 × 10 5 kg/s. The mass flow of the ammonia and hydrogen fuels is distributed among the double ring-shaped burner as (1) supplying 100% NH 3 through the inner ring-shaped burner, i.e., m ̇ NH 3 = 9.99 × 10 6 kg/s and (2) 100% hydrogen (i.e., m ̇ H 2 = 1.0 × 10 8 kg/s) via the outer ring-shaped burner. Figure 10 illustrates the time evolution of the fluctuations of the acoustic properties (acoustic pressure and velocity) and thermodynamic properties such as temperature T ( t ) and heat of reaction Q ( t ). The mean heat of reaction Q ¯ is approximately 260 W. It can be seen that self-sustained intermittent pulsating oscillations are successfully produced. Furthermore, there are two dominant modes. One is at approximately 280 Hz, and the other mode is at 2.0–4.0 Hz. The higher frequency corresponds to the longitudinal acoustic mode-shape of the combustor.50–52 However, such extremely low-frequency pulsations at 2.0–4.0 Hz are most likely due to the pulsating fuel flow resulting from the force convection and shear layer separation through the double ring-shaped burner. It may also result from the thermal-diffusive instability.53 It excites the flame sheet to move toward and away from the double ring-shaped burner periodically.

FIG. 10.

Time evolution of (a) the envelope of the acoustic pressure fluctuations, (b) the envelope of the acoustic velocity fluctuations, (c) the envelope of the temperature fluctuations, (d) time evolution of the heat of reaction fluctuations Q ( t ), as m ̇ f = 1.0 × 10 5 kg/s and m ̇ H 2 = 2.0 × 10 8 kg/s.

FIG. 10.

Time evolution of (a) the envelope of the acoustic pressure fluctuations, (b) the envelope of the acoustic velocity fluctuations, (c) the envelope of the temperature fluctuations, (d) time evolution of the heat of reaction fluctuations Q ( t ), as m ̇ f = 1.0 × 10 5 kg/s and m ̇ H 2 = 2.0 × 10 8 kg/s.

Close modal

It is worth noting that further reduction of the mass flow rate of the blended hydrogen is shown to lead to an unsuccessful generation of such pulsating oscillations54 or unsuccessful ignition. Therefore, we confirm that the minimum hydrogen blended to enhance ammonia combustion in the presence of self-excited pulsating oscillations55–58 is 0.1% in terms of the total mass flow rate m ̇ f.

Unlike the conventional hydrocarbon fuels ( C x H y) or fossil fuels, ammonia and hydrogen are non-carbon fuels, which exhibit a larger mass diffusivity feature. If the ammonia and hydrogen were burnt with pure oxygen, then the products of the complete combustion are described below.
(3)
In practical applications, ammonia and hydrogen are burnt with the air, which consists of both oxygen and nitrogen. Then, it is expected that there are more complicated chemical reaction mechanism and more complex products. For example, the combustion product may compose of NO and NO 2.59,60
(4)

Figures 11 and 12 below compare the profiles of the reactants and combustion product species of NO , NO 2, and H 2 O in both axial (see Fig. 11) and radial (see Fig. 12) directions, as the total fuel flow rate m ̇ f is set to 1.0 × 10 5 kg/s, and 0.1% hydrogen is blended with ammonia. Figures 11(a)–11(c) show the variation of oxygen, hydrogen, and ammonia distribution at three different timings, i.e., t = 1 , t = 3, and t = 5 s. Along the axial direction, the hydrogen and ammonia are consumed along the axial direction. There is a sudden change near the heat exchangers, which could lead to the effect of “blocking.” Complete depletion of the hydrogen and ammonia is achieved after the heat exchanger zone. As for the oxygen is concerned, there is excessive oxygen flowing out of the combustor. It reveals that the present combustion system is operated at lean combustion condition. Figures 11(d) and 11(e) depicts the combustion product variations along the axial direction. It can be seen that there is an increase in NO 2 emission and a decrease in NO at the heat exchangers location. This means that the combustion is more complete. Finally, at different combustion times, i.e., t = 1, 3, and 5 s, the reactants and combustion products are different in term of the number of moles. This reveals that the periodic thermodynamic oscillations61 contribute to such emission difference. This finding is consistent with the radial observation as shown in Figs. 12(d)–12(f).

FIG. 11.

Variation of the reactants of O 2 (a), H 2 (b), and NH 3 (c) and combustion products of NO 2 (d) , NO (e), and H 2 O (f) along the axial direction at different times, as m ̇ f = 1.0 × 10 5 kg/s and ω ̇ H 2 = 0.1 %.

FIG. 11.

Variation of the reactants of O 2 (a), H 2 (b), and NH 3 (c) and combustion products of NO 2 (d) , NO (e), and H 2 O (f) along the axial direction at different times, as m ̇ f = 1.0 × 10 5 kg/s and ω ̇ H 2 = 0.1 %.

Close modal
FIG. 12.

Variation of the reactants of O 2 (a), H 2 (b), and NH 3 (c) and combustion products of NO 2 (d) , NO (e), and H 2 O (f) along the radial direction at different times, as m ̇ f = 1.0 × 10 5 kg/s and ω ̇ H 2 = 0.1 %. NM denotes the number of moles of each species.

FIG. 12.

Variation of the reactants of O 2 (a), H 2 (b), and NH 3 (c) and combustion products of NO 2 (d) , NO (e), and H 2 O (f) along the radial direction at different times, as m ̇ f = 1.0 × 10 5 kg/s and ω ̇ H 2 = 0.1 %. NM denotes the number of moles of each species.

Close modal

Figures 12(a)–12(c) show the variation of the reactants in the radial direction, as t is set to three different values. It can be seen that there are non-uniform profiles of the reactant distribution in the radial direction. The reactants are depleted more near the axis of the combustor than the regions near the combustor wall. Same non-uniform profiles of the combustion products are observed from Figs. 12(d)–12(f). This is most likely due to the periodic pressure and flow oscillations. Such pulsating fluctuations62–64 propagate along the combustor and the double ring-shaped fuel tubes. They in turn affect the local molar fractions of the mixture species such as H 2 , H 2 O , NH 3, NO, and NO 2. In general, burning ammonia and/or hydrogen does not produce any CO X but NO X. The combustion product of NO X is also undesirable. However, it is widely applied to use selective catalytic reduction (SCR) process to reduce NO X emission from burning ammonia.65 Unlike the conventional hydrocarbon fuels,66,67 less-flammable ammonia with a high ignition temperature is an alternative but promising carbon-free fuel. It may contribute to the significant reduction of the greenhouse gas of CO 2 emission from thermoacoustic combustors.68,69 In order to neutralize CO 2 emission, next-generation gas turbines70–72 are expected to burn NH 3 completely or the mixture of blending NH 3 with a small amount of H 2 or conventional hydrocarbon fuels.73,74 In practice, it is more interesting to know the minimum hydrogen being blended with the ammonia13,75–77 and the NO X emission. This motivated the present work.

In this work, 2D numerical unsteady RANS (Reynolds-averaged Navier–Stokes) simulations are conducted on an acoustically open–open combustor with double ring-shaped fuel injectors implemented. Here, premixed ammonia–hydrogen-burnt flames are confined near the combustor inlet. The turbulence model of k–ϵ and combustion model of eddy dissipation concept (EDC) are applied. The model is fist validated by using the theoretical and experimental data available in the literature. It is then used to study the effects of (1) the mass flow rate m ̇ f of the mixture consisting of NH 3 and H 2; (2) the mass fraction w ̇ H 2 of H 2 blended with NH 3 with respect to the total fuel mixture of NH 3 and H 2. It is found that nonlinear self-excited pulsating combustion oscillations could be generated, depending on the mass flow rate m ̇ f of the fuel mixture. Furthermore, the amplitude and frequency of the dominant mode is found to strongly depend on the mass fraction w ̇ H 2 of the hydrogen. In addition, as w ̇ H 2 is reduced from 88% to 25%, intermittent combustion oscillations are observed. A pulsating mode at extremely low frequency of approximately 2.0 4.0 Hz is found to occur. Such periodic mode is shown to related to the thermo-diffusive instability of the NH 3. When the total mass flow rate m ̇ f of NH 3 and H 2 is small, for example, 1.0 × 10 6 kg/s, no apparent limit cycle oscillations are observed. This is most likely due to fact that the heat-to-sound energy conversion is not as intensified as the acoustic losses. In addition, the minimum hydrogen being blended with ammonia is determined by varying the mass fraction of the hydrogen with respect to the ammonia while maintaining the total fuel flow rate m ̇ f to be constant. It is shown to be 0.1%, i.e., w ̇ H 2 min = 0.1 %. Note that the minimum hydrogen is determined, as the pulsating oscillations are successfully generated.

Combustion products/emissions from the premixed ammonia–hydrogen flames are examined in the presence of such large-amplitude pulsating oscillations. It is observed that NO emission is reduced, and the product of H 2 O is increased along the axial direction. However, NO 2 emission is almost unchanged. Along the radial direction, the reactants such as NH 3 and H 2 are consumed/depleted more and faster near the centerline of the combustor in comparison with the regions near the combustor inner wall. The combustion products are not uniformly distributed along the radial direction. The presence of the pulsating oscillations leads to the combustion products to be varied “periodically” along both the axial and radial directions.

Finally, with the generation mechanism of self-excited pulsating combustion oscillations better understood, temperature-variable heat exchangers as an “acoustic absorber” are applied and numerically examined for controlling and enhancing the combustion stability behaviors. Such control approach is proven to be “slowly responding” and somehow dampen the intermittent pulsating combustion oscillations (at both 2.0 4.0 and 250 Hz), as its surface temperature is set to TH = 1300 K. However, it is found to be insufficient to completely attenuate such pulsating oscillations. Furthermore, increasing the surface temperature by another 200–300 K and varying its axial location may be applicable. This may lead to a further attenuation of the pulsating oscillations. However, this may lead to excessive additional heat input (thermal energy input) to the combustion system to increase its average temperature. This may prevent its application in practical engines.

This work is supported by University of Canterbury with a Grant No. 452DISDZ. This financial support is gratefully acknowledged.

The authors have no conflicts to disclose.

Yiheng Guan: Investigation (equal); Software (equal); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal). Dan Zhao: Conceptualization (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Writing – original draft (equal); Writing – review & editing (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
A.
Valera-Medina
,
H.
Xiao
,
M.
Owen-Jones
,
W. I. F.
David
, and
P. J.
Bowen
, “
Ammonia for power
,”
Prog. Energy Combust. Sci.
69
(
1
),
63
102
(
2018
).
2.
A.
Valera-Medina
,
F.
Amer-Hatem
,
A. K.
Azad
,
I. C.
Dedoussi
,
M.
de Joannon
,
R. X.
Fernandes
,
P.
Glarborg
,
H.
Hashemi
,
X.
He
,
S.
Mashruk
,
J.
McGowan
,
C.
Mounaim-Rouselle
,
A.
Ortiz-Prado
,
A.
Ortiz-Valera
,
I.
Rossetti
,
B.
Shu
,
M.
Yehia
,
H.
Xiao
, and
M.
Costa
, “
Review on ammonia as a potential fuel: From synthesis to economics
,”
Energy Fuels
35
(
1
),
6964
7029
(
2021
).
3.
H.
Kobayashi
,
A.
Hayakawa
,
K. D. K. A.
Somarathne
, and
E. C.
Okafor
, “
Science and technology of ammonia combustion
,”
Proc. Combust. Inst.
37
(
1
),
109
133
(
2019
).
4.
M.-C.
Chiong
,
C. T.
Chong
,
J.-H.
Ng
,
S.
Mashruk
,
W. W. F.
Chong
,
N. A.
Samiran
,
G. R.
Mong
, and
A.
Valera-Medina
, “
Advancements of combustion technologies in the ammonia-fuelled engines
,”
Energy Convers. Manage.
244
(
1
),
114460
(
2021
).
5.
C.
Duynslaegher
,
H.
Jenmart
, and
J.
Vandooren
, “
Ammonia combustion at elevated pressure and temperature conditions
,”
Fuel
89
(
1
),
3540
3545
(
2010
).
6.
C.
Duynslaegher
,
F.
Contino
,
J.
Vandooren
, and
H.
Jeanmart
, “
Modeling of ammonia combustion at low pressure
,”
Combust. Flame
159
(
9
),
2799
2805
(
2012
).
7.
W. S.
Chai
,
Y.
Bao
,
P.
Jin
,
G.
Tang
, and
L.
Zhou
, “
A review on ammonia, ammonia-hydrogen and ammonia-methane fuels
,”
Renewable Sustainable Energy Rev.
147
(
1
),
111254
(
2021
).
8.
F. J.
Verkamp
,
M. C.
Hardin
, and
J. R.
Williams
, “
Ammonia combustion properties and performance in gas-turbine burners
,”
Symp. (Int.) Combust.
11
(
1
),
985
992
(
1967
).
9.
A. M.
Elbaz
,
S.
Wang
,
T. F.
Guiberti
, and
W. L.
Roberts
, “
Review on the recent advances on ammonia combustion from the fundamentals to the applications
,”
Fuel Commun.
10
(
1
),
100053
(
2022
).
10.
S.
Ni
,
D.
Zhao
,
W. W.
Wu
, and
Y.
Guan
, “
NOx emission reduction reaction of ammonia-hydrogen with self-sustained pulsating oscillations
,”
Therm. Sci. Eng. Prog.
19
(
1
),
100615
(
2020
).
11.
A.
Hayakawa
,
Y.
Arakawa
,
R.
Mimoto
,
K. D. K. A.
Somarathne
,
T.
Kudo
, and
H.
Kobayashi
, “
Experimental investigation of stabilization and emission characteristics of ammonia/air premixed flames in a swirl combustor
,”
Int. J. Hydrogen Energy
42
(
19
),
14010
14018
(
2017
).
12.
E. C.
Okafor
,
M.
Tsukamoto
,
A.
Hayakawa
,
K. D. K. A.
Somarathne
,
T.
Kudo
,
T.
Tsujimura
, and
H.
Kobayashi
, “
Influence of wall heat loss on the emission characteristics of premixed ammonia-air swirling flames interacting with the combustor wall
,”
Proc. Combust. Inst.
38
(
4
),
5139
5146
(
2021
).
13.
A.
Katoch
,
T. F.
Guiberti
,
D. V.
de Campos
, and
D. A.
Lacoste
, “
Dual-fuel, dual-swirl burner for the mitigation of thermoacoustic instabilities in turbulent ammonia-hydrogen flames
,”
Combust. Flame
246
(
1
),
112392
(
2022
).
14.
B.
Mei
,
J.
Zhang
,
X.
Shi
,
Z.
Xi
, and
Y.
Li
, “
Enhancement of ammonia combustion with partial fuel cracking strategy: Laminar flame propagation and kinetic modeling investigation of NH3/H2/N2/air mixtures up to 10 atm
,”
Combust. Flame
231
(
1
),
111472
(
2021
).
15.
K.
Bioche
,
L.
Bricteux
,
A.
Bertolino
,
A.
Parente
, and
J.
Blondeau
, “
Large eddy simulation of rich ammonia/hydrogen/air combustion in a gas turbine burner
,”
Int. J. Hydrogen Energy
46
(
79
),
39548
39562
(
2021
).
16.
T.
Cai
and
D.
Zhao
, “
Enhancing and assessing ammonia-air combustion performance by blending with dimethyl ether
,”
Renewable Sustainable Energy Rev.
156
(
1
),
112003
(
2022
).
17.
H.
Xiao
,
A.
Valera-Medina
, and
P. J.
Bowen
, “
Modeling combustion of ammonia/hydrogen fuel blends under gas turbine conditions
,”
Energy Fuels
31
(
8
),
8631
8642
(
2017
).
18.
H.
Lee
and
M. J.
Lee
, “
Recent advances in ammonia combustion technology in thermal power generation system for carbon emission reduction
,”
Energies
14
(
18
),
5604
5629
(
2021
).
19.
D.
Wang
,
C.
Ji
,
S.
Wang
,
Z.
Wang
,
J.
Yang
, and
Q.
Zhao
, “
Numerical study on the premixed oxygen-enriched ammonia combustion
,”
Energy Fuels
34
(
12
),
16903
16917
(
2020
).
20.
A.
Mohammadpour
,
K.
Mazaheri
,
A.
Alipoor
, and
M.
Mohammadpour
, “
Effect of oxidizer dilution on thermophotovoltaic performance of premixed ammonia/hydrogen/air combustion in radial micro-combustor
,”
Fuel
331
(
2
),
125874
(
2023
).
21.
S.
Oh
,
C.
Park
,
J.
Oh
,
S.
Kim
,
Y.
Kim
,
Y.
Choi
, and
C.
Kim
, “
Combustion, emissions, and performance of natural gas ammonia dual fuel spark ignited engine at full load condition
,”
Energy
258
(
1
),
124837
(
2022
).
22.
J.
Frost
,
A.
Tall
,
A. M.
Sheriff
,
A.
Schönborn
, and
P.
Hellier
, “
An experimental and modelling study of dual fuel aqueous ammonia and diesel combustion in a single cylinder compression ignition engine
,”
Int. J. Hydrogen Energy
46
(
71
),
35495
35510
(
2021
).
23.
J.
Li
,
H.
Huang
,
N.
Kobayashi
,
Z.
He
, and
Y.
Nagai
, “
Study on using hydrogen and ammonia as fuels: Combustion characteristics and NOx formation
,”
Int. J. Energy Res.
38
(
9
),
1214
1223
(
2014
).
24.
J. G.
Lee
and
D. A.
Santavicca
, “
Experimental diagnostics for the study of combustion instabilities in lean premixed combustors
,”
J. Propul. Power
19
(
5
),
737
750
(
2003
).
25.
K. C.
Schadow
and
E.
Gutmark
, “
Combustion instability related to vortex shedding in dump combustors and their passive control
,”
Prog. Energy Combust. Sci.
18
(
2
),
117
132
(
1992
).
26.
W.
Li
,
D.
Zhao
,
X.
Chen
,
Y.
Sun
,
S.
Ni
,
D.
Guan
, and
B.
Wang
, “
Numerical investigations on solid-fueled ramjet inlet thermodynamic properties effects on generating self-sustained combustion instability
,”
Aerosp. Sci. Technol.
119
,
107097
(
2021
).
27.
Y. C.
Yu
,
J. C.
Sisco
,
S.
Rosen
,
A.
Madhav
, and
W. E.
Anderson
, “
Spontaneous longitudinal combustion instability in a continuously-variable resonance combustor
,”
J. Propul. Power
28
(
5
),
876
887
(
2012
).
28.
K. K.
Venkataraman
,
L. H.
Preston
,
D. W.
Simons
,
B. J.
Lee
,
J. G.
Lee
, and
D. A.
Santavicca
, “
Mechanism of combustion instability in a lean premixed dump combustor
,”
J. Propul. Power
15
(
6
),
909
918
(
1999
).
29.
V.
Yang
and
F. E. C.
Culick
, “
Analysis of low frequency combustion instabilities in a laboratory ramjet combustor
,”
Combust. Sci. Technol.
45
(
1–2
),
1
25
(
1986
).
30.
K. H.
Yu
,
A.
Trouvé
, and
J. W.
Daily
, “
Low-frequency pressure oscillations in a model ramjet combustor
,”
J. Fluid Mech.
232
(
1
),
47
72
(
1991
).
31.
F.
Qin
,
G.
He
,
P.
Liu
, and
J.
Li
, in
45th AIAA/ASME Joint Propulsion Conference and Exhibit
,
2009
.
32.
W.
Li
,
D.
Zhao
,
X.
Chen
,
L.
Zhu
, and
S.
Ni
, “
Numerical investigation of inlet thermodynamic conditions on solid fuel ramjet performances
,”
Int. J. Aerosp. Eng.
2021
,
8868288
.
33.
J.
Peng
,
Z.
Cao
,
X.
Yu
,
S.
Yang
,
Y.
Yu
,
H.
Ren
,
Y.
Ma
,
S.
Zhang
,
S.
Chen
, and
Y.
Zhao
, “
Analysis of combustion instability of hydrogen fueled scramjet combustor on high-speed OH-PLIF measurements and dynamic mode decomposition
,”
Int. J. Hydrogen Energy
45
(
23
),
13108
13118
(
2020
).
34.
Y.
Guan
,
D.
Zhao
,
S.
Becker
,
J.
Xu
,
M.
Shahsavari
, and
J.
Schluter
, “
Entropy generation and CO2 emission in presence of pulsating oscillations in a bifurcating thermoacoustic combustor with a Helmholtz resonator at off-design conditions
,”
Aerosp. Sci. Technol.
136
(
1
),
108204
(
2023
).
35.
S.
Nakaya
,
H.
Yamana
, and
M.
Tsue
, “
Experimental investigation of ethylene/air combustion instability in a model scramjet combustor using image-based methods
,”
Proc. Combust. Inst.
38
(
3
),
3869
3880
(
2021
).
36.
X.
Fu
,
Z.
Guo
, and
F.
Yang
, “
Combustion instability mode transition in a pilot bluff-body stabilized combustor
,”
J. Propul. Power
32
(
1
),
1
12
(
2016
).
37.
Y.
Cheng
,
T.
Jin
,
K.
Luo
,
Z.
Li
,
H.
Wang
, and
J.
Fan
, “
Large eddy simulations of spray combustion instability in an aero-engine combustor at elevated temperature and pressure
,”
Aerosp. Sci. Technol.
108
(
1
),
106329
(
2021
).
38.
G.
Zhao
,
M.
Sun
,
J.
Wu
,
X.
Cui
, and
H.
Wang
, “
Investigation of flame flashback phenomenon in a supersonic crossflow with ethylene injection upstream of cavity flameholder
,”
Aerosp. Sci. Technol.
87
(
1
),
190
206
(
2019
).
39.
S.
Ducruix
,
T.
Schuller
,
D.
Durox
, and
S.
Candel
, “
Combustion dynamics and instabilities: Elementary coupling and driving mechanisms
,”
J. Propul. Power
19
(
5
),
722
731
(
2003
).
40.
J.
O'Connor
,
V.
Acharya
, and
T.
Lieuwen
, “
Transverse combustion instabilities: Acoustic, fluid mechanic, and flame processes
,”
Prog. Energy Combust. Sci.
49
(
1
),
1
39
(
2015
).
41.
D.
Zhao
and
X. Y.
Li
, “
A review of acoustic dampers applied to combustion chambers in aerospace industry
,”
Prog. Aerosp. Sci.
74
(
1
),
114
130
(
2015
).
42.
D.
Zhao
and
X. Y.
Li
, “
A review of active control approaches in stabilizing combustion systems in aerospace industry
,”
Prog. Aerosp. Sci.
97
(
1
),
35
60
(
2018
).
43.
W.
Li
and
D.
Zhao
, “
Proper orthogonal and dynamic mode decomposition analyses of nonlinear combustion instabilities in a solid-fuel ramjet combustor
,”
Ther. Sci. Eng. Prog.
27
(
1
),
101147
(
2022
).
44.
O.
Musa
,
X.
Chen
, and
C.
Zhou
, “
Experimental and numerical investigation on the ignition and combustion stability in solid fuel ramjet with swirling flow
,”
Acta Astronaut.
137
(
1
),
157
167
(
2017
).
45.
D.
Zhao
and
Y.
Guan
, “
Characterizing modal exponential growth behaviors of self-excited transverse and longitudinal thermoacoustic instabilities
,”
Phys. Fluids
34
(
1
),
024109
(
2022
).
46.
R. L.
Raun
,
M. W.
Beckstead
,
J. C.
Finlinson
, and
K. P.
Brooks
, “
A review of Rijke tubes, Rijke burners and related devices
,”
Prog. Energy Combust. Sci.
19
(
4
),
313
364
(
1993
).
47.
J. W. S. B.
Rayleigh
,
The Theory of Sound
(
Dover
,
New York
,
1945
), Vol.
II
.
48.
D.
Zhao
and
Z. H.
Chow
, “
Thermoacoustic instability of a laminar premixed flame in Rijke tube with a hydrodynamic region
,”
J. Sound Vib.
332
(
14
),
3419
3437
(
2013
).
49.
D.
Zhao
, “
Transient growth of flow disturbances in triggering a Rijke tube combustion instability
,”
Combust. Flame
159
(
6
),
2126
2137
(
2012
).
50.
L.
Luo
,
T.
Li
,
J.
Deng
,
R.
Zhao
, and
J.
Wang
, “
Experimental investigation on self-excited thermoacoustic instability in a Rijke tube
,”
Appl. Sci.
12
(
15
),
8046
(
2022
).
51.
T.
Liu
,
J.
Li
,
S.
Zhu
, and
L.
Yang
, “
Determination of the heat conduction transfer function within the thermoacoustic instability limit cycle in a Rijke tube
,”
Appl. Therm. Eng.
206
(
1
),
118084
(
2022
).
52.
Z.
Cao
,
Y.
Xiao
,
W.
Song
, and
W.
Li
, “
Experimental study on thermoacoustic instability and emission characteristics of premixed swirl flame in a longitudinal combustor
,”
Int. J. Therm. Sci.
172
,
107354
(
2022
).
53.
D.
Zhao
and
E.
Gutmark
, “
Mitigating self-excited flame pulsating and thermoacoustic oscillations using perforated liners
,”
Sci. Bull.
64
(
13
),
941
952
(
2019
).
54.
E.
Flores-Montoya
,
V.
Muntean
,
M.
Sánchez-Sanz
, and
D.
Martínez-Ruiz
, “
Non-adiabatic modulation of premixed-flame thermoacoustic frequencies in slender tubes
,”
J. Fluid Mech.
933
(
1
),
A50
(
2022
).
55.
L.
Miniero
,
G. A.
Mensah
,
C.
Bourquard
, and
N.
Noiray
, “
Failure of thermoacoustic instability control due to periodic hot gas ingestion in Helmholtz dampers
,”
J. Sound Vib.
548
(
1
),
117544
(
2023
).
56.
Q.
Zeng
,
C.
Hu
, and
J.
Zhu
, “
Prediction of acoustic pressure of thermoacoustic combustion instability based on Elman neural network
,”
J. Low Freq. Noise Vib. Act. Control
(published online 2023).
57.
X.
Liu
,
H.
Zhou
,
Y.
Lai
, and
Y.
Zhang
, “
Equivalence ratio independence and dependence ranges of system responses for a nonlinear thermoacoustic oscillation in a Rijke tube
,”
J. Sound Vib.
547
(
1
),
117545
(
2023
).
58.
Z. X.
Chen
,
N.
Swaminathan
,
M.
Mazur
, and
N. A.
Worth
, “
Numerical investigation of azimuthal thermoacoustic instability in a gas turbine model combustor
,”
Fuel
339
(
1
),
127405
(
2023
).
59.
B.
Hu
and
H.
Zhou
, “
Experimental investigation of combustion instability and NOX emissions control by adjustable swirl jet in lean premixed flames
,”
J. Energy Inst
106
(
1
),
101149
(
2023
).
60.
P. W.
Agostinelli
,
D.
Laera
,
I.
Chterev
, and
I.
Boxx
, “
Large eddy simulations of mean pressure and H2 addition effects on the stabilization and dynamics of a partially-premixed swirled-stabilized methane flame
,”
Combust. Flame
249
(
1
),
112592
(
2023
).
61.
X.
Wang
and
X.
Han
, “
Transitions of thermoacoustic modes and flame dynamics in a centrally-staged swirl combustor
,”
Energy
263
,
125813
(
2023
).
62.
F.
Li
,
L.
Xu
, and
L.
Yang
, “
Experimental characterization of self-excited combustion pulsation in a thermoacoustic combustor
,”
Fuel
334
(
2
),
126423
(
2023
).
63.
Y.
Sun
,
D.
Zhao
, and
C.
Ji
, “
Large eddy simulations of self-excited thermoacoustic instability in a premixed swirling combustor with an outlet nozzle
,”
Phys. Fluids
34
(
1
),
044112
(
2022
).
64.
M.
Jafari
and
A. J.
Jaworski
, “
On the nonlinear behaviour of oscillatory flow in a high pressure amplitude standing-wave thermoacoustic heat engine
,”
Int. J. Heat Mass Transfer
201
(
2
),
123595
(
2023
).
65.
T.
Cai
,
D.
Zhao
, and
E.
Gutmark
, “
Overview of fundamental kinetic mechanisms and emission mitigation in ammonia combustion
,”
Chem. Eng. J.
458
(
1
),
141391
(
2023
).
66.
R. S.
Bhavi
,
I.
Pavithran
,
A.
Roy
, and
R. I.
Sujith
, “
Abrupt transitions in turbulent thermoacoustic systems
,”
J. Sound Vib.
547
(
1
),
117478
(
2023
).
67.
N. A.
Baraiya
,
V.
Ramanan
,
B.
Nagarajan
,
C. S.
Vegad
, and
S. R.
Chakravarthy
, “
Dynamic mode decomposition of syngas (H2/CO) flame during transition to high-frequency instability in turbulent combustor
,”
Energy
263
(
1
),
125998
(
2023
).
68.
Z.
Qin
,
X.
Wang
,
X.
Han
, and
Y.
Lin
, “
Pre-trained combustion model and transfer learning in thermoacoustic instability
,”
Phys. Fluids
35
(
1
),
037117
(
2023
).
69.
D.
Zhao
, “
System identification and impulse response of a methane-fuelled swirling thermoacoustic combustor
,”
Therm. Sci. Eng. Prog.
37
(
1
),
101575
(
2023
).
70.
K.
Guo
,
L.
Ni
,
Y.
Ren
,
W.
Lin
,
Y.
Tong
, and
W.
Nie
, “
Experimental investigation on effects of injection distribution on combustion instability in a model rocket combustor
,”
Phys. Fluids
34
(
1
),
011702
(
2022
).
71.
M.
Balasubramaniyan
,
A.
Kushwaha
,
Y.
Guan
,
J.
Feng
,
P.
Liu
,
V.
Gupta
, and
L. K. B.
Li
, “
Global hydrodynamic instability and blowoff dynamics of a bluff-body stabilized lean-premixed flame
,”
Phys. Fluids
33
(
1
),
034103
(
2021
).
72.
A. I.
Krikunova
,
K. Y.
Arefyev
,
A. S.
Saveliev
,
A. D.
Kossov
, and
G. A.
Cheshko
, “
Inverted conical methane/air flame shape transformation under acoustic excitation
,”
Phys. Fluids
33
(
1
),
053610
(
2021
).
73.
Y.
Fang
,
Y.
Yang
,
K.
Hu
,
G.
Wang
,
J.
Li
, and
Y.
Zheng
, “
Experimental study on self-excited thermoacoustic instabilities and intermittent switching of azimuthal and longitudinal modes in an annular combustor
,”
Phys. Fluids
33
(
1
),
084104
(
2021
).
74.
X.
Wang
,
X.
Han
,
H.
Song
,
D.
Yang
, and
C.-J.
Sung
, “
Multi-bifurcation behaviors of stability regimes in a centrally staged swirl burner
,”
Phys. Fluids
33
(
1
),
095121
(
2021
).
75.
J.
Li
,
R.
Zhang
,
J.
Pan
,
H.
Wei
,
G.
Shu
, and
L.
Chen
, “
Ammonia and hydrogen blending effects on combustion stabilities in optical SI engines
,”
Energy Convers. Manage.
280
(
1
),
116827
(
2023
).
76.
T.
Cai
and
D.
Zhao
, “
Removal and mechanism analysis of NOx emissions in carbon-free ammonia combustion systems with a secondary fuel injection
,”
Fuel
344
(
1
),
128088
(
2023
).
77.
H.
Zhao
and
D.
Zhao
, “
No emission and enhanced thermal performances studies on counter-flow double-channel hydrogen/ammonia-fuelled microcombustors with oval-shaped internal threads
,”
Fuel
341
(
1
),
127665
(
2023
).
Published open access through an agreement with University of Canterbury