Phase equilibrium calculations with quasi-Newton methods

被引:21
作者
Nichita, Dan Vladimir [1 ]
Petitfrere, Martin [1 ,2 ]
机构
[1] Univ Pau & Pays Adour, CNRS, UMR 5150, Lab Fluides Complexes & Leurs Reservoirs, F-64013 Pau, France
[2] Ctr Sci & Tech Jean Feger, F-64018 Pau, France
关键词
Phase split; Phase stability; Convergence; Scaling; Condition number; Quasi-Newton method; ISOTHERMAL FLASH PROBLEM; EQUATION-OF-STATE; BOUND-CONSTRAINED OPTIMIZATION; LIQUID-EQUILIBRIA; REDUCTION METHOD; SUCCESSIVE SUBSTITUTION; MULTICOMPONENT MIXTURES; CONVERGENCE CONDITIONS; STABILITY ANALYSIS; CUBIC EQUATIONS;
D O I
10.1016/j.fluid.2015.07.035
中图分类号
O414.1 [热力学];
学科分类号
摘要
The phase split problem, formulated as an unconstrained minimization of the Gibbs free energy, is commonly solved by the second-order Newton method, preceded by a number of first-order successive substitutions. For difficult problems, the convergence radius of the Newton method may be small and a high number of successive substitution iterations may be required before the switch, or in repeated switch-backs. An interesting alternative is given by the quasi-Newton methods, representing a good compromise between complexity and convergence speed. The quasi-Newton BFGS method exhibits a super-linear convergence rate (in some cases without step length control) and a rank two update of the Hessian matrix guarantees a hereditary positive definiteness. In this work, a scaling methodology is proposed for finding the appropriate change of variables for phase equilibrium problems; applied to the two-phase split problem, the resulting change of variables leads to a Hessian matrix of the form H = I + D + ND, where I is the identity matrix, D is a diagonal matrix with elements vanishing at the solution, and ND is an effective low-rank matrix. The results of numerical experiments carried out on several test cases show that the BFGS method using the proposed variables is more robust and efficient than previous implementations (from the literature and open source codes). A two-parameter cubic equation of state was used in this work, but any equation of state can be used. The quasi-Newton methods are particularly suited for thermodynamic models for which the Hessian matrix is difficult or costly to obtain. (C) 2015 Elsevier B.V. All rights reserved.
引用
收藏
页码:194 / 208
页数:15
相关论文
共 50 条
[1]   THE ISOTHERMAL FLASH PROBLEM - NEW METHODS FOR PHASE SPLIT CALCULATIONS [J].
AMMAR, MN ;
RENON, H .
AICHE JOURNAL, 1987, 33 (06) :926-939
[2]  
[Anonymous], SPE ANN TECHN C EXH
[4]   A LIMITED MEMORY ALGORITHM FOR BOUND CONSTRAINED OPTIMIZATION [J].
BYRD, RH ;
LU, PH ;
NOCEDAL, J ;
ZHU, CY .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1995, 16 (05) :1190-1208
[5]   QUASI-NEWTON METHODS, MOTIVATION AND THEORY [J].
DENNIS, JE ;
MORE, JJ .
SIAM REVIEW, 1977, 19 (01) :46-89
[6]  
FLETCHER R, 1994, NATO ADV SCI INST SE, V434, P109
[7]  
Fletcher R., 2000, PRACTICAL METHODS OP
[8]   Vapor-liquid equilibria of nitrogen-hydrocarbon systems using the PC-SAFT equation of state [J].
García-Sánchez, F ;
Eliosa-Jiménez, G ;
Silva-Oliver, G ;
Vázquez-Román, R .
FLUID PHASE EQUILIBRIA, 2004, 217 (02) :241-253
[9]   Modeling of microemulsion phase diagrams from excess Gibbs energy models [J].
García-Sánchez, F ;
Eliosa-Jiménez, G ;
Salas-Padrón, A ;
Hernández-Garduza, O ;
Apam-Martínez, D .
CHEMICAL ENGINEERING JOURNAL, 2001, 84 (03) :257-274
[10]   Modeling of multiphase liquid equilibria for multicomponent mixtures [J].
GarciaSanchez, F ;
Schwartzentruber, J ;
Ammar, MN ;
Renon, H .
FLUID PHASE EQUILIBRIA, 1996, 121 (1-2) :207-225