  Thermoacoustic instabilities observed in turbulent combustion systems have
disastrous consequences and are notoriously challenging to model, predict and
control. Here, we introduce a mean-field model of thermoacoustic transitions,
where the nonlinear flame response is modeled as the amplitude weighted
response of an ensemble of phase oscillators constrained to collectively evolve
at the rhythm of acoustic fluctuations. Starting from the acoustic wave
equation coupled with the phase oscillators, we derive the evolution equations
for the amplitude and phase and obtain the limit cycle solution. We show that
the model captures abrupt and continuous transition to thermoacoustic
instability observed in disparate combustors. We obtain quantitative insights
into the model by estimating the model parameters from the experimental data
using parameter optimisation. Importantly, our approach provides an explanation
of spatiotemporal synchronization and pattern-formation underlying the
transition to thermoacoustic instability while encapsulating the statistical
properties of desynchronization, chimeras, and global phase synchronization. We
further show using the model that continuous and abrupt transitions to limit
cycle oscillations in turbulent combustors corresponds to synchronization
transitions of \textit{second-order} and \textit{first-order}, respectively.
The present formulation provides a highly interpretable model of thermoacoustic
transitions: changes in empirical bifurcation parameters which lead to limit
cycle oscillations amounts to an increase in the coupling strength of the phase
oscillators, promoting global phase synchronization. The generality of the
model in capturing different types of transitions and states of
pattern-formation highlights the possibility of extending the present model to
a broad range of fluid-dynamical phenomena beyond thermoacoustics.
