南瓜丶奇迹师上周一个做相变储热的朋友发来一段C++代码，说他的仿真程序在处理相变材料充放电过程时，温度曲线在相变点附近总是像抽风一样上下震荡，不管怎么加密时间步长都压不下去。我让他别急着改代码，先把"潜热处理"那一段的算法逻辑用伪代码完整写出来看看。半天后他把伪代码发过来，问题一眼就看明白了：液相分数的更新把当前时刻温度和上一时刻温度混在一起用，所谓的内部迭代根本没有真正迭代起来。这已经不是第一次有人栽在这种问题上——相变潜热处理本身就是一个强非线性、多物理场耦合的数值问题，物理逻辑一旦在代码和公式之间"翻译"错了，后面所有的调参都是在浪费时间。
这篇文章就从那次debug聊起。我会把"伪代码示意相变潜热处理"这件事拆开讲透：相变潜热的物理本质到底是什么、焓-孔隙率法这套主流数值框架在计算层面经历了哪几步、如何在纸面上提前规避数值振荡，以及伪代码作为一种算法设计语言，具体该怎么组织才不至于把读者绕晕。不管你是做CFD仿真的工程师、学传热学的研究生，还是刚接触数值编程、想搞懂"算法的伪代码怎么写"的开发者，这篇文章都可以直接当参考。
1. 为什么相变潜热问题要先写伪代码——一次debug给我的教训
先复盘一下朋友那个案例。他的程序本身写得很规整，类封装、参数列表、注释都不算差，但问题恰恰出在"代码太具体"上：几百行的循环里夹杂着物理更新、边界处理、数值求解和输出逻辑，时间步进的外循环和相变源项的内迭代搅在一起，我很难一眼看出他到底用的是哪种时间离散、哪种迭代策略。后来让他把伪代码写出来，压缩到三十行左右，所有问题立刻暴露了——f_l（液相分数）在第4步用新温度更新，但第2步组装系数矩阵时用的还是更新前的旧值，中间没有任何迭代路径能把这个滞后的信息补回来。
这就是伪代码的第一个价值：它强迫你只保留算法的骨架，砍掉语言层面的枝叶。当你只写T_new ← TDMA(A, b)而不是去纠结std::vector的迭代器怎么传参时，判断一个算法是否逻辑自洽会容易得多。尤其对相变潜热这种问题，它不像普通的稳态导热，写个线性方程求解器就能完事；相变过程自带移动边界、潜热释放和强源项非线性，每一步时间推进里至少嵌套了三层物理过程：
- 温度场决定液相分数（相变是否发生、发生多少）；
- 液相分数反过来通过源项影响温度场；
- 潜热释放又会改变温度场的分布，进而改变下一时刻的液相分数。
这种相互耦合如果不用伪代码把"谁先算、谁后算、循环几层、收敛判断用什么量"的规定清楚，直接在真实代码里调，很容易陷入"改一个参数、另一个地方崩掉"的泥潭。伪代码相当于把算法从"机器要执行的指令序列"翻译回"人脑能推演的物理-数学逻辑"，在动手编码之前先拿它在纸面上做一次推演。
还有一个常被忽略的好处：伪代码是团队评审的最佳载体。真实代码里有大量的工程细节，别人很难在半小时内抓住核心算法；但一份写得好的伪代码，组里任何懂物理的人都能直接review。我身边很多做储能仿真和金属凝固模拟的团队，现在都要求论文复现或模块开发之前必须提交一份伪代码评审稿，目的就是让算法逻辑在投入大量编码时间之前，先经过一轮"物理正确性"检查。
所以这里我给出的建议是：对于相变潜热处理这类涉及多场耦合的问题，伪代码不是给老师交作业用的形式主义，而是一个真正的设计工具。下文所有内容都会围绕"如何设计一份能落地的相变潜热处理伪代码"展开。
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 相变潜热的物理与数值本质：不搞懂这两点，伪代码无从下手
如果说伪代码是表达手段，那么物理本质就是你要表达的内容。如果对相变过程中的能量方程理解不到位，写出来的伪代码再漂亮也是空中楼阁。
2.1 潜热到底"潜"在哪里
先做一个生活类比的铺垫。烧一壶水，你会发现一个非常有意思的现象：水从20℃加热到100℃只需要几分钟，但从100℃开始沸腾到整壶水烧干，却要花十几分钟，而且整个沸腾阶段温度一直钉在100℃不动。你持续输入的热量去哪儿了？答案是它们被水分子"吸收"去完成液—气相变，这个过程中温度不升高，能量以分子间势能的形式储存起来，这就是潜热（Latent Heat）。固液相变同理：冰在0℃融化时，温度也不会上升，直到所有冰都变成水，温度才会继续往上走。
在数值模拟里，这个问题比"温度不变"要复杂得多。对于一个纯物质，相变发生在单一温度点，固液界面是数学上的一个间断面；但对于合金、石蜡、盐类水合物等工程常用的相变材料，相变发生在一个温度区间（固相线温度 T_s 到液相线温度 T_l 之间），在这个区间内材料处于"糊状区"（mushy zone），固体骨架和液体共存。数值模拟必须回答三个问题：每个网格现在处于什么相态？相变过程中能量怎么分配？相变前沿在哪、以什么速度移动？
这三个问题在数学上被归结为著名的Stefan问题——一个带移动边界的热传导方程。移动边界的麻烦之处在于它不是一个固定的几何位置，而是由温度场本身决定的一条等温线（或等温面），并且边界上需要满足能量守恒：界面吸收/释放的潜热等于两侧热流之差。如果直接追踪这条界面，网格会不断变形或重构，计算量非常大，而且在三维复杂几何里几乎不可行。
2.2 等效热容法的局限
早期工程上常用等效热容法（Apparent Heat Capacity Method）回避移动边界追踪。思路很简单：既然相变在温度区间 T_s~T_l 内发生，那我就在这个区间内把比热容人为放大，把潜热"摊到"这个温度区间里：
c_eff = c_p + L / (T_l - T_s)（当 T_s < T < T_l）
这样方程形式上还是一个普通的导热方程，只是物性参数分区赋值。这个方法实现简单，但有两个致命弱点。第