## Strange Nonchaotic Attractor in an Unforced Turbulent Reactive Flow System

Thonti Beeraiah, 1, 2 Shruti Tandon, 1, 2 Premraj Durairaj, 3 and R. I. Sujith 1, 2

2

1 Department of Aerospace Engineering, Indian Institute of Technology Madras, Chennai 600 036, India Centre of Excellence for Studying Critical Transition in Complex Systems,

Indian Institute of Technology Madras, Chennai 600 036, India

3 Institute of Systems Science, Huaqiao University, Xiamen, 361021, China

(Dated: August 5, 2024)

We discover strange nonchaotic attractor (SNA) through experiments in an unforced system comprising turbulent reactive flow. While models suggest SNAs are common in dynamical systems, experimental observations are primarily limited to systems with external forcing. We observe SNA prior to the emergence of periodic oscillations from chaotic fluctuations. In complex systems, self-organization can lead to order, and inherent nonlinearity can induce chaos. The occurrence of SNA, which is nonchaotic yet nonperiodic in one such complex system, is intriguing.

Introduction -Strange nonchaotic attractor (SNA) is a non-trivial dynamical state that has attracted sustained interest in the nonlinear dynamics community [1-6]. SNAs exhibit intricate fractal geometries akin to those observed in chaotic systems. However, unlike chaotic systems, SNAs do not exhibit exponential sensitivity to initial conditions and are hence characterized as nonchaotic. Typically, a positive largest Lyapunov exponent indicates sensitivity to initial conditions; for SNAs, this exponent is nonpositive [1].

We discover the state of SNA in a system comprising turbulent reactive flow. Our finding is particularly intriguing given that our system is not externally forced and is rather a selforganized complex system owing to a web of inter-subsystem interactions [7, 8]. Turbulent flows are ubiquitous in nature and consist of eddies ranging from very small scale to very large scale [9]. These eddies interact nonlinearly, making turbulent flow a quintessential complex system [10]. Introducing an exothermic reaction into a turbulent flow adds to the intricacy and multifaceted nature of the system.

Turbulent reactive flows occur in natural systems, including stars, which are thermonuclear reactors composed of plasma, bounded by gravity. These stars undergo turbulent mixing [11], and in rotating stars, this turbulent mixing generates a magnetic field which, in turn, establishes a nonlinear interaction between turbulent mixing, rotation, and magnetic field. Such stars may exhibit chaotic fluctuations [12] or unstable pulsation modes [11]. Turbulent reactive flows also exist in human-made systems, such as combustors in gas turbine and rocket engines. Turbulent reactive flows in combustors comprise subsystems, namely the hydrodynamic field, heat release (flame) field, and acoustic field that interact in a nonlinear manner, increasing the overall complexity [13]. During the stable operation in such systems, the acoustic pressure fluctuations ( p ′ ) are low-amplitude and aperiodic in nature and are characterized as high-dimensional chaos [14, 15]. During this state, the subsystems are desynchronized [16, 17]. In contrast, during unstable operation, high-amplitude periodic acoustic pressure oscillations occur, characterized as limit cycle oscillations (order). Order emerges from the chaotic fluctuations when positive feedback is established between the heat release rate and the acoustic field, leading to their synchronization [7]. In these systems, the presence of chaotic and ordered states, and the transition from chaotic to ordered states via a state of intermittency-where high amplitude periodic epochs are interspersed with low amplitude aperiodic fluctuations-is well established [18, 19]. However, the discovery of the state of SNA, which is nonchaotic yet nonperiodic, in such a complex system is surprising.

FIG. 1. (a) Schematic of the experimental setup, (b) bluff body mounted on a shaft fixed at 30 mm from the backward-facing step, and (c) variation of root-mean-square of acoustic pressure fluctuations ( p ′ ) as a function of Re during the transition, with the highlighted points indicating various dynamical states present.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01131v1-figure-9ca93fca-fa99-4995-a27d-6d04e95d8482.jpg)

