Завышение температуры адиабационного пламени (газообразный водород)

Я работаю над кодом, который вычисляет адиабатическую температуру пламени с учетом сбалансированного уравнения и соответствующих термодинамических свойств, и я начинаю с самой простой реакции горения, которую только могу придумать: стехиометрического сжигания чистого водорода с чистым кислородом.

$2H_2+O_2 \Rightarrow 2H_2O$

Я использую эти термодинамические значения:

$h_{f,H_2} = 0.0 \text{ kJ/kmol}$

$h_{f,O_2} = 0.0 \text{ kJ/kmol}$

$h_{f,H_2O} = -241820 \text{ kJ/kmol}$

$c_p = 1.864$кДж/кг-К при 300К для$H_2O_{(g)}$

$c_p = 3.217$кДж/кг-К при 4000К для$H_2O_{(g)}$

Если уравнение сгорания, которое я написал, основано на кмоль, я получаю общее изменение энтальпии 483640 кДж при сгорании. Если это тепло идет на нагрев продукта (водяного пара), я получаю повышение температуры в диапазоне от 4173K до 7201K (конечные температуры 4471K и 7499K) в зависимости от используемой вами удельной теплотворной способности.

$\frac{483640\text{ kJ}}{(3.217 \text{ kJ/kg-K})(36.03\text{ kg})} = 4173\text{ K}$

$\frac{483640\text{ kJ}}{(1.864 \text{ kJ/kg-K})(36.03\text{ kg})} = 7201\text{ K}$

Когда я запускаю код, который фактически учитывает изменения удельной теплоемкости на протяжении всего процесса нагрева (интерполяция на основе таблицы на каждом шаге), я получаю 5024К. Даже мое самое низкое значение, предполагающее постоянную удельную теплоемкость при самом высоком значении, дает повышение температуры выше, чем число, указанное в Википедии, ~ 3500K. Я делаю какую-то ошибку или есть какой-то фактор, который отличает реальную вещь от идеализации? Спасибо.