热力学基础
热力学定律、熵、热力循环与㶲分析
热力学是研究能量、能量转换及其基本极限的科学。对于从事内燃机、燃气轮机、暖通空调系统或任何能量转换装置CAE仿真的工程师而言，热力学设定了物理上可实现的硬性上限——没有任何仿真能够违背第二定律。深入理解这些原理，能让你更准确地解读仿真结果，并设计出更高效的系统。
1. 热力学四大定律
教授，我听说热力学有四条定律，编号是0、1、2、3。为什么从零开始数？难道有人忘记编了一条？
这个问题问得很好。第零定律其实是在第一定律和第二定律之后才被明确提出的，但物理学家意识到它在逻辑上必须先于其他定律——所以给它编了第零号。它建立的是热平衡的传递性：如果物体A和C处于热平衡，物体B和C也处于热平衡，那么A和B必然也处于热平衡。这正是"温度"这个物理量能够一致定义和测量的基础。没有第零定律，温度计就没有意义了。
1.1 第零定律——热平衡
若 $A \sim C$ 且 $B \sim C$（$\sim$ 表示热平衡），则 $A \sim B$。这定义了温度作为等价类的数学含义，也是温度计能正常工作的物理基础。
1.2 第一定律——能量守恒
能量既不能凭空产生，也不能凭空消失，只能在不同形式之间转换。对于封闭系统：
对于开放系统（控制体），稳态流动能量方程为：
其中 $h = u + pv$ 为比焓，$\dot{W}_s$ 为轴功率。在大多数工程设备中，动能和势能项通常可以忽略，喷嘴和泵除外。
在燃气轮机的CFD仿真中，第一定律具体是怎么体现的？求解器好像并不是直接写 $Q-W$ 的……
在CFD中，第一定律变成了能量输运方程。求解器对每个网格单元求解总能量 $E = u + V^2/2$（内能加动能）。能量通过单元面的通量——包括压力做功项——加上粘性耗散和热传导，都被纳入计算。当你在Fluent或OpenFOAM中设置绝热壁面时，本质上就是令热通量 $q = 0$，这就是第一定律的边界条件。后处理看到的温度场，正是在质量方程、动量方程约束下求解能量方程的结果。
1.3 第二定律——不可逆性与熵
任何自发过程都使宇宙的总熵增加：
等号仅在可逆（理想化）过程中成立。该定律决定了热量不能自发地从低温流向高温，也决定了任何热机的效率都不可能达到100%。
1.4 第三定律——绝对零度
完美晶体在绝对零度时熵为零：
这为熵的计算提供了绝对基准，同时也解释了为什么绝对零度在实践中不可达到——移除最后一个量子的热能需要无穷大的功。
| 定律 | 表述 | CAE意义 |
|---|---|---|
| 第零定律 | 热平衡具有传递性 | 温度边界条件的物理基础 |
| 第一定律 | 能量守恒 | CFD/FEM中的能量方程 |
| 第二定律 | 宇宙熵总是增大 | 设定效率上限；确定过程方向 |
| 第三定律 | $T \to 0$ K 时熵趋于零 | 低温材料属性的绝对基准 |
2. 熵：第二定律深度解析
熵经常被解释为"无序度"，但这个说法感觉很模糊。在发动机仿真里，我应该怎么具体地理解熵产？
把"无序"这个说法忘掉吧——那是科普简化。在工程上，熵产与损失的功直接成正比。系统中的每一个不可逆性——流体粘性摩擦、有限温差下的传热、两种流体的混合——都会产生熵，而产生的熵就代表永久失去的有用功。在发动机里，熵产的大头来自：燃烧（不可逆化学反应）、缸壁漏热、排气放空，以及节流损失。第二定律分析能精确告诉你哪个因素对效率损失贡献最大。
2.1 克劳修斯不等式与熵
对任何循环，克劳修斯不等式为：
由此引出熵作为状态函数的定义。对于可逆过程：
对任意实际（不可逆）过程，从状态1到状态2：
2.2 体积熵产率
流体中的体积熵产率 $\dot{s}_{\text{gen}}$ 为：
第一项来自粘性耗散（速度梯度），第二项来自温度梯度下的热传导，二者均恒为非负值，从局部验证了第二定律。现代CFD后处理软件可以计算 $\dot{s}_{\text{gen}}$ 的场云图，直观显示损失集中在哪里，是强大的设计优化工具。
2.3 $T$–$s$ 图
$T$–$s$ 图是工程师分析循环的核心工具。在 $T$–$s$ 图上：
- 过程曲线下的面积 = 热量：$q = \int T\,ds$
- 可逆绝热（等熵）过程为垂直线（$ds = 0$）
- 循环热效率 = 封闭面积 / 吸热面积
- 不可逆性使过程曲线向右偏移，减小了封闭面积
3. 理想气体与实际气体模型
理想气体假设在什么情况下会失效？我在做700 bar高压储氢的仿真，同事说需要用实际气体模型。
直觉很对，这时候确实要用实际气体。理想气体假设分子体积为零、分子间无相互作用力。氢气在700 bar（车用储氢罐的典型压力）下，这两个假设都严重失效——压缩因子Z偏离1高达20%~40%。这时需要Peng-Robinson方程或Redlich-Kwong方程这类实际气体状态方程。经验准则：对比压力 $P_r = P/P_c > 0.1$ 或对比温度 $T_r = T/T_c < 2$ 时，实际气体效应显著。对比一下常温常压下的空气：$T_r \approx 10$，$P_r \approx 0.001$，理想气体完全适用。
3.1 理想气体方程
其中 $\bar{R} = 8.314$ J/(mol·K) 为通用气体常数，$M$ 为摩尔质量。对于空气：$R_{\text{空气}} = 287$ J/(kg·K)。
理想气体的内能和焓仅是温度的函数：
3.2 等熵关系
对于理想气体的可逆绝热过程：
这组关系式决定了压气机和透平的等熵效率。对于压比为 $\pi$ 的轴流压气机级：
3.3 van der Waals 与 Peng-Robinson 方程
van der Waals方程考虑了分子体积（$b$）和分子间吸引力（$a$）：
工业CFD软件（ANSYS Fluent、STAR-CCM+实际气体模型）普遍采用精度更高的Peng-Robinson方程：
4. 热力循环
我在用GT-Power仿真汽车发动机，软件把奥托循环称为"理想"参考循环。但实际发动机的P-V图跟奥托循环完全不一样，研究这个理想循环有什么意义？
理想循环给你两样东西：效率的理论上限（让你知道实际发动机离理想还有多远），以及快速参数分析用的解析公式。假设你的GT-Power仿真给出38%的热效率，而奥托理论在相同压缩比下预测57%，这19个百分点的差距精确告诉你不可逆性——壁面漏热、摩擦、非瞬时燃烧——总共损失了多少。没有理想循环作参照，你只是看到一个数字，毫无上下文。
4.1 卡诺循环——理论极限
卡诺循环在温度为 $T_H$ 和 $T_L$ 的两个热源间工作，实现了任何热机可能达到的最高效率：
对于 $T_H = 1600$ K、$T_L = 300$ K 的现代燃气轮机：$\eta_{\text{卡诺}} = 81.25\%$。实际燃气轮机效率仅达40~45%，差距代表了热力学分析可以帮助量化和减少的不可逆损失。
卡诺循环在 $T$–$s$ 图上是完美的矩形：两条等温线（水平）和两条等熵线（垂直）。
4.2 奥托循环——火花点火发动机
理想奥托循环模拟汽油发动机，由四个内部可逆过程组成：等熵压缩、等容加热、等熵膨胀、等容散热。
其中 $r_c = V_{\text{下止点}}/V_{\text{上止点}}$ 为压缩比。当 $r_c = 10$，$\gamma = 1.4$ 时：$\eta_{\text{奥托}} = 60.2\%$。实际发动机因热损失、摩擦和不完全燃烧，热效率仅为30~38%。
| 过程 | 类型 | 说明 |
|---|---|---|
| 1→2 | 等熵 | 活塞压缩 |
| 2→3 | 等容 | 燃烧（等容加热） |
| 3→4 | 等熵 | 膨胀做功冲程 |
| 4→1 | 等容 | 排气放热 |
4.3 狄塞尔循环——压燃式发动机
狄塞尔循环中，加热在等压条件下进行（活塞运动时喷油），而非等容。截止比 $r_c$ 表征等压燃烧的程度：
其中 $r = V_1/V_2$ 为压缩比（柴油机通常为14~22）。相同压缩比时，$\eta_{\text{狄塞尔}} < \eta_{\text{奥托}}$，但柴油机压缩比更高，因此实际中往往效率更佳。
4.4 布雷顿循环——燃气轮机
布雷顿循环是喷气发动机和工业燃气轮机的理论基础。过程：等熵压缩、等压加热（燃烧）、等熵膨胀（透平）、等压散热（排气）。
其中 $\Pi = P_3/P_2$ 为压比。当 $\Pi = 20$，$\gamma = 1.4$ 时：$\eta_{\text{布雷顿}} = 57.5\%$。
联合循环电站效率能超过60%，这怎么可能超过布雷顿循环的理论上限？
布雷顿循环把500~600℃的高温尾气直接排掉——这是巨大的浪费。联合循环电站在后面加了一个余热锅炉（HRSG）和朗肯蒸汽循环来回收这部分废热。整个系统工作在燃烧温度（约1600℃）和凝汽器温度（约30℃）之间，更接近卡诺效率。整体效率公式为：$\eta_{\text{联合}} = \eta_{\text{布雷顿}} + (1 - \eta_{\text{布雷顿}})\eta_{\text{朗肯}}$。现代联合循环电站能达到62~65%，不是魔法，是叠加两个循环充分利用全部温差的结果。
5. 㶲分析
㶲（也称有效能或可用能）是系统与环境死态（$T_0, P_0$）达到平衡过程中可获取的最大有用功。能量守恒但㶲因不可逆性而被破坏。
5.1 比㶲
$(h - h_0) - T_0(s - s_0)$ 项为热力学（流动）㶲。对于温度为 $T_R$ 的热源传热量 $Q$，其携带的㶲为：
这就是卡诺系数——1000℃热源的热量有77%可转化为功；相对20℃环境，50℃热源的热量仅有8%的㶲价值。
5.2 㶲损
Gouy-Stodola定理将㶲损与熵产直接联系：
这极为有用：凡是CFD仿真显示熵产率高的地方——激波、分离区、高梯度燃烧区——恰恰就是㶲损最大、性能损失最严重的地方。
5.3 第二定律效率
| 部件 | 典型 $\eta_{\text{II}}$ | 主要损失机制 |
|---|---|---|
| 现代轴流压气机级 | 88~92% | 叶片附面层粘性摩擦、叶顶间隙 |
| 燃烧室 | 65~75% | 化学反应不可逆性、散热 |
| 透平级 | 88~93% | 叶顶间隙、二次流、冷却气 |
| 工业换热器 | 50~80% | 有限温差传热 |
| 节流阀 | 0%（纯㶲损） | 无功回收的压降 |
6. CAE工程应用
我在建立燃气轮机燃烧室的CFD模型，怎么确认能量方程的物理设置是正确的？
好问题。首先检查焓的表述方式——燃烧CFD应该以总焓（静焓加动能）作为输运量，而不是静温。燃烧释放的热量必须以物种消耗速率为基础，以源项形式加入。仿真完成后做全局能量守恒：对所有边界面的热通量积分，检查其是否等于燃料热值输入减去输出功。超过1~2%的不平衡说明设置有问题。另外要验证 $c_p(T)$ 曲线的准确性——1000 K以上的燃烧产物，$c_p$ 随温度显著增大，用常数值会带来实质性误差。
6.1 可压缩流动：滞止参数
对于高速流动，CAE工程师使用滞止（总）参数来考虑动能的贡献：
滞止温度在绝热流动中守恒（即使有摩擦）。滞止压力仅在等熵流动中守恒——任何不可逆性都会降低 $P_0$。因此，总压损失是压气机、透平和管道系统中流动不可逆性的直接度量。
6.2 发动机一维热力学仿真
GT-Power、Ricardo WAVE、AVL BOOST等工具将内燃机建模为由管道连接的控制体网络。每个气缸用开放系统第一定律建模：
Wiebe函数模拟燃烧放热率（已燃质量分数）：
其中 $\theta$ 为曲轴转角，SI发动机典型参数为 $a = 5$，$m = 2$。
6.3 CFD热力学一致性检验
- 全局能量守恒： 所有壁面热通量之和 + 输出功率 = 燃料热值输入，容差 < 1%。
- 熵审计： 计算域内 $\int \dot{s}_{\text{gen}}\,dV$，验证其为正值且与第二定律预测一致。
- 绝热壁温： 高速可压缩流动中，绝热壁温 $T_{aw} = T_\infty(1 + r\frac{\gamma-1}{2}M^2)$，湍流时回收系数 $r \approx 0.85$，用此验证壁面边界条件设置。