saturation_pressure Function

public function saturation_pressure(model, n, t, kind, p0, y0, max_iters)

Uses

Saturation pressure calculation function.

Calculates the saturation pressure of a multicomponent mixture with a given molar composition n. It is possible to calculate:

  • Bubble point: kind="bubble"
  • Dew point: kind="dew"
  • Liquid-Liquid point: kind="liquid-liquid"

Arguments

Type IntentOptional Attributes Name
class(ArModel), intent(in), target :: model
real(kind=pr), intent(in) :: n(:)

Composition vector [moles / molar fraction]

real(kind=pr), intent(in) :: t

Temperature [K]

character(len=*), intent(in) :: kind

[bubble|dew|liquid-liquid]

real(kind=pr), intent(in), optional :: p0

Initial pressure [bar]

real(kind=pr), intent(in), optional :: y0(:)

Initial composition

integer, intent(in), optional :: max_iters

Maximum number of iterations

Return Value type(EquilibriumState)


Variables

Type Visibility Attributes Name Initial
real(kind=pr), public :: dlnphi_dp_y(size(n))
real(kind=pr), public :: dlnphi_dp_z(size(n))
real(kind=pr), public :: f
integer, public :: i
character(len=50), public :: incipient
integer, public :: iterations
integer, public :: its
real(kind=pr), public :: k(size(n))
real(kind=pr), public :: lnfug_y(size(n))
real(kind=pr), public :: lnfug_z(size(n))
real(kind=pr), public :: lnk(size(n))
character(len=50), public :: main
real(kind=pr), public :: p
real(kind=pr), public :: step
real(kind=pr), public :: vy
real(kind=pr), public :: vz
real(kind=pr), public :: y(size(n))
real(kind=pr), public :: z(size(n))