Abstract
In this work, we present a method for the calculation of two-phase equilibria in multireaction systems. The procedure uses a transformed composition variable in the orthogonal derivatives of the Gibbs function and the tangent plane equation to form a system of non-linear equations. We solve this system with a Newton-Raphson method and our initialization procedure uses results from the reactive stability analysis and the reactive equal area rule. With this initialization strategy, convergence occurs with only a few iterations in the numerical method. Several examples with multiple chemical reactions demonstrate the performance of our approach.
Users
Please
log in to take part in the discussion (add own reviews or comments).