热扩散率α是什么？

热扩散率α = λ/(ρc)，即导热系数除以体积热容。α越大（如铜：约1.2×10⁻⁴ m²/s），材料对温度变化的响应越快；α越小（如玻璃棉：约4×10⁻⁷ m²/s），隔热性能越好。

显式格式的稳定性条件是什么？

傅里叶数 Fo = α·Δt/Δx² 必须满足 Fo ≤ 0.5，否则数值解会发散。本模拟器会自动计算稳定的时间步长Δt。当热扩散率α或网格间距变化时，Δt会自动调整。

对流边界条件如何工作？

对流边界条件应用牛顿冷却定律：q = h(T_wall − T_∞)。对流换热系数h越大（强制对流：25–250 W/m²K），壁面温度越快趋近于环境温度T_∞；自然对流通常h ≈ 5–25 W/m²K。

稳态与非稳态热传导有什么区别？

稳态时∂T/∂t=0，温度场不随时间变化。对于固定温度边界，稳态解为线性分布；对于绝热边界，最终趋于均匀温度。非稳态分析用于捕捉热冲击、电子器件温度尖峰和过程控制等瞬态响应。

› 一维非稳态热传导模拟器返回

热分析

一维非稳态热传导模拟器

基于有限差分法实时模拟一维非稳态热传导。调节初始条件与边界条件，直观观察温度分布随时间的演变过程。

材料与参数

材料预设

热扩散率 α 1.2×10⁻⁵

初始条件

初始温度分布

边界条件

左端边界条件

左端温度 T_L 0°C

右端边界条件

右端温度 T_R 0°C

对流系数 h 10 W/m²K

中心温度 T_mid

—

最高温度 T_max

—

时刻 t

0 s

傅里叶数 Fo

—

控制方程：$\dfrac{\partial T}{\partial t}= \alpha \dfrac{\partial^2 T}{\partial x^2}$
显式格式：$T_i^{n+1}= T_i^n + Fo(T_{i+1}^n - 2T_i^n + T_{i-1}^n)$
稳定条件：$Fo = \alpha\Delta t/\Delta x^2 \leq 0.5$

网格数: 80 Δt: — s 长度 L: 1.0 m

温度分布 T(x, t)

中心点瞬态响应 T_mid(t)

什么是一维非稳态热传导

🧑‍🎓

“非稳态热传导”是什么？听起来好复杂。

🎓

简单来说，就是温度会随时间变化的热传递过程。比如，你刚把一根冰冷的铁棍一端放进火里，整根铁棍的温度并不会瞬间变热，而是热量从热端慢慢“扩散”到冷端，这个过程就是非稳态的。试着在模拟器里把初始条件选为“阶跃”，左边热右边冷，然后点击播放，你就能直观看到热量是如何从左向右传播的了。

🧑‍🎓

诶，真的吗？那旁边那个“热扩散率”的滑块是干嘛的？我调大了会怎样？

🎓

热扩散率就像材料的“导热速度”。你把它调大，相当于把铁棍换成了铜棒，热量会跑得飞快，温度分布很快就变得均匀了。在实际工程中，铜的散热片就是利用其高热扩散率快速散热的。你试试看，在模拟器里先用一个小的热扩散率（比如代表隔热材料）运行，再换成一个大的（比如代表金属），观察温度曲线变平的速度，差别会非常明显！

🧑‍🎓

原来如此！那“边界条件”里“对流”是什么意思？和“固定温度”有啥不同？

🎓

“固定温度”就像把物体端部始终泡在冰水或沸水里，温度不变。“对流”则更贴近现实，比如一根高温铁棒暴露在空气中，它的端部会与空气进行热交换，温度会逐渐接近空气温度。在模拟器里，你把右端边界条件从“绝热”改成“对流”，并设置一个较低的环境温度，你会看到热量不仅向右传导，还会从右端表面“散失”到空气中，整个降温过程会快很多。这正是汽车发动机散热器翅片的工作原理。

物理模型与关键公式

本模拟器求解的核心是一维非稳态热传导的控制方程，也称为傅里叶热传导定律的瞬态形式：

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