Typically, SNAs are found in dynamical systems driven by external quasiperiodic [16, 20-25] or periodic forcing [26, 27]. SNAs have also been observed in experimental systems with periodic forcing [28, 29]. Furthermore, a few studies using mathematical models have demonstrated that neither quasiperiodic force nor periodic force is essential for the occurrence of SNA [30, 31]. However, to date, only two physical systems have been reported to exhibit SNAs without being subjected to external force: a pulsating star KIC 5520878 network [32] and a laminar thermoacoustic system [33].

In this letter, we present the experimental discovery of state of SNA in a highly complex turbulent reactive flow system. In the experiments, we increase the Reynolds number ( Re ) as a bifurcation parameter in a quasi-static manner and observe a transition from lowamplitude aperiodic fluctuations to high-amplitude periodic oscillations. We show that the state of SNA occurs prior to the onset of high-amplitude periodic oscillations. While the computation of Lyapunov exponents is generally useful for determining the nature of dynamical states, ranging from regular to chaotic, it is often unreliable for noise-contaminated experimental data [34]. Therefore, we employ Fast Fourier Transform (FFT) and phase space reconstruction to distinguish different dynamical states. We then apply the 0 -1 test [35] and correlation dimension test [36] to characterize these states. Finally, to confirm the presence of SNA, we use singular continuous spectrum analysis [37].

Experimental setup -Weperform experiments in a lab-scale backward-facing step turbulent combustor, using a bluff body as a flame stabilizer, as depicted in figure 1a. The bluff body is fixed at a location of 30 mm from the backward-facing step throughout the

FIG. 2. Time series (a-e), Amplitude spectrum (f-j), and reconstructed phase space (k-o) of the five dynamical states. During the transition, periodicity emerges, as reflected in the amplitude spectrum and reconstructed phase space.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01131v1-figure-d876164b-7a9f-4d5f-927e-571ff2c6c18b.jpg)

transition. The setup comprises three main components: a plenum chamber, a burner, and a combustion chamber. The combustion chamber has a square cross-section measuring 90 mm × 90 mm and a length of 1100 mm. The plenum chamber connects to the burner, which further leads to the combustion chamber. Gaseous fuel is injected into the air in the burner section, creating a partially premixed fuel-air mixture by the time it reaches the combustion chamber, where all reactions occur. To vary the Reynolds number ( Re ), we fix the fuel flow rate and increase the air flow rate using mass flow controllers (Alicat Scientific, MCR). The uncertainty in the measurement of the flow rate is ± (0.8% of reading + 0.2% of full scale), resulting in a maximum uncertainty of ± 2 5% in . Re . We flush-mount a piezoelectric pressure transducer (PCB103B02) on the combustion chamber at 40 mm from the backward-facing step to measure acoustic pressure fluctuations ( p ′ ). The sensor has an uncertainty of ± 0 15 . Pa. During the transition, we measure the acoustic pressure oscillations ( p ′ ) at a sampling rate of 20 kHz.

Results - In our experiments, we vary Re from 28 000 to 44 000 and observe a transition , , from low-amplitude aperiodic fluctuations to high-amplitude periodic oscillations (Fig. 1c). We primarily use Fast Fourier Transform and phase space reconstruction of p ′ to distinguish the dynamical states present, labeled as states I to V (Fig. 2). As the transition occurs from state I to V, we observe an increase in the amplitude of fluctuations in p ′ (Fig. 2a-

e). Concurrently, the periodicity of the signal increases, culminating in the emergence of a single prominent peak in the amplitude spectrum during state V (Fig. 2j) compared to the broad peak during state I (Fig. 2f). To further characterize these states, we apply 0 -1 test and correlation dimension test. Finally, we employ the singular continuous spectrum test to confirm the SNA, which we identify through the preceding analyses.

a. 0 -1 test Originally, Gottwald and Melbourne [35] proposed this method to discriminate between regular and chaotic time series. The initial step of the 0 -1 test is to transform the original signal into a new space where the underlying dynamics becomes more apparent. For the same purpose we calculate the translation variables ( p n ( n ) and q n ( n )), as defined by,

