# CHARACTERIZATION OF LAMINAR ULTRA-LEAN H₂-AIR ANNULAR FLAMES USING SPECTRAL PROPER ORTHOGONAL DECOMPOSITION

N. L. Junqueira¹,², L. F. Figueira da Silva² and V. Morgenthaler¹

nicole.lopes-monteiro-de-barros-junqueira@ensma.fr

¹ Ansys, France

² Institut Pprime CNRS - ISAE-ENSMA and Université de Poitiers, France

## Abstract

Lean hydrogen-air flames have the benefit of consuming less fuel and reaching lower temperatures than stoichiometric flames, but they exhibit intrinsic instabilities that affect flame dynamics and propagation. The goal of this study is to analyze the spatio-frequential structure of ultra-lean, burner stabilized, hydrogen-air flames in two distinct situations. The first corresponds to the presence of intrinsic instabilities only, whereas, in the second, longitudinal pulsations related to the burner acoustic arise. The analysis uses measured time resolved schlieren data and the spectral proper orthogonal decomposition technique (SPOD). The SPOD is a reduction technique for isolating patterns and structures in complex flows, commonly applied in aeroacoustics and jet studies to represent flow structures. The results for the intrinsically unstable flame show a monotonic energy decay with frequency, with an exponential decay of around 19 modes per energy decade. The spatial distribution remains similar across frequencies, depicting intercalated wave packets across the burner, resembling the intrinsic instability's transversal movement, as observed in schlieren images. For the longitudinal pulsation case, the energy spectrum reveals multiple peaks, with the first mode at the harmonic frequencies exhibiting significantly higher energy. From the second mode on, an exponential decrease of about 18 modes per energy decade is observed. The spatial distribution varies with frequency, where the first mode at harmonics displays longitudinally intercalated wave packets, whereas the second mode appears to capture a superposition from both instabilities. The SPOD analysis highlights the distinct interactions between intrinsic instabilities and longitudinal pulsations due to burner acoustic modes across different frequencies and spatial modes.

## Introduction

Hydrogen, a source of potentially sustainable alternative fuel, is being increasingly studied [1]. The behavior of the stable premixed laminar H₂-air flames when the Lewis number controlling the combustion is close to unity, i.e., mixtures close to stoichiometry, is well known [2]. However, burner stabilized, ultra-lean H₂-air mixtures are still to be thoroughly examined. Those mixtures are interesting for a few aspects: they consume less fuel and generate lower NOx emissions than stoichiometric ones, and temperature of the burned gases is smaller, i.e., they impose less thermal stresses on the combustion chambers.

However, lean premixed hydrogen/air flames exhibit intrinsic flame instabilities that influence flame dynamics and propagation [3, 4]. One of these instabilities, present in all premixed flames, is the Darrieus-Landau instability, that arises from density variations across the flame sheet [4]. Another intrinsic instability is the thermodiffusive one, controlled by the hydrogen's low Lewis number [5, 6], that is associated to a substantial difference between the diffusivities of temperature and fuel, thus leading to differential diffusion effects within the flame front. These instabilities play a crucial role in shaping the behavior and characteristics of lean premixed hydrogen-air flames [7].

Several recent works, based on detailed numerical simulations, have analyzed the of flame front structure, flame speed enhancements, and the corresponding influence of parameters such as temperature, pressure, and equivalence ratio [6–11]. These studies also showed that thermodiffusive instabilities not only change the flame speed but also introduce intrinsic length scales
in the flame [12, 13], which could impact the design and optimization of hydrogen combustion practical systems.

Recent application-oriented experimental studies also explore lean hydrogen-air flames, focusing on flame structure, thermo-acoustic instabilities, and propagation [14–16], thus contributing to the development of hydrogen-based combustion systems. In particular, the combustion instabilities have been studied to reveal the oscillatory modes and frequencies that influence flame behavior and efficiency in lean hydrogen-air combustion [16]. However, the dynamics of burner stabilized, unstable ultra-lean laminar hydrogen-air flames has not yet been thoroughly experimentally examined. While existing experimental studies have advanced understanding of these instabilities, further insights can be gained by examining the dominant energy modes and spatial-temporal structures within these flames, which is the goal of this study.

