1 历史与发展
1.1 热传导理论的早期萌芽
关于热在物体中传播的认识，早期主要来自经验观察。人们很早就注意到，金属比木材更容易传热，厚壁容器比薄壁容器更能延缓温度变化。这些经验为后来的定量研究奠定了基础。随着经典力学与微积分的发展，热现象开始被视为可用数学描述的连续过程。
1.2 傅里叶的贡献
热传导方程的系统化建立，通常与傅里叶的工作密切相关。他将热流与温度梯度联系起来，并用偏微分方程刻画热量在空间中的扩散。傅里叶不仅提出了导热的基本定律，还推动了傅里叶级数在热问题中的应用，使热传导从经验学科转变为可分析、可计算的数学物理问题。
1.3 现代偏微分方程框架下的扩展
在现代数学中，热传导方程被纳入更一般的偏微分方程理论，研究重点已从单一解析表达式扩展到解的存在性、唯一性、正则性与稳定性等性质。与此同时，随着计算技术的发展，热传导模型被广泛用于复杂几何、非均匀材料和多物理耦合问题中，形成了更为完整的传热分析框架。
2 基本概念
2.1 温度场与热流密度
温度场是指物体内部各点温度随空间和时间变化的分布。热流密度则描述单位时间内通过单位面积的热量，反映热传递的方向和强弱。两者之间存在直接关联：温度差越明显，通常热流也越大。
2.2 热传导与热扩散
热传导是热量借助物质内部微观粒子碰撞和能量交换而传播的过程。热扩散强调这种传播会使温度差逐渐减小，最终趋向均匀。由此，热传导方程常被视为一种典型的扩散型方程，具有明显的“抹平”空间不均匀性的特征。
2.3 热导率与材料参数
热导率是衡量材料导热能力的重要物理量。不同材料的热导率差异很大，金属通常较高，绝热材料则较低。除热导率外，密度、比热容等参数也会影响温度变化速度，因此热传导方程中的材料参数通常共同决定系统的动态响应。
2.4 热源项与能量守恒
热源项用于表示物体内部或外部对系统持续输入或抽取的热量，例如电热、化学反应放热或辐射吸收等。将热源纳入模型后，热传导方程能够更真实地反映实际过程。其建立通常以能量守恒为基础，即一个区域内温度变化由热流输入输出与内部热生成共同决定。
3 数学表达形式
3.1 一维热传导方程
一维热传导方程通常写为温度对时间的一阶导数等于温度对空间二阶导数的某个系数倍。它适用于细长杆、薄片近似或温度只沿一个方向变化的情形，是理解热扩散机制的最基本模型。
3.2 二维热传导方程
二维形式用于描述平面区域中的温度演化，例如薄板、平板结构或横截面上的温度分布。此时热量可以在两个空间方向上传播，方程中包含两个方向上的二阶空间导数，能够更细致地反映实际几何的影响。
3.3 三维热传导方程
三维热传导方程适用于一般实体中的导热过程。它考虑三个空间方向上的温度变化，是最完整的经典导热模型之一。对于复杂机械构件、岩体或生物组织，三维描述往往更符合实际情况。
3.4 含热源项的热传导方程
当系统内部存在热生成或热吸收时，方程右端会加入热源项。这样不仅可描述纯扩散过程，还能反映持续加热、化学反应或局部能量注入对温度场的影响。该类方程在工程分析中十分常见。
3.5 非均匀介质中的热传导方程
在非均匀介质中，热导率和其他材料参数可能随空间位置变化。此时热传导方程的系数不再是常数，而成为位置相关函数。由于材料性质不均一，温度传播速度和方向会出现局部差异，模型也更接近真实材料。
4 方程的推导
4.1 基于能量守恒定律的推导
热传导方程的经典推导从任意控制体出发，比较其内部能量变化与边界热流及热源输入。若将能量增加率与净热流和内部热生成联系起来，再结合材料的热学参数，就可得到描述温度演化的微分方程。这种方法突出守恒思想，逻辑清晰。
4.2 基于傅里叶定律的推导
傅里叶定律指出，热流与温度梯度成正比，方向指向温度降低的一侧。将这一经验定律与能量守恒结合，便可推出热传导方程。该推导方式最能体现“温差驱动热流”的物理本质，也是教科书中最常见的来源之一。
4.3 微元体分析方法
微元体分析将空间中极小体积单元作为研究对象，通过比较各方向上的热流进出，建立局部平衡关系。由于微元足够小，可以用导数刻画变量变化，从而得到偏微分方程。这种方法在处理多维问题和复杂边界时尤为直观。
4.4 各向同性与各向异性材料下的形式
在各向同性材料中，导热性质与方向无关，因此方程形式相对简洁，系数通常为标量。对于各向异性材料，热量在不同方向上的传播能力不同，导热系数需用张量表示，方程结构也随之复杂。这类差异在晶体材料和复合材料中尤为明显。
5 边界条件与初始条件
5.1 初始条件
初始条件给出某一时刻，通常是起始时刻，整个区域内的温度分布。它决定了热传导问题的时间起点，是求解非稳态问题不可缺少的组成部分。没有初始条件，温度演化往往无法唯一确定。
5.2 定温边界条件
定温边界条件规定边界上的温度随时间保持已知值。它常用于与大型恒温环境接触的情况，如物体表面被保持在固定温度的热浴中。数学上，这类条件称为第一类边界条件。
5.3 定热流边界条件
定热流边界条件直接指定边界上的热通量，表示单位面积上的热输入或输出已知。它适合描述加热器、绝热近似下的局部热交换或指定功率输入的边界。此类条件对应第二类边界条件。
5.4 对流边界条件
对流边界条件用于描述固体表面与周围流体之间的热交换。边界热流与表面温度和环境温度之差成正比，常见于空气冷却、水冷或自然散热过程。它也被称为第三类边界条件。
5.5 混合边界条件
混合边界条件是多种边界类型的组合，例如不同边界段采用不同约束，或在同一边界上同时涉及温度与热流关系。这种条件在实际工程中较为普遍，尤其适用于几何形状复杂、换热机制多样的系统。
6 解析解方法
6.1 分离变量法
分离变量法假设解可以写成时间函数与空间函数的乘积，再将原方程拆分为若干常微分方程。该方法适用于边界条件较规则、区域几何较简单的问题，常用于矩形区域、圆柱区域或球对称问题的求解。
6.2 傅里叶级数法
傅里叶级数法通过将初始温度分布展开为正弦、余弦或更一般的三角函数级数，构造满足边界条件的解。它特别适合周期性或定边界问题，能够将复杂初值转化为一系列可单独演化的模态。
6.3 格林函数法
格林函数法把点热源响应作为基本解，再通过叠加原理处理一般热源和初值。它适用于线性问题，能够有效表示边界与源项的影响。由于具有较强的构造性，这一方法在理论分析和积分表示中都很重要。
6.4 特征函数展开法
特征函数展开法利用边值问题的特征函数系构造解，将温度场表示为若干基函数的线性组合。每个特征函数对应一个衰减模式，便于分析系统的长期行为。该方法与傅里叶级数法关系密切，但适用范围更广。
6.5 相似解与自相似变量
相似解方法通过引入组合变量，将偏微分方程化为常微分方程。对于某些具有尺度不变性的热传导问题，温度分布会呈现自相似结构，即不同时间下的曲线可通过适当缩放相互对应。这类解常用于半无限域和热波前问题。
7 数值求解方法
7.1 有限差分法
有限差分法用离散网格上的差商近似导数，将连续问题转化为代数方程组。该方法实现简单，适合规则网格和基础热传导模型，是早期数值模拟中的重要工具。
7.1.1 显式格式
显式格式直接利用当前时刻的已知温度计算下一时刻的值，形式直观、编程简便。其缺点是稳定性条件较严格，时间步长通常需要足够小，否则解可能发散。
7.1.2 隐式格式
隐式格式在更新时同时涉及未知时刻的多个网格点，因此通常需要求解线性方程组。它的稳定性较好，适合较大时间步长，但计算代价相对更高。
7.1.3 Crank-Nicolson格式
Crank-Nicolson格式结合了显式与隐式思想，在时间方向上取平均处理，兼具较高精度和较好稳定性。它常被视为热传导数值求解中的经典格式，尤其适用于平滑解问题。
7.2 有限元法
有限元法通过将计算区域划分为若干单元，在局部基函数上建立近似解。它对复杂几何和非均匀材料具有较强适应性，已成为工程传热分析中最常用的方法之一。
7.3 有限体积法
有限体积法以守恒控制体为基础，直接离散热量通量，因此在保持能量守恒方面表现突出。该方法广泛应用于流体与传热耦合问题，也常见于计算工程软件中。
7.4 稳定性与收敛性分析
数值方法是否可靠，取决于稳定性与收敛性。稳定性反映误差是否被放大，收敛性则说明网格加密后近似解是否趋近真实解。对于热传导问题，这两项分析通常与时间步长、空间网格和离散格式密切相关。
7.5 计算误差与网格选择
离散化会引入截断误差、舍入误差和建模误差。合理选择网格大小与时间步长，既能提高精度，又能控制计算成本。实际应用中常需在局部加密与整体效率之间取得平衡。
8 物理性质与数学性质
8.1 扩散性与平滑性
热传导方程最显著的特征是扩散性，即初始的不均匀温度会随时间逐渐被“摊平”。从数学上看，解通常具有平滑效应，随着时间推进，温度场会比初值更光滑。
8.2 最大值原理
最大值原理说明，在适当条件下，热传导方程的解在内部不会超过边界或初始数据给定的极值范围。这一性质对判断解的合理性和建立先验估计非常有用。
8.3 唯一性定理
唯一性定理表明，在给定初始条件和边界条件后，热传导问题的解通常是唯一的。也就是说，若有两个满足相同条件的解，它们实际上应当一致。这为模型的物理可解释性提供了保障。
8.4 热核与基本解
热核是热传导方程在无界空间中的基本传播函数，也可理解为点热源随时间扩散后的响应。它在构造一般解、研究卷积表示以及分析局部行为方面具有核心作用。
8.5 长时间行为与稳态极限
随着时间增加，许多热传导系统会逐渐趋向稳态，即温度不再随时间变化。若存在稳定热源或边界维持不变，系统往往收敛到对应的稳态分布；若没有持续输入，温度差则可能最终消失，整体接近均匀状态。
9 特殊情形与推广
9.1 非线性热传导方程
当热导率依赖温度、材料参数随状态变化或存在强耦合效应时，热传导方程会呈现非线性特征。非线性问题通常难以获得解析解，但更能反映高温、强辐射或复杂材料中的实际行为。
9.2 含相变的热传导模型
在材料发生熔化、凝固等相变时，除了温度变化，还涉及潜热吸收或释放。这类模型需要同时描述温度场与相界面移动，因此比普通热传导问题更复杂，常见于铸造、冷冻和热储能研究。
9.3 多孔介质中的热传导
多孔介质内部同时存在固体骨架与孔隙流体，热量传播往往受到孔结构、饱和度和流体运动影响。相应模型需要考虑有效热导率及局部热交换，在土壤、岩石和过滤材料中应用广泛。
9.4 各向异性导热问题
各向异性导热指材料在不同方向上的传热能力不同，常见于层状材料、晶体或复合结构。此时温度等值线会呈现方向性拉伸或扭曲，分析时需采用张量形式的导热系数。
9.5 分数阶热传导方程
分数阶热传导方程引入非整数阶导数，用于描述具有记忆效应或异常扩散特征的系统。与经典模型相比，它能够更好地刻画某些复杂介质中的慢扩散或非局域传热现象，属于近年来较活跃的研究方向。
10 应用领域
10.1 材料科学
在材料科学中，热传导方程用于分析材料的导热性能、热稳定性和温度响应。它可帮助研究热处理过程、微结构变化及复合材料中的热分布，对材料设计具有指导意义。
10.2 土木与建筑工程
在土木与建筑工程中，热传导模型用于评估墙体保温、屋顶隔热、混凝土养护及火灾升温等问题。通过计算温度场，可以优化结构安全性与能源效率。
10.3 机械与制造工程
机械与制造领域常利用热传导方程研究切削发热、焊接温度场、铸件冷却和电子器件散热。温度分布直接影响材料性能、加工精度和设备寿命，因此该方程具有重要工程价值。
10.4 地球物理与环境模拟
在地球物理和环境研究中，热传导用于描述地壳温度分布、地下热流以及地表热交换。它也可用于模拟土壤温度、冻土层变化和污染热扩散等环境过程。
10.5 生物热传递
生物热传递研究人体或生物组织中的热量传播，涉及血液灌注、代谢产热和外部加热等因素。热传导方程在此基础上可构成更复杂的生物热模型，用于医学热疗、组织冷却和体表热分析。
11 相关概念与对比
11.1 波动方程
波动方程描述的是振动或扰动以有限速度传播的过程，解通常呈现波峰和波谷的传播特征。与热传导方程不同，它不具有同样明显的平滑效应，且信息传播机制更偏向惯性驱动。
11.2 拉普拉斯方程
拉普拉斯方程是稳态情形下的典型椭圆型方程，用于描述时间不变的温度分布、电势场等。若热传导系统达到稳态，原来的热方程通常会退化为拉普拉斯方程或其带源形式。
11.3 对流扩散方程
对流扩散方程在热扩散基础上加入了流体运动导致的对流项，用于描述温度随流场迁移和扩散的过程。它比纯热传导方程更适合流动介质中的换热问题。
11.4 传热学中的其他控制方程
除热传导方程外，传热学还涉及辐射传热方程、对流换热关联式以及与流体动力学耦合的能量方程等。不同方程适用于不同物理机制，实际工程常需综合多种模型共同分析。