建筑用PCM（相变化蓄热材）的热模拟
建筑用PCM（相变化蓄热材）热的理论基础
PCM是什么
老师，我听说建筑中使用PCM，墙壁会融化吗？好像很可怕...
墙壁本身不会融化。PCM（Phase Change Material = 相变化蓄热材）是指从固体变为液体时吸收大量热量的材料。在建筑中，石蜡的一种被微胶囊化成直径5～50μm的颗粒，然后混合到石膏板或砂浆中。
微胶囊？也就是说石蜡被壳包裹住了，只在壳内部溶解和固化对吧？
完全正确。例如BASF的Micronal DS系列是将融点23°C左右的石蜡用三聚氰胺树脂壳包裹。白天当室温超过23°C时，墙内的PCM开始熔化，将太阳热作为潜热蓄存。到晚上室温下降时，它再次固化并放热。通过这个"蓄热→放热"循环，能够将室温变动抑制2～4°C。
原来如此，普通墙壁也能蓄热啊？PCM的优势是什么？
很好的问题。普通混凝土的蓄热是显热—— 通过温度上升来蓄存。这样蓄热密度很低。水的显热是每升高1°C时4.18 kJ/kg，但石蜡C18的溶解潜热约为244 kJ/kg。相比之下，同样质量的PCM仅温度变化1～2°C范围内，就能蓄存等于混凝土升温50°C才能达到的热量。而且PCM在舒适温度附近工作，直接有利于冷暖房负荷的削减。
PCM蓄放热中的全焓变化由固相显热＋潜热＋液相显热表示。
其中 $m$ 是质量，$c_s, c_l$ 是固相和液相比热，$T_m$ 是融点，$L_f$ 是溶解潜热，$T_1, T_2$ 是初始和最终温度。下表示出建筑用石蜡系PCM的代表物性值。
| 物性 | 值 | 单位 |
|---|---|---|
| 融点 $T_m$ | 21～26 | °C |
| 溶解潜热 $L_f$ | 170～250 | kJ/kg |
| 密度 $\rho$（固体） | 850～950 | kg/m³ |
| 热传导率 $k$ | 0.15～0.25 | W/(m·K) |
| 比热 $c_p$（固/液） | 1.8～2.4 | kJ/(kg·K) |
Stefan问题与移动边界
PCM溶解或固化时，固液边界会移动吧？这怎么定式化啊？
这就是Stefan问题。1891年奥地利物理学家Jozef Stefan首次定式化了这个移动边界问题。对固相和液相区域分别求解热传导方程，在边界面 $x = s(t)$ 处满足能量保存条件（Stefan条件）。一维情况下如下所示。
固相区域 ($0 < x < s(t)$):
液相区域 ($s(t) < x < L$):
Stefan条件（界面 $x = s(t)$ 处的能量平衡）:
界面位置 $s(t)$ 是未知数，还随着时间移动，这能分析求解吗？
半无限固体单侧加热这样的特殊条件有Neumann解 $s(t) = 2\lambda\sqrt{\alpha t}$ 的分析解。但像建筑壁体这样的有限厚度、两面温度变动、多材料层的问题不可能有分析解。必须用数值方法，这时候焓法就派上用场了。
焓法的支配方程
焓法是求解Stefan问题的方便方法吗？什么原理呢？
Stefan问题的麻烦之处在于"界面位置追踪"。焓法是一个天才的想法，用体积焓 $H$ 代替温度 $T$ 作为因变量。由于相变化前后 $H(T)$ 函数会在潜热处跳跃，可以用单一方程统一处理固相、液相和浆状区域（mushy zone）。不需要显式追踪界面。
焓法的支配方程可以写成：
其中体积焓 $H(T)$ 定义为温度的函数：
建筑用PCM中，相变化不是在单一温度，而是从 $T_s$ 到 $T_l$ 的温度范围（通常宽度2～4°C）。这个范围称为浆状区域，固相和液相共存。
明白了！用温度写方程的话会在界面处不连续，很困难。但用焓来写就能光滑处理，对吧。
正是。当浆状区宽度 $\Delta T = T_l - T_s$ 趋向零时，退化为理想的"尖锐界面"Stefan问题。实际的建筑用PCM具有 $\Delta T = 2\sim4$°C的有限宽度，所以焓法其实是最自然的定式化。
PCM壁体的传热机制
实际的建筑壁体不只是PCM层，还有保温材料、混凝土等各种层对吧？整体怎么建模？
是的。典型的PCM壁体结构例如"外装材（10mm）→保温材（50mm）→含PCM石膏板（12.5mm）→空气层→内装仕上"这样的多层结构。各层的一维非稳态热传导方程联立求解。
多层壁体各层 $i$ 的方程：
层间的连接条件（温度的连续性和热流保存）：
外壳和内壳的边界条件包括对流（$h_\mathrm{ext}$, $h_\mathrm{int}$）、日射吸收（$\alpha_s I_\mathrm{sol}$）、长波辐射（$\varepsilon \sigma (T^4 - T_\mathrm{sky}^4)$）等。
只有PCM层 $c_p(T)$ 会因温度依赖而成为非线性，其他层保持通常的热传导方程对吧。
对。所以计算的关键就是"PCM层是否足够精细地离散化"。保温材和混凝土可以用粗网格，但PCM层在相变温度宽度 $\Delta T$ 内焓急速变化，需要细致划分。
Stefan问题的历史
Stefan问题得名于1891年奥地利物理学家Jozef Stefan（1835-1893），他曾是玻尔兹曼的老师。他研究的是"北冰洋冰层如何增长"的问题。已知冰面温度和海水温度时，冰的厚度 $s(t)$ 与 $\sqrt{t}$ 成正比的Neumann解，直到今天仍被用作PCM模拟的验证基准问题。建筑领域的本格应用始于1980年代，由美国能源部（DOE）的Doris Boehm等人的研究开启。
建筑用PCM（相变化蓄热材）热的数值计算方法
焓法的离散化
焓法在计算机上怎么实现，怎样离散化？
建筑BES（Building Energy Simulation）工具主要用有限差分法（FDM）在厚度方向离散化壁体。以EnergyPlus的CondFD算法为例，PCM层分成 $N$ 个节点，在各节点 $j$ 求解以下隐式差分式：
$H$ 和 $T$ 都是未知数，用 $H(T)$ 的关系式连接，这样成为非线性吧？
说对了。在各时间步需要用逆函数 $T = H^{-1}(H)$ 从 $H^{n+1}$ 求出 $T^{n+1}$。具体上，以前一时间步的温度为初值估计，在 $H(T)$ 曲线上迭代更新温度。浆状区内 $T = T_s + (H - H_s)\Delta T / (\rho L_f)$ 能明示求解，通常1～3次迭代就收敛。
视比热法与比较
还听说过"视比热法"，与焓法有什么不同？
视比热法（Apparent Heat Capacity法）是将潜热表示为温度依存的等效比热 $c_\mathrm{eff}(T)$：
用视比热法就能用通常的热传导方程 $\rho c_\mathrm{eff}(T) \partial T / \partial t = \nabla \cdot (k \nabla T)$ 来解，易于在既有求解器上实现。但有大的陷阱——时间步粗的话，相变温度范围会被"跳过"，潜热可能从计算中完全消失。
那太恐怖了！就是说一步从 $T_s$ 跳到 $T_l$ 以上的话，潜热的能量就没了？
完全正确。这是视比热法最大的弱点。相反，焓法直接追踪 $H$，即便时间步稍粗，也保证了潜热的保存。实务上用法如下。
| 方法 | 优点 | 缺点 | 推荐场景 |
|---|---|---|---|
| 焓法 | 能量保存严格、稳健 | 需要 $H \to T$ 的逆变换 | BES工具（EnergyPlus等） |
| 视比热法 | 易于在既有求解器上实现 | 时间步依存、潜热遗漏风险 | CFD/FEM（COMSOL、Fluent等） |
时间步长与网格划分
PCM壁体分析适当的时间步长与网格细度大概是多少？
这个挺关键的。PCM的热扩散率 $\alpha = k/(\rho c_p) \approx 8 \times 10^{-8}$ m²/s，非常小。12.5mm厚的PCM石膏板的特性时间约 $\tau = L^2/\alpha \approx 2400$ 秒（40分钟）。基于此：
- 时间步长：1～3分钟（年间模拟实用上是3分钟，精密评价用1分钟）
- 空间划分：PCM层至少4份（$\Delta x \le 3$ mm），推荐6～8份
- Fourier数约束：$\mathrm{Fo} = \alpha \Delta t / \Delta x^2 \le 0.5$（显式格式。隐式格式无约束）
年间模拟的话 365天 $\times$ 1440分钟 = 525,600步... 3分钟间隔也要175,200步。计算量如何？
1D壁体的有限差分计算代价微小。EnergyPlus的1模型年间模拟，含PCM壁体也数十秒～数分钟完成。多整栋建筑约10分钟。COMSOL的2D/3D详细分析是另外一个量级，但除非要看墙体内自然对流或界面近处的细节，通常1D-BES足够。
非线性收敛处理
PCM层的非线性会对收敛造成问题吗？
$c_\mathrm{eff}(T)$ 在融点附近急速变化，用Newton-Raphson法求解时Jacobian条件数恶化。实务对策有3个。
- 加宽浆状区宽度：把 $\Delta T$ 设为实测值（通常2～4°C），而不是人为地缩小到0.1°C等
- 采用焓法：比视比热法更稳健
- 引入弛豫系数：$T^{n+1}_\mathrm{new} = \omega T^{n+1}_\mathrm{calc} + (1-\omega) T^n$，$\omega = 0.5\sim0.8$
建筑用PCM（相变化蓄热材）热的实务应用
EnergyPlus中的PCM设置
EnergyPlus中怎样建模PCM壁体？具体步骤是什么？
EnergyPlus（v24.1以后）用 MaterialProperty:PhaseChange 对象定义温度-焓曲线。步骤如下。
- 定义壁体构成：
Construction对象定义多层墙。PCM层作为常规Material输入基本物性 - 补充PCM物性：
MaterialProperty:PhaseChange指定8～16个温度-焓数值对 - 启用CondFD：
HeatBalanceAlgorithm设为ConductionFiniteDifference（默认的CTF算法不支持PCM） - 设置时间步：
Timestep设为20以上（= 3分钟间隔）
温度-焓曲线的输入，具体是什么数值？
以BASF Micronal DS 5040 X（融点23°C、潜热110 kJ/kg）为例：
| 温度 [°C] | 焓 [J/kg] | 说明 |
|---|---|---|
| -20 | 0 | 基准点 |
| 15 | 31,500 | 固相显热积累 |
| 20 | 40,500 | 相变开始前 |
| 21 | 55,000 | 浆状区前期 |
| 23 | 120,000 | 融点（潜热峰） |
| 25 | 155,000 | 浆状区后期 |
| 26 | 162,000 | 相变完成 |
| 50 | 210,000 | 液相显热积累 |
重要的是这条曲线必须单调递增。焓有减少的点会导致EnergyPlus不预警地输出异常结果。另外，CondFD的空间节点数由 CondFD:SpaceDiscretizationConstant 控制，默认值3对PCM层太粗。PCM含有层设为0.5以下，细致划分。
TRNSYS中的PCM模型
TRNSYS中处理PCM的方法与EnergyPlus有区别吗？
TRNSYS用Type 399（Kuznik等开发）或Type 1270（TREFFTZ模块）建模PCM壁体。特色是支持滞后。TRNSYS的PCM模块能设定融解和凝固的不同焓曲线。
滞后，是说溶解和固化的温度不同对吧？
完全正确。实际PCM中凝固比融解低1～3°C开始（过冷）。不考虑这个会把年间蓄放热量高估最多15%。特别是中间期（春秋）完整的融解-凝固循环无法成立的日子影响很大。TRNSYS能再现这个效应，精密评价时TRNSYS有利。
COMSOL/Ansys详细分析
1D-BES不能解什么问题？
需要2D/3D分析主要是3种情况。
- 热桥评价：窗框、阳台接合部、钢骨柱周边的PCM局部行为
- 散蓄热体内自然对流：地板下蓄热槽或PCM板内液相对流为主导的情况
- 微胶囊级伝热分析：胶囊-基质界面的局部温度分布、充填率最优化
COMSOL的"Heat Transfer with Phase Change"界面直接提供基于视比热法的PCM模型。Gaussian型的 $c_\mathrm{eff}(T)$ 分布自动按浆状区宽度生成，设置轻松。Ansys Fluent用Solidification/Melting模型，VOF法追踪界面。但建筑用途的年间模拟太重了，主要用于理解详细现象和验证简易模型。
年间能源评价工作流
评价PCM壁体的省能效果时，具体怎样进行？
典型工作流如下。
- 基础模型构建：制作无PCM的建筑模型，算出年间冷暖房负荷
- PCM参数研究：融点（18～28°C）、含有率（10～30wt%）、设置位置（内壁/外壁/天井）作为参数系统变化
- 气象数据选定：ASHRAE TMY3或EA气象数据。至少3个都市（酷热·温暖·寒冷）对比
- 灵敏度分析：PCM物性（特别是潜热和融点）±10%变动对冷暖房负荷的响应确认
- 经济评价：ΔkWh × 电力单价 vs. PCM追加费用 → 回收年数计算
融点的最优值随气候变吗？
变化很大。简言之，冷房负荷削减目标时融点约为冷房设定温度 - 1～2°C（日本住宅的话24～26°C），暖房负荷削减时为暖房设定温度 + 1～2°C（20～22°C）。两者同时期望时融点23°C前后为折中。东京等冷暖房都需要的气候，研究着有2种融点的级联PCM，但成本倍增，实用化还远。
| 气候区分 | 城市例 | 推荐融点 | 主要效果 | 冷房负荷削减 |
|---|---|---|---|---|
| 亚热带 | 那霸 | 25～27°C | 冷房峰值削减 | 15～25% |
| 温暖湿润 | 东京 | 23～25°C | 冷房+中间期舒适 | 10～20% |
| 温带 | 札幌 | 20～22°C | 暖房负荷削减 | 5～15% |
| 干燥高温 | 利雅得 | 26～28°C | 冷房大幅削减 | 20～30% |
沙特阿拉伯PCM壁体实证
国王阿卜杜拉科技大学（KAUST）在2015年安装了混有Micronal DS 5040 X（融点23°C）的混凝土板，进行了2年实测。在沙漠外气温超过45°C的环境中，室内最高气温降低3～5K，空调负荷削减12%。与EnergyPlus PCM模块对比，预测精度达到±0.8°C。国内竹中工务店在2018年对都内办公楼PCM天井板进行了实证试验，报告了夏季峰值电力削减8%。
建筑用PCM（相变化蓄热材）热的软件比较
BES（建筑能源模拟）工具
建筑PCM模拟能用什么工具？
大致分为两个系统：BES（Building Energy Simulation）工具和CFD/FEM工具。先看BES系。
| 工具 | PCM模型 | 方法 | 滞后 | 许可 |
|---|---|---|---|---|
| EnergyPlus | MaterialProperty:PhaseChange | 焓法（CondFD） | 非对应 | 免费（DOE开发） |
| TRNSYS | Type 399 / Type 1270 | 焓法（FDM） | 对应 | 商用 |
| ESP-r | Special material module | 视比热法 | 非对应 | 免费（OSS） |
| WUFI Plus | PCM material property | 焓法 | 对应 | 商用 |
CFD/FEM工具
CFD/FEM系呢？
| 工具 | PCM模型 | 方法 | 自然对流 | 用途 |
|---|---|---|---|---|
| COMSOL Multiphysics | Heat Transfer with Phase Change | 视比热法（Gaussian） | 对应 | 详细壁体分析、热桥 |
| Ansys Fluent | Solidification/Melting | 焓-多孔介质法 | 对应 | 散蓄热体、蓄热槽 |
| OpenFOAM | chtMultiRegionFoam + PCM扩展 | 视比热法 | 对应 | 研究用途 |
| Ansys Mechanical | Phase Change material property | 焓法 | 非对应 | 热应力连成 |
工具选择指南
最终该用哪个工具？判断基准教我。
按目的分别使用为正确答案。
- 年间省能评价 → EnergyPlus（免费、经验丰富、1D壁体足够）
- 含滞后的精密评价 → TRNSYS（商用但PCM模型高机能）
- 热桥的2D/3D分析 → COMSOL（多物理对应）
- 蓄热槽内流动分析 → Ansys Fluent 或 OpenFOAM
实务上常见EnergyPlus年间评价 → 问题部位用COMSOL详细分析的两阶段方法。
建筑用PCM主要厂商
BASF（德国）的Micronal DS系列微胶囊型产品全球最大市占。Rubitherm（德国）以石蜡系散蓄热体产品对研究机构有优势。BioPCM（Phase Change Energy Solutions公司、美国）推出大豆油衍生的生物基PCM。国内DNP大日本印刷从事PCM膜片。价格从石蜡系每kg 500～2,000円、无机盐水化物系200～800円，与断熱材相比10～30倍成本差是普及的壁垒。但考虑电力峰值削减带来的需要费削减，办公楼可能在5～10年内回收。
建筑用PCM（相变化蓄热材）热的先进研究
级联PCM与复合融点设计
PCM的先进技术往什么方向发展？
注目的是级联PCM。不同融点的多种PCM在壁体厚度方向分层堆积。例如外壁侧融点28°C的PCM，内壁侧融点21°C的PCM，夏季外侧吸收日射热，冬季内侧保持室内热。年间省能效果比单一融点PCM提高15～25%有报告。
有意思！但模拟肯定变复杂了。
原理上很简单，各层分别赋予不同的 $H(T)$ 曲线就行。EnergyPlus也能对应，各PCM层定义分别的 MaterialProperty:PhaseChange 。问题在最优化—— 融点组合、层厚比、设置位置要网罗探索，这时候遗传算法或机器学习的最优设计在研究中。
机器学习优化配置
机器学习怎样用于PCM最优化？
典型方法是进行数百～数千案例的 EnergyPlus年间模拟，建训练数据，构造代理模型（Gaussian Process Regression或随机森林）。输入变量是PCM融点、潜热、含有率、壁体内位置、气候条件。输出是年间冷暖房负荷削减量。代理模型完成后，用贝叶斯最优化瞬间探索最优参数。
最新研究有用PINN（Physics-Informed Neural Network）让神经网络学习焓方程本身的方法。这样不用数值解法，可能实时预测温度分布。不过还在研究阶段，实务降级还远。
纳米粒子强化PCM
PCM的弱点是热传导率低对吧？改善方法有吗？
很好的着眼点。石蜡的 $k = 0.2$ W/(m·K) 远低于铝（237 W/(m·K)）的千分之一。所以纳米粒子强化PCM（NePCM）在研究。碳纳米管、石墨烯纳米板、Al₂O₃纳米粒子添加1～5wt%，热传导率能提高50～300%。但有代价——纳米粒子添加会减少潜热5～15%，粘度增大充填性恶化。
模拟上，用Maxwell-Garnett理论或Bruggeman式推估NePCM的有效热传导率，作为 $k_\mathrm{eff}$ 输入。
这里 $\phi$ 是纳米粒子体积分率，$k_\mathrm{np}$ 是纳米粒子热传导率。
建筑用PCM（相变化蓄热材）热的故障应对
收敛失败与温度振荡
老师，PCM模拟常见的故障是什么？学长说"加PCM的话计算就乱"...
最多的是融点附近的温度振荡。视比热法时 $\Delta T$（浆状区宽）设太窄，$c_\mathrm{eff}$ 的峰值过锐造成数值不稳定。具体上温度每步在融点挟着振荡。
| 症状 | 原因 | 对策 |
|---|---|---|
| 融点附近温度每步振荡 | 浆状区太窄 / 时间步太粗 | $\Delta T$ 设为2°C以上，缩短时间步 |
| 潜热无效（PCM无结果） | 用CTF算法（EnergyPlus） | HeatBalanceAlgorithm改为CondFD |
| 能量平衡误差大 | 空间分割太粗 | 增加PCM层分割数（最少4份） |
| 年间蓄热量 $\ne$ 放热量 | 完全融解-凝固循环未成立 | 确认融点是否在室温范围内 |
| COMSOL"ill-conditioned matrix" | $c_\mathrm{eff}$ 急峻峰值 | 加宽Gaussian平滑宽度 |
滞后现象的忽视引起的误差
刚才提过滞后，忽视的话误差有多大？
Brij等人在Journal of Building Performance Simulation（2019年）的报告，不考虑滞后会把年间冷暖房负荷削减效果高估8～15%。特别是中间期（春秋）有问题。日较差小的日子常发生融解启动但凝固无法完成的"不完全循环"。滞后时凝固开始温度低于融解温度，到早晨PCM还部分为液体——蓄热余力消失的状态迎接新的一天。
对策是用支持滞后的工具（TRNSYS、WUFI Plus），或EnergyPlus采用保守推估（用凝固侧焓曲线）吧？
完全对。用EnergyPlus时实务处理法是用凝固侧焓曲线（保守推估）或用融解·凝固平均曲线（折中推估）。总之要从数据表里确认是否有融解和凝固两条曲线，向厂商确认是第一步。
过冷却的处理
"过冷却"是PCM温度低于融点还不固化的现象吧？这对模拟有影响？
石蜡系PCM过冷几乎没有（1°C以下）所以实用上没问题。问题在无机盐水化物系（例：CaCl₂·6H₂O、Na₂SO₄·10H₂O），可显示3～10°C过冷却。过冷大时，夜间室温降到融点以下也凝固不启动，蓄存的潜热无法放出。
模拟上的对处是，设定凝固开始温度 $T_\mathrm{nucleation}$ 低于融点。在凝固侧重新定焓曲线为从 $T_\mathrm{nucleation}$ 的放热。但核生成是概率现象，决定论模拟无法完全再现。用无机盐系PCM时选择添加了核生成剂、过冷却抑制在1°C以下的材料为实务上更合理。
那就是说材料选定阶段就要把过冷却问题解决，比在模拟上强行应对更确实吧。
完全赞同。建筑用PCM实务的成功关键是"取得正确的DSC数据"、"与厂商一致物性值"这种模拟前的模拟工作。再高级的求解器，入力物性数据质量低也是垃圾进垃圾出。PCM模拟是材料数据的质量直接决定结果信頼性的分野。
价值
更详细
错误