import math
# 常数
R = 208.12 # J/kg.K，气体常数
M = 0.039944 # kg/mol，氩气的分子量
Pcr = 5 # MPa，临界压力
roucr = 536 # kg/m3，临界密度
tb = -185.86 # ℃，沸点
tm = -189.37 # ℃，熔点
h0 = -330 # J/kg 温度 315.56℃，压力 0.1013MPa 时氩气的焓
s0 = 2189.5 # J/kg.K 温度 315.56℃，压力 0.1013MPa 时氩气的熵
T0 = 273.15 # K
p0 = 100000 # Pa
# 比容计算，m3/kg
def f_vpT_Ar(V, P, T):
A0 = 81.033
a = 0.000583312
B0 = 0.000984966
b = 0
c = 1500.88
f = ((R * T * (1 - c / (V * math.pow(T, 3)))) / (math.pow(V, 2))) * (V + B0 * (1 - b / V))
f = f - (1 - a / V) * A0 / math.pow(V, 2)
return f - P
# 比容，m3/kg
def V_PT_Ar(P, T):
vi = [0.01, 2]
WuCha = 0.00001
if f_vpT_Ar(vi[0], P, T) != 0 and f_vpT_Ar(vi[1], P, T) != 0:
vi[2] = vi[1] - f_vpT_Ar(vi[1], P, T) * (vi[1] - vi[0]) / (f_vpT_Ar(vi[1], P, T) -
f_vpT_Ar(vi[0], P, T))
if f_vpT_Ar(vi[1], P, T) != 0 and f_vpT_Ar(vi[2], P, T) != 0:
vi[3] = vi[2] - f_vpT_Ar(vi[2], P, T) * (vi[2] - vi[0]) / (f_vpT_Ar(vi[2], P, T) -
f_vpT_Ar(vi[0], P, T))
i = 3
while abs(f_vpT_Ar(vi[i], P, T)) >= 0.00001:
vi.append(vi[i] - f_vpT_Ar(vi[i], P, T) * (vi[i] - vi[0]) / (f_vpT_Ar(vi[i], P, T) -
资源评论
