氮气物性计算-源代码-Python-20240307.docx-完整代码
81 浏览量
2024-03-07
17:09:42
上传
评论
收藏 17KB DOCX 举报
import math
# 常数
R = 296.78 # J/kg.K，气体常数
M = 0.028016 # kg/mol，氮气的分子量
Pcr = 3.396 # MPa，临界压力
roucr = 304 # kg/m3，临界密度
tb = -195.8 # ℃，沸点
tm = -210 # ℃，熔点
T0 = 273.15 # K
p0 = 100000 # Pa
# 比容计算，m3/kg
def f(V, P, T):
A0 = 173.571
a = 0.000934109
B0 = 0.00177511
b = -0.000246645
c = 1499.14
f = ((R * T * (1 - c / (V * T ** 3))) / (V ** 2)) * (V + B0 * (1 - b / V))
f -= (1 - a / V) * A0 / V ** 2
f -= P
return f
# 比容，m3/kg
def V_N2(P, T):
vi = [0.01, 4]
WuCha = 0.00001
if f(vi[0], P, T) != 0 and f(vi[1], P, T) != 0:
vi[2] = vi[1] - f(vi[1], P, T) * (vi[1] - vi[0]) / (f(vi[1], P, T) - f(vi[0], P, T))
if f(vi[1], P, T) != 0 and f(vi[2], P, T) != 0:
vi[3] = vi[2] - f(vi[2], P, T) * (vi[2] - vi[0]) / (f(vi[2], P, T) - f(vi[0], P, T))
i = 3
while abs(f(vi[i], P, T)) >= 0.00001:
vi.append(vi[i] - f(vi[i], P, T) * (vi[i] - vi[0]) / (f(vi[i], P, T) - f(vi[0], P, T)))
i += 1
if i >= 10000:
资源评论
