2D热扩散模拟器 返回
交互式模拟器

2D热扩散模拟器

点击/拖拽放置热源。实时体验傅里叶热方程。
自由更改材料、边界条件和速度,观察温度场的时间演变。

100×100网格有限差分法 Jet彩色图 实时计算 温度历史绘图
热场
100°C
图例
0°C
绘图模式
模拟速度
速度
x
笔刷大小
材料预设
α = 1.20 (相对值)
边界条件
播放控制
操作
温度场快照
实时统计
计算结果
最高温度 (°C)
最低温度 (°C)
平均温度 (°C)
0
模拟时间 (s)
彩色标尺
0°C25°C50°C75°C100°C
探针点温度历史(鼠标位置)
位置: (50, 50) — °C
探针
理论与数值方法

关于2D热扩散模拟器

🙋
我用模拟器点击放置热源后,看到颜色慢慢扩散开来,特别有趣!但这个扩散速度是怎么决定的呢?
🎓
就是由"热扩散率 α"这个数值决定的。你试试在左边的控制面板里点击不同的材料。选铜(α≈112)的话,颜色一下子就散开了,而选木材(α≈0.13)的话,热根本扩散不出去。铜的热扩散率比木材高约860倍!这就是为什么锅的把手要用木头或塑料做——它们的热扩散率很低,不容易烫手。
🙋
哇,我理解了!还有"边界条件"下面的诺依曼和迪利克雷,这两个词听起来很复杂。有什么区别吗?
🎓
其实很简单。迪利克雷条件就是"固定边界的温度"。比如说,把物体的四面都贴在温度0°C的冰上,边界就始终是0°C。这样热源的热就会逐渐逸出,最终达到稳定状态。诺依曼条件是"边界绝热",也就是说热无法从边界逃逸。用模拟器试试——诺依曼下,热源放在中心的话,整个区域会慢慢升温;迪利克雷下,边界会吸收热,温度不会无限上升。
🙋
我放了一个热源和一个冷源,然后等了一会儿,结果看到一条直线。这条线代表什么?
🎓
那就是"定常状态"!在这个状态下,温度分布不再变化。从数学上讲,一维的稳态热传导方程 d²T/dx² = 0 的解就是直线 T = ax + b。所以如果你的热源和冷源固定不动,放足够长的时间,温度分布最终会是直线。这条线的斜率由热源的温度和冷源的温度决定。这个原理在建筑保温、芯片散热等实际工程中都很重要。

常见问题解答

不完全是。热传导遵循傅里叶定律 $q = -k \nabla T$,描述的是热的传递过程本身;热扩散是温度场的时间演变 $\partial T/\partial t = \alpha \nabla^2 T$,其中 $\alpha = k/(\rho c_p)$ 是热扩散率。热扩散率不仅取决于热导率 $k$,还取决于密度 $\rho$ 和比热 $c_p$,所以"导热快"和"温度变化快"是两个不同的概念。
本模拟器采用显式有限差分法(FTCS法)求解,使用100×100的网格。稳定性条件为 $\alpha \Delta t / (\Delta x)^2 \leq 0.25$(二维)。空间分解能足以理解定性现象,但对于工程精密计算需要更细的网格和隐式求解法。
红色笔刷画的"热源"模拟的是那个位置的温度被强行固定在100°C(迪利克雷边界条件)。物理上相当于一个温度为100°C的恒温热源与物体接触。蓝色"冷源"表示温度固定在0°C的冷却装置。用橡皮擦除后,该点会自由演变,不再被强制固定。
本工具使用等间距网格的有限差分法(FDM),而实际CAE软件采用有限元法(FEM)或有限体积法(FVM),支持任意复杂的几何形状和非均匀网格。此外,CAE软件可以耦合流场、对流、辐射等多物理场,还能处理时间相关的瞬态问题。但基础的热传导原理(傅里叶定律)是一样的,这个工具帮助理解基本概念。

实际应用举例

电子产品散热设计:CPU和功率芯片的发热必须快速传导到基板和散热器上。用本模拟器放一个代表芯片的热源,观察热如何向周围扩散,可以直观理解为什么高功率器件需要铜或铝的散热路径,而不能用木材或陶瓷。

建筑保温与结露防止:冬季,室内暖气通过墙体向室外散热。当墙体内部温度降至露点,就会结露。用模拟器在墙的两侧分别设置热源(室内20°C)和冷源(室外-5°C),可以看到温度梯度如何分布,进而判断是否需要加保温层。

焊接和热处理:焊接时,焊缝附近会形成一个热影响区(HAZ)。用模拟器在一小区域放热源,改变材料的热扩散率(钢、铝、铜),可以观察到热影响范围的差异。这直接影响焊接工艺参数(焊接电流、速度等)的选择。

食品加工:罐头杀菌、真空包装食品加热时,需要保证容器中心达到足够高的温度。容器的中心升温慢,因为热只能通过周围的食品(低热扩散率)传导。这个模拟器可以帮助理解为什么大罐头需要更长加热时间。

常见误区和注意事项

误区1:热导率高 = 温度容易升高 正确的理解应该是:热导率高指的是热传递快,但物体的温度升高快慢还要看比热。金属的比热往往很小,所以加热时温度上升快;而水的比热大,同样加热反而温度升高慢。这就是热扩散率 $\alpha = k/(\rho c_p)$ 的含义。

误区2:模拟器上的"速度"就是真实物理时间 模拟器上方的速度滑块只是改变每帧的计算步数,屏幕右下角的"t = xxx s"才是真实模拟时间。如果想知道实际的热扩散需要多长时间,必须查看这个时间显示。

误区3:2D模拟可以代替3D** 2D模拟适用于薄板的面内传热,但厚度方向的热流被忽视了。对于球体、厚圆柱等,3D模拟才是必需的。此模拟器纯粹用于理论学习。

注意事项:模型假设 本模拟器假设材料各向同性、热性质不随温度变化、无热源项(除了笔刷放置的点)、无对流和辐射。现实中,这些假设往往不完全成立,特别是在大温差或高温条件下。

使用指南

  1. 用笔刷滑块设置笔刷大小(1~8格子),在温度图上点击或拖动放置热源(固定100°C)或冷源(固定0°C),观察颜色向外扩散
  2. 用速度滑块调节时间演变速度(1~20倍),控制傅里叶热方程 ∂T/∂t=α∇²T 的计算速度
  3. 在材料选项中切换铜(α=11.6相对)、钢(α=1.2相对)、木材(α=0.13相对),观察热扩散率差异,设定边界条件(绝热/恒温)后开始计算

具体计算示例

盘面以均匀 25°C 开始。在中央用笔刷大小 3 画一个"热源"(固定 100°C),以速度 5× 运行:铜(相对 α=11.6)热量迅速向四周扩散,短时间内整个盘面变热;改用木材(相对 α=0.13)时,由于 α 约小 90 倍,热量滞留在局部、扩散缓慢。材料导热性的差异一目了然(真实热扩散率:铜≈116×10⁻⁶ m²/s,木材≈0.13×10⁻⁶ m²/s)。

工程应用注意事项