Procedures of the core calculations of CubicEoS mixing rules.
This module holds all the basic math to use mixing rules in other codes. Keeping it simple and accesible.
bi = [0.2, 0.3]
lij = reshape([0.0, 0.2, 0.2, 0], [2,2])
! Calculate B parameter with Quadratric Mixing Rules.
call bmix_qmr(n, bi, lij, b, dbi, dbij)
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=pr), | intent(in) | :: | n(:) | |||
| real(kind=pr), | intent(in) | :: | bi(:) | |||
| real(kind=pr), | intent(out) | :: | b | |||
| real(kind=pr), | intent(out) | :: | dbi(:) | |||
| real(kind=pr), | intent(out) | :: | dbij(:,:) |
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=pr), | intent(in) | :: | n(:) | |||
| real(kind=pr), | intent(in) | :: | bi(:) | |||
| real(kind=pr), | intent(in) | :: | lij(:,:) | |||
| real(kind=pr), | intent(out) | :: | b | |||
| real(kind=pr), | intent(out) | :: | dbi(:) | |||
| real(kind=pr), | intent(out) | :: | dbij(:,:) |
RKPR parameter mixing rule.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=pr), | intent(in) | :: | n(:) | |||
| real(kind=pr), | intent(in) | :: | d1i(:) | |||
| real(kind=pr), | intent(out) | :: | D1 | |||
| real(kind=pr), | intent(out) | :: | dD1i(:) | |||
| real(kind=pr), | intent(out) | :: | dD1ij(:,:) |
Infinite pressure limit parameter
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=pr), | intent(in) | :: | d1 | |||
| real(kind=pr), | intent(in) | :: | dd1i(:) | |||
| real(kind=pr), | intent(in) | :: | dd1ij(:,:) | |||
| real(kind=pr), | intent(out) | :: | L | |||
| real(kind=pr), | intent(out) | :: | dLi(:) | |||
| real(kind=pr), | intent(out) | :: | dLij(:,:) |