A method for computing multicomponent chemical equilibria based on equilibrium constants

The equilibrium activities and concentrations of N chemical species in a multicomponent system may be calculated given M independent equilibrium constants relating these species and N-M mass balances. These N equations are solved by Newton-Raphson iteration. Where initial concentration estimates are poor or the number of species large, an additional curve crawler technique introduces components by stepwise infinitesimal increments.