Usage

from 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