The paper deals with development and application the numerical model for solution of processes at combustion chamber of the thermal power plant boiler. Mathematical simulation is based on solution of physical and chemical processes occuring at burning pulverized coal in the furnace model. Three-dimensional flows, heat and mass transfer, chemical kinetics of the processes, effects of thermal radiation are considered. Obtained results give quantitative information on velocity distributions, temperature and concentration profiles of the components, the amount of combustion products including harmful substances. The numerical model becomes a tool for investigation and design of combustion chambers with high-efficiency and reliable operation of boiler at thermal power plants.