Thank you John. I followed your advice and got more reasonable change in enthalpies.
But now I need to simulate the combustor in the gas turbine. None of the equilibrate() combinations below adequately simulates the combustion in the combustor. since temperature, entropy, internal energy and volume are not constant. The equilibrate() holds the pair in the bracket constant.
gas1.TP
gas1.TD
gas1.HP
gas1.UV
gas1.SP
gas1.SV
If I equilibrate (HP) since the combustion in a gas turbine is at constant enthalpy and constant pressure to get the adiabetic flame temperature and the equilibrate again at constant pressure and constant temperature at the adiabetic temperature to obtain the new enthalpy. Will this be a good representation of the combustion in a gas turbine combustor?
I have attached the updated model in Jupyter Notebook.
Regards,
Kofi.