kij_exp_tdep Subroutine

public subroutine kij_exp_tdep(self, T, a, dadt, dadt2, aij, daijdt, daijdt2)

Uses

- hyperdual_mod

kij_exp_tdep

Combining rule that uses temperature dependant $k_{ij}$ values. With the following expression: $k_{ij}(T) = k_{ij}^0 + k_{ij}^\infty \exp\left(\frac{-T}{T^*}\right)$

[ a_{ij} = \sqrt{a_i a_j} (1 - k_{ij}) ]

Type Bound

Arguments

Type	Intent	Optional	Attributes		Name
class(QMRTD),	intent(in)			::	self
real(kind=pr),	intent(in)			::	T	Temperature [K]
real(kind=pr),	intent(in)			::	a(:)	Pure components attractive parameters (\a_i)
real(kind=pr),	intent(in)			::	dadt(:)	$\frac{da_i}{dT}$
real(kind=pr),	intent(in)			::	dadt2(:)	$\frac{d^2a_i}{dT^2}$
real(kind=pr),	intent(out)			::	aij(:,:)	$a_{ij}$ Matrix
real(kind=pr),	intent(out)			::	daijdt(:,:)	$\frac{da_{ij}{dT}$
real(kind=pr),	intent(out)			::	daijdt2(:,:)	$\frac{d^2a_{ij}{dT^2}$

Variables

Type	Visibility	Attributes		Name		Initial
type(hyperdual),	public		::	T_hd
real(kind=pr),	public		::	Tstar(size(a),size(a))
type(hyperdual),	public		::	a_hd(size(a))
type(hyperdual),	public		::	aij_hd(size(a),size(a))
integer,	public		::	i
integer,	public		::	j
real(kind=pr),	public		::	k0(size(a),size(a))
type(hyperdual),	public		::	kij_hd(size(a),size(a))
real(kind=pr),	public		::	kinf(size(a),size(a))
integer,	public		::	nc