This paper presents a computer model for studying the processes of ignition and combustion of air suspended aluminum particles, based on a system of differential equations of the gas phase balance and the probability density function for temperature and particle size. The model is implemented as software for a two-dimensional non-steady state case. A new implicit difference scheme consisting in two - and three-point approximations for the first and second derivatives is proposed and implemented. A set of calculations is performed and the results of the reacting mixture characteristics distribution, the areas of the maximum heat release and their relationship with the most air suspension probable state are presented.