lnphi_vt Subroutine

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

Calculate natural logarithm of fugacity coefficent.

Calculate the natural logarithm of the fugacity coefficient and its derivatives given volume and temperature. The routine gives the possibility to calculate the pressure and it’s derivatives.

Examples

eos = PengRobinson76(Tc, Pc, w)

n = [1.0_pr, 1.0_pr]
T = 300.0_pr
V = 1.0_pr

call eos%lnphi_vt(&
   n, V, T, lnPhi=lnPhi, &
   dlnPhidP=dlnPhidP, dlnPhidT=dlnPhidT, dlnPhidn=dlnPhidn &
   )

Type Bound

ArModel

Arguments

Type IntentOptional 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))

vector
real(kind=pr), intent(out), optional :: dlnPhidP(size(n))

Presssure derivative
real(kind=pr), intent(out), optional :: dlnPhidT(size(n))

Temp derivative
real(kind=pr), intent(out), optional :: dlnPhidn(size(n),size(n))

compositional derivative
real(kind=pr), intent(out), optional :: dPdV

real(kind=pr), intent(out), optional :: dPdT

real(kind=pr), intent(out), optional :: dPdn(:)


Variables

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