预设场景
数值参数
热扩散率 α10 mm²/s
网格间距 Δx10 mm
总长度 L200 mm
边界条件
左端温度 T_L200 °C
右端温度 T_R / T_∞25 °C
初始温度25 °C
计算结果
—
傅里叶数 r
—
Fo = αt/L²
—
T_max (°C)
0
步数 n
显式有限差分法：$T_i^{n+1}= T_i^n + r(T_{i+1}^n - 2T_i^n + T_{i-1}^n)$，稳定性条件 $r = \alpha\Delta t/\Delta x^2 \leq 0.5$。实时动画演示一维热传导过程。
这是模拟的核心——显式有限差分格式。它告诉我们，下一个时间步某一点的温度，由当前时刻该点及其左右相邻点的温度决定。
$$T_i^{n+1}= T_i^n + r(T_{i+1}^n - 2T_i^n + T_{i-1}^n)$$这里，$T_i^n$ 是第 $i$ 个网格点在第 $n$ 个时间步的温度。$r$ 是一个关键的无量纲数，叫做网格傅里叶数。
网格傅里叶数 $r$ 直接决定了计算的稳定性和物理真实性。它必须小于等于0.5，否则解会发散。
$$r = \frac{\alpha \Delta t}{\Delta x^2} \leq 0.5$$$\alpha$ 是材料的热扩散率（m²/s），$\Delta t$ 是时间步长（s），$\Delta x$ 是空间步长（m）。工程师通过调整 $\Delta t$ 和 $\Delta x$ 来控制 $r$。
汽车与航空航天：模拟发动机活塞、涡轮叶片在剧烈燃烧和冷却循环中的温度分布。使用显式有限差分法可以快速评估不同冷却通道设计的效果，防止材料因热应力而失效。
电子设备散热：分析手机CPU或显卡芯片的热量如何通过散热片和风扇散逸。Robin边界条件在这里至关重要，因为它精确描述了芯片表面与冷却气流之间的对流换热。
材料加工与焊接：预测焊接过程中，热量在金属工件中的传播以及随后的冷却速率。这直接关系到焊缝区域的微观组织和机械性能，是保证焊接质量的关键模拟。
建筑节能与地热：计算墙壁在不同季节的传热过程，或评估地埋管换热器的长期热性能。通过模拟温度场的瞬态变化，可以优化保温材料厚度和地源热泵的设计方案。
首先，请勿混淆“热扩散率α”与“热导率k”。虽然在模拟器中直接调整的是热扩散率α，但实际工程中材料参数表通常提供的是热导率k。两者关系为 $\alpha = k / (\rho c_p)$。例如，铜的热导率较高（约400 W/mK），因此α值也较大，热量扩散迅速。反之，泡沫塑料热导率较低，α值较小，热量容易积聚。设置参数时，必须同时考虑密度ρ和比热$c_p$。
其次，“网格尺寸Δx与时间步长Δt不能独立选择”。稳定性条件 $r \le 0.5$ 可通过减小Δt或增大Δx来满足。例如，将Δx从1mm细化至0.5mm虽能提高空间分辨率，但需将Δt缩减至原值的1/4，否则r值会过大。这将导致模拟相同物理时间所需的计算步数激增，大幅延长计算时间。实际工程中需时刻权衡“精度”与“计算成本”。
最后，边界条件中“罗宾条件”的对流换热系数h随工况变化显著。本工具虽可调整h=10或100等值，但实际应用中自然对流（空气中约5〜25 W/m²K）、强制对流（若有风扇约25〜250 W/m²K）与水冷（约500〜10000 W/m²K以上）的数值量级截然不同。若设置偏离实际的值，模拟结果将完全无法反映真实情况，因此务必通过文献或实验数据仔细校验。
这一维热传导模拟器背后的核心思想，正是描述“扩散现象”的偏微分方程本身。这意味着完全相同的数学模型也出现在众多其他工程领域。
例如物质扩散（化学工程）。半导体制造中掺杂剂的扩散、电池内部锂离子的迁移，均可用与热扩散方程形式完全相同的“菲克第二定律”描述——仅需将温度T替换为浓度C，热扩散率α替换为扩散系数D。边界条件也高度相似，如“表面浓度固定（狄利克雷条件）”或“表面化学反应（类罗宾条件）”。
另一领域是地下水流动力学（水文与土木工程）。含水层中的地下水流动，在不可压缩流条件下可建模为“拉普拉斯方程”或“扩散型方程”。此时温度对应水头（水压高度），热扩散率对应渗透系数。钻井抽水的状态恰似设置热源的边界条件。
此外，结构力学分支（材料工程）中亦有应用。无论是梁挠度的基本方程，还是具时间依赖性的蠕变现象简化模型，尽管物理量不同，但“某量的变化正比于该量空间曲率（二阶微分）”的结构是共通的。通过本模拟器所体会的“稳定性条件”与“边界条件影响”，将成为学习这些领域数值模拟的坚实基础。
熟悉本工具后，建议学习“隐式解法”的原理与意义。显式解法虽计算简便，但存在严格的时间步长限制（r≤0.5）。隐式解法通过求解包含未来时刻(n+1)相邻节点温度的方程组，理论上无论取多大Δt计算都不会发散，但单步计算量更大。实际工程中，常对热传导缓慢区域采用大Δt隐式解法、对变化剧烈区域采用小Δt显式解法的“显隐混合解法”等高级技术。
数学背景方面，推荐理解“偏微分方程分类”。本文涉及的热传导方程属“抛物型”，其特点是初始影响随时间平滑衰减。相比之下，波动方程（双曲型）与静电场方程（椭圆型）性质截然不同。这种差异决定了适用的数值解法与边界条件设置的根本逻辑。
实践层面的下一步，可尝试二维热传导模拟。其核心思想与一维完全一致：某点下一时刻温度由“当前自身温度”与“上下左右四点温度”决定。公式形如 $T_{i,j}^{n+1} = T_{i,j}^n + r(T_{i+1,j}^n + T_{i-1,j}^n + T_{i,j+1}^n + T_{i,j-1}^n - 4T_{i,j}^n)$ 。二维模拟能处理L型散热片效率、多热源干涉等更具实用性与趣味性的问题。