Fully compressible simulations of the impact of acoustic waves on the dynamics of laminar premixed flames for engine-relevant conditions

Thermo-acoustic instabilities remain problematic in the design of propulsion systems such as gas turbine engines, rocket motors, and ramjets. They arise from the constructive interaction of heat release rate and acoustic pressure oscillations, and can result in increased noise and mechanical fatigue. In the present work, we are concerned with the flame response to the thermodynamic fluctuations that accompany an incident acoustic wave. The objective is to investigate the flame dynamics under engine-relevant conditions using high-fidelity numerical simulations and detailed chemical kinetics. The focus is placed on the combustion of hydrogen and n-heptane, as they are both of practical interest and behave very differently when subjected to acoustic waves. We extract the phase and gain of the unsteady heat release response, which are directly related to the Rayleigh criterion and thus the stability of the system. We highlight the differences between results obtained using the fully compressible Navier-Stokes equations and the low Mach number approximation. The two simulation frameworks agree very well for acoustic wavelengths much larger than the flame thickness. However, they differ significantly at high frequencies. The gain erroneously reaches a plateau under the low Mach number approximation, while it decays to zero using the fully compressible framework. This difference is attributed to the spatial variations in the acoustic pressure, which are not captured by the low Mach number approximation.     

Ethanol turbulent spray flame response to gas velocity modulation

A numerical investigation of the interaction between a spray flame and an acoustic forcing of the velocity field is presented in this paper. In combustion systems, a thermoacoustic instability is the result of a process of coupling between oscillations in heat released and acoustic waves. When liquid fuels are used, the atomisation and the evaporation process also undergo the effects of such instabilities, and the computational fluid dynamics of these complex phenomena becomes a challenging task. In this paper, an acoustic perturbation is applied to the mass flow of the gas phase at the inlet and its effect on the evaporating fuel spray and on the flame front is investigated with unsteady Reynolds averaged Navier-Stokes numerical simulations. Two flames are simulated: a partially premixed ethanol/air spray flame and a premixed pre-vaporised ethanol/air flame, with and without acoustic forcing. The frequencies used to perturb the flames are 200 and 2500 Hz, which are representative for two different regimes. Those regimes are classified based on the Strouhal number St = (D/U)f_f: at 200 Hz, St = 0.07, and at 2500 Hz, St = 0.8. The exposure of the flame to a 200 Hz signal results in a stretching of the flame which causes gas field fluctuations, a delay of the evaporation and an increase of the reaction rate. The coupling between the flame and the flow excitation is such that the flame breaks up periodically. At 2500 Hz, the evaporation rate increases but the response of the gas field is weak and the flame is more stable. The presence of droplets does not play a crucial role at 2500 Hz, as shown by a comparison of the discrete flame function in the case of spray and pre-vaporised flame. At low Strouhal number, the forced response of the pre-vaporised flame is much higher compared to that of the spray flame.     

Combustion dynamics of inverted conical flames

An inverted conical flame anchored on a central bluff-body in an unconfined burner configuration features a distinctive acoustic response. This configuration typifies more complex situations in which the thermo-acoustic instability is driven by the interaction of a flame with a convective vorticity mode. The axisymmetric geometry investigated in this article features a shear region between the reactive jet and the surrounding atmosphere. It exhibits self-sustained oscillations for certain operating conditions involving a powerful flame collapse phenomenon with sudden annihilation of flame surface area. This is caused by a strong interaction between the flame and vortices created in the outer jet shear layer, a process which determines the amplitude of heat release fluctuation and its time delay with respect to incident velocity perturbations. This process also generates an acoustic field that excites the burner and synchronizes the vortex shedding mechanism. The transfer functions between the velocity signal at the burner outlet and heat release are obtained experimentally for a set of flow velocities fluctuations levels. It is found that heat release fluctuations are a strong function of the incoming velocity perturbation amplitude and that the time delay between these two quantities is mainly determined by the convection of the large scale vortices formed in the jet shear layer. A model is formulated, which suitably describes the observed instabilities.     

The Rayleigh integral is always positive in steadily operated combustors

