matlab水蒸气焓值计算_利用python对燃料电池进行科学计算：定压比热容、焓、熵与自由能...-优快云博客

本文链接：https://blog.youkuaiyun.com/weixin_33037713/article/details/112157027

本文介绍了使用Python进行燃料电池相关物质的定压比热容、焓、熵和自由能的计算。作者分享了计算方法，包括从经验公式获取氢、氧、水蒸气的数据，以及对液态水使用线性插值法处理。文章还探讨了焓和熵的温度变化关系，并提供了吉布斯自由能的计算方式。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

大家好，断更一年后的我又回来了。这次回来打算开一个天坑，写写自己正在学的东西，一来是当写笔记，二来是跟大家交流一下。

开篇建议大家配合以下几本书看：

《fuel cell system explained》重点推荐，首选

《PEM燃料电池理论与实践》学看中文的话，推荐这本

《质子交换膜燃料电池建模与MATLAB仿真》有中文，有英文，但是小错漏较多，看的时候注意

《燃料电池系统建模与优化控制》不太推荐的一本，不太适合当作入门读物

之后的内容基本都来自上述四本书，还有少量来自网络，来自网络内容我都会尽量指出来源，方便大家进一步了解。

第一篇我们先来搞清楚，氢、氧、水蒸气与液态水的定压比热容、焓、熵与自由能.

1. 定压比热容

氢、氧、水蒸气的定压比热容均可通过经验公式来计算。

至于上面的常数的具体数值，书里让我自行到NIST查找，可是我找不到，求好心人指导。

因此我只能从百度文库找到了。

原文写的范围是300~1000K，我对过273-300K几个点，发现计算值和查表值差别也不大，我就将就这用了。

import numpy as np
import math

#object index
objects = np.array(['h2', 'o2', 'vapor'])
# Row 0:Hydrogen, 1:Oxygen, 2:vapor
coefficience_data = np.array([[3.057, 2.677*10**-3, -5.81*10**-6, 5.521*10**-9, -1.812*10**-12],
                              [3.626, -1.878*10**-3, 7.055*10**-6, -6.764*10**-9, 2.156*10**-12],