Usage 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21from property_utils.units import KILO_JOULE, MOL from chemical_utils.substances import METHANE, WATER, CARBON_MONOXIDE, HYDROGEN2 from chemical_utils.reactions import r, WATER_GAS_SHIFT reforming = r(METHANE + WATER, CARBON_MONOXIDE + 3 * HYDROGEN2) total_reaction = reforming + WATER_GAS_SHIFT print(total_reaction) print( "Standard enthalpy change: ", total_reaction.standard_enthalpy_change.to_unit(KILO_JOULE / MOL), ) water_gas_shift = total_reaction - reforming print(water_gas_shift) print( "Standard Gibbs energy change: ", water_gas_shift.standard_gibbs_energy_change.to_unit(KILO_JOULE / MOL), ) Result: 2H2O + CH4 + CO -> CO + 4H2 + CO2 Standard enthalpy change: 164.638 kJ / mol H2O + CO -> H2 + CO2 Standard Gibbs energy change: -28.63 kJ / mol