The Rayleigh index has been used for decades by a large number of researchers as an indicator to determine if a flame is driving or damping thermoacoustic interaction mechanisms. The use of the Rayleigh criterion has found applications in rocket combustors, gas turbine combustion technology and basic combustion research. The global Rayleigh index or integral is obtained by integrating the product of heat release rate and pressure fluctuations over space and time. Depending on the phase between pressure oscillations and heat release rate response, the oscillations can be enhanced or damped. It is commonly assumed in literature that the sign of the Rayleigh index from steady state data can be used to determine if the thermoacoustic feedback loop is stabilizing or destabilizing. However, we show in this paper that under fairly general conditions, a correctly measured Rayleigh index is always positive if evaluated from statistically stationary data. This proves to be true even if the heat release rate response to pressure fluctuations is in phase opposition to those pressure fluctuations. This is shown in a straightforward manner by substituting the wave equation with a heat release rate source term into the Rayleigh index. This was verified experimentally on a fully premixed combustion system by measuring the flame chemiluminescence using a photo multiplier and pressure fluctuations using a microphone placed sufficiently close to the flame to ensure acoustic compactness for the frequency range of interest. A large range of operating conditions have been tested, spanning linearly stable and unstable stationary thermoacoustic states, respectively corresponding to resonance or a limit cycle driven by the inherent stochastic forcing from the turbulent combustion noise. The experimental results corroborated the analytic finding: the Rayleigh index is found to be positive for all frequencies and all operating conditions.     

氢、碳氢燃料对向扩散火焰

对氢、正烷烃碳氢燃料与氧的对向扩散火焰,其中正烷烃包含了在工业用燃料中广泛应用的C_nH_2n+2正烷烃CH₄～C1₆H₃₄,对这些燃料的火焰结构进行了分析和比较,系统地分析了压力和拉伸率对火焰行为和热释放率等的影响,其中包含了2115个组分8157个可逆反应．研究结果表明,所有燃料的火焰厚度和热释放率与压力和拉伸率的乘积的平方根成线性关系．在相同工况下,氢的火焰厚度总是大于所有的碳氢燃料,而CH₄～C1₆H₃₄所有的碳氢燃料在相同工况下总是具有几乎相同的燃烧温度分布、燃烧产物分布、火焰厚度和热释放率,该结果表明由这些碳氢燃料组成的混合燃料具有同样的火焰特性.     

Phase synchronization and collective interaction of multiple flamelets in a laboratory scale spray combustor

In this paper, we investigate the coupled behvior of the acoustic field in the confinement and the unsteady flame dynamics in a laboratory scale spray combustor. We study this interaction during the intermittency route to thermoacoustic instability when the location of the flame is varied inside the combustor. As the flame location is changed, the synchronization properties of the coupled acoustic pressure and heat release rate signals change from desynchronized aperiodicity (combustion noise) to phase synchronized periodicity (thermoacoustic instability) through intermittent phase synchronization (intermittency). We also characterize the collective interaction between the multiple flamelets anchored at the flame holder and the acoustic field in the system, during different dynamical states observed in the combustor operation. When the signals are desynchronized, we notice that the flamelets exhibit a steady combustion without the exhibition of a prominent feedback with the acoustic field. In a state of intermittent phase synchronization, we observe the existence of a short-term coupling between the heat release rate and the acoustic field. We notice that the onset of collective synchronization in the oscillations of multiple flamelets and the acoustic field leads to the simultaneous emergence of periodicity in the global dynamics of the system. This collective periodicity in both the subsystems causes enhancement of oscillations during epochs of amplitude growth in the intermittency signal. On the contrary, the weakening of the coupling induces suppression of periodic oscillations during epochs of amplitude decay in the intermittency signal. During phase synchronization, we notice a sustained synchronized movement of all flamelets with the periodicity of the acoustic cycle in the system.     

Mechanism of combustion dynamics in a backward-facing step stabilized premixed flame

