MaximaLabs
All unit operations

Gibbs reactor simulation

1+ inlets -> 1 outlet (free-energy min)

Governing equations

The exact equations the solver works for a gibbs reactor — the same math shown in the app's "Theory" panel, not a black box.

minni GRT=ini[ΔGf,i(T)RT+lnnintot+lnPP]\min_{n_i}\ \frac{G}{RT} = \sum_i n_i\left[\frac{\Delta G^\circ_{f,i}(T)}{RT} + \ln\frac{n_i}{n_{tot}} + \ln\frac{P}{P^\circ}\right]
s.t.iniaij=bj (atoms of element j),ni0\text{s.t.}\quad \sum_i n_i\,a_{ij} = b_j \ (\text{atoms of element } j),\quad n_i \ge 0
GG
total Gibbs free energy of the mixture [J] — minimized to find equilibrium
nin_i
moles of species i [mol/s] (the unknowns)
ntotn_{tot}
total moles [mol/s]
ΔGf,i(T)\Delta G^\circ_{f,i}(T)
standard Gibbs energy of formation of species i at T [J/mol]
PP
pressure [Pa]
PP^\circ
reference pressure (1 atm)
aija_{ij}
atoms of element j in one molecule of species i
bjb_j
total atoms of element j (conserved from the feed)
RR
gas constant

Parameters

temperature [K, required]; optional product_species [list of candidate product component ids the reactor may form, e.g. ['co','co2','h2'] for CH4/H2O reforming], pressure_drop [Pa]. Finds the equilibrium composition by Gibbs free-energy minimization under element conservation (Aspen RGibbs equivalent) — no reaction network needed. Every considered species must be in the curated cited formation set (H2/N2/O2/H2O/CO2/CO/CH4/ethane/propane/ethylene/propylene/acetylene/methanol/ethanol/NH3/H2S/SO2/NO/NO2/benzene/toluene). Isothermal

Example flowsheets that use it

Stop fighting legacy software. Build your first flowsheet in 60 seconds.