Iron powder can serve as an energy carrier and help generate heat and power in a clean and sustainable manner. In this paper, we analyze one of the challenges related to iron combustion: the mode of flame propagation. Depending on the conditions, an iron dust flame can propagate in either the continuous mode, with neighboring particles burning simultaneously, or in a discrete mode, with each particle burning separately, i.e. as a percolating wave. We use particle-resolved simulations including the heterogeneous reaction of iron with oxygen and detailed transport models. Simulations are carried out with multiple particles in a single row with hot gas at the inlet for initiating combustion under weak convection effects. Continuous and discrete flame propagation modes are identified and the particle distance, at which the transition from the continuous to the discrete propagation mode occurs, is calculated for different particle diameters. It is found that the transition from one mode to another occurs close to stoichiometry and the maximum flame reaction front speed is within the range of experiments with similarly sized particles. The flame propagation mode transitions to the discrete mode for larger particles at slightly lower particle distances than smaller particles due to higher thermal inertia. The discreteness parameter is modified to take into account the effect of heating from burnt particles on the unburnt ones. This results in an overall adequate prediction of the numerically found transitional distances under weak convection effects.Novelty and significance: As a carbon-free energy carrier, there is a great requirement for exploring the fundamental and application aspects of iron combustion. One critical issue for the heterogeneous combustion of iron is the mode of flame propagation. While previous studies have relied on simplified models and simplifying assumptions, in this work, we carry out detailed boundary layer resolved simulations of multiple particles with varying distances between the particles. Beyond a critical particle spacing, the mode of flame propagation transitions from the continuous mode where multiple particles burn simultaneously, to the discrete mode where the particles burn standalone. The scaling of such a transition is improved with the help of appropriate temperature based properties and the length scale. Our work is significant in understanding the fundamentals of flame propagation in iron clusters and can help in the development of Euler–Lagrangian models where tens of thousands of particles need to be simulated.