其中，$T$ 是温度，$t$ 是时间，$x$ 是空间位置。$\alpha$ 就是热扩散率，它由材料本身决定：$\alpha = \lambda / (\rho c_p)$，$\lambda$是导热系数，$\rho$是密度，$c_p$是比热容。这个方程描述了温度随时间的变化率与温度在空间上的曲率成正比。

对于对流边界条件，我们使用牛顿冷却定律来描述物体表面与流体之间的热交换：

$$-\lambda \frac{\partial T}{\partial x}= h (T_s - T_{\infty})$$

这里，等式左边是物体内部传导到表面的热量，$h$ 是对流换热系数，$T_s$是表面温度，$T_{\infty}$是环境流体温度。这个条件决定了边界处热量流出的快慢。

现实世界中的应用

电子设备散热：CPU芯片产生的热量通过导热材料传递到散热鳍片，这个过程就是典型的一维（或扩展为多维）非稳态热传导。工程师使用CAE模拟来优化散热片材料和结构，防止芯片过热降频。

建筑材料热工性能评估：评估墙壁在夏季室外高温和冬季室内供暖下的温度响应。模拟不同隔热材料（热扩散率不同）和空气层（对流边界）对室内温度延迟和衰减的影响，从而设计出更节能的建筑围护结构。

热处理工艺：在金属淬火过程中，炽热的工件被浸入冷却液（对流边界），工件内部温度急剧变化并产生应力。通过模拟非稳态温度场，可以预测组织转变和变形，优化冷却速度以避免零件开裂。

地源热泵系统设计：地下埋管与周围土壤进行热交换。土壤中的热传导是非稳态过程，其热扩散率直接影响换热效率。模拟不同土壤类型和运行模式下的长期温度场，是设计高效可持续地热系统的关键。

常见误解与注意事项

首先，请勿混淆“热扩散率”与“热导率”。模拟器中直接调整的“热扩散率α”决定了温度变化的“速度”。而对流条件中涉及的“热导率λ”则表征热量“传递的难易程度”本身。例如，隔热材料玻璃棉的热导率极小（不易传热），因此其热扩散率也较小，导致温度变化缓慢。若不理解这一区别，在查阅材料数据表设置参数时容易产生混淆。

其次，“对流”条件中的传热系数h取值非常关键。使用默认值进行探索虽无不可，但在实际应用中必须调研合适的数值。例如，自然对流（空气）通常为5〜25 W/(m²·K)，强制对流（风扇驱动空气）为25〜250 W/(m²·K)，水冷则高达500〜10000 W/(m²·K)，数量级差异显著。若在此处统一采用“暂设100”，将导致结果严重偏离现实。

最后，请始终意识到一维模型的局限性。本工具虽适用于棒材或厚板“厚度方向”的分析，但实际热量是在三维空间扩散的。例如，手机电路板上芯片产生的热量会沿基板平面（二维）扩散。一维模型应被视为理解“最主要传热路径”或“截面代表性行为”的“第一近似”，这一点需谨记于心。

进阶学习建议

熟悉本模拟器后，建议下一步接触“离散化”与“数值解法”的概念。计算机无法直接求解连续的微分方程，其方法是将物体分割为大量微小单元（网格），并将时间划分为细小步长进行近似计算，这种方法称为“有限差分法”。例如，将热传导方程的时间微分近似为 $$\frac{\partial T}{\partial t} \approx \frac{T_{\text{new}} - T_{\text{old}}}{\Delta t}$$ 是理解该方法的起点。掌握这一思路后，便能更深入地探讨模拟结果的精度与计算时间问题。

在数学层面，学习偏微分方程的分类可拓宽视野。本文涉及的是“抛物型”偏微分方程。与之相对，还有描述振动的“双曲型”与表征静态场的“椭圆型”。各类方程的解法和性质各异，了解其区别将为熟练运用更复杂的CAE软件奠定坚实的理论基础。

关于实践方向的后续主题，建议进阶至“二维热传导模拟”。从一维扩展到二维后，可处理热量在平面上的扩散行为、不同材料接触的“界面条件”等更贴近现实的复杂现象。理解这些内容后，解读商用专业CAE软件输出结果的能力也将大幅提升。