To isolate patterns and structures within complex, turbulent flows, in the field of aero-acoustics and jets, the spectral proper orthogonal decomposition (SPOD) has been applied to a wide range of purposes [17–20]. This technique empirically obtains hidden structures of the data and it is particularly effective for analyzing data where coherent, energy-dominant structures evolve in both time and space [21]. It finds an optimal orthogonal basis, with modes that vary in both space and time [22]. These modes are orthogonal under a space-time inner product, meaning that they can be ranked in terms of energy. The SPOD thus combines temporal and spatial information, delivering a set of orthogonal modes ranked by energy at discrete frequencies. Accordingly, this study aims to characterize the spectral proper orthogonal decomposition structures of the instantaneous schlieren of the unstable lean premixed hydrogen-air laminar flames stabilized on an annular burner for the equivalence ratios between the lean blowoff and the flashback limits. This will enable us to unveil the spectral and modal behavior of these ultra-lean flames, in cases where thermodiffusive instability arises only, and, where flame pulsations associated to the burner acoustics modes are also observed.

## Methodology

In this section it is presented the developed burner, the experimental method and the spectral proper orthogonal decomposition.

### *Experimental methodology*

The burner developed for this study, shown in Fig. 1 is based on the “ANÔ” burner, that was used to perform a stochastic characterization of unstable lean hydrogen-air flames [5]. The new slit burner was designed to provide a smaller pressure drop throughout the inlet system, and to enable the stabilization of lower Reynolds number flames. The central part of the new annular-slit burner is open, thus allowing ambient air to flow freely. As it is shown in Fig. 1a, the mixture of hydrogen and air occurs in a ellipsoidal-shaped settling chamber connected by ten 26 cm length tubes of 6 mm inner diameter to the bottom of the burner. The burner cross section is illustrated in Fig. 1b: a porous plate, which homogenizes the mixture and stabilizes the flow is located upstream to the burner slit at a distance of 80 mm. The flame is located outside the burner, attached to two discs that form a 2 mm annular slit, with a thickness of 10 mm. The inner radius of the slit is 39 mm and the outer radius is 41 mm. The resulting perimeter to slit ratio is around 120, which allows us to observe the nearly two-dimensional flame instabilities. The material of the burner is stainless steel and the outer disk surface temperature is monitored by a thermocouple.

To explore the ultra-lean laminar hydrogen-air flames, three mass flow controllers (Brooks models SLA5850 and SLA5851) were used, with flow meter ranges of 20, 50 and 100 nlpm. In this burner, the blow-off limit is around the equivalence ratio of $\phi$ =0.26, whereas the first point at which the flame is anchored on both sides of the slit is $\phi$ =0.28, and the flashback limit is around the equivalence ratio of 0.40. The obtained Reynolds number, defined using the fuel-air mixture flow rate and the burner hydraulic diameter, ranges from 480 to 860, and the equivalence ratio ranges from 0.28 to 0.40.
<figure><img src="image_1.png" /><figcaption>**Figure 1:** Annular slit burner scheme.</figcaption></figure>

The experimental technique that is the focus of this work is the dual-lens schlieren imaging [9], that allow us to observe the temperature gradient along the flame. In this setup, a LED lamp, the light source, is located 58 mm from the burner center. The light path on this schlieren setup is: the light passes thought the flame, it is focused by a 100 mm diameter, 750 mm focal length plano-convex lens to a knife edge of 45° that blocks the non-deviated light. Then, the deviated-light is focused by a bi-convex lens of 95 mm diameter and 150 mm focal distance to the camera, Photron FASTCAM Mini AX200 CMOS detector. The camera operated with two temporal resolutions, depending on the considered flame regime: 2,000 and 13,500 frames per second, with a shutter speed of 4 µs, and a spatial resolution of 10 px/mm. For each flow condition a set of 10,000 images was obtained. The raw images, sized 896 × 560 px in the (x, y) directions, were processed using a pixel-wise adaptive low-pass filter with a 5 × 5 pixel neighborhood to estimate local image statistics, thus reducing Gaussian white noise while preserving image features across each frame of the sequence. Then, a set of instantaneous frames is used to characterize the flame through the spectral proper orthogonal decomposition.

Two successive sets of experiments were carried out: the first aimed to determine the flash-back and blow-off limits, and the second one to characterize the flame regimes as a function of the equivalence ratio and Reynolds number.

