Abstract:Abstract We present an overview of novel numerical methods for chemical equilibrium and kinetic calculations for complex non-ideal multiphase systems. The methods we present for equilibrium calculations are based either on Gibbs energy minimization (GEM) calculations or on solving the system of extended law of mass-action (xLMA) equations. In both methods, no a posteriori phase stability tests, and thus no tentative addition or removal of phases during or at the end of the calculations, are necessary. All potentially stable phases are considered from the beginning of the calculation, and stability indices are immediately available at the end of the computation to determine which phases are actually stable at equilibrium. Both GEM and xLMA equilibrium methods are tailored for computationally demanding applications that require many rapid local equilibrium calculations, such as reactive transport modeling. The numerical method for chemical kinetic calculations we present supports both closed and open systems, and it considers a partial equilibrium simplification for fast reactions. The method employs an implicit integration scheme that improves stability and speed when solving the often stiff differential equations in kinetic calculations. As such, it requires compositional derivatives of the reaction rates to assemble the Jacobian matrix of the resultant implicit algebraic equations that are solved at every time step. We present a detailed procedure to calculate these derivatives, and we show how the partial equilibrium assumption affects their computation. These numerical methods have been implemented in Reaktoro (reaktoro.org), an open-source software for modeling chemically reactive systems. We finish with a discussion on the comparison of these methods with others in the literature.

An Algorithm for Simultaneous Chemical and Phase-Equilibrium Calculation

A Modified NRTL Equation for the Calculation of Phase Equilibrium of Polymer Solutions

Measurement of Phase Diagrams for New Aqueous Two-Phase Systems and Prediction by A Generalized Multicomponent Osmotic Virial Equation

A Simultaneous Method for Two‐ and Three‐liquid–phase Stability Determination

Simultaneous Predictions of Chemical and Phase Equilibria in Systems with an Esterification Reaction Using PC-SAFT

Phase Field Simulation for Non-Isothermal Solidification of Multicomponent Alloys Coupled with Thermodynamics Database

One-step Algorithm for Multiphase Equilibrium Calculation Using State Equations

Calculation of Phase Equilibria Based on the Levenberg-Marquardt Method

Convex Envelope Method for determining liquid multi-phase equilibria in systems with arbitrary number of components

Equation-Oriented Approach for Handling the Perturbed-Chain SAFT Equation of State in Simulation and Optimization of Polymerization Processes

Elemental Composition Method for Computation and Analysis of Simultaneous Chemical and Phase Equilibrium

A study on SAFT-CP for phase equilibria of fluid mixtures with a component near its critical point

A New Algorithm for LLV Equilibrium Calculation

Unified approach for distillation column simulation and optimization based on equation-oriented algorithms

An initial value insensitive method for phase equilibrium calculation: Constrained backtracking search optimization algorithm

The Concept of the Estimation of Phase Diagrams (An Optimised Set of Simplified Equations to Estimate Equilibrium Liquidus and Solidus Temperatures, Partition Ratios, and Liquidus Slopes for Quick Access to Equilibrium Data in Solidification Software) Part I: Binary Equilibrium Phase Diagrams

A simple chemical equilibrium algorithm applied for single and multiple reaction systems

Python and MATLAB implementations of a fast numerical method for the solution of multiple chemical equilibria

Phase Equilibrium Calculation of Mixtures: Comparison of SAFT, Modified SAFT, and BACK EOS for Supercritical CO2-C2H5OH System

Surrogate Equations of State for Equation Oriented Optimization of Polymerization Processes

An overview of computational methods for chemical equilibrium and kinetic calculations for geochemical and reactive transport modeling