$$p _ { n } ( n ) = \sum _ { i = 1 } ^ { n } p ^ { \prime } ( i ) \cos ( i c ),$$

$$q _ { n } ( n ) = \sum _ { i = 1 } ^ { n } p ^ { \prime } ( i ) \sin ( i c ), n = 1, 2, 3, \dots, & & ( 1 b )$$

where, p ′ ( ) i denotes the acoustic pressure fluctuations at instant i and c is a parameter selected arbitrarily within π/ 5 to 4 π/ 5. The dynamics of variables p n ( n ) and q n ( n ) are closely related to the discrete variable n , which is significantly smaller than the length of the signal ( N ) (typically n = N/ 10). To quantify the evolution of the trajectory in the [ p n ( n , q ) n ( n )] plane as n increases, the mean square displacement D n ( ) is utilized, defined as:

$$D ( n ) = \frac { 1 } { N } ( \sum _ { i = 1 } ^ { N } [ p _ { n } ( i + n ) - p _ { n } ( i ) ] ^ { 2 } + [ q _ { n } ( i + n ) - q _ { n } ( i ) ] ^ { 2 } ).$$

To address convergence issues with D n ( ), a modified mean square displacement M n ( ) [38, 39] is proposed, which is calculated as:

$$M ( n ) = D ( n ) - V _ { o s c } ( c, n ),$$

where, V osc ( c, n ) = ∑ n i =1 p ′ ( ) i 1 -cos( ic ) 1 -cos( c ) , is an oscillatory term. The asymptotic growth rate K c is determined through linear regression fit of M n ( ) as:

$$K _ { c } = \lim _ { n \to \infty } \frac { \log M ( n ) } { \log n }.$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01131v1-figure-b13f1ae6-41d3-4f62-b0bf-2c66a621da74.jpg)

n

FIG. 3. Modified mean square displacement ( M n ( )) (top row) and asymptotic growth rate ( K c ) (bottom row) during the dynamical states of (a, f) I, (b, g) II, (c, h) III, (d, i) IV, and (e, j) V. The value of K c decreases gradually from 1 to 0 during the transition from state I to V, indicating the emergence of periodicity.

The dynamical states are distinguished using the variation of M n ( ) and magnitude of K c . The value of M n ( ) grows linearly for chaotic dynamics owing to sensitivity to initial conditions while it remains bounded for periodic dynamics [38, 39]. This behavior causes the value of K c to fall within the range of 1 to 0, where a value of K c = 1 suggests chaotic dynamics and K c = 0 implies periodic behavior. For the dynamical state of SNA, the value of K c is distributed between 0 and 1, reflecting its nature as a state that is neither chaotic nor periodic.

In figure 3, we plot M n ( ) and K c for the dynamical states I-V. For state I, the value of M n ( ) increases linearly on average(Fig. 3a), and the value of K c is fairly lying close to 1 (Fig. 3f) depicting a chaotic nature of the state. During state II also, the value of M n ( ) increases linearly on average with slight fluctuations (Fig. 3b) and the values of K c distributed just below 1 (Fig. 3g), showing a chaotic nature.

For state III, the value of M n ( ) either increases or decreases in a bounded range on average as evident in figure 3c and the value of K c is distributed between 1 and 0 (Fig. 3h), depicting that the state III as neither chaotic nor periodic. During state IV, the value of M n ( ) becomes bounded and the average value remains constant with fluctuations present (Fig. 3d). The value of K c distributed just above 0 (Fig. 3i), showing a periodic behavior. During state V also, the value of M n ( ) bounded and the average value remains constant with fluctuations present (Fig. 3e), the values of K c are fairly lying close to 0 (Fig. 3j), delineating a clear periodic behavior, compared to the previous state.

