pt_F_three_phases Subroutine

private subroutine pt_F_three_phases(model, z, Xvars, ns, S, F, df)

Uses

  • proc~~pt_f_three_phases~~UsesGraph proc~pt_f_three_phases pt_F_three_phases iso_fortran_env iso_fortran_env proc~pt_f_three_phases->iso_fortran_env

Function to solve at each point of a three phase envelope.

The vector of variables X corresponds to:

While the equations are:

Arguments

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

Vector of variables
integer, intent(in) :: ns

Number of specification
real(kind=pr), intent(in) :: S

Specification value
real(kind=pr), intent(out) :: F(size(Xvars))

Vector of functions valuated
real(kind=pr), intent(out) :: df(size(Xvars),size(Xvars))

Jacobian matrix

Calls

proc~~pt_f_three_phases~~CallsGraph proc~pt_f_three_phases pt_F_three_phases proc~lnphi_pt ArModel%lnphi_pt proc~pt_f_three_phases->proc~lnphi_pt proc~lnphi_vt ArModel%lnphi_vt proc~lnphi_pt->proc~lnphi_vt proc~volume~3 ArModel%volume proc~lnphi_pt->proc~volume~3 residual_helmholtz residual_helmholtz proc~lnphi_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~~pt_f_three_phases~~CalledByGraph proc~pt_f_three_phases pt_F_three_phases proc~solve_point~5 solve_point proc~solve_point~5->proc~pt_f_three_phases proc~pt_envelope_3ph pt_envelope_3ph proc~pt_envelope_3ph->proc~solve_point~5

Variables

Type Visibility Attributes Name Initial
real(kind=pr), private :: Kx((Size(Xvars)-3)/2)
real(kind=pr), private :: Ky((Size(Xvars)-3)/2)
real(kind=pr), private :: P
real(kind=pr), private :: T
real(kind=pr), private :: Vw
real(kind=pr), private :: Vx
real(kind=pr), private :: Vy
real(kind=pr), private :: beta
real(kind=pr), private, dimension((Size(Xvars)-3)/2, (Size(Xvars)-3)/2) :: dlnphi_dn_w
real(kind=pr), private, dimension((Size(Xvars)-3)/2, (Size(Xvars)-3)/2) :: dlnphi_dn_x
real(kind=pr), private, dimension((Size(Xvars)-3)/2, (Size(Xvars)-3)/2) :: dlnphi_dn_y
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: dlnphi_dp_w
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: dlnphi_dp_x
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: dlnphi_dp_y
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: dlnphi_dt_w
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: dlnphi_dt_x
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: dlnphi_dt_y
real(kind=pr), private :: dwdKx((Size(Xvars)-3)/2)
real(kind=pr), private :: dwdKy((Size(Xvars)-3)/2)
real(kind=pr), private :: dwdb((Size(Xvars)-3)/2)
real(kind=pr), private :: dxdKx((Size(Xvars)-3)/2)
real(kind=pr), private :: dxdKy((Size(Xvars)-3)/2)
real(kind=pr), private :: dydKx((Size(Xvars)-3)/2)
real(kind=pr), private :: dydKy((Size(Xvars)-3)/2)
integer, private :: i
integer, private :: j
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: lnphi_w
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: lnphi_x
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: lnphi_y
integer, private :: nc
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: w
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: x
real(kind=pr), private, dimension((Size(Xvars)-3)/2) :: y