Pythonで気体物性を取得する

▶︎ 気体の物性を扱えるライブラリは以前はあまりなかったのですが便利なものも出てきています。

ThermodynamicProperty.py

import pyromat as pm
import numpy as np
import matplotlib.pyplot as plt

T   = np.linspace(300.,5000.,201)
air = pm.get('ig.air')

plt.xlabel('Temperature[K]')
plt.ylabel('Cp[-]')
plt.plot(T, air.cp(T))
plt.show()

plt.xlabel('Temperature[K]')
plt.ylabel('Cv[-]')
plt.plot(T, air.cv(T))
plt.show()

plt.xlabel('Temperature[K]')
plt.ylabel('gamma=Cp/Cv[-]')
plt.ylim([1.0,1.5])
plt.plot(T, air.cp(T)/air.cv(T))
plt.show()