yaeos__models_ar_cubic_implementations Module

Implemented Cubic Equations of State.

PengRobinson76
PengRobinson78
SoaveRedlichKwong
RKPR

Uses

Functions

public function PSRK(tc, pc, w, molecules, c1, c2, c3) result(model)

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
type(Groups),	intent(in)			::	molecules(:)
real(kind=pr),	intent(in),	optional		::	c1(:)
real(kind=pr),	intent(in),	optional		::	c2(:)
real(kind=pr),	intent(in),	optional		::	c3(:)

Return Value type(CubicEoS)

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

PengRobinson76.

Arguments

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

Return Value type(CubicEoS)

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

PengRobinson78.

Arguments

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

Return Value type(CubicEoS)

public function RKPR(tc, pc, w, zc, kij, lij, delta_1, k) result(model)

RKPR Equation of State

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)			::	zc(:)	Critical compressibility
real(kind=pr),	intent(in),	optional		::	kij(:,:)	k_{ij} matrix
real(kind=pr),	intent(in),	optional		::	lij(:,:)	l_{ij} matrix
real(kind=pr),	intent(in),	optional		::	delta_1(:)
real(kind=pr),	intent(in),	optional		::	k(:)

Return Value type(CubicEoS)

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

SoaveRedlichKwong.

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)

Subroutines

public subroutine refit_rkpr_k(model, component)

Refit the $k$ parameter of the RKPR EoS to match the acentric factor

Arguments

Type	Intent	Optional	Attributes		Name
type(CubicEoS),	intent(inout)			::	model	The model to be refitted
integer,	intent(in)			::	component	Component index to refit

private subroutine get_OMa_OMb(del1, OMa, OMb)

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=pr),	intent(in)			::	del1(:)
real(kind=pr),	intent(out)			::	OMa(size(del1))
real(kind=pr),	intent(out)			::	OMb(size(del1))