Abstract. The reduction of NOx-formation in biomass fired grate furnaces requires the development of numerical models. To represent the variety in scales and physical processes playing a role in the conversion, newly developed submodels are required. Here, a submodel for the reverse combustion process in the solid fuel layer on the grate is described. The submodel is shown to give good predictions for the velocity of the combustion front as well as for the spatial profiles of porosity, oxygen mass fraction and temperature. These predictions are essential input for NOx-calculations.