Combustion dynamics leading to thermoacoustic instability in a rearward-facing step stabilized premixed flame is experimentally examined with the objective of investigating the fluid dynamic mechanism that drives heat release rate fluctuations, and how it couples with the acoustic field. The field is probed visually, using linear photodiode arrays that capture the spatiotemporal distribution of CH* and OH*; an equivalence ratio monitor; and a number of pressure sensors. Results show resonance between the acoustic quarter wave mode of the combustion tunnel and a fluid dynamic mode of the wake. Under unstable conditions, the flame is convoluted around a large vortex that extends several step heights downstream. During a typical cycle, while the velocity is decreasing, the vortex grows, and the flame extends downstream around its outer edge. As the velocity reaches its minimum, becoming mostly negative, the vortex reaches its maximum size, and the flame collides with the upper wall; its leading edge folds, trapping reactants pockets, and its trailing edge propagates far upstream of the step. In the next phase, while the velocity is increasing, the heat release grows rapidly as trapped reactant’ pockets are consumed by flames converging towards their centers, and the upstream flame is dislodged back downstream. The heat release rate reaches its maximum halfway into the velocity rise period, leading the maximum velocity by about 90°. In this quarter-wave mode, the pressure leads the velocity by 90° as well, that is, it is in phase with the heat release rate. Numerical modeling results support this mechanism. Equivalence ratio contribution to the instability mechanism is shown to be minor, i.e., heat release dynamics are governed by the cyclical formation of the wake vortex and its interaction with the flame.     

Impact of symmetry breaking on the Flame Transfer Function of a laminar premixed flame

This work presents a numerical study of the acoustic response of a laminar flame with tunable asymmetry. A V-shaped premixed flame is stabilised in the wake of a cylindrical flame holder that can be rotated. The configuration is symmetric when the flame holder is fixed but increasing its rotation rate breaks the symmetry of the flow. This configuration is submitted to acoustic forcing to measure the effect of rotation of the flame holder on the Flame Transfer Functions. It appears that the asymmetry of the two flame branches changes their respective time delays, resulting in interference in the global unsteady heat release rate fluctuations. Consequently, the Flame Transfer Function exhibits dips and bumps, which are studied via laminar Direct Numerical Simulation. Potential applications for the control of combustion instabilities are discussed.     

Dynamics of spray and swirling flame under acoustic oscillations : A joint experimental and LES investigation

The dynamics of spray swirling flames is investigated by combining experiments on a single sector generic combustor and large eddy simulations of the same configuration. Measurements and calculations correspond to a self-sustained limit cycle operation where combustion coupled by an axial quarter wave acoustic mode induces large amplitude oscillations of pressure in the system. A detailed analysis of the mechanisms controlling the process is carried out first by comparing the measured and calculated spray and flame dynamics. Considering in a second stage that the spray and flame are compact with respect to the acoustic wavelength the analysis can be simplified by defining state variables that are obtained by taking averages over the combustor cross section and representing the behavior of these average quantities as a function of the axial coordinate and time. This reveals a first region in which essentially convective processes prevail. The convective heat release rate then couples further downstream with the pressure field giving rise to positive Rayleigh source terms which feed energy in the axial acoustic mode. In the convective region, the swirl number features oscillations around its mean value with an impact on the flow aerodynamics and flame radial displacement. Fluctuations in the fuel flow rate are initiated at the injector exhaust and likewise convected downstream. The total mass flow rate that exhibits strong convective disturbances is dominated further downstream by the acoustic motion. This information provides new insights on the convective-acoustic coupling that controls the heat release rate disturbances and reveals the time delays governing the combustion oscillation process.     

Premixed flame response to equivalence ratio perturbations

This paper studies the heat-release oscillation response of premixed flames to oscillations in reactant stream fuel/air ratio. Prior analyses have studied this problem in the linear regime and have shown that heat release dynamics are controlled by the superposition of three processes: flame speed, heat of reaction, and flame surface area oscillations. Each contribution has somewhat different dynamics, leading to complex frequency and mean fuel/air ratio dependencies. The present work extends these analyses to include stretch and non quasi-steady effects on the linear flame dynamics, as well as analysis of nonlinearities in flame response characteristics. Because the flame response is controlled by a superposition of multiple processes, each with a highly nonlinear dependence upon fuel/air ratio, the results are quite rich and the key nonlinearity mechanism varies with mean fuel/air ratio, frequency, and amplitude of excitation. In the quasi-steady framework, two key mechanisms leading to heat-release saturation have been identified. The first of these is the flame-kinematic mechanism, previously studied in the context of premixed flame response to flow oscillations and recently highlighted by Birbaud et al. (Combustion and Flame 154 (2008), 356–367). This mechanism arises due to fluctuations in flame position associated with the oscillations in flame speed. The second mechanism is due to the intrinsically nonlinear dependence of flame speed and mixture heat of reaction upon fuel/air ratio oscillations. This second mechanism is particularly dominant at perturbation amplitudes that cause the instantaneous stoichiometry to oscillate between lean and rich values, thereby causing non-monotonic variation of local flame speed and heat of reaction with equivalence ratio.     

