Sobreestimación de la temperatura de la llama de adiabación (gas hidrógeno)

Estoy trabajando en un código que calcula la temperatura de llama adiabática dada una ecuación balanceada y las propiedades termodinámicas relevantes y estoy comenzando con la reacción de combustión más simple que se me ocurre: combustión estequiométrica de hidrógeno puro con oxígeno puro.

$2H_2+O_2 \Rightarrow 2H_2O$

Estoy usando estos valores termodinámicos:

$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$kJ/kg-K a 300K para$H_2O_{(g)}$

$c_p = 3.217$kJ/kg-K a 4000K para$H_2O_{(g)}$

Si la ecuación de combustión que escribí está basada en kmol, obtengo un cambio total en la entalpía de 483640 kJ con la combustión. Si ese calor se usa para calentar el producto (vapor de agua), obtengo aumentos de temperatura que van desde 4173K a 7201K (temperaturas finales de 4471K y 7499K) según el valor de calor específico que use.

$\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}$

Cuando ejecuto el código que realmente considera los cambios en el calor específico durante todo el proceso de calentamiento (interpolando según una tabla en cada paso), obtengo 5024K. Incluso mi valor más bajo, suponiendo un calor específico constante en el valor más alto, da un aumento de temperatura superior al número que figura en Wikipedia, ~3500K. ¿Estoy cometiendo algún tipo de error o hay algún tipo de factor que hace que lo real sea diferente de la idealización? Gracias.