Critical state properties

Print the critical state properties for the fluids for which Cantera has built-in liquid/vapor equations of state.

Tags: Python thermodynamics multiphase non-ideal fluid

Critical State Properties
               Fluid      Tc [K]     Pc [Pa]          Zc
               water        647.3    2.209E+07      0.2333
            nitrogen        126.2      3.4E+06      0.2891
             methane        190.6    4.599E+06      0.2904
            hydrogen        32.94    1.284E+06      0.3013
              oxygen        154.6    5.043E+06      0.2879
      carbon dioxide        304.2    7.384E+06      0.2769
             heptane        537.7     2.62E+06      0.2972
             hfc134a        374.2    4.059E+06        0.26

import cantera as ct

fluids = {'water': ct.Water(),
          'nitrogen': ct.Nitrogen(),
          'methane': ct.Methane(),
          'hydrogen': ct.Hydrogen(),
          'oxygen': ct.Oxygen(),
          'carbon dioxide': ct.CarbonDioxide(),
          'heptane': ct.Heptane(),
          'hfc134a': ct.Hfc134a()
          }

print('Critical State Properties')
print('%20s  %10s  %10s  %10s' % ('Fluid','Tc [K]', 'Pc [Pa]', 'Zc'))
for name in fluids:
    f = fluids[name]
    tc = f.critical_temperature
    pc = f.critical_pressure
    rc = f.critical_density
    mw = f.mean_molecular_weight
    zc = pc * mw / (rc * ct.gas_constant * tc)
    print('%20s   %10.4g   %10.4G  %10.4G' % (name, tc, pc, zc))