Primal-Dual Interior-Point Method for an Optimization Problem Related to the Modeling of Atmospheric Organic Aerosols

Abstract

A mathematical model for the computation of the phase equilibrium related to atmospheric organic aerosols is presented. The phase equilibrium is given by the global minimum of the Gibbs free energy for a system that involves water and organic components. This minimization problem is equivalent to the determination of the convex hull of the corresponding molar Gibbs free energy function. A geometrical notion of phase simplex related to the convex hull is introduced to characterize mathematically the phases at equilibrium. A primal-dual interior-point algorithm for the efficient solution of the phase equilibrium problem is presented. A novel initialization of the algorithm, based on the properties of the phase simplex, is proposed to ensure the convergence to a global minimum of the Gibbs free energy. For a finite termination of the interior-point method, an active phase identification procedure is incorporated. Numerical results show the robustness and efficiency of the approach for the prediction of liquid-liquid equilibrium in multicomponent mixtures.