PNG 520
Phase Relations in Reservoir Engineering

Solution Algorithms for VLE Problems

PrintPrint

We have essentially reduced the typical VLE problem to that of solving a system of non-linear algebraic equations. The objective function is the nucleus of the VLE calculation. Rather than solve the material balance equations in their elemental form, the use of an objective function will allow for a more stable iteration scheme. We have seen that three objective functions are available for this calculation:

F y ( α g )= i=1 n z i K i 1+ α g ( K i 1) 1=0 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.7a)
 
F x ( α g )= i=1 n z i 1+ α g ( K i 1) 1=0 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.7b)
 
F( α g )= i=1 n z i ( K i 1) 1+ α g ( K i 1) =0 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.7c)
 

Expression (17.7c) is the Rachford-Rice Objective Function; we have proven that its numerical solution is cleaner than that of equations (17.7a) and (17.7b). We also discussed that the key problem in solving for “ α g This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers. ” is the need for values for the equilibrium ratios (Ki’s). With the aid of equations (17.5) and (17.6) these problems can be solved simultaneously by iteration techniques.

In the problem we have posed, we have (3n+1) unknowns:

K i  i=1, 2, ... n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
 
y i  i=1, 2, ... n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
 
x i  i=1, 2, ... n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
 

And, α g This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.

And we also have (3n+1) equations:

K i = ϕ i L (P,T, x i ) ϕ i V (P,T, y i ) i=1,2,...n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.8a)
 
y i = K i x i   i=1,2,...n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.8b)
 
x i = z i 1+ α g ( K i 1)   i=1,2,...n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.8c)
 
i=1 n z i ( K i 1) 1+ α g ( K i 1) =0 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.8d)
 

where P, T, and zi are the values that we know. See how this is only a restatement of the original set of equations that contain the equilibrium conditions, material balances, and molar fraction constraints:

f i V = f i L   i=1, 2, ... n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.9a)
 
y i = K i, x i   i=1, 2, ... n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.9b)
 
y i α g + x i (1 α g )= z i   i=1, 2, ... n This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.9c)
 
i=1 n c x i =1  or   i=1 n c y i =1 This equation is not rendering properly due to an incompatible browser. See Technical Requirements in the Orientation for a list of compatible browsers.
(17.9d)
 

Since we have the same number of equations as unknowns, this system is determined and has a unique solution. To solve this system of simultaneous non-linear equations, there are two categories of solution techniques:

  • Newton-type methods,
  • Substitution-type methods.

Newton-type methods for more than one unknown require finding (n x n) elements of a Jacobian matrix at every iteration step. This could be computationally expensive. Additionally, the Newton Raphson method requires a very good initial guess (i.e., close to the actual solution) to guarantee convergence to the true values. This is not always possible, especially at the start of the procedure.

The most popular method, and the easiest to implement, is the substitution type. However, substitution type methods can be quite slow for some conditions of interest. In these sort of cases, we either switch to a Newton-Raphson solver (that has a much better rate of convergence) or implement an acceleration procedure to the substitution method itself.