Writing script to model combustion equilibrium of hydrogen-air flame and reporting major/minor species equilibriums along with associated adiabatic flame temperatures

• Flame temperature & equilibrium concentrations as a function of equivalence ratio​

• Effect of increase in pressure

• Effect of calorically perfect gas assumption

• Effect of adding neon diluent

• Comparison of NOx emission concentration based on all above changes

Comparison and validation of written script to available commercial programs (GasEq & NASA CEA)

Equilibrium calculation script required use of multiple iteration methods (bisection method and solve/linsolve)

• Developed code included unique numerical integration method to account for changing specific heats of gasses over temperature range​

Written script closely matches commercial code calculations for adiabatic flame temperature

Major species concentration calculations all fell within 5% of commercial code outputs

Minor species concentrations all within same order of magnitude

See below for full report of code characteristics and algorithms