### Spectral proper orthogonal decomposition

Figure 2 shows a scheme of the SPOD algorithm [22]. It begins by dividing the temporal dataset, $Sc(t)$), into overlapping time-domain blocks ($Sc^{(B_i)}$). For each of these blocks, the Direct Fourier Transform (DFT) is applied to convert the data to the frequency domain ($\hat{Sc}^{(B_i)}$). After obtaining the frequency representation of each block, the data is reorganized by frequency. For each frequency, Proper Orthogonal Decomposition (POD) is then applied. This procedure extracts the dominant coherent structures for each frequency by identifying the modes that contribute the most energy, $\hat{Sc}_{fk}$. Such an algorithm has already been presented and used for different applications, for additional details refer to [21–24].

From the SPOD, the energy spectrum that gives the decay of energy as a function of frequency for each mode is obtained. Here, the term “energy” refers to the quantity of information that a mode captures from the data, indicating its importance in the overall representation of the data set. In other words, it is the magnitude of a mode’s contribution to describing the structure or variability of the data. Overall, in a POD decomposition, the first few modes have a higher energy, and as the number of modes increases, this energy decreases. As a result, it is possible to evaluate whether a low-rank approximation to the problem is feasible, since the most rele-
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /><figcaption>**Figure 2:** Spectral proper orthogonal decomposition (SPOD) scheme, adapted from [24].</figcaption></figure>

vant information in the system (in this sense) has a higher energy. This characteristic comes from the singular value decomposition used in the POD step, in which the eigenvalues are organized in descending order so that the information provided by the first modes retains the main characteristics of the system analyzed. However, because SPOD requires the convergence of second-order space-time statistics, large datasets and long time series are necessary for accurate results [22].

In this work, the SPOD algorithm is applied to different regions of the instantaneous images, which are the complete flame, or, the center of the burner. Those regions of interest have an area of 896×210 px (ROI 1) and 200×210 px (ROI 2), respectively. Those regions are shown overlaid to the schlieren image in Fig 1b, where the blue rectangle shows the ROI 1, and, the orange rectangle is the ROI 2. In both cases, the vertical height (y-direction in Fig. 1b) of the image is the same and spans from the burner surface to around 10 mm above the tip of the flame, for the highest flame height obtained in the set of experiments. Each processed temporal dataset, that represents a Reynolds number and equivalence ratio (Re, φ) pair, is formed by 8,200 snapshots of the instantaneous schlieren images for the ROI 2 and 6,000 snapshots for the ROI 1. Regarding the SPOD parameters, the default spectral estimation parameters are used: the Hamming temporal window; 50% of the length of the temporal window as the number of segments overlapped; uniform spatial inner product weight. Therefore, for each (Re, φ) point at the ROI 2, there are a total of 31 temporal blocks, where each block has a total of 512 snapshots, and 256 of which are overlapping for the SPOD applied to the central region of the burner. For the full burner region, ROI 1, there are 22 temporal blocks with 528 snapshots, where 264 of those are overlapping. Note that the number of temporal blocks limits that of the SPOD modes.

## Results and Discussion

The first part of this section presents the particularities of the newly developed burner regarding the flame regimes and its flame structures, that were using the classical schlieren imaging technique. This is followed by the spectral proper orthogonal decomposition analysis, where is presented and discussed the energy and the wave packet of the different flame structures found.

### Flame regimes and structures

In this burner two characteristic flame structures are observed: the first is when only the intrinsic instabilities of the lean flames arise, i.e., the thermo-diffusive and the Darieus-Landau. The second flame structure is when there is a longitudinal pulsation caused by the coupling of the natural instability of the flame heat release rate with the acoustic modes of the burner. For the sake of simplicity, here, we refer to these cases as the intrinsic instability and longitudinal
pulsation cases, respectively. In Fig. 3, a regime diagram is showing, as a function of the Reynolds number and the equivalence ratio, the instabilities boundaries, as well as, the blow-off and flashback limits. It is possible to notice that the region where the longitudinal pulsations occurs lies between the equivalence ratios of 0.30 and 0.35. Moreover, it was observed that, in this $\phi$ range, the flame emits audible sound for the lower Reynolds numbers, even if the longitudinal pulsation is noticed in the schlieren for larger Re. Since there seems to be an acoustic feedback loop, acoustic measurements are not the subject of current work. Therefore, only the analysis of the different flame structures obtained is presented here.

