17.5: The Accelerated Successive Substitution Method (ASSM)

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$ $$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\kernel}{\mathrm{null}\,}$$ $$\newcommand{\range}{\mathrm{range}\,}$$ $$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$ $$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}[1]{\| #1 \|}$$ $$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\kernel}{\mathrm{null}\,}$$ $$\newcommand{\range}{\mathrm{range}\,}$$ $$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$ $$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}[1]{\| #1 \|}$$ $$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$

When the system is close to the critical point and fugacities are strongly composition-dependent, a slowing-down of the convergence rate of the SSM (Successive Substitution Method) is to be expected. In an attempt to avoid slow convergence problems, some methods have been proposed. Among the most popular are the Minimum Variable Newton Raphson (MVNR) Method and the Accelerated and Stabilized Successive Substitution Method (ASSM).

The ASSM is basically an accelerated version of the SSM procedure, and thus follows a similar theory. Such procedure is implemented to accelerate the calculation of Ki-values, especially in the region close to critical point where the use of the SSM alone will not be efficient. The ASSM technique was presented by Rinses et al. (1981) and consists of the following steps:

1. Use the SSM technique to initiate the updating of the Ki-values the first time.
2. Check all following criteria at every step during iterations using the SSM: (17.13a)
agnew - agold | < 0.1 (17.13b)
(17.13c)
0 < agnew < 1(17.13d)

These criteria show that you have sufficient proximity to the conditions to ensure the efficiency of the method. Rri is the ratio of liquid fugacity to gas fugacity of the i-th component and ‘αg’ is molar gas fraction of the two-phase system.

(17.14)

3. If the system satisfies ALL above criteria, the iteration technique is then switched from the SSM to the ASSM. Otherwise, SSM is used for the update of the Ki-values. The following expressions are used to update Ki-values in ASSM:

Ki new= Kiold Rri λi

where λi =[ ( Rriold -1) / (Rriold - Rrinew ) ]

In some cases, using a constant acceleration value of λi =2 is good enough.

4. Once all the criteria in step (2) are satisfied, skip step (2) for the subsequent iterations and use the ASSM technique to update Ki-values until convergence is attained, unless it does not give acceptable new estimates (as stated next).
5. When ASSM is used, it must always be tested to show that it leads to an improved solution (i.e., that it brings fugacity ratios closer to unity). If not, it must be rejected and switched back to SSM.

Even though we are outlining Risnes et al.’s version of the accelerated Successive Substitution Method, there are several other published algorithms whose main purpose have also been to accelerate the successive substitution method. Fussel and Yanosik (1978), Michelsen (1982), and Mehra et al. (1983) are examples of such attempts. Risnes et al. version is the easiest and most straightforward to implement, but it is subjected to limitations.