pressure_equality_V_beta_xy Subroutine

public subroutine pressure_equality_V_beta_xy(model, T, V, beta, x, y, Vx, Vy, P)

Solve pressure equality between two phases at a given temperature, total volume, vapor molar fractions and compositions.

Arguments

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

Temperature [K]

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

Total volume [L/mol]

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

Molar fraction of light-phase

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

Molar fractions of heavy-phase

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

Molar fractions of light-phase

real(kind=pr), intent(inout) :: Vx

Heavy-phase molar volume [L/mol]

real(kind=pr), intent(inout) :: Vy

Light-Phase molar volume [L/mol]

real(kind=pr), intent(out) :: P

Pressure [bar]


Variables

Type Visibility Attributes Name Initial
real(kind=pr), public :: Bx

Liquid phase covolume

real(kind=pr), public :: Px
real(kind=pr), public :: Py
real(kind=pr), public :: dPxdV
real(kind=pr), public :: dPydV
real(kind=pr), public :: dVydVx

Derivative of Vy wrt Vx

real(kind=pr), public :: dh

dh/

real(kind=pr), public :: h

Pressure equality

integer, public :: its
real(kind=pr), public :: stepv