Class SystemPrLeeKeslerEos

All Implemented Interfaces:
Serializable, Cloneable, SystemInterface

public class SystemPrLeeKeslerEos extends SystemPrEos
Thermodynamic system combining the Peng-Robinson EOS with the Lee-Kesler BWR method for enthalpy/entropy/Cp departures (PR-LK).

This property package uses two separate models:

  1. VLE (phase equilibrium): standard PR76 alpha function applied to all acentric factors:
      m = 0.37464 + 1.54226 \omega - 0.26992 \omega^2
    
  2. Enthalpy / Entropy / Cp departures: Lee-Kesler BWR correlation instead of PR departure functions — hence the "PR-LK" label.

The key difference from SystemPrEos1978 is that the PR76 alpha applies to all components including heavy pseudo-fractions (ω > 0.49), whereas PR1978 uses a modified polynomial for those components. For typical natural-gas and oil systems the vapour fraction predicted by this class will differ from the PR1978 formulation for components with ω > 0.49.

Note on enthalpy: the Lee-Kesler BWR enthalpy override is currently not implemented; H/S/Cp departures are calculated using the standard PR EOS departure functions. A future enhancement should override getEnthalpy in the phase classes to use LK BWR tables.

Author:
Even Solbraa
See Also:
  • Field Details

    • serialVersionUID

      private static final long serialVersionUID
      Serialization version UID.
      See Also:
  • Constructor Details

    • SystemPrLeeKeslerEos

      public SystemPrLeeKeslerEos()
      Default constructor creating a PR-LK system at 298.15 K and 1 bara.
    • SystemPrLeeKeslerEos

      public SystemPrLeeKeslerEos(double T, double P)
      Constructs a PR-LK system at the specified temperature and pressure.
      Parameters:
      T - the temperature in Kelvin
      P - the pressure in bara (absolute pressure)
    • SystemPrLeeKeslerEos

      public SystemPrLeeKeslerEos(double T, double P, boolean checkForSolids)
      Constructs a PR-LK system with optional solid phase calculations.
      Parameters:
      T - the temperature in Kelvin
      P - the pressure in bara (absolute pressure)
      checkForSolids - set true to enable solid phase calculations
  • Method Details