Предлагается вычислительная схема для совместного численного решения уравнений энергии, переноса энергии излучения и уравнений модели горения. Для получения дискретных аналогов уравнений энергии и модели горения используется метод сплайнколлокации в сочетании с методом конечных разностей. Радиационный теплообмен рассматривается в рамках P1 – приближения метода сферических гармоник. Учитывается переменность теплофизических свойств продуктов сгорания.