什么是二维热传导FDM模拟
🙋
这个模拟器里显示的彩色云图是什么？看起来像一张温度地图。
🎓
简单来说，这就是一块平板的温度分布图。比如一块手机芯片的散热片，不同颜色代表不同温度。红色是高温，蓝色是低温。你试着拖动“上边界温度”和“下边界温度”的滑块，让它们一个高一个低，就能立刻看到热量是怎么从上往下“流”的，等温线也会跟着变密。
🙋
诶，真的吗？那中间那个“中心热源Q”是干嘛的？我把它调大，云图中间就出现一个红点。
🎓
对的！这就像在平板中心放了一个发热的电子元件，比如CPU。Q越大，发热越厉害。在实际工程中，我们最怕这种局部“热点”。你试试把热源调高，同时把四周的边界温度都设低，模拟强制散热。你会看到热量从中心红点向四周扩散的路径，等温线会变成一圈圈的同心圆。
🙋
原来如此！那下面那条“中心截面温度曲线”又怎么看呢？为什么有时候是弯的，有时候是直的？
🎓
这条曲线是从平板中间“切一刀”的温度剖面。如果是直的，说明温度均匀变化，比如只有两边温度不同。如果是弯的，尤其是中间鼓起来，就说明有内部热源在加热。工程现场常见的是，我们用它来快速判断散热是否有效。你操作一下：先把热源Q设为零，左右边界设成不同温度，曲线是斜直线；然后打开热源，曲线中间就会立刻隆起，直观显示热源的影响。
物理模型与关键公式
这个模拟器求解的核心是二维稳态热传导的泊松方程。它描述了在没有时间变化的情况下，物体内部温度分布必须满足的平衡关系：热扩散与内部产热相抵消。
$$\frac{\partial^2 T}{\partial x^2}+\frac{\partial^2 T}{\partial y^2}+\frac{Q}{k}=0$$
这里，$T$是温度(°C)，$Q$是内部热源强度(W/m³)，$k$是材料的导热系数(W/m·K)。当$Q=0$时，方程退化为拉普拉斯方程，表示没有内热源的纯导热。
为了在计算机上求解，我们使用有限差分法（FDM）进行离散。高斯-赛德尔迭代公式用于计算每个网格点$(i, j)$的温度，它用周围四个邻居点的温度平均值加上热源贡献来更新当前点温度。
$$T_{i,j}^{(new)}=\frac{T_{i+1,j}^{(old)}+T_{i-1,j}^{(new)}+T_{i,j+1}^{(old)}+T_{i,j-1}^{(new)}+\frac{Q\,\Delta x^2}{k}}{4}$$
$T_{i,j}$是网格点温度，$\Delta x$是网格间距。上标$(new)$和$(old)$代表迭代中的新值与旧值。这个公式的妙处在于，它一边算一边就用上新值，所以收敛比雅可比法快得多。
现实世界中的应用
电子设备散热设计：比如手机主板或CPU散热器。工程师用此类模拟确定芯片（热源）的位置、散热片的形状和风扇（对应边界冷却条项）的布局，防止过热导致性能下降或损坏。
建筑节能与热舒适性分析：分析墙体或玻璃幕墙的温度分布。将室外高温和室内空调温度设为边界条件，评估隔热材料的性能（导热系数k），优化能源消耗。
核反应堆堆芯监控：核燃料棒是强大的内部热源（Q极大）。模拟其温度分布对于防止堆芯熔毁至关重要，需要结合冷却剂（边界条项）的流动来确保安全。
材料加工与热处理：例如大型金属部件的淬火或焊接。模拟工件在加热炉（边界加热）和内部相变潜热（等效为热源）共同作用下的温度场，以控制材料微观结构和避免热应力开裂。
常见误解与注意事项
首先需要注意，本工具中设置的“内部热源 Q”是简化为单位面积发热量的。实际工程中通常采用单位体积发热量 [W/m³]，因此在分析厚度关键的板材时，需要将此值除以厚度等进行换算。其次，边界条件设置错误是常见问题。若本打算将所有边界设为“绝热”，却不慎将某处设为固定温度，会导致热量向该处聚集的非现实分布。例如，若将四条边全部设为“绝热”滑块，当存在内部热源Q时，整体温度将无限上升（现实中不可能发生）。这是因为不满足热平衡条件（发热量=散热量），实际设备设计中必须考虑散热路径。
最后，需要理解网格尺寸的依赖性。本工具虽固定为20×20网格，但实际上网格数量会显著影响结果精度。例如，在等温线密集区域（温度梯度陡峭处），网格过粗会导致分布比实际情况更平缓。反之，网格过细会使计算量爆炸式增长。实际工程中，确定必要且足够分辨率的“网格敏感性分析”是不可或缺的步骤。