[デジタル開発]気体の物性をPythonで自動取得

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