saturation_temperature Function

public function saturation_temperature(model, n, p, kind, t0, y0, max_iters)

Uses

- yaeos__math_continuation
- stdlib_optval

Saturation temperature calculation function.

Calculates the saturation pressure of a multicomponent mixture with a given molar composition n. It is possible to calculate:

Bubble point: kind="bubble"
Dew point: kind="dew"
Liquid-Liquid point: kind="liquid-liquid"

Arguments

Type	Intent	Optional	Attributes		Name
class(ArModel),	intent(in),		target	::	model
real(kind=pr),	intent(in)			::	n(:)	Composition vector [moles / molar fraction]
real(kind=pr),	intent(in)			::	p	Pressure [bar]
character(len=*),	intent(in)			::	kind	[bubble\|dew\|liquid-liquid]
real(kind=pr),	intent(in),	optional		::	t0	Initial temperature [K]
real(kind=pr),	intent(in),	optional		::	y0(:)	Initial composition
integer,	intent(in),	optional		::	max_iters	Maximum number of iterations

Return Value type(EquilibriumState)

Variables

Type	Visibility	Attributes		Name		Initial
real(kind=pr),	public		::	dlnphi_dt_y(size(n))
real(kind=pr),	public		::	dlnphi_dt_z(size(n))
real(kind=pr),	public		::	f
character(len=50),	public		::	incipient
logical,	public		::	is_incipient(size(n))
integer,	public		::	iterations
integer,	public		::	its
real(kind=pr),	public		::	k(size(n))
real(kind=pr),	public		::	lnfug_y(size(n))
real(kind=pr),	public		::	lnfug_z(size(n))
real(kind=pr),	public		::	lnk(size(n))
character(len=50),	public		::	main
real(kind=pr),	public		::	step
real(kind=pr),	public		::	t
real(kind=pr),	public		::	vy
real(kind=pr),	public		::	vz
real(kind=pr),	public		::	y(size(n))
real(kind=pr),	public		::	z(size(n))