modelCriticalPointEoS
Information
The critical constants in fluidConstants must be the critical point of the EoS itself,
where (dp/dd)@T = (d2p/dd2)@T = 0, because they are used for all region checks (two-phase or not).
They can differ from published critical values and from the reducing parameters of the EoS
(helmholtzCoefficients.T_reducing, d_reducing) and of the ancillary equations
(ancillaryCoefficients.T_crit, p_crit, d_crit).
This model determines the critical point of the EoS by two nested bisections: for a given temperature, the inflection
density with (d2p/dd2)@T = 0, and the temperature at which (dp/dd)@T at the inflection density is zero.
The results are written to criticalpoint.txt, and an assert checks that they agree with the critical constants
(relative tolerance 1e-9 for temperature and pressure, 1e-7 for density).
Components
| Type | Name | Default | Description |
|---|---|---|---|
| Medium.Temperature | T_c | critical temperature of the EoS | |
| Medium.Density | d_c | critical density of the EoS | |
| Medium.AbsolutePressure | p_c | critical pressure of the EoS |
Contents
| Name | Description |
|---|---|
 Medium | |
| inflectionDensityprotected | density where (d2p/dd2)@T = 0, by bisection |
| criticalPointprotected | temperature where the minimum of (dp/dd)@T, at the inflection density, is zero |