参数设置
材料预设
边界条件
初始条件预设
0°C
100°C
Click: 熱源添加 (Click＆Drag連続追加)
理论与主要公式
中心差分离散化格式：
$$T_{i,j}^{n+1}= T_{i,j}^n + r\!\left(T_{i+1,j}^n + T_{i-1,j}^n + T_{i,j+1}^n + T_{i,j-1}^n - 4T_{i,j}^n\right)$$
60×60网格，每帧计算多个时间步（由速度倍率决定）。边界条件支持Dirichlet（温度固定）和Neumann（绝热，∂T/∂n=0）。
控制方程
二维非稳态热传导方程（傅里叶定律）：
$$\frac{\partial T}{\partial t}= \alpha\!\left(\frac{\partial^2 T}{\partial x^2}+ \frac{\partial^2 T}{\partial y^2}\right)$$
稳定条件（库朗数）: $r = \alpha\Delta t/\Delta x^2 \leq 0.25$
什么是二维瞬态热传导
🙋
这个模拟器里，温度像水波纹一样扩散开，好神奇！这背后的“二维瞬态热传导”到底是什么？
🎓
简单来说，就是热量在一个平面上，随着时间推移，从热的地方传到冷的地方的动态过程。比如你刚把一块烧红的铁板放在金属工作台上，热量就会从接触点向整个台面扩散。在我们的模拟器里，你点击或拖动鼠标添加的热源，就相当于那块“红铁板”，你可以实时看到热量（用颜色表示）如何传播开来。
🙋
诶，真的吗？那为什么不同材料扩散的速度不一样？我注意到模拟器里可以选“材料预设”。
🎓
问得好！这关键就在于一个叫“热扩散率”的参数，它决定了材料导热的本领。铝的热扩散率比钢大得多，所以铝传热更快。你试着在模拟器里把材料从“钢”切换到“铝”，保持其他条件不变，你会发现热量“跑”得更快了。这就是为什么电脑CPU散热片多用铝做的原因。
🙋
原来如此！那如果我疯狂调快“模拟速度”滑块，动画会不会出错或者“爆炸”啊？
🎓
哈哈，你问到点子上了！在实际的数值计算中，确实有“爆炸”的可能，这叫做数值发散。为了保证稳定，计算必须满足一个叫“库朗条件”的限制。在我们的模拟器里，当你把模拟速度调得过高，导致时间步长太大时，画面可能会开始剧烈抖动或出现不合理的色块，这就是在模拟数值不稳定。你可以试试看，这是一个理解CAE仿真底层限制的绝佳机会。
物理模型与关键公式
模拟器求解的核心是二维非稳态热传导方程，它描述了TemperatureT随时间t和空间(x, y)的变化关系：
$$\frac{\partial T}{\partial t}= \alpha\!\left(\frac{\partial^2 T}{\partial x^2}+ \frac{\partial^2 T}{\partial y^2}\right)$$
其中，$T$是温度(°C)，$t$是时间(s)，$\alpha$是热扩散率(m²/s)，它由材料的导热系数、密度和比热容决定。$\frac{\partial^2 T}{\partial x^2}+ \frac{\partial^2 T}{\partial y^2}$ 是温度的拉普拉斯算子，衡量空间上的温度不均匀程度。
模拟器使用“显式有限差分法”将上面的偏微分方程离散化，得到可以直接在网格上计算的递推公式：
$$T_{i,j}^{n+1}= T_{i,j}^n + r\!\left(T_{i+1,j}^n + T_{i-1,j}^n + T_{i,j+1}^n + T_{i,j-1}^n - 4T_{i,j}^n\right)$$
这里，$T_{i,j}^n$ 表示第n个时间步、网格点(i, j)的温度。$r = \alpha \Delta t / (\Delta x)^2$ 被称为库朗数。为了保证计算稳定，必须满足 $r \le 0.25$。这就是为什么“模拟速度”（影响$\Delta t$）不能无限调快的原因。
现实世界中的应用
电子设备散热设计：比如智能手机或电脑主板上的芯片（CPU/GPU）是热源，工程师使用此类模拟来设计散热片、热管和风扇的布局，防止元件过热降频或损坏。
建筑节能与暖通空调（HVAC）：模拟太阳辐射、室内热源（如人体、电器）对建筑墙体温度场的影响，用于优化保温材料厚度、窗户设计以及空调系统的送风策略，实现节能。
材料加工与焊接：在激光焊接或热处理过程中，极高的热量在极短时间内输入工件。模拟热传导可以预测材料的熔化区、热影响区大小，从而控制工艺参数，保证接头质量。
地热资源评估：分析地下岩层的瞬态热传导过程，帮助评估地热田的热储量以及热量提取的可持续性，为地热电站的选址和设计提供依据。
常见误解与注意事项
首先，请勿混淆“仿真速度”与“物理时间流逝速度”。即使通过滑块加快动画播放速度，热量实际传递的速率并不会改变。这仅仅是可视化所需的播放速率。物理时间的流逝由“时间步长Δt”和“计算步数”决定。例如，若设定Δt=0.1秒并计算1000步，则相当于观察100秒后的状态。
其次，“热扩散率α”与“热导率λ”是不同的物理量。α表征“热量扩散的难易程度”，λ表征“热量传递的难易程度”。二者关系式为 $\alpha = \lambda / (\rho c)$，其中涉及密度ρ与比热c。隔热材料虽然λ较小，但有时其ρ或c也较小，因此α可能意外地较大。在分析材料的“热特性”时，请养成将这三个物性参数（λ, ρ, c）作为整体考量的习惯。
最后，本工具的“发散”现象并非真实爆炸。当库朗数r超过0.25时计算会崩溃，但这是数值解法（显式解法）的固有局限。在实际工程中，细化网格往往需要极端缩小Δt，导致计算量剧增。正因如此，业界才开发了隐式解法、不完全隐式解法等其他稳定算法。
具体计算示例
以100×100 mm铜板为例，初始温度20°C，左上角加热源设置100°C。取Δx=2 mm，α=113×10⁻⁶ m²/s，Δt=0.1 s，则Fo=0.28（满足稳定性条件）。中心差分格式：T(i,j,n+1)=T(i,j,n)+Fo[T(i+1,j,n)+T(i-1,j,n)+T(i,j+1,n)+T(i,j-1,n)-4T(i,j,n)]。经过5分钟仿真，热源周围50 mm范围内温度可升至约60°C，T平均约32°C。