FIG. 4. (a-e) Log-log plots showing the variation of correlation sum (C( )) with ϵ ϵ for the states I to V. The plots are shown for the embedding dimensions ranging from 5( · ), 6 ( ⋆ ) and 7 ( □ ). The solid line represents the power law fit for the straight-line portion of this plot, indicating the value of the correlation dimension ( D ). The value of D decreases as the system transitions from state I to state V, indicating the decrease of complexity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01131v1-figure-07237860-28ab-4905-8821-7025cad62e8e.jpg)

b. Correlation Dimension The correlation dimension is a topological measure of the fractal dimension of an attractor [40]. We compute the correlation dimension using the algorithm proposed by Grassberger and Procaccia [36]. The initial step is to calculate the delayed vectors s n = { p ′ i ( m ) } N -( m -1) τ i =1 from the time series of acoustic pressure fluctuations ( p ′ ). Here, N is the length of the time series, m is the embedding dimension, and τ is the time delay for which the average mutual information of the signal is minimum [41]. Then, we calculate the correlation sum as follows:

$$C ( \epsilon ) = \frac { 2 } { N ( N - 1 ) } \sum _ { j = 1 } ^ { N } \sum _ { k = j + 1 } ^ { N } \Theta ( \epsilon - | | s _ { j } - s _ { k } | | )$$

where Θ is the Heaviside function defined as:

$$\text{where $\Theta$ is the Heaviside} \\ \Theta ( x ) = \begin{cases} 1, & x > 0, \\ 0, & x \leq 0, \end{cases}$$

and ϵ is the threshold distance. The sum counts the number of pairs ( s , s j k ), whose distance is smaller than ϵ . In the limit of infinite data ( N →∞ ) and ϵ → 0, C ϵ ( ) scales as a power law i.e., C ϵ ( ) ∝ ϵ D . We can define the correlation dimension ( D ) from the correlation sum as follows:

$$d ( N, \epsilon ) = \frac { \partial ( \ln C ( \epsilon, N ) ) } { \partial ( \ln ( \epsilon ) ) },$$

$$D = \lim _ { N \to \infty } d ( N, \epsilon ).$$

The value of D can be used to estimate the dynamical state of a system. In general, the value of D is noninteger for the strange attractor, and it is an integer and approximately one for periodic attractors [42].

It is sufficient to choose a value of m greater than the box dimension to estimate the value of D [43]. To ensure the convergence of c ϵ ( ), we compute it for embedding dimensions ranging from 5 -7. In figure 4, we plot the variation of C ϵ ( ) with increasing ϵ at three embedding dimensions in a logarithmic plot for the states I to V. We apply a power law fit to the straight-line portion to obtain the correlation dimension ( D ).

For state I, the value of D is 4 13 (Fig. 4a). . This high fractional value of D indicates the necessity of a high-dimensional structure to describe the underlying chaotic dynamics. For the state II, the value of D is 3 87 (Fig. 4b), reflecting the need for a high-dimensional . structure to capture the dynamics of epochs of high-amplitude periodic oscillations amidst epochs of low-amplitude aperiodic fluctuations (Fig. 2b). The dynamics in the phase space indicate that the system switches between low-amplitude aperiodic fluctuations and highamplitude periodic oscillations (Fig. 2l), characteristics of the state of intermittency.

For state III, the reconstructed phase space exhibits complex folding and stretching (Fig. 2m), requiring a high value of dimension to fully describe the dynamics. The observed value of D of 2 89 (Fig. 4c) suggests a complex structure, although less intricate than chaos . and intermittency.

For state IV, the reconstructed phase space is found to be a thick ring (Fig. 2n). The trajectory during a few acoustic cycles deviates from the mean, resulting in a large thickness. This is reflected by a value of D of 1 62 (Fig. 4d), indicating the need for an additional . effective dimension to capture the noisy behavior. These noisy periodic oscillations with amplitude modulations are also called noisy limit cycle oscillations.