<figure><img src="image_4.png" /><figcaption>**Figure 3:** Flame regimes diagram. The lines red and blue delimit the flashback and blow-off limits of the burner. The green region shows the flame regimes where it is observed the coupling of two instabilities: intrinsic and longitudinal pulsation. Symbols: □- the intrinsically unstable flame; ○- the longitudinal pulsation case.</figcaption></figure>

Figure 4 shows two examples of representative schlieren images, one for the intrinsic instability and another when there is a longitudinal pulsation. A typical mean schlieren image of the annular flame is shown in Figs. 4a and 4b. For both these regimes, it is possible to observe the boundaries between the fresh and burned gases, and between those the ambient air, as well as the open flame tip, characteristic of these low Lewis number, around 0.28 and 0.55, flames. These mean images are similar to those analyzed in the previous burner version [5].

<figure><img src="image_5.png" /><figcaption>**Figure 4:** Schlieren images for two cases with $Re = 645$ and different equivalence ratios: intrinsic instability (left); longitudinal pulsation (right).</figcaption></figure>

Regarding the instantaneous images, it is possible to differentiate the two regimes. Note that these instantaneous images are obtained by dividing each frame by the schlieren mean, thus corresponding to the fluctuating part of the temperature gradient. For the case where only the intrinsic instabilities appear, Fig. 4c shows the flame surface wrinkles along the transversal direction of the burner. These wrinkles move rather slowly along the transversal direction. A
rather different instantaneous structure is observed for the case where there is a longitudinal pulsation, Fig. 4d. Indeed, the flame tip opens and closes at a certain frequency, causing a pocket of fresh gases to detach and, perhaps, burn downstream of the average flame tip. Since these flame structures are typical of those found throughout the (Re, φ) ranges given in Fig. 3, for the sake of brevity, only these two cases will be discussed here: the intrinsic instability (Re =645, φ =0.29) and the longitudinal pulsation regime (Re =645, φ =0.31). This Reynolds number value was chosen to lie in the middle of the experimental range given in Fig 3.

Intrinsic instability case SPOD

In order to examine the measured flame structures via SPOD, Fig. 5a shows the SPOD energy spectrum for the intrinsic instability case considering the full image, i.e., ROI 1. It can be seen, for each frequency range obtained through the direct Fourier transform, the energy contained in each SPOD mode. The first mode, highlighted in red, always has the highest energy, which decreases monotonically with increasing frequency. The subsequent modes are colored from black to light gray, thus emphasizing the importance of their energy contribution. For the lowest frequency in the spectra, the energy is around 10¹ and for the highest frequency, the energy drops to the order of 10⁻⁴. This figure suggest that the energy decay rate beyond f ≈ 200 Hz is constant (E*SPOD* ∝ f⁻ⁿ), with 6 ≤ 2n ≤ 7. Beyond approximately 800 Hz the flat spectrum indicates no information remains in the dataset.

<figure><img src="image_6.png" /><figcaption>**Figure 5:** SPOD of the intrinsic instability case (Re = 645, φ = 0.29), ROI 2: (a) Energy spectrum, the first mode is the red line and the following modes go from black to light gray; (b) Energy content (%), normalized by the sum, as a function of the SPOD modes for the arbitrary frequency (367 Hz).</figcaption></figure>

To analyze the percentage of the modal energy distribution at a given frequency, Fig. 5b
shows the SPOD energy normalized by the sum as a function of each mode for the arbitrary
frequency of 367 Hz. The obtained 22 modes are depicted, where the first three modes have
an energy higher than 10¹ and the last two modes have an energy around 10⁰, thus the majority
of these modes are within an energy range of one order of magnitude. Furthermore, it can be
seen that the energy difference between the first and second modes is smaller than 2%. Overall,
the energy decays by one decade approximately every 19 modes. This slow mode decay rate
suggests that, in order to build a reduced order model of intrinsically unstable flames, a large
number of modes would have to be accounted for.

The spatial distribution of the modes, also known as wave packets [17, 18], is illustrated
in Figs. 6 and 7, where we can distinguish the different wave amplitudes. To understand how this distribution behaves along the burner, Fig. 6a shows the wave packet of the first mode for different frequencies obtained from 6,000 frames of instantaneous schlieren. Here, colored regions are the amplitudes of the wave, where blue denotes a negative amplitude and red a positive amplitude.

