fugacity_vt Subroutine

public subroutine fugacity_vt(eos, n, V, T, P, lnPhi, dlnPhidP, dlnPhidT, dlnPhidn, dPdV, dPdT, dPdn)

Calculate fugacity coefficent given volume and temperature.

Note

While the natural output variable is $ln \phi_i P$ . The calculated derivatives will be the derivatives of the fugacity coefficient $ln \phi_i$

Type Bound

Arguments

Type	Intent	Optional	Attributes		Name
class(ArModel)				::	eos	Model
real(kind=pr),	intent(in)			::	n(:)	Mixture mole numbers
real(kind=pr),	intent(in)			::	V	Volume [L]
real(kind=pr),	intent(in)			::	T	Temperature [K]
real(kind=pr),	intent(out),	optional		::	P	Pressure [bar]
real(kind=pr),	intent(out),	optional		::	lnPhi(size(n))	$\ln(\phi_i*P)$ vector
real(kind=pr),	intent(out),	optional		::	dlnPhidP(size(n))	$ln(phi_i)$ Presssure derivative
real(kind=pr),	intent(out),	optional		::	dlnPhidT(size(n))	$ln(phi_i)$ Temp derivative
real(kind=pr),	intent(out),	optional		::	dlnPhidn(size(n),size(n))	$ln(phi_i)$ compositional derivative
real(kind=pr),	intent(out),	optional		::	dPdV	$\frac{dP}{dV}$
real(kind=pr),	intent(out),	optional		::	dPdT	$\frac{dP}{dT}$
real(kind=pr),	intent(out),	optional		::	dPdn(:)	$\frac{dP}{dn_i}$

Variables

Type	Visibility	Attributes		Name		Initial
real(kind=pr),	public		::	Ar
real(kind=pr),	public		::	ArTV
real(kind=pr),	public,	dimension(size(n))	::	ArTn
real(kind=pr),	public		::	ArV
real(kind=pr),	public		::	ArV2
real(kind=pr),	public,	dimension(size(n))	::	ArVn
real(kind=pr),	public,	dimension(size(n))	::	Arn
real(kind=pr),	public		::	Arn2(size(n),size(n))
real(kind=pr),	public		::	P_in
real(kind=pr),	public		::	RT
real(kind=pr),	public		::	Z
real(kind=pr),	public		::	dPdT_in
real(kind=pr),	public		::	dPdV_in
real(kind=pr),	public		::	dPdn_in(size(n))
integer,	public		::	i
integer,	public		::	j
integer,	public		::	nc
real(kind=pr),	public		::	totn