The incineration of municipal solid waste (MSW) contributes significantly to the presence of heavy metals in urban area aerosols. It is thus important to ascertain the quantities and chemical forms of the heavy metals (HM) that are emitted from the incineration plant stacks. The behaviour of HM, which depends strongly on the thermal and chemical environments, was investigated herein with a modelling approach, consisting of several parts. First, a refuse bed combustion model was developed for simulating on-grate MSW incineration. It describes most of the physico-chemical and thermal phenomena occurring during, waste combustion. Second, results from the bed model were taken as boundary conditions to perform 3D simulations of the post-combustion zone and of the boiler. The case studied was of the Strasbourg incineration plant. Finally, the local thermal conditions and the local elementary compositions of gas and solid phases obtained from these simulations were used to carry out thermodynamic calculations of the speciation of HM at each point in the incinerator. The results for four metals (Cd, Zn, Pb, Cr) are presented, discussed and compared to available data. Predicted species are in agreement with observations for volatile metals, except lead, whose volatilization seems overestimated.