SoaveRedlichKwong Function

public function SoaveRedlichKwong(tc, pc, w, kij, lij) result(model)

SoaveRedlichKwong.

Using the critical constants setup the parameters to use the SoaveRedlichKwong Equation of State

• $$\alpha(T_r) = (1 + k (1 - \sqrt{T_r}))^2$$
• $$k = 0.48 + 1.574 \omega - 0.175 \omega^2$$
• $$a_c = 0.427480 R^2 * T_c^2/P_c$$
• $$b = 0.086640 R T_c/P_c$$
• $$\delta_1 = 1$$
• $$\delta_2 = 0$$

There is also the optional posibility to include the k_{ij} and l_{ij} matrices. Using by default Classic Van der Waals mixing rules.

After setting up the model, it is possible to redefine either the mixing rule or the alpha function using a different derived type defined outside the function.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=pr),	intent(in)			::	tc(:)	Critical temperature [K]
real(kind=pr),	intent(in)			::	pc(:)	Critical pressure [bar]
real(kind=pr),	intent(in)			::	w(:)	Acentric factor
real(kind=pr),	intent(in),	optional		::	kij(:,:)	$k_{ij}$ matrix
real(kind=pr),	intent(in),	optional		::	lij(:,:)	$l_{ij}$ matrix

Return Value type(CubicEoS)

Variables

Type	Visibility	Attributes		Name		Initial
type(AlphaSoave),	private		::	alpha
type(Substances),	private		::	composition
integer,	private		::	i
type(QMR),	private		::	mixrule
integer,	private		::	nc