%cat 0Source_Citation.txt
%matplotlib inline
# %matplotlib notebook # for interactive
For high dpi displays.
%config InlineBackend.figure_format = 'retina'
This example compares pressure calculated from pytheos
and original publication for the platinum scale by Yokoo 2009.
import matplotlib.pyplot as plt
import numpy as np
from uncertainties import unumpy as unp
import pytheos as eos
eta = np.linspace(1., 0.60, 21)
print(eta)
yokoo_pt = eos.platinum.Yokoo2009()
yokoo_pt.print_equations()
yokoo_pt.print_equations()
yokoo_pt.print_parameters()
v0 = 60.37930856339099
yokoo_pt.three_r
v = v0 * (eta)
temp = 3000.
p = yokoo_pt.cal_p(v, temp * np.ones_like(v))
print('for T = ', temp)
for eta_i, p_i in zip(eta, p):
print("{0: .3f} {1: .2f}".format(eta_i, p_i))
It is alarming that even 300 K isotherm does not match with table value. The difference is 1%.
v = yokoo_pt.cal_v(p, temp * np.ones_like(p), min_strain=0.6)
print(1.-(v/v0))