detect_critical Subroutine

public subroutine detect_critical(nc, np, point, kinds_x, kind_w, binary_stop, Xold, X, dXdS, ns, dS, S, found_critical, Xc)

Uses

  • proc~~detect_critical~3~~UsesGraph proc~detect_critical~3 detect_critical module~yaeos__math yaeos__math proc~detect_critical~3->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

detect_critical

Detect if the system is close to a critical point.

Description

When the system is close to a critical point, the values are close to zero, since the composition of the incipient phase and the phase are similar (equal in the critical point). This can be used to detect if the system is close to a critical point and force a jump above it.

Arguments

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

Number of components in the mixture.
integer, intent(in) :: np

Number of main phases.
integer, intent(in) :: point

Point number in the phase boundary.
character(len=14), intent(inout) :: kinds_x(np)

Kinds of the main phases.
character(len=14), intent(inout) :: kind_w

Kind of the incipient phase.
logical, intent(in) :: binary_stop

If true, stop at the critical point if its a binary system.
real(kind=pr), intent(in) :: Xold(:)

Old vector of variables.
real(kind=pr), intent(inout) :: X(:)

Vector of variables.
real(kind=pr), intent(inout) :: dXdS(:)

Sensitivity of the variables wrt the specification.
integer, intent(inout) :: ns

Number of the specified variable.
real(kind=pr), intent(inout) :: dS

Step size of the specification for the next point.
real(kind=pr), intent(inout) :: S

Specification value.
logical, intent(out) :: found_critical

If true, a critical point was found.
real(kind=pr) :: Xc(size(X))

Vector of variables at the critical point.

Calls

proc~~detect_critical~3~~CallsGraph proc~detect_critical~3 detect_critical proc~interpol interpol proc~detect_critical~3->proc~interpol

Called by

proc~~detect_critical~3~~CalledByGraph proc~detect_critical~3 detect_critical proc~pt_envelope pt_envelope proc~pt_envelope->proc~detect_critical~3 proc~px_envelope px_envelope proc~px_envelope->proc~detect_critical~3 proc~tx_envelope tx_envelope proc~tx_envelope->proc~detect_critical~3

Variables

Type Visibility Attributes Name Initial
real(kind=pr), public :: P
real(kind=pr), public :: T
real(kind=pr), public :: Xnew(size(X))
real(kind=pr), public :: a
integer, public :: i
character(len=14), public :: incipient_kind
integer, public :: lb
real(kind=pr), public :: limit
real(kind=pr), public :: lnK(nc)
real(kind=pr), public :: lnKold(nc)
integer, public :: ncomp
integer, public :: ub