mirror of
https://github.com/CoolProp/CoolProp.git
synced 2026-01-21 03:48:08 -05:00
41 lines
1.3 KiB
Python
41 lines
1.3 KiB
Python
import matplotlib.pyplot as plt
|
|
import CoolProp.CoolProp as CP
|
|
import itertools, numpy as np
|
|
|
|
variables = ['T','P','H','S','U','D']
|
|
|
|
for fluid in ['Water','CO2','n-Propane']:
|
|
T = np.linspace(CP.Props(fluid,'Tmin'),CP.Props(fluid,'Tcrit'),1000)
|
|
|
|
fig = plt.figure(1, figsize=(10, 10), dpi=100)
|
|
for i, types in enumerate(itertools.combinations(variables,2)):
|
|
ax = plt.subplot(5, 3, i+1)
|
|
|
|
types = list(types)
|
|
|
|
if types[0] in ['T','P'] and types != ['T','P']:
|
|
types[0],types[1] = types[1],types[0]
|
|
|
|
xL = CP.Props(types[0],'T',T,'Q',0,fluid)
|
|
yL = CP.Props(types[1],'T',T,'Q',0,fluid)
|
|
xV = CP.Props(types[0],'T',T,'Q',1,fluid)
|
|
yV = CP.Props(types[1],'T',T,'Q',1,fluid)
|
|
xc = CP.Props(types[0],'T',CP.Props(fluid,'Tcrit'),'Q',1,fluid)
|
|
yc = CP.Props(types[1],'T',CP.Props(fluid,'Tcrit'),'Q',1,fluid)
|
|
|
|
ax.plot(xL,yL,'k')
|
|
ax.plot(xV,yV,'k')
|
|
ax.plot(xc,yc,'o')
|
|
|
|
if types[0] in ['P','D']:
|
|
ax.set_xscale('log')
|
|
if types[1] in ['P','D']:
|
|
ax.set_yscale('log')
|
|
|
|
plt.title(' '.join(types))
|
|
plt.xlabel(types[0])
|
|
plt.ylabel(types[1])
|
|
|
|
plt.tight_layout()
|
|
plt.savefig('AllFlash'+fluid+'.pdf')
|
|
plt.show() |