Comparison of reduced and conventional two-phase flash calculations

Seyhan Emre Gorucu, Russell Taylor Johns

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

Phase-equilibrium calculations become computationally intensive in compositional simulation as the number of components and phases increases. Reduced methods were developed to address this problem, where the binary-interaction-parameter (BIP) matrix is approximated either by spectral decomposition (SD), as performed by Hendriks and van Bergen (1992), or with the twoparameter BIP formula of Li and Johns (2006). Several authors have recently stated that the SD method-And by reference all reduced methods-is not as fast as previously reported in the literature. In this paper we present the first study that compares all eight reduced and conventional methods published to date by use of optimized code and compilers. The results show that the SD method and its variants are not as fast as other reduced methods, and can be slower than the conventional approach when fewer than 10 components are used. These conclusions confirm the findings of recently published papers. The reason for the slow speed is the requirement that the code must allow for a variable number of eigenvalues. We show that the reduced method of Li and Johns (2006) and its variants, however, are faster because the number of reduced parameters is fixed to six, which is independent of the number of components. Speed up in flash calculations for their formula is achieved for all fluids studied when more than six components are used. For example, for 10-component fluids, a speed up of 2-3 in the computational time for Newton-Raphson (NR) iterations is obtained compared with the conventional method modeled after minimization of Gibbs energy. The reduced method modeled after the linearized approach of Nichita and Graciaa (2011), which uses the twoparameter BIP formula of Li and Johns (2006), is also demonstrated to have a significantly larger radius of convergence than other reduced and conventional methods for five fluids studied.

Original languageEnglish (US)
Pages (from-to)294-305
Number of pages12
JournalSPE Journal
Volume20
Issue number2
DOIs
StatePublished - Apr 1 2015

Fingerprint

Decomposition
Fluids
Gibbs free energy
Phase equilibria
decomposition
fluid
comparison
calculation
method
eigenvalue
phase equilibrium
matrix
parameter
simulation
energy
speed

All Science Journal Classification (ASJC) codes

  • Energy Engineering and Power Technology
  • Geotechnical Engineering and Engineering Geology

Cite this

@article{445dc261e7094229b7de05246e0931d0,
title = "Comparison of reduced and conventional two-phase flash calculations",
abstract = "Phase-equilibrium calculations become computationally intensive in compositional simulation as the number of components and phases increases. Reduced methods were developed to address this problem, where the binary-interaction-parameter (BIP) matrix is approximated either by spectral decomposition (SD), as performed by Hendriks and van Bergen (1992), or with the twoparameter BIP formula of Li and Johns (2006). Several authors have recently stated that the SD method-And by reference all reduced methods-is not as fast as previously reported in the literature. In this paper we present the first study that compares all eight reduced and conventional methods published to date by use of optimized code and compilers. The results show that the SD method and its variants are not as fast as other reduced methods, and can be slower than the conventional approach when fewer than 10 components are used. These conclusions confirm the findings of recently published papers. The reason for the slow speed is the requirement that the code must allow for a variable number of eigenvalues. We show that the reduced method of Li and Johns (2006) and its variants, however, are faster because the number of reduced parameters is fixed to six, which is independent of the number of components. Speed up in flash calculations for their formula is achieved for all fluids studied when more than six components are used. For example, for 10-component fluids, a speed up of 2-3 in the computational time for Newton-Raphson (NR) iterations is obtained compared with the conventional method modeled after minimization of Gibbs energy. The reduced method modeled after the linearized approach of Nichita and Graciaa (2011), which uses the twoparameter BIP formula of Li and Johns (2006), is also demonstrated to have a significantly larger radius of convergence than other reduced and conventional methods for five fluids studied.",
author = "Gorucu, {Seyhan Emre} and Johns, {Russell Taylor}",
year = "2015",
month = "4",
day = "1",
doi = "10.2118/163577-PA",
language = "English (US)",
volume = "20",
pages = "294--305",
journal = "SPE Journal",
issn = "1086-055X",
publisher = "Society of Petroleum Engineers (SPE)",
number = "2",

}

Comparison of reduced and conventional two-phase flash calculations. / Gorucu, Seyhan Emre; Johns, Russell Taylor.

In: SPE Journal, Vol. 20, No. 2, 01.04.2015, p. 294-305.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Comparison of reduced and conventional two-phase flash calculations

AU - Gorucu, Seyhan Emre

AU - Johns, Russell Taylor

PY - 2015/4/1

Y1 - 2015/4/1

N2 - Phase-equilibrium calculations become computationally intensive in compositional simulation as the number of components and phases increases. Reduced methods were developed to address this problem, where the binary-interaction-parameter (BIP) matrix is approximated either by spectral decomposition (SD), as performed by Hendriks and van Bergen (1992), or with the twoparameter BIP formula of Li and Johns (2006). Several authors have recently stated that the SD method-And by reference all reduced methods-is not as fast as previously reported in the literature. In this paper we present the first study that compares all eight reduced and conventional methods published to date by use of optimized code and compilers. The results show that the SD method and its variants are not as fast as other reduced methods, and can be slower than the conventional approach when fewer than 10 components are used. These conclusions confirm the findings of recently published papers. The reason for the slow speed is the requirement that the code must allow for a variable number of eigenvalues. We show that the reduced method of Li and Johns (2006) and its variants, however, are faster because the number of reduced parameters is fixed to six, which is independent of the number of components. Speed up in flash calculations for their formula is achieved for all fluids studied when more than six components are used. For example, for 10-component fluids, a speed up of 2-3 in the computational time for Newton-Raphson (NR) iterations is obtained compared with the conventional method modeled after minimization of Gibbs energy. The reduced method modeled after the linearized approach of Nichita and Graciaa (2011), which uses the twoparameter BIP formula of Li and Johns (2006), is also demonstrated to have a significantly larger radius of convergence than other reduced and conventional methods for five fluids studied.

AB - Phase-equilibrium calculations become computationally intensive in compositional simulation as the number of components and phases increases. Reduced methods were developed to address this problem, where the binary-interaction-parameter (BIP) matrix is approximated either by spectral decomposition (SD), as performed by Hendriks and van Bergen (1992), or with the twoparameter BIP formula of Li and Johns (2006). Several authors have recently stated that the SD method-And by reference all reduced methods-is not as fast as previously reported in the literature. In this paper we present the first study that compares all eight reduced and conventional methods published to date by use of optimized code and compilers. The results show that the SD method and its variants are not as fast as other reduced methods, and can be slower than the conventional approach when fewer than 10 components are used. These conclusions confirm the findings of recently published papers. The reason for the slow speed is the requirement that the code must allow for a variable number of eigenvalues. We show that the reduced method of Li and Johns (2006) and its variants, however, are faster because the number of reduced parameters is fixed to six, which is independent of the number of components. Speed up in flash calculations for their formula is achieved for all fluids studied when more than six components are used. For example, for 10-component fluids, a speed up of 2-3 in the computational time for Newton-Raphson (NR) iterations is obtained compared with the conventional method modeled after minimization of Gibbs energy. The reduced method modeled after the linearized approach of Nichita and Graciaa (2011), which uses the twoparameter BIP formula of Li and Johns (2006), is also demonstrated to have a significantly larger radius of convergence than other reduced and conventional methods for five fluids studied.

UR - http://www.scopus.com/inward/record.url?scp=84928491979&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84928491979&partnerID=8YFLogxK

U2 - 10.2118/163577-PA

DO - 10.2118/163577-PA

M3 - Article

VL - 20

SP - 294

EP - 305

JO - SPE Journal

JF - SPE Journal

SN - 1086-055X

IS - 2

ER -