We present a numerical method for chemical equilibrium calculations based on a Gibbs energy minimization approach.
The method is capable of determining the stable phases at equilibrium without techniques that are dependent on their types (e.g., pure minerals).
The chemical equilibrium algorithm is applied in reactive transport simulations to demonstrate its effective use in intensive applications.
The equilibrium calculations are shown to converge within a few iterations (1–5) with a warm-start approach in reactive transport simulations.