For state V, the reconstructed phase space exhibits a thin ring structure (Fig. 2o), indicating periodic oscillations with fewer amplitude modulations. The value of D is observed to be 1 08 (Fig. 4e), which is very close to the integer value 1. . These clean periodic oscillations essentially clean limit cycle oscillations.

The results of 0 -1 test show that state III is nonchaotic and nonperiodic. Additionally, the value of D falls between those of the chaotic and periodic states suggesting fractal characteristics. These findings imply that state III may be a strange nonchaotic attractor (SNA), which is a nonchaotic fractal. To confirm this, we perform a singular continuous

FIG. 5. Singular continuous spectrum analysis of acoustic pressure oscillations ( p ′ ) for state III. (a) log 10 | X α,N ( ) | 2 vs log 10 N . We have | X α,N ( ) | 2 ∼ N 1 49 . . The subplot shows the corresponding trajectory (Re X , Im X ), which is a fractal. (b) fractal dimension of 1 9 calculated . using boxcounting method following a power law with increasing box size.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01131v1-figure-5fb7424f-5f61-4b10-bf79-6113b5e20af6.jpg)

spectrum analysis.

- c. Singular continuous spectrum In dynamical systems, the power spectrum is continuous during chaotic dynamics and discrete during regular dynamics. A continuous spectrum displays a broadband distribution, while sharp peaks characterize a discrete spectrum. The singular continuous spectrum is a mixture of both continuous and discrete spectrums. For the time series of acoustic pressure fluctuations p ′ , the singular continuous spectrum is defined as [37, 44],

$$X ( \alpha, N ) = \sum _ { k = 1 } ^ { N } p ^ { \prime } ( k ) e ^ { 2 \pi i k \alpha }.$$

Here, α = √ 5 -1 2 , the golden mean ratio, and N is the length of the time series of p ′ . X α,N ( ) is a complex value. For chaotic signals, the power of the signal | X α,N ( ) | 2 is proportional to N , indicating a continuous (broadband) spectrum. For periodic signals, the power of signal | X α,N ( ) | 2 is proportional to N 2 , illustrating a discrete (sharp) spectrum. For the state of SNA, the spectrum exhibits a singular continuous spectrum. For SNA the power of signal | X α,N ( ) | 2 is proportional to N γ and the value of γ ranges between 1 and 2.

We apply singular continuous spectrum analysis for the state III. We observe that the value of γ is 1 49 (Fig. 5a) and the trajectory of . X α,N ( ) on the complex plane is a fractal as shown in the subplot of figure 5a. To confirm the fractal structure of the complex plane, we apply the box-counting dimension and calculate the fractal dimension of the structure.

- d. Box-counting dimension The concept of dimension is fundamental to fractal geometry, with self-similarity being a key characteristic of fractals [45]. There are various

definitions of fractal dimension (FD) [46], each suitable for specific contexts. Among the methods available to estimate FD, the box-counting method is frequently used due to its simplicity [47]. The FD can be calculated using the method as follows:

$$F D = - \lim _ { \epsilon \rightarrow 0 } \frac { \log ( N ( \epsilon ) ) } { \log ( \epsilon ) },$$

where N ϵ ( ) is the number of boxes of size ϵ required to completely cover the structure. The FD is estimated from the least square linear fit of log-log plot of N ϵ ( ) against ϵ .

We apply the box-counting method on the image of the structure in the complex plane shown in the subplot of figure 5b and observe that the FD of the image is 1 9 (Fig. 5c). . This fractional value of FD indicates that the trajectory (Re X , Im X ) in the complex plane is fractal in nature, confirming the presence of SNA.

In summary, we report the state of SNA in a self-organized system comprising turbulent reactive flow. We observe the state of SNA prior to the emergence of periodic oscillations from chaos. We distinguish the dynamical states present during the transition using tools from nonlinear dynamics, such as phase space reconstruction, 0 -1 test, and correlation dimension test. We substantiate the existence of the state of SNA with singular continuous spectrum analysis. The presence of SNAs offers the advantage of robustness to noise over chaos in synchronization phenomena owing to their insensitivity to initial conditions. This advantage of SNA and the possibility of SNAs in self-organized complex systems opens a new avenue for research. Further, exploration of spatiotemporal dynamics during the state of SNA could provide deeper insights into synchronization and other nonlinear phenomena.