Response of non-premixed flames to bulk flow perturbations

This paper describes the dynamics of non-premixed flames responding to bulk velocity fluctuations, and compares the dynamics of the flame sheet position and spatially integrated heat release to that of a premixed flame. The space–time dynamics of the non-premixed flame sheet in the fast chemistry limit is described by the stoichiometric mixture fraction surface, extracted from the solution of the

-equation. This procedure has some analogies to premixed flames, where the premixed flame sheet location is extracted from the G = 0 surface of the solution of the G-equation. A key difference between the premixed and non-premixed flame dynamics, however, is the fact that the non-premixed flame sheet dynamics are a function of the disturbance field everywhere, and not just at the reaction sheet, as in the premixed flame problem. A second key difference is that the non-premixed flame does not propagate and so flame wrinkles are convected downstream at the axial flow velocity, while wrinkles in premixed flames convect downstream at a vector sum of the flame speed and axial velocity. With the exception of the flame wrinkle propagation speed, however, we show that that the solutions for the space–time dynamics of the premixed and non-premixed reaction sheets in high velocity axial flows are quite similar. In contrast, there are important differences in their spatially integrated unsteady heat release dynamics. Premixed flame heat release fluctuations are dominated by area fluctuations, while non-premixed flames are dominated by mass burning rate fluctuations. At low Strouhal numbers, the resultant sensitivity of both flames to flow disturbances is the same, but the non-premixed flame response rolls off slower with frequency. Hence, this analysis suggests that non-premixed flames are more sensitive to flow perturbations than premixed flames at O(1) Strouhal numbers.     

Interaction of a Gaussian acoustic wave with a turbulent non-premixed flame

A Direct Numerical Simulation (DNS) of a turbulent non-premixed flame interacting with a Gaussian acoustic wave is carried out in this work. This numerical simulation takes into account detailed transport phenomena including the Soret effect as well as complete chemical kinetics on a very fine mesh. Turbulent non-premixed flame calculations are carried out both with and without an acoustic wave and results are recorded at the same time. By a simple difference it is then possible to obtain the influence of the acoustic wave/turbulent flame interaction. Using an extension of the non-linear Rayleigh criterion to a system with many species and elementary reactions, the obtained results can be further analysed. The initially planar acoustic wave develops strong perturbations along its transverse direction because of the interaction process, even at very early times. The amplitude of the pressure perturbation presents locally high positive as well as negative values, demonstrating the importance of focussing/defocussing effects and local amplification (resp. damping) phenomena. In the same way, the heat release rate is locally modified (either increased or decreased) during the interaction process. Finally, the presented Rayleigh criterion is used to identify regions where local amplification (respectively damping) takes place. Both amplification and damping zones coexist directly close to each other inside the reaction zone. The observed, resulting global effect is thus based on an average over highly varying local conditions within the flame front, leading to a smoothing effect. The complexity of the coupling procedure leading to this global wave amplification or damping is demonstrated by the present analysis.     

Impact of swirl and bluff-body on the transfer function of premixed flames