<figure><img src="image_7.png" /><figcaption>**Figure 6:** Spatial distribution of the first SPOD mode (wave packet), obtained from the full burner (ROI 1). (a) Intrinsic instability case, for the frequency 70 Hz. (b) Longitudinal pulsation case, for the fundamental frequency (366 Hz). The colorbar goes from [-10<sup>-2</sup>, 10<sup>-2</sup>].</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 7:** Spatial distribution of the first and second SPOD modes ($m_1, m_2$) for three arbitrary frequencies ($f_1, f_2, f_3$) = (70, 207 and 367 Hz) in ROI 2. The frequency increases from left to right and the modes increases from top to bottom. The colorbar is identical to Fig. 6.</figcaption></figure>

In Fig. 6a, we can see that the wave packet has a distinctive pattern across the burner. The blue and red zones are intercalated transversely and are regularly arranged, being closely spaced farther away from the center due to the burner curvature effect. This wave packet structure is a reminiscent of the transversal structure of the instability across the burner. Indeed, compared to Fig. 4c, this wave packet resembles the thermodiffusive instability observed in the schlieren. Another aspect to note is the intensity of the colors on the SPOD spatial distribution: the more intense the color, means a higher absolute value of the amplitude. Accordingly, the modulus of the amplitude is higher at the sides of the burner than in the center, which indicates that there is a stronger signal at the sides of the burner, as in the schlieren. This is related to the larger light detection angle at the sides, compared to the center.

Since the SPOD behavior is similar across the burner, the central region (ROI 2) is analyzed in further detail. Figure 7 shows the wave packet of the first and second modes for three
arbitrary frequencies chosen to represent the lower, intermediate and higher frequency ranges of Fig. 5a. Regarding mode 1 behavior, note that a similar pattern to Fig. 6a is seen for all the frequencies, Figs. 7a - 7c. Indeed, departing from the burner surface, the highest absolute amplitude values are near the center of the ROI 2, and decrease further downstream. Moreover, the transversal wavelength of the wave packet seems to be the same for these frequencies. Concerning now mode 2, the overall wave packet patterns are very similar for the each given frequencies, Figs. 7d - 7f, when compared to those of the first mode. The main difference regarding the maximum amplitude occurs in the transverse direction also.

*Longitudinal pulsation case*

Figure 8a shows the SPOD energy spectrum for the case where the longitudinal pulsation occurs. As in the previous case, the energy spectrum of the SPOD modes has a higher energy in the first mode and decreases with the number of modes. However, one remarkable difference between the spectra is that, in Fig. 8a, several peaks are clearly seen. The frequencies of these peaks are the system's harmonics, with the fundamental frequency at 366 Hz, the second harmonic at 732 Hz, the third at 1098 Hz, and so on. The energy difference between the first and the second modes is higher, thus ensuring that these frequencies are significant for the system studied. Note that, for the frequencies between the harmonics, the differences between the energy of the modes are lower and, similar to that of Fig. 5a. This is also the case for the second to the last modes in the harmonics frequencies. Again, the region of high-frequency beyond 4,000 Hz has very low constant energy, bearing is no information, which is corroborated by the spatial distribution of the modes on these frequencies, not shown here for the sake of brevity.

<figure><img src="image_9.png" /><figcaption>**Figure 8:** SPOD of the longitudinal pulsation case ($Re = 645, \phi = 0.31$), ROI 2: (a) Energy spectrum, the first mode is the red line and the following modes goes from black to light gray. The dots mark the frequencies that are analyzed, the green dots are the harmonics and the pink dots are the frequencies between the harmonics that have a local minimum of energy. (b) Energy content (%), normalized by the sum, as a function of the SPOD modes for the fundamental frequency (366 Hz).</figcaption></figure>

To further examine the energy distribution in the harmonics, Fig. 8b shows the SPOD energy as a function of each mode for the fundamental frequency, 366 Hz. The first remarkable difference between this case and the previous, shown in Fig. 5b, is the energy content in the first mode, that is much higher than the second mode, with energies around 49% and 8% respectively. In total, 31 modes are shown and, besides the first mode, the energy of the modes is exponentially decreasing. Furthermore, beyond the second mode, the energy decays by one
decade approximately every 18 modes, which is similar to what was observed in the previous case.

