一维非定常热传导模拟器 返回
热分析

一维 非定常 热传导模拟器 — 初值边界问题与显式格式

有限差分法实现一维非定常热传导实时动画。改变初值条件和边界条件,体验温度分布的变化。

材料·参数

初值条件

边界条件

计算结果
中点温度 T_mid
最高温度 T_max
时刻 t
0 s
傅里叶数 Fo
热扩散率 α
左端温度
右端温度
对流 h
网格数: 80 Δt: s 长度 L: 1.0 m
温度分布 T(x, t)
中点的温度响应 T_mid(t)
瞬间响应
理论·主要公式

$$q = -k \frac{dT}{dx}$$

傅里叶热传导定律:热流量 $q$(W/m²)是温度梯度乘以热导率 $k$(W/mK)。

$$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2}$$

一维热扩散方程:$\alpha = k/(\rho c_p)$ 是热扩散率(m²/s)。

$$T(x,\infty) = T_L + \frac{T_R - T_L}{L}x$$

定常解:定常状态下温度呈线性分布(边界条件为常数时)。

一维非定常热传导简介

🙋
什么是「非定常」热传导?与普通热传导有什么区别吗?
🎓
简单来说,非定常热传导是指温度随时间变化的热传导。例如,从冰箱中取出一根钢棒放在室温下,刚开始棒的中心很冷,但随着时间推移,整个棒都会接近室温。这个「温度随时间不断变化的过程」就是非定常分析。在这个模拟器中,你可以把材料从「铜」换成「混凝土」,会看到升温速度差异很大,这就能直观体验材料性质的影响。
🙋
那么「边界条件」该怎么选择呢?「绝热」是什么状态?
🎓
边界条件是用来指定端部的状态。例如,将左端设为「恒定温度」(高温),右端设为「绝热」,就能模拟像保温瓶那样热量不会从右端逃逸的情况。反之,如果右端选择「对流」,就能看到空气或水冷却的效果。实际产品设计中,常见的组合是发动机部件(恒定高温)与外界空气(对流)的混合。你可以在模拟器中试着改变「左端BC」和「右端BC」的组合,观察温度分布如何变化。
🙋
当选择「初值温度分布」为「中央高温」时,为什么温度会逐渐下降?热量是向端部逃散吗?
🎓
完全正确!只有中央高温的状态在自然界中是不稳定的。热量总是从高温区流向低温区,所以中央的热量会向两端「扩散」。这就是「热扩散」现象。如果材料选为「玻璃棉」,这个扩散过程会非常缓慢,你可以立即看到绝热材料的效果。通过模拟器的动画,你能清楚地看到热像波一样逐渐向外扩散的过程。

常见问题

这是因为时间步长Δt过大,不满足有限差分法的数值稳定条件(CFL条件)。请减小Δt或增大空间步长Δx,使得αΔt/(Δx)² ≤ 0.5 左右。
初值条件是指t=0时刻整个物体的温度分布(例如均匀20°C)。边界条件是指物体两端的热进出方式,可以选择恒定温度(如左端100°C)、绝热或对流等。初学者建议从简单的恒定温度开始,更容易理解其中规律。
α越大,热传播越快。例如铜的α≈1.1×10⁻⁴ m²/s,温度变化传播非常迅速;而木材的α≈1×10⁻⁷ m²/s,传播极为缓慢。α增加10倍,相同时间内温度到达的距离约增加3倍。
热传递系数h [W/m²K]表示物体表面与周围流体之间热交换的容易程度。典型值:静止空气中5~25、风扇强制对流10~100、水中100~1000。h值越大,表面温度越快接近周围温度。

实际应用

电子设备热设计:智能手机和PC的CPU是热源。模拟基板上热的传播(可用一维近似)以及筐体表面的散热(对流边界条件),能防止过热。这是设计高效冷却系统的基础。

建筑隔热评估:可以分析墙体内部温度分布随时间的变化。面对外界气温变化(非定常问题),评估室内温度的稳定程度,验证隔热材料的效果。

铸造与热处理工艺:金属铸入铸模或进行焙烧热处理时,需要预测制品内部温度随时间如何变化。材料的热扩散率$\alpha$是决定冷却速度的关键参数。

地热利用:地下一定深度保持基本恒温。通过分析埋入地中管道周围温度分布的时间变化,可优化地温热泵系统的设计。

常见误解和注意事项

首先,不要混淆「热扩散率」与「热导率」。模拟器中直接调整的「热扩散率α」决定了温度变化的「快慢」。而对流条件中出现的「热导率λ」是热的「流动难易程度」本身。例如,玻璃棉的热导率极小(隔热性好),所以热扩散率也很小,导致温度变化传播缓慢。不理解这个区别,在查阅材料数据表设置参数时容易出错。

其次,对流条件的热传递系数h值非常敏感。虽然用默认值可以体验效果,但实际应用中必须查找合适的值。例如,自然对流(空气)约5~25 W/(m²·K),强制对流(风扇吹风)约25~250 W/(m²·K),水冷却约500~10000 W/(m²·K),数量级差异极大。如果随意设h=100,会与现实相去甚远。

最后,要牢记一维模型的局限性。本工具最适合分析棒材或厚板的「厚度方向」热传导,但现实中热总是在三维空间中传播。例如,手机芯片产生的热不仅在厚度方向传导,也会在基板平面内(二维)扩散。一维模型只是「最主要热流路径」或「某个断面的代表性行为」的第一近似,不要误认为这就是完整的物理图景。

使用指南

  1. 设定材料预设或热扩散率α(m²/s):从钢0.12×10⁻⁴、铝0.97×10⁻⁴、混凝土0.3×10⁻⁶等中选择
  2. 指定左端(x=0)和右端(x=L)的边界条件:恒定温度TL、TR,或通过对流系数h(W/m²K)模拟自然对流条件
  3. 输入初值温度分布和时间步长,点击模拟开始按钮;显式格式(FTCS差分格式)对温度场u(x,t)的时间发展进行动画显示

具体计算示例

厚度L=0.1m的不锈钢板(α=4.2×10⁻⁶m²/s),初始温度100°C。左端固定在0°C(TL=0),右端固定在100°C(TR=100),取Δt=1秒、Δx=0.01m。计算结果显示,中点x=0.05m的温度从初始50°C上升到约120秒后的62°C。加入对流条件h=25W/m²K与周围环境(Ta=20°C)接触时,放热作用导致温度上升抑制约25%。

实务中的注意事项