The frequency response of three lean methane/air flames submitted to flowrate perturbations is analyzed for flames featuring the same equivalence ratio and thermal power, but a different stabilization mechanism. The first flame is stabilized by a central bluff body without swirl, the second one by the same bluff body with the addition of swirl and the last one only by swirl without central insert. In the two last cases, the swirl level is roughly the same. These three flames feature different shapes and heat release distributions, but their Flame Transfer Function (FTF) feature about the same phase lag at low frequencies. The gain of the FTF also shows the same behavior for the flame stabilized by the central insert without swirl and the one fully aerodynamically stabilized by swirl. Shedding of vortical structures from the injector nozzle that grow and rollup the flame tip controls the FTF of these flames. The flame stabilized by the swirler-plus-bluff-body system features a peculiar response with a large drop of the FTF gain around a frequency at which large swirl number oscillations are observed. Velocity measurements in cold flow conditions reveal a strong reduction of the size of the vortical structures shed from the injector lip at this forcing condition. The flame stabilized aerodynamically only by swirl and the one stabilized by the bluff body without swirl do not exhibit any FTF gain drop at low frequencies. In the former case, large swirl number oscillations are still identified, but large vortical structures shed from the nozzle also persist at the same forcing frequency in the cold flow response. These different flame responses are found to be intimately related to the dynamics of the internal recirculation region, which response strongly differs depending upon the injector used to stabilize the flame.     

Nonlinear hydrodynamic and thermoacoustic oscillations of a bluff-body stabilised turbulent premixed flame

Turbulent premixed flames often experience thermoacoustic instabilities when the combustion heat release rate is in phase with acoustic pressure fluctuations. Linear methods often assume a priori that oscillations are periodic and occur at a dominant frequency with a fixed amplitude. Such assumptions are not made when using nonlinear analysis. When an oscillation is fully saturated, nonlinear analysis can serve as a useful avenue to reveal flame behaviour far more elaborate than period-one limit cycles, including quasi-periodicity and chaos in hydrodynamically or thermoacoustically self-excited system. In this paper, the behaviour of a bluff-body stabilised turbulent premixed propane/air flame in a model jet-engine afterburner configuration is investigated using computational fluid dynamics. For the frequencies of interest in this investigation, an unsteady Reynolds-averaged Navier–Stokes approach is found to be appropriate. Combustion is represented using a modified laminar flamelet approach with an algebraic closure for the flame surface density. The results are validated by comparison with existing experimental data and with large eddy simulation, and the observed self-excited oscillations in pressure and heat release are studied using methods derived from dynamical systems theory. A systematic analysis is carried out by increasing the equivalence ratio of the reactant stream supplied to the premixed flame. A strong variation in the global flame structure is observed. The flame exhibits a self-excited hydrodynamic oscillation at low equivalence ratios, becomes steady as the equivalence ratio is increased to intermediate values, and again exhibits a self-excited thermoacoustic oscillation at higher equivalence ratios. Rich nonlinear behaviour is observed and the investigation demonstrates that turbulent premixed flames can exhibit complex dynamical behaviour including quasiperiodicity, limit cycles and period-two limit cycles due to the interactions of various physical mechanisms. This has implications in selecting the operating conditions for such flames and for devising proper control strategies for the avoidance of thermoacoustic instability.     

Dynamic response of turbulent swirling flames to acoustic perturbations

The dynamic response of a turbulent, perfectly premixed flame, stabilized by means of an aerodynamic flameholder, to an upstream acoustic perturbation of the approaching flow is investigated by means of experimental and analytical tools, and simulated through a large eddy simulation of the reacting flow. It is found that the main contribution to the unsteady heat release rate is due to the fluctuation in area of the flame front, which in turn is strongly influenced by the corresponding response of the flow field to the acoustic perturbation. Numerical data show that perturbing a swirling flow that undergoes vortex breakdown results in a strong displacement of the breakdown position along its axis, while its outer part only weakly responds to the perturbation. This results in a translational motion of the flame's anchoring point, which ultimately leads to an unsteady variation of the flame area and, therefore, of the amount of heat released. This unsteady heat release mechanism can be described in a way similar to that used for characterizing the dynamic behaviour of ducted flames, stabilized by means of a bluff-body flameholder; differently from these models, however, the anchoring point of the flame can now fluctuate freely in space, and the time delay of the system is no longer identified with the travelling time of a perturbation of the flame element along it, but is now associated with the oscillation of the breakdown position. Controlling the interaction between breakdown and acoustics should allow for obtaining optimal flame dynamics, so as to limit and possibly avoid the occurrence of strong pulsation peaks whenever the combustion device is operated in an acoustically closed system.     

Combustion driven oscillations

