Abstract A can-annular combustor consists of a set of nominally identical cans, in which the flames burn in an essentially isolated manner. However, adjacent cans are able to communicate acoustically, which provides dynamic coupling of the entire can-annular arrangement. Recently, it was shown that the acoustic coupling is not negligible and can cause clustering of eigenfrequencies. In this study, we present a low-order modeling framework for self-excited thermoacoustic oscillations in generic can-annular combustors consisting of N identical cans. The dynamics of the flames are modeled with the nonlinear G-equation; the acoustic model accounts for plane acoustic waves inside the cans and can-to-can communication. The latter is enabled through a coupling boundary condition that is based on conservation of mass and a Rayleigh conductivity. For weak coupling between adjacent cans, the thermoacoustic feedback cycle shows clusters of linearly unstable modes of different azimuthal order, which are close in frequency and growth rate. Their interaction in the nonlinear regime is investigated using time-domain simulations. Two simulations for generic can-annular combustors consisting of 4 and 6 cans with weak acoustic coupling are discussed in this study. We observe a strong interaction between the modes, which can cause long transition times and allows modes that do not dominate the system dynamics in the linear regime to be dominant in the nonlinear regime. While the N = 6 case converges to a periodic oscillation pattern with one dominant frequency, the N = 4 case converges to a quasi-periodic oscillation involving modes of different azimuthal order. Moreover, we observe a synchronization of these modes. These results raise the questions whether it is possible to predict which mode(s) will dominate the system in the saturated state and under which conditions synchronization of clustered modes can occur.