相变化材料（PCM）热解析模拟器

此模拟的核心是伴有相变的一维热传导方程（Stefan问题）。使用焓法可以在不显式追踪移动固液界面的情况下进行计算。

$$ \rho \frac{\partial H}{\partial t}= \frac{\partial}{\partial x}\left( k(T) \frac{\partial T}{\partial x}\right) $$

其中，$\rho$是密度[kg/m³]，$H$是焓[J/kg]，$t$是时间[s]，$k(T)$是温度相关的热导率[W/mK]，$T$是温度[K或℃]，$x$是位置[m]。焓$H$是将温度$T$和相状态（固相率）联系起来的变量，其中潜热效应已包含。

为对比，还显示简单一维半无限体Stefan问题的解析解。这是常数温度热源加热、初始温度为融点的固体时熔融前线位置$s(t)$的表达式。

$$ s(t) = 2 \lambda \sqrt{\alpha_l t} $$

其中，$s(t)$是熔融前线位置[m]，$\alpha_l$是液相热扩散率[m²/s]，$t$是时间[s]，$\lambda$是从超越方程确定的常数。从这个公式可以看出，前线位置与时间$t$的平方根成正比。这就是为什么在模拟中前线前进的速度随时间逐渐减缓。

建筑·建材(被动空调)：在墙壁或天花板中嵌入PCM，白天抑制室温上升，夜间释放热量，减少空调负荷。用模拟器调节石蜡参数可以找到最佳融点(如25℃左右)。

太阳能发电·蓄热系统：集中日光加热高温熔融盐(硝酸盐)，利用其潜热和显热驱动涡轮。在工具中选择\"硝酸盐熔融盐\"预设，设置高\"表面温度\"，可理解电厂蓄热体行为。

电子器件热管理：在发热但断续的电子部件(如CPU)中集成PCM，吸收过渡发热峰值，缓解温度上升。选择热导率高的PCM(增大k_s、k_l)是性能提升的关键。

食品低温运输·冷链保存：在冷箱中使用冰或其他有机PCM，融化时释放冷热，长时间维持恒定温度。用工具选\"冰\"并计算\"半域熔融时间\"可估算所需的冷却时间。

使用这个模拟器时，初学者容易陷入的一个大误区是\"潜热越大性能就越好\"。确实融解潜热L_f越大，蓄热容量越大，但如果热导率k低，热量无法深入材料，无法充分利用这种容量。例如，石蜡潜热大但热导率低(约0.2 W/mK)。在工具中对比k=0.1和k=0.5的情况，潜热相同但热传播速度截然不同，半域熔融时间大幅改变。实际应用中需要用翅片或金属框架来改善热导率。

另一个问题是边界条件的现实性。该工具采用\"表面温度恒定\"的简化条件。但在实际应用(如建筑墙体)中，外气温随日射变化，对流传热也不可忽视。该工具的固定\"表面温度 T_s\"模拟是第一步，真实设计还需考虑\"第三类边界条件(对流)\"。试着快速改变\"表面温度\"观察现象，体会与现实的差异。

最后是一维模型的局限性。此计算基于板状或无限长圆柱等一维热流为主的几何。实际蓄热装置有三维热传播和自然对流的影响。用工具模拟\"冰\"融化时，融化方式与真实冰块不同就反映了这一点。该工具适合基本原理理解和参数灵敏度分析，详细设计需要更高阶的仿真。