We thank S. Anand, S. Thilagaraj, G. Sudha, M. Raghunathan, S. Sudarsanan, Ramesh S Bhavi, and Anaswara Bhaskaran for their help during the experiments. T.B. acknowledges the research assistantship from the Ministry of Human Resource Development, India, and the Indian Institute of Technology Madras. S.T. acknowledges the support from Prime Minister's Research Fellowship, Govt. of India. R.I.S. acknowledges the funding from the Science and Engineering Research Board (SERB) of the Department of Science and Technology (DST) through a J. C. Bose Fellowship (No. JCB 2018 000034 SSC) and from the IOE initiative / / /

(No. SP22231222CPETWOCTSHOC).

- [1] C. Grebogi, E. Ott, S. Pelikan, and J. A. Yorke, Physica D: Nonlinear Phenomena 13 , 261 (1984).
- [2] M. Ding, C. Grebogi, and E. Ott, Physical Review A 39 , 2593 (1989).
- [3] U. Feudel, J. Kurths, and A. S. Pikovsky, Physica D: Nonlinear Phenomena 88 , 176 (1995).
- [4] F. J. Romeiras and E. Ott, Physical Review A 35 , 4404 (1987).
- [5] J. Heagy and W. Ditto, Journal of Nonlinear Science 1 , 423 (1991).
- [6] T. Zhou, F. Moss, and A. Bulsara, Physical Review A 45 , 5394 (1992).
- [7] R. I. Sujith and S. A. Pawar, Thermoacoustic Instability (Springer, 2021).
- [8] S. Tandon and R. I. Sujith, Journal of Fluid Mechanics 966 , A9 (2023).
- [9] H. Tennekes and J. L. Lumley, A First Course in Turbulence , MIT Press (M.I.T.Pr., 1972).
- [10] R. T. Pierrehumbert, Journal of Fluid Mechanics 938 , F1 (2022).
- [11] W. D. Arnett and C. Meakin, Reports on Progress in Physics 79 , 102901 (2016).
- [12] J. R. Buchler, T. Serre, Z. Koll´th, a and J. Mattei, Phys. Rev. Lett. 74 , 842 (1995).
- [13] B.-T. Chu and L. S. G. Kovasznay, Journal of Fluid Mechanics 3 , 494-514 (1958).
- [14] V. Nair, G. Thampi, S. Karuppusamy, S. Gopalan, and R. I. Sujith, International Journal of Spray and Combustion Dynamics 5 , 273 (2013).
- [15] J. Tony, E. Gopalakrishnan, E. Sreelekha, and R. I. Sujith, Physical Review E 92 , 062902 (2015).
- [16] S. A. Pawar, A. Seshadri, V. R. Unni, and R. I. Sujith, Journal of Fluid Mechanics 827 , 664 (2017).
- [17] S. Mondal, V. R. Unni, and R. I. Sujith, Journal of Fluid Mechanics 811 , 659 (2017).
- [18] V. Nair, G. Thampi, and R. I. Sujith, Journal of Fluid Mechanics 756 , 470 (2014).
- [19] C. Aoki, H. Gotoda, S. Yoshida, and S. Tachibana, Journal of Applied Physics 127 , 224903 (2020).
- [20] M. Ding and J. S. Kelso, International Journal of Bifurcation and Chaos 4 , 553 (1994).
- [21] A. Bondeson, E. Ott, and T. M. Antonsen Jr, Physical Review Letters 55 , 2103 (1985).
- [22] J. Heagy and S. Hammel, Physica D: Nonlinear Phenomena 70 , 140 (1994).
- [23] A. Prasad, V. Mehra, and R. Ramaswamy, Physical Review Letters 79 , 4127 (1997).