Regarding the spatial distribution of the SPOD modes, Fig. 6b shows the wave packet of the first mode for a frequency of 366 Hz in ROI 1. It is possible to observe that, in most of the domain, the amplitude is nearly zero, having the local maxima and local minimum concentrated on both sides of the image, the lateral part of the burner. Due to this prevalence of the flame lateral region, in contrast with the intrinsic instability case previously presented, Fig. 6b does not show a transversal pattern across the burner for the fundamental frequency, 366 Hz. To analyze in further detail the pattern of the longitudinal pulsation, the following analysis considers the center of the burner (ROI 2) only.

In order to evidence the effect of the longitudinal pulsation on the wave packets, these will be examined for the frequencies marked by dots in Fig. 8a. Accordingly, six frequencies will be analyzed for the case with longitudinal pulsation: the three harmonics; a low frequency (220 Hz) and the two frequencies between the harmonics that have the lowest energy (586 and 952 Hz). Accordingly, Figs. 9 and 10 show the spatial distribution of the three first SPOD modes for these different frequencies.

<figure><img src="image_10.png" /><figcaption>**Figure 9:** Spatial distribution of the first and second SPOD modes ($m_1, m_2$) for the fundamental frequency ($f_1$ = 366 Hz) and the second and third harmonics ($f_2, f_3$) = (732 Hz and 1098 Hz) in ROI 2. The frequency increases from left to right and the modes increases from top to bottom. The colorbar is identical to Fig. 6.</figcaption></figure>

For the first mode (Figs. 9a - 9c) there is a pattern along the *y*-direction, which intermingles zones with positive amplitude (red) with zones that have a negative amplitude (blue). This can be seen for the first mode of all harmonic frequencies, in Fig. 8a. The main difference between these frequencies wave packets is the wavelength, i.e., the distance between the longitudinal red regions. Indeed, higher frequencies correspond to smaller wavelengths, i.e., around 18 mm, 9.7 mm, and, 6.5 mm for the fundamental, first and second harmonic, respectively. Another difference is the position of the highest amplitude: for the fundamental frequency, it is found close to the burner surface, between 0 and 5 mm in the *y*-direction. However, as the frequency increases, this local maximum position is now located progressively further downstream: for the second and third harmonics it is located around 13 and 15 mm in the *y*-direction, respectively.

The first mode of the fundamental frequency is related to the instantaneous schlieren shown in Fig. 4d. Comparing both images, along the *y*-direction, the wave packet structure is a con-
sequence of the pulsating flame behavior observed in the schlieren, in which pockets of fresh gases periodically separate from the flame. The acoustic of the burner system seems to dominate the flame movement, since the first mode is the most energetic.

Even if the second SPOD mode has less energy than the first, it is interesting to analyze the corresponding spatial distribution. Accordingly, Figs. 9d - 9f show the spatial distribution of the wave packets for the first three harmonics of the second mode. These figures show that the wave packet has distinctive structure which seems to be a superposition of longitudinal and transversal patterns. Moreover, it is possible to notice that the regions connecting the positive extreme has an inclination with respect to the x-axis of the burner. This inclination has a tendency to decrease as the frequency increases: for to the fundamental frequency, 366 Hz, this inclination is approximately 72°, and it decreases to 64° and 52° for the second and third harmonics, respectively. Therefore, this wave packet structure of the second mode seems to be influenced both by the longitudinal pulsations, which wavelength decreases with frequency, and by the intrinsic instabilities, which wavelength is independent of the frequency.

Regarding the frequencies that are not the harmonics, Fig. 10 shows the wave packet for first and second modes of the frequencies 220, 586 and 952 Hz. For the frequency 220 Hz, Fig. 10a and 10d show that the wave packet of the first and second modes resembles that of the intrinsic instability case for the frequencies 70, 207 and 367 Hz, which has a transversal pattern. This suggest that the longitudinal pulsations do not have an impact on the wave packets on the low frequencies. Concerning now the frequencies between the harmonics, the first mode, Figs. 10b and 10c, resembles the second mode of the wave packets of the harmonics, i.e., an inclination is seen, which means that, even though the frequency is not an harmonic, the pulsations of the system still influence the wave packet. Concerning the second mode, Figs. 10d - 10f, the wave packets resemble their respective first mode wave packet.

