look_for_cep Subroutine

private subroutine look_for_cep(model, z0, zi, Pc, Vc, Tc, a, u, found, CEP)

Uses

  • proc~~look_for_cep~~UsesGraph proc~look_for_cep look_for_cep module~yaeos__math yaeos__math proc~look_for_cep->module~yaeos__math module~yaeos__constants yaeos__constants module~yaeos__math->module~yaeos__constants module~yaeos__math_continuation yaeos__math_continuation module~yaeos__math->module~yaeos__math_continuation module~yaeos__math_linalg yaeos__math_linalg module~yaeos__math->module~yaeos__math_linalg iso_fortran_env iso_fortran_env module~yaeos__constants->iso_fortran_env module~yaeos__math_continuation->module~yaeos__constants module~yaeos__math_continuation->module~yaeos__math_linalg module~yaeos__math_linalg->module~yaeos__constants

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
real(kind=pr), intent(inout) :: u(:)

Eigen-vector
logical, intent(out) :: found

Found a Critical End Point
type(EquilibriumState), intent(out) :: CEP

Critical End Point

Calls

proc~~look_for_cep~~CallsGraph proc~look_for_cep look_for_cep proc~df_cep df_cep proc~look_for_cep->proc~df_cep proc~f_cep F_cep proc~look_for_cep->proc~f_cep proc~get_a get_a proc~look_for_cep->proc~get_a proc~pressure ArModel%pressure proc~look_for_cep->proc~pressure proc~set_a set_a proc~look_for_cep->proc~set_a proc~solve_system solve_system proc~look_for_cep->proc~solve_system proc~stability_check stability_check proc~look_for_cep->proc~stability_check proc~df_cep->proc~f_cep proc~df_cep->proc~get_a proc~f_cep->proc~get_a proc~lambda1 lambda1 proc~f_cep->proc~lambda1 proc~lnfug_vt ArModel%lnfug_vt proc~f_cep->proc~lnfug_vt residual_helmholtz residual_helmholtz proc~pressure->residual_helmholtz none~dgesv dgesv proc~solve_system->none~dgesv proc~stability_check->proc~lnfug_vt proc~volume~3 ArModel%volume proc~stability_check->proc~volume~3 proc~lambda1->proc~get_a proc~lambda1->proc~lnfug_vt proc~eigen eigen proc~lambda1->proc~eigen 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 none~dsyev dsyev proc~eigen->none~dsyev

Called by

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

Variables

Type Visibility Attributes Name Initial
real(kind=pr), private :: Fcep(size(z0)+4)
real(kind=pr), private :: Fcep_new(size(z0)+4)
real(kind=pr), private :: V_cep
real(kind=pr), private :: Xcep(size(z0)+4)
real(kind=pr), private :: alpha
real(kind=pr), private :: dFcep(size(z0)+4,size(z0)+4)
real(kind=pr), private :: dXcep(size(z0)+4)
real(kind=pr), private :: damp
integer, private :: inits
integer, private :: its
integer, private :: nc
real(kind=pr), private :: y_cep(size(z0))