1 基本概念
1.1 传热方程的定义
传热方程是描述热量在空间中随时间传播的偏微分方程，用于刻画温度随位置与时间的变化规律。它通常以温度场为未知量，也可与热流密度、材料参数和热源项共同构成完整模型。由于热量传播受介质性质影响，传热方程能够反映导热、扩散以及局部能量积累等现象。
1.2 物理背景与热传导现象
在物理上，传热方程主要来源于热传导过程。温度较高区域会向温度较低区域传递能量，形成热流。若材料内部存在空间温差，热量会逐渐重新分布，最终趋向更均匀的状态。该过程不仅存在于固体中，也可用于描述液体、气体或多孔介质中的温度演化。
1.3 与扩散方程的关系
传热方程在数学结构上与扩散方程高度相似，二者都属于描述某种量从高浓度区向低浓度区扩散的模型。若将温度视为待研究的“扩散量”，则热传导可看作扩散现象的一种具体表现。因此，在很多教材与研究中，传热方程也常被归入扩散型方程的范畴。
1.4 方程类型与数学特征
传热方程通常属于抛物型偏微分方程，具有明显的时间演化特征。与椭圆型方程相比，它强调过程而非静态分布；与双曲型方程相比，它更突出平滑化效应。其数学特征一般包括初值问题、边值问题、唯一性、稳定性以及长时间衰减等内容。
2 数学表达形式
2.1 一维传热方程
2.1.1 标准形式
一维传热方程常写为 \[ \frac{\partial u}{\partial t}= \alpha \frac{\partial^2 u}{\partial x^2}, \] 其中 \(u(x,t)\) 表示温度，\(\alpha\) 为热扩散系数。该形式适用于沿单一空间方向变化的细长杆、薄片近似或理想化导热模型。
2.1.2 初始条件
为确定具体解，通常需要给定初始温度分布，例如 \[ u(x,0)=f(x). \] 初始条件表明在起始时刻各位置的温度状态，是传热问题中不可缺少的信息。
2.1.3 边界条件
一维模型还需配合边界条件，如定温、定热流或对流边界等。边界条件决定了系统与外界的热交换方式，对解的形式与演化结果具有直接影响。
2.2 二维与三维传热方程
2.2.1 坐标系下的表达
在二维或三维空间中，温度通常依赖多个坐标变量。二维情形可写为 \[ \frac{\partial u}{\partial t}=\alpha\left(\frac{\partial^2 u}{\partial x^2}+\frac{\partial^2 u}{\partial y^2}\right), \] 三维情形则进一步加入 \(z\) 方向项。该类方程适用于板材、块体以及复杂空间结构中的热分析。
2.2.2 拉普拉斯算子形式
利用拉普拉斯算子 \(\Delta\)，传热方程可统一写成 \[ \frac{\partial u}{\partial t}=\alpha \Delta u. \] 这种写法简洁而通用，便于在不同坐标系下进行分析，如直角坐标、柱坐标和球坐标等。
2.3 含热源项的传热方程
2.3.1 均匀热源
当介质内部存在均匀分布的热源时，方程右端可加入常数项或仅依赖时间的源项。此时温度变化不仅来自扩散，还来自内部持续供热，常用于电加热、化学反应放热等场景。
2.3.2 非均匀热源
若热源随位置变化，则需引入空间相关的源函数。非均匀热源会导致温度场出现局部峰值或不对称分布，常见于材料内部功率密度不均、局部发热元件或复杂边界加热问题。
2.4 非线性传热方程
2.4.1 温度依赖系数
在某些材料中，导热系数、比热容或密度会随温度变化，从而使方程呈现非线性形式。此类模型更接近真实材料行为，但求解难度也明显增加。
2.4.2 非线性边界项
当边界热交换关系不再是线性时，也会产生非线性传热方程。例如辐射散热常引入温度的高次幂项，使边界条件与温度之间形成非线性耦合。
3 数学推导
3.1 能量守恒原理
3.1.1 微元体分析
传热方程的推导通常从一个微小体积单元出发。比较单位时间内流入、流出以及内部产生的热量，就能建立温度变化与热通量之间的关系。
3.1.2 热流守恒关系
若微元体内热量的净变化等于内部能量增加，则可得到守恒方程。这个关系是从物理过程走向数学表达的核心步骤。
3.2 傅里叶定律
3.2.1 热流与温度梯度
傅里叶定律指出，热流密度与温度梯度成正比，方向指向温度降低的一侧。它是建立传热方程的基本经验定律，也是导热模型的重要基础。
3.2.2 传热系数的引入
在具体介质中，比例系数通常与材料导热性能相关。通过引入热扩散系数或导热系数，可把物理材料参数纳入方程，使模型具备定量描述能力。
3.3 从守恒律到偏微分方程
3.3.1 体积积分形式
先对微元体或控制体积进行积分，可得到热量平衡的积分表达式。该形式直观反映“输入—输出—积累”的关系，便于物理解释。
3.3.2 微分形式推导
在体积趋于无穷小的极限下，积分关系可转化为局部微分方程，最终得到传热方程的标准形式。此过程体现了连续介质假设与局部平衡思想。
4 解析解方法
4.1 分离变量法
4.1.1 特征函数展开
分离变量法通常将解写成空间部分与时间部分的乘积，再将空间解展开为特征函数组合。这种方法适用于规则区域和线性边界条件问题。
4.1.2 本征值问题
在分离变量后，空间部分往往转化为本征值问题。本征值和本征函数决定了解的模态结构，是解析解的重要组成部分。
4.2 傅里叶级数法
4.2.1 正弦级数
当边界条件为定温或零边界值时，正弦级数常用于展开初始温度分布。该方法形式简洁，适合区间问题。
4.2.2 余弦级数
在对称问题或零法向热流边界条件下，余弦级数更为自然。它能够较好地表达偶对称的温度分布。
4.3 格林函数法
4.3.1 点热源响应
格林函数可看作系统对单位点热源的响应。通过叠加原理，可由点源响应构造一般热源下的解。
4.3.2 边界影响处理
格林函数法不仅能描述源项，还可纳入边界条件的影响。通过适当构造满足边界要求的核函数，可得到更广泛的解析表达。
4.4 相似解与变量变换
4.4.1 自相似变量
某些传热问题具有尺度不变性，可通过引入自相似变量将偏微分方程化为常微分方程。这类方法常用于半无限域、热前沿传播等情形。
4.4.2 尺度分析
尺度分析用于判断不同物理量之间的量级关系，并识别主导项。它可帮助发现合适的无量纲变量与相似结构。
5 数值解方法
5.1 有限差分法
5.1.1 显式格式
显式格式直接用当前时刻的已知值计算下一时刻温度，形式简单、实现方便，但稳定性往往受时间步长限制。
5.1.2 隐式格式
隐式格式在未知时刻同时出现多个网格点变量，通常需要解线性方程组。其稳定性较好，适合较长时间模拟。
5.1.3 Crank-Nicolson格式
Crank-Nicolson格式兼具显式与隐式思想，通常具有较高精度和较好的稳定性。它在工程计算中应用广泛。
5.2 有限元法
5.2.1 弱形式
有限元法先将传热方程写成弱形式，再选择试探函数空间进行近似。该方法特别适合复杂几何和非均匀材料。
5.2.2 单元离散
在离散阶段，计算区域被划分为若干单元，局部近似在单元上拼接形成整体解。网格质量会直接影响结果精度。
5.3 有限体积法
5.3.1 通量守恒
有限体积法以控制体为核心，通过离散通量保证局部守恒特性，适合强调物理守恒的传热计算。
5.3.2 网格控制体
每个控制体对应一个离散平衡方程。该方法在流体与热耦合模拟中尤为常见，也便于处理复杂边界。
5.4 稳定性与收敛性
5.4.1 条件稳定
部分离散格式仅在满足特定步长和网格条件时才稳定。稳定性分析可避免数值振荡和发散。
5.4.2 数值耗散
数值耗散会使温度峰值被人为削弱，影响传播细节。合理选择离散格式有助于减少这类误差。
6 边界与初始条件
6.1 初始温度分布
初始条件描述系统在起始时刻的温度状态，是时间演化的出发点。不同初始分布会导致后续热扩散过程显著不同。
6.2 定温边界条件
定温边界条件规定边界温度为已知函数或常数，适用于与恒温环境接触的场景。它能直接限制边界处的温度变化。
6.3 定热流边界条件
定热流边界条件规定边界上的热通量值，常用于给定加热功率或绝热近似问题。若热流为零，则可视为绝热边界。
6.4 对流边界条件
对流边界条件把边界热交换与环境温度联系起来，通常反映表面与外界空气或流体之间的换热关系。它在工程中非常常见。
6.5 混合边界条件
混合边界条件同时包含温度与热流信息，或在不同边界上采用不同类型条件。此类设定更贴近实际装置中的复杂热交换情形。
7 典型解与性质
7.1 热扩散特征
7.1.1 温度平滑化
传热方程的一个重要特征是平滑化效应。初始温度分布中的尖峰和不连续性会随着时间推移逐渐变缓。
7.1.2 长时间行为
在无持续热源或特定边界约束下，系统往往趋向稳定状态或整体降温。长时间行为通常与边界条件和总能量有关。
7.2 最大值原理
最大值原理说明，在适当条件下，温度的内部极值受边界与初始数据控制。它是分析解的基本工具，也有助于判断数值结果是否合理。
7.3 唯一性与存在性
在常见条件下，传热方程的解通常具有唯一性，即给定初始与边界信息后解被确定。存在性则保证问题确实有可求解的数学结果。
7.4 能量衰减性质
随着时间推进，系统内部的温差通常减小，某些能量范数会表现出衰减趋势。这一性质反映了热扩散趋于均匀化的本质。
8 应用领域
8.1 工程热分析
8.1.1 金属材料冷却
在金属加工中，传热方程常用于分析淬火、退火和冷却过程中的温度变化，以评估组织转变和应力分布。
8.1.2 电子器件散热
电子器件工作时会产生局部热量，传热模型可用于设计散热片、热界面材料和冷却方案，防止温度过高影响性能。
8.2 地球与环境科学
8.2.1 地温变化
传热方程可用于研究土壤和岩层中的地温随季节与深度变化的规律，适合分析地下热环境。
8.2.2 大气与海洋温度模拟
在环境模拟中，传热思想可与其他过程共同描述大气和海洋中的温度演化。该类模型常用于研究热交换与季节性变化。
8.3 生物医学建模
8.3.1 组织热疗
在热疗技术中，传热方程用于模拟组织内部温度上升与扩散，以辅助控制治疗区域和热损伤范围。
8.3.2 生物体温分布
传热模型还可用于分析人体或动物体内温度分布，研究局部供血、散热和外部环境对体温的影响。
8.4 日常与实验场景
8.4.1 食物加热
烹饪和食品加工中，传热方程可解释加热过程中的温度渗透和熟化速度，帮助理解表面先热、内部后热的现象。
8.4.2 实验室温控系统
在恒温箱、反应器和精密仪器中，传热方程可用于设计温控系统，实现更稳定的热管理。
9 相关扩展
9.1 反传热问题
9.1.1 参数识别
反传热问题关注从观测温度反推出材料参数、热源分布或边界条件，常用于实验数据分析与工程诊断。
9.1.2 病态性分析
这类问题通常具有较强病态性，对噪声敏感，微小测量误差可能导致推断结果明显波动，因此需要正则化等方法辅助处理。
9.2 耦合传热模型
9.2.1 热弹耦合
热弹耦合模型将温度变化与材料变形联系起来，用于分析热胀冷缩导致的应力和位移问题。
9.2.2 流固耦合
在流体流动与固体传热同时发生时，二者会相互影响。流固耦合模型可用于发动机冷却、换热器和复杂设备分析。
9.3 传热与扩散的统一视角
9.3.1 数学类比
从数学上看，传热方程与许多扩散模型共享相同的结构特征，如梯度驱动、守恒形式和抛物型演化规律。
9.3.2 物理对应关系
在物理层面，热量传播、粒子扩散和某些概率演化过程都表现出类似的“由集中走向均匀”的趋势，因此可在统一框架下理解。