Calculates the critical line between two mixtures using the continuation method. The two mixtures compositions are restricted to the relation between them, by a parameter , which represents the molar fraction of the second fluid with respect to the whole mixture.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(ArModel), | intent(in) | :: | model |
Equation of state model |
||
| real(kind=pr), | intent(in) | :: | a0 |
Initial value |
||
| real(kind=pr), | intent(in) | :: | z0(:) |
Molar fractions of the first fluid |
||
| real(kind=pr), | intent(in) | :: | zi(:) |
Molar fractions of the second fluid |
||
| integer, | intent(in) | :: | ns0 |
Position of the specification variable |
||
| real(kind=pr), | intent(in) | :: | S0 |
Specified value |
||
| real(kind=pr), | intent(in) | :: | dS0 |
Initial step size |
||
| real(kind=pr), | intent(in), | optional | :: | v0 |
Initial volume [L/mol] |
|
| real(kind=pr), | intent(in), | optional | :: | t0 |
Initial temperature [K] |
|
| real(kind=pr), | intent(in), | optional | :: | p0 |
Initial pressure [bar] |
|
| integer, | intent(in), | optional | :: | max_points |
Maximum number of points |
|
| real(kind=pr), | intent(in), | optional | :: | maxP |
Maximum pressure |
|
| type(EquilibriumState), | intent(in), | optional | :: | first_point | ||
| logical, | optional | :: | stability_analysis |
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| real(kind=pr), | private | :: | P | ||||
| real(kind=pr), | private | :: | S | ||||
| real(kind=pr), | private | :: | T | ||||
| real(kind=pr), | private | :: | V | ||||
| real(kind=pr), | private | :: | V_cep | ||||
| real(kind=pr), | private | :: | X0(4) | ||||
| real(kind=pr), | private, | allocatable | :: | XS(:,:) |
Full set of solved points |
||
| real(kind=pr), | private, | allocatable | :: | XS_i(:) |
Full set of solved points |
||
| real(kind=pr), | private | :: | a | ||||
| logical, | private | :: | found_cep | ||||
| integer, | private | :: | i | ||||
| real(kind=pr), | private | :: | max_P | ||||
| integer, | private | :: | npoints | ||||
| integer, | private | :: | ns | ||||
| logical, | private | :: | stab_anal | ||||
| real(kind=pr), | private | :: | u(size(z0)) |
eigen-vector |
|||
| real(kind=pr), | private | :: | u_new(size(z0)) |
eigen-vector |
|||
| real(kind=pr), | private | :: | y_cep(size(z0)) | ||||
| real(kind=pr), | private | :: | z(size(z0)) |