Combustion driven oscillations can occur when a turbulent flame is enclosed in a tube or cavity. Interaction between heat fluctuations and the internal standing wave field at one of the natural frequencies of the air column produces strong organ pipe tones. The sound power emitted by this thermal-acoustic interaction depends on the impedance either side of the combustion zone and on a transfer function defining the response of the flame to sound wave disturbances. If this power exceeds the rate at which energy is dissipated at the cavity boundaries then there is a growth of the internal pressure field and an increase in the radiated sound. Plane wave theory is used to calculate the flame transfer function and adjacent impedances for a simple gas fired tube assembly. The predicted instability frequencies are then compared with experimental data. The results indicate that the flame transfer function plays a dominant role in determining the acoustic stability of the cavity and that insufficient data is available for accurately predicting unsteady flame front behaviour.     

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

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

Computational investigations of the coupling between transient flame dynamics and thermo-acoustic instability in a self-excited resonance combustor

Combustion instability due to thermo-acoustic interactions is a critical combustion problem that requires a thorough understanding because of its adverse impact on stable and reliable operation of combustors in high-speed propulsion devices like gas turbines and rockets. This work conducts computational investigations of the coupling between the transient flame dynamics such as the ignition delay and local extinction and the thermo-acoustic instability developed in a self-excited resonance combustor to gain deep insights into the mechanisms of thermo-acoustic instability. A 2D modelling framework that employs different flamelet models (the steady flamelet model and the flamelet/progress variable approach) is developed to enable the examination of the effect of the transient flame dynamics caused by the strong coupling of the turbulent mixing and finite-rate chemical kinetics on the occurrence of thermo-acoustic instability. The models are validated by using the available experimental data for the pressure signal. Parametric studies are performed to examine the effect of the occurrence of the transient flame dynamics, the effect of artificial amplification of the Damköhler number, and the effect of neglecting mixture fraction fluctuations on the predictions of the thermo-acoustic instability. The parametric studies reveal that the occurrence of transient flame dynamics has a strong influence on the onset of the thermo-acoustic instability. Further analysis is then conducted to localise the effect of a particular flame dynamic event, the ignition delay, on the thermo-acoustic instability. The reverse effect of the occurrence of the thermo-acoustic instability on the transient flame dynamics in the combustor is also investigated by examining the temporal evolution of the local flame events in conjunction with the pressure wave propagation. The above observed two-way coupling between the transient flame dynamics (the ignition delay) and the thermo-acoustic instability provides a plausible mechanism of the self-excited and sustained thermo-acoustic instability observed in the combustor despite the fact that the results are obtained from 2D simulations. The same analysis is expected to be extensible to fully 3D simulations.     

Turbulent premixed hydrogen/air flame-wall interaction with heterogeneous surface reactions

A premixed ${\mathrm{H}}_2$/air flame impinging onto a flat surface in statistically stationary state is studied for both reactive and inert wall cases to gain insights into the effects of the heterogeneous surface reactions. Direct numerical simulation (DNS) results with detailed gas-phase chemistry and surface adsorption and desorption mechanisms indicate differences in the flame front topology and near-wall flame dynamics between these two cases. In the reactive surface case, gas-phase free radicals are inclined to be adsorbed with much reduced near-wall concentration. Consequently, the gas-phase heat release rate (HRR) close to the wall decreases as well because of the low availability of free radicals. However, extra heat released from the reactive surface partially compensates for such difference. Moreover, wall reactions will intensify the turbulent fluctuations of the wall temperature and wall heat flux, while for these two cases the mean temperature profile along the flame propagating direction remains similar.     

Interplay of noise and coupling in heterogeneous ensembles of phase oscillators

We study the effects of noise on the collective dynamics of an ensemble of coupled phase oscillators whose natural frequencies are all identical, but whose coupling strengths are not the same all over the ensemble. The intensity of noise can also be heterogeneous, representing diversity in the individual responses to external fluctuations. We show that the desynchronization transition induced by noise may be completely suppressed, even for arbitrarily large noise intensities, is the distribution of coupling strengths decays slowly enough for large couplings. Equivalently, if the response to noise of a sufficiently large fraction of the ensemble is weak enough, desynchronization cannot occur. The two effects combine with each other when the response to noise and the coupling strength of each oscillator are correlated. This combination is quantitatively characterized and illustrated with explicit examples.