Cantera 内置液体/蒸汽状态方程 的流体的临界状态属性

使用python进行编程

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))

得到输出结果

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