legacy_thermo_properties Module


Uses


Subroutines

public subroutine PUREFUG_CALC(nc, icomp, T, P, V, phi)

Fugacity of a pure component

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nc
integer, intent(in) :: icomp
real(kind=pr), intent(in) :: T
real(kind=pr), intent(in) :: P
real(kind=pr), intent(in) :: V
real(kind=pr), intent(out) :: phi

public subroutine TERMO(nc, mtyp, indic, t, p, rn, v, PHILOG, DLPHIP, DLPHIT, FUGN)

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nc

Number of components
integer, intent(in) :: mtyp

Type of root desired (-1 vapor, 1 liquid, 0 lower Gr)
integer, intent(in) :: indic

Desired element, this should be setted with optionals
real(kind=pr), intent(in) :: t

Temperature [K]
real(kind=pr), intent(in) :: p

Pressure [bar]
real(kind=pr), intent(in) :: rn(nc)

Mixture mole numbers
real(kind=pr), intent(out) :: v

Volume [L]
real(kind=pr), intent(out) :: PHILOG(nc)

ln(phi*p) vector
real(kind=pr), intent(out), optional :: DLPHIP(nc)

ln(phi) Presssure derivative
real(kind=pr), intent(out), optional :: DLPHIT(nc)

ln(phi) Temp derivative
real(kind=pr), intent(out), optional :: FUGN(nc,nc)

ln(phi) compositional derivative

public recursive subroutine VCALC(ITYP, nc, ntemp, rn, T, P, V)

ROUTINE FOR CALCULATION OF VOLUME, GIVEN PRESSURE

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: ITYP

TYPE OF ROOT DESIRED (-1 vapor, 1 liquid, 0 lower Gibbs energy phase)
integer, intent(in) :: nc

NO. OF COMPONENTS
integer, intent(in) :: ntemp

1 if T-derivatives are required
real(kind=pr), intent(in) :: rn(nc)

FEED MOELS
real(kind=pr), intent(in) :: T

TEMPERATURE
real(kind=pr), intent(in) :: P

PRESURE
real(kind=pr), intent(out) :: V

VOLUME

public subroutine zTVTERMO(nc, indic, t, rn, v, p, dpv, PHILOG, DLPHIP, DLPHIT, FUGN)

Calculation of lnphiP and derivatives rn mixture mole numbers (input) t temperature (k) (input) v volume (L) (input) p pressure (bar) (output) PHILOG vector of ln(phi(i)P) (output) 0 < INDIC < 5 DLPHIT t-derivative of ln(phi(i)) (const P, n) (output) 0 < INDIC = 2 or 4 DLPHIP P-derivative of ln(phi(i)) (const T, n) (output) 0 < INDIC < 5 FUGN comp-derivative of ln(phi(i)) (const t & P)(output) 2 < INDIC

Read more…

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nc
integer, intent(in) :: indic
real(kind=pr), intent(in) :: t
real(kind=pr), intent(in) :: rn(nc)
real(kind=pr), intent(in) :: v
real(kind=pr), intent(out) :: p
real(kind=pr), intent(out) :: dpv
real(kind=pr), intent(out) :: PHILOG(nc)
real(kind=pr), intent(out) :: DLPHIP(nc)
real(kind=pr), intent(out) :: DLPHIT(nc)
real(kind=pr), intent(out) :: FUGN(nc,nc)