1. 课程导入与热传导问题界定

热传导是自然界与工程系统中最常见的能量传递方式之一。从电子芯片的散热设计、航空发动机涡轮叶片的冷却，到建筑墙体的保温性能评估、锂电池包的热失控防护，工程师都需要回答同一个问题：温度在空间和时间上究竟如何分布？这正是“温度场”要描述的对象。本课程围绕热传导温度场的数值模拟展开，目标是让学习者从物理定律出发，建立可计算的离散模型，并能够判断计算结果是否可信。

一、热传导基本定律：傅里叶定律

热传导的定量描述建立在傅里叶定律之上。一维情形下，热流密度与温度梯度成正比：

q = -k · dT/dx

其中 q 为热流密度（W/m²），k 为材料导热系数（W/(m·K)），dT/dx 为温度梯度。负号表示热量沿温度降低的方向传递。三维一般形式为 q = -k∇T。这里必须明确符号约定：本课程统一取温度升高方向为梯度正方向，热流方向与之相反，后续所有公式与图示均遵循这一约定。

傅里叶定律只描述“某一点、某一时刻”的热流，要得到整个区域的温度分布，还需结合能量守恒，导出热传导微分方程。以常物性、无内热源的三维瞬态问题为例：

ρ·c_p·∂T/∂t = k·(∂²T/∂x² + ∂²T/∂y² + ∂²T/∂z²)

其中 ρ 为密度，c_p 为比热容，α = k/(ρc_p) 为热扩散率，表征温度扰动在材料中传播的快慢。方程本身并不难写，难的是边界条件与初始条件：第一类边界给定边界温度，第二类给定热流密度，第三类给定对流换热（-k∂T/∂n = h(T_s - T_f)）。工程问题的复杂性往往不在方程，而在边界。

二、温度场的概念与解析解的局限

温度场是温度在空间与时间上的分布函数 T(x, y, z, t)。按时间特性可分为稳态温度场（∂T/∂t = 0）与瞬态温度场。只有极少数几何简单、边界规则的问题（如无限大平板、长圆柱、球体的一维稳态导热）能获得解析解。一旦几何含复杂曲面、材料多层且物性随温度变化、边界为对流与辐射耦合，解析方法基本失效。这正是数值模拟的用武之地。

三、数值模拟在工程中的定位

数值模拟不是替代理论分析与实验，而是与二者构成三角验证关系。理论给出控制方程与量纲关系，实验提供真实边界与验证数据，数值模拟则在两者之间填补“复杂几何+复杂边界”的计算空白。工程中它主要承担四类任务：方案筛选（比较不同散热结构）、参数敏感性分析（导热系数、对流系数的影响）、极端工况预测（瞬态升温峰值）、以及实验难以测量的内部场重构。

需要清醒认识的是：数值模拟给出的是离散近似解，其可信度取决于网格质量、离散格式、边界条件准确性和收敛判据。一个常见的工程误区是“算出彩色云图就认为任务完成”，而忽略了网格无关性验证与能量守恒校核。

四、课程目标与学习路径

本课程的学习路径按“物理—离散—求解—验证—应用”递进组织。第一阶段建立热传导方程与定解条件的完整认识，明确符号与坐标系约定；第二阶段掌握有限差分与有限体积的离散思想，理解网格、时间步长与稳定性条件（如显式格式的 Fo ≤ 0.5）之间的关系；第三阶段进入代数方程求解与边界条件实现；第四阶段学习结果后处理与可信度评估；最后通过典型工程案例完成从问题定义到结论输出的闭环。

以一个具体教学示例说明本课程的典型工作方式（以下参数为教学假设，非实测）：一块 200 mm × 100 mm 的铝板，k = 237 W/(m·K)，ρ = 2700 kg/m³，c_p = 900 J/(kg·K)，初始温度 25 ℃，左边界恒温 100 ℃，其余边界对流换热 h = 10 W/(m²·K)、环境温度 25 ℃。若采用显式有限差分，网格 Δx = 5 mm，则热扩散率 α ≈ 9.75×10⁻⁵ m²/s，稳定性要求 Δt ≤ Δx²/(4α) ≈ 0.064 s。学习者应能自行完成这一估算，并理解为何网格加密会显著压缩允许时间步长。

学完本章，你应当能够用自己的语言回答三个问题：傅里叶定律中的负号意味着什么；为什么复杂工程问题必须借助数值方法；以及一个数值结果需要经过哪些校核才能被采信。这三个问题构成了后续所有章节的判断基准。