stability_check Subroutine

private subroutine stability_check(model, z0, zi, Pc, Vc, Tc, a, unstable, y_other, V_other)

Uses

  • proc~~stability_check~~UsesGraph proc~stability_check stability_check module~yaeos__equilibria_stability yaeos__equilibria_stability proc~stability_check->module~yaeos__equilibria_stability module~yaeos__constants yaeos__constants module~yaeos__equilibria_stability->module~yaeos__constants module~yaeos__models_ar yaeos__models_ar module~yaeos__equilibria_stability->module~yaeos__models_ar module~yaeos__models_base yaeos__models_base module~yaeos__equilibria_stability->module~yaeos__models_base module~yaeos__models_ge yaeos__models_ge module~yaeos__equilibria_stability->module~yaeos__models_ge iso_fortran_env iso_fortran_env module~yaeos__constants->iso_fortran_env module~yaeos__models_ar->module~yaeos__constants module~yaeos__models_ar->module~yaeos__models_base module~yaeos__substance yaeos__substance module~yaeos__models_base->module~yaeos__substance module~yaeos__models_ge->module~yaeos__constants module~yaeos__models_ge->module~yaeos__models_base module~yaeos__substance->module~yaeos__constants

stability_check

Description

Check the stability of a point in the critical line. The stability is determined by tpd analysis.

Arguments

Type IntentOptional Attributes Name
class(ArModel), intent(in) :: model

Equation of state model
real(kind=pr), intent(in) :: z0(:)

Molar fractions of the first fluid
real(kind=pr), intent(in) :: zi(:)

Molar fractions of the second fluid
real(kind=pr), intent(in) :: Pc

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

Volume [L/mol]
real(kind=pr), intent(in) :: Tc

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

Molar fraction of the second fluid
logical, intent(out) :: unstable

Stability of the point)
real(kind=pr), intent(out) :: y_other(:)

Molar fractions of the second fluid
real(kind=pr), intent(out) :: V_other

Volume [L/mol]

Calls

proc~~stability_check~~CallsGraph proc~stability_check stability_check proc~lnfug_vt ArModel%lnfug_vt proc~stability_check->proc~lnfug_vt proc~volume~3 ArModel%volume proc~stability_check->proc~volume~3 residual_helmholtz residual_helmholtz proc~lnfug_vt->residual_helmholtz get_v0 get_v0 proc~volume~3->get_v0 interface~newton newton proc~volume~3->interface~newton proc~newton_1d newton_1d interface~newton->proc~newton_1d

Called by

proc~~stability_check~~CalledByGraph proc~stability_check stability_check proc~look_for_cep look_for_cep proc~look_for_cep->proc~stability_check proc~critical_line critical_line proc~critical_line->proc~look_for_cep

Variables

Type Visibility Attributes Name Initial
real(kind=pr), private :: P
real(kind=pr), private :: dy
logical, private :: first
real(kind=pr), private :: fug_y(2)
real(kind=pr), private :: fug_z(2)
integer, private :: istab
integer, private :: istab0
real(kind=pr), private :: mintpd
logical, private :: possible
real(kind=pr), private :: tpd
real(kind=pr), private :: y(2)
real(kind=pr), private :: z(2)