- [24] X. Wang, M. Zhan, C.-H. Lai, and Y.-C. Lai, Physical Review Letters 92 , 074102 (2004).
- [25] A. Y. Jalnine and S. Kuznetsov, Technical Physics 52 , 401 (2007).
- [26] V. Anishchenko, T. Vadivasova, and O. Sosnovtseva, Physical Review E 54 , 3231 (1996).
- [27] Z. LIU, International Journal of Bifurcation and Chaos 11 , 225 (2001).
- [28] Y. Guan, M. Murugesan, and L. K. Li, Chaos: An Interdisciplinary Journal of Nonlinear Science 28 (2018), 10.1063/1.5026252.
- [29] A. Chithra and I. Raja Mohamed, Nonlinear Dynamics 105 , 3615 (2021).
- [30] S. S. Negi and R. Ramaswamy, Pramana 56 , 47 (2001).
- [31] S. Kumarasamy, A. Srinivasan, M. Ramasamy, and K. Rajagopal, Chaos: An Interdisciplinary Journal of Nonlinear Science 32 (2022), 10.1063/5.0089373.
- [32] J. F. Lindner, V. Kohar, B. Kia, M. Hippke, J. G. Learned, and W. L. Ditto, Physical Review Letters 114 , 054101 (2015).
- [33] D. Premraj, S. A. Pawar, L. Kabiraj, and R. I. Sujith, Europhysics Letters 128 , 54005 (2020).
- [34] P. Bryant, R. Brown, and H. D. I. Abarbanel, Phys. Rev. Lett. 65 , 1523 (1990).
- [35] G. A. Gottwald and I. Melbourne, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 460 , 603 (2004).
- [36] P. Grassberger and I. Procaccia, Physical Review Letters 50 , 346 (1983).
- [37] A. S. Pikovsky and U. Feudel, Journal of Physics A: Mathematical and General 27 , 5209 (1994).
- [38] G. A. Gottwald and I. Melbourne, SIAM Journal on Applied Dynamical Systems 8 , 129 (2009).
- [39] G. A. Gottwald and I. Melbourne, Nonlinearity 22 , 1367 (2009).
- [40] H. Kantz and T. Schreiber, Nonlinear Time Series Analysis , 2nd ed. (Cambridge University Press, 2003).
- [41] A. M. Fraser and H. L. Swinney, Physical Review A 33 , 1134 (1986).
- [42] C. J. McMahon, J. P. Toomey, and D. M. Kane, PLoS One 12 , e0181559 (2017).
- [43] T. Sauer and J. A. Yorke, International Journal of Bifurcation and Chaos 03 , 737 (1993).
- [44] T. Yal¸ cınkaya and Y.-C. Lai, Phys. Rev. E 56 , 1623 (1997).
- [45] B. B. Mandelbrot, The Fractal Geometry of Nature , Vol. 1 (WH Freeman New York, 1982).
- [46] K. Falconer, 'Alternative definitions of dimension,' in Fractal Geometry (John Wiley and Sons, Ltd, 2003) Chap. 3.
- [47] H.-O. Peitgen, H. J¨rgens, u D. Saupe, and M. J. Feigenbaum, Chaos and Fractals: New

Frontiers of Science , Vol. 106 (Springer, 2004).
Beeraiah Thonti Shruti Tandon Premraj Durairaj R. I. Sujith
Strange Nonchaotic Attractor in an Unforced Turbulent Reactive Flow System
We discover strange nonchaotic attractor (SNA) through experiments in an
unforced system comprising turbulent reactive flow. While models suggest SNAs
are common in dynamical systems, experimental observations are primarily
limited to systems with external forcing. We observe SNA prior to the emergence
of periodic oscillations from chaotic fluctuations. In complex systems,
self-organization can lead to order, and inherent nonlinearity can induce
chaos. The occurrence of SNA, which is nonchaotic yet nonperiodic in one such
complex system, is intriguing.