<figure><img src="image_11.png" /><figcaption>**Figure 10:** Spatial distribution of the first and second SPOD modes ($m_1, m_2$) for the three intermediate frequencies ($f_1, f_2, f_3$) = (220, 586 and 952 Hz) in ROI 2. The frequency increases from left to right and the modes increases from top to bottom. The colorbar is identical to Fig. 6.</figcaption></figure>

## Conclusions

The SPOD characterization analysis of ultra-lean H₂-air annular flames provides observations into the energy distribution and spatial structures of intrinsic instabilities and longitudinal pulsations in the burner system. In particular, when intrinsic instabilities are present only, the
schlieren fluctuations energy spectrum decays monotonically with frequency. The SPOD mode
energy was found to decrease rather slowly; i.e., an exponential decrease over 19 modes per
energy decade. Both the frequency and mode energy behaviors suggest that it would be hard
to propose a low rank representation to this intrinsically unstable flame regime. The spatial
distribution analysis reveals intercalated wave packet structures that remain similar across fre-
quencies, resembling the transverse movement of instability along the burner. These findings
underline the persistence and spatial organization of intrinsic instabilities across modes and fre-
quencies.

In contrast, the longitudinal pulsation case analysis shows a distinct energy spectrum with pronounced peaks, i.e., harmonic frequencies, with significant energy content. The first mode in the fundamental frequency carries substantially higher energy than subsequent modes, whereas the decay rate from the second mode on resembles the decay observed in intrinsic instability case, with a 18 mode-per-decade energy decay. Regarding the spatial distribution of the SPOD modes at the system’s harmonics, the first mode exhibits a wave packet that is intercalated longitudinally, with a wavelength that decreases with the frequency. On the other hand, the second mode appears to be a superposition of longitudinal and transverse patterns. This spatial arrangement emphasizes the complex coupling between flame heat release and burner acoustics, highlighting how both intrinsic instabilities and longitudinal pulsation due to the acoustic modes of the burner interact in distinctive behavior across different frequencies and spatial modes.

Acknowledgements

For the purpose of Open Access, a CC-BY public copyright license has been applied by the
authors to the present document and will be applied to all subsequent versions up to the Author
Accepted Manuscript arising from this submission. The authors are grateful for the experimen-
tal support provided by V. Montassier, A. Claverie and N. Papin.

**References**

[1] Shadidi, B., Najafi, G., Yusaf, T., “A review of hydrogen as a fuel in internal combustion engines”, *Energies* 14(19) (2021).

[2] Konnov, A.A., Mohammad, A., Kishore, V.R., Kim, N.I., Prathap, C., Kumar, S., “A comprehensive review of measurements and data analysis of laminar burning velocities for various fuel+air mixtures”, *Progress in Energy and Combustion Science* 68: 197–267 (2018).

[3] Clavin, P., Searby, G., *Combustion Waves and Fronts in Flows: Flames, Shocks, Detona-tions, Ablation Fronts and Explosion of Stars*, Cambridge University Press (2016).

[4] Matalon, M., “The Darrieus–Landau instability of premixed flames”, *Fluid Dynamics Research* 50(5): 051412 (2018).

[5] Figueira da Silva, L.F., “Stochastic characterisation of unstable lean hydrogen–air annular premixed flames”, *Combustion Theory and Modelling* 28(3): 317–343 (2024).

[6] Berger, L., Grinberg, M., Jürgens, B., Lapenna, P.E., Creta, F., Attili, A., Pitsch, H., “Flame fingers and interactions of hydrodynamic and thermodiffusive instabilities in laminar lean hydrogen flames”, *Proceedings of the Combustion Institute* 39(2): 1525–1534 (2023).

[7] Zirwes, T., Zhang, F., Kaiser, T.L., Oberleithner, K., Stein, O.T., Bockhorn, H., Kronen-
burg, A., “The role of thermodiffusion and dimensionality in the formation of cellular
instabilities in hydrogen flames”, *Proceedings of the Combustion Institute* 40(1): 105665
(2024).

[8] Berger, L., Attili, A., Pitsch, H., “Intrinsic instabilities in premixed hydrogen flames: Para-metric variation of pressure, equivalence ratio, and temperature. part I - dispersion relations in the linear regime”, *Combustion and Flame* 240: 111935 (2022).

