public subroutine rachford_rice(z, K, beta, rr, drrdb)

rachford_rice

Description

Rachford-Rice equation for a two phase system. This equation is used to calculate the $\beta$ value that satisfies the mass balance between two phases.

$$rr(\beta) = \sum_i \frac{z_i(K_i - 1)}{1 + \beta(K_i - 1)}]$$

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=pr),	intent(in)			::	z(:)	Mole fractions vector
real(kind=pr),	intent(in)			::	K(:)	K-factors
real(kind=pr),	intent(in)			::	beta	$\beta$ value
real(kind=pr),	intent(out)			::	rr	Rachford-Rice function value
real(kind=pr),	intent(out)			::	drrdb	Derivative of the Rachford-Rice function

Variables

Type	Visibility	Attributes		Name		Initial
real(kind=pr),	public		::	denom(size(z))