[9] Berger, L., Attili, A., Pitsch, H., “Intrinsic instabilities in premixed hydrogen flames: para-
metric variation of pressure, equivalence ratio, and temperature. part II – non-linear regime and flame speed enhancement", *Combustion and Flame* 240: 111936 (2022).

[10] Berger, L., Attili, A., Pitsch, H., "Synergistic interactions of thermodiffusive instabilities and turbulence in lean hydrogen flames", *Combustion and Flame* 244: 112254 (2022).

[11] Berger, L., Kleinheinz, K., Attili, A., Pitsch, H., "Characteristic patterns of thermodiffusively unstable premixed lean hydrogen flames", *Proceedings of the Combustion Institute* 37(2): 1879-1886 (2019).

[12] Wen, X., Berger, L., Cai, L., Parente, A., Pitsch, H., "Thermodiffusively unstable laminar hydrogen flame in a sufficiently large 3D computational domain – part I: Characteristic patterns", *Combustion and Flame* 263: 113278 (2024).

[13] Böttler, H., Kaddar, D., Karpowski, T.J.P., Ferraro, F., Scholtissek, A., Nicolai, H., Hasse, C., "Can flamelet manifolds capture the interactions of thermo-diffusive instabilities and turbulence in lean hydrogen flames?—an a-priori analysis", *International Journal of Hydrogen Energy* 56: 1397-1407 (2024).

[14] Welch, C., Erhard, J., Shi, H., Dreizler, A., Böhm, B., "An experimental investigation of lean hydrogen flame instabilities in spark-ignition engines", *Proceedings of the Combustion Institute* 40(1): 105391 (2024).

[15] Ezenwajiaku, C., Balachandran, R., Ducci, A., Picciani, M., Talibi, M., "Experimental characterisation of the dynamics of partially premixed hydrogen flames in a lean direct injection (LDI) combustor", *Proceedings of the ASME Turbo Expo 2023: Turbo-machinery Technical Conference and Exposition* 3A: Combustion, Fuels, and Emissions: V03AT04A073 (2023).

[16] Kang, H., Lee, M., Kim, K.T., "Measurements of self-excited instabilities and nitrogen oxides emissions in a multi-element lean-premixed hydrogen/methane/air flame ensemble", *Proceedings of the Combustion Institute* 39(4): 4721-4729 (2023).

[17] Padilla-Montero, I., Rodríguez, D., Jaunet, V., Jordan, P., "Eduction of coherent structures from schlieren images of twin jets using spod informed with momentum potential theory in the spectral domain", *Theoretical and Computational Fluid Dynamics* 38: 375-401 (2024).

[18] Kaplan, O., Jordan, P., Cavalieri, A.V., Brès, G.A., "Dynamics of turbulent boundary layers exciting wavepackets in subsonic jets", in *25th AIAA/CEAS Aeroacoustics Conference* (2019).

[19] Fiore, M., Parisot-Dupuis, H., Etchebarne, B., Gojon, R., "Spectral proper orthogonal decomposition of coupled hydrodynamic and acoustic fields: Application to impinging jet configurations (draft)", *Computers & Fluids* 241: 105484 (2022).

[20] Muthichur, N., Vempati, C., Hemchandra, S., Samanta, A., "Reduced-order models of aeroacoustic sources for sound radiated in twin subsonic jets", *Journal of Fluid Mechanics* 972: A33 (2023).

[21] Sieber, M., Paschereit, C.O., Oberleithner, K., "Spectral proper orthogonal decomposition", *Journal of Fluid Mechanics* 792: 798-828 (2016).

[22] Schmidt, O.T., Colonius, T., "Guide to spectral proper orthogonal decomposition", *AIAA Journal* 58(3): 1023-1033 (2020).

[23] Schmidt, O.T., "Spectral proper orthogonal decomposition using multitaper estimates", *Theoretical and Computational Fluid Dynamics* 36(5): 741-754 (2022).

[24] Towne, A., Schmidt, O.T., Colonius, T., "Spectral proper orthogonal decomposition and its relationship to dynamic mode decomposition and resolvent analysis", *Journal of Fluid Mechanics* 847: 821-867 (2018).