1. 课程导入与热电耦合建模全景

热电器件（如半导体制冷片、温差发电模块）的工作本质，是热场与电场在同一块材料内部相互交织、彼此制约的过程。电流通过热电材料时既携带电荷，也携带热量；温度梯度反过来又驱动载流子扩散，产生电动势。因此，任何只考虑单一物理场的模型都无法准确预测器件性能。本课程的目标，是带领大家从物理效应出发，逐步构建一个可计算、可验证的热电耦合模型，并最终形成从工程需求到模型交付的完整工作链条。

三种基本效应：耦合的物理源头

热电耦合的数学描述建立在三个基本效应之上，理解它们的符号约定与物理含义，是后续建模的前提。

• 塞贝克效应（Seebeck）：温度梯度驱动载流子扩散，在材料两端产生开路电压。定义 S = −ΔV/ΔT，单位为 V/K。约定温度升高方向为 ΔT 正方向，P 型与 N 型材料的塞贝克系数符号相反，这正是热电偶对可以串联升压的原因。
• 珀尔帖效应（Peltier）：电流通过异质结时在界面吸收或释放热量，热流 Q = Π·I，其中 Π = S·T 为珀尔帖系数。电流方向决定吸热还是放热，这是制冷片冷端与热端形成的直接原因。
• 汤姆孙效应（Thomson）：在同一均质材料内部，电流与温度梯度共存时产生的可逆吸放热，热流密度正比于 J·∇T，系数为汤姆孙系数 τ = T·dS/dT。它体现了塞贝克系数随温度变化带来的附加效应。

三者并非独立现象，而是同一组耦合本构方程在不同条件下的表现。工程建模时，通常以塞贝克系数 S、电导率 σ、热导率 κ 三个参数作为材料输入，珀尔帖与汤姆孙效应由它们与温度共同导出，无需单独设定。

多物理场耦合框架

从控制方程看，热电耦合由两个守恒方程联立构成：电荷守恒 ∇·J = 0，能量守恒 ∇·q = J·E − ∇·(ΠJ)。本构关系为 J = σ(E − S∇T)，q = ΠJ − κ∇T。可以看到，温度场通过 S∇T 项进入电流方程，电流场又通过焦耳热与珀尔帖热进入能量方程，形成双向耦合。下图展示了这一耦合框架以及从需求到交付的整体流程。

从工程需求到模型交付

建模不是从打开软件开始的，而是从明确工程指标开始的。以一个典型制冷应用为例：需求可能是“在热端 50 ℃、输入电流 3 A 条件下，冷端温度不高于 5 ℃，制冷量不低于 20 W”。这一句话已经隐含了边界条件、载荷和验收指标，后续所有建模动作都围绕它展开。

第一步是几何与材料准备。几何上需要建立 P 型与 N 型热电偶对、金属互连层、陶瓷基板的装配关系，注意热电偶对之间的间距与串联方式。材料上需要给出 S(T)、σ(T)、κ(T) 的温度曲线，通常来自厂商数据表或文献，建模时以分段线性插值录入，避免用室温常数近似全温区。

第二步是边界与载荷设定。电流载荷施加在串联回路端部，冷热端分别设置对流换热系数与环境温度，侧表面按绝热处理。这里最容易出错的是接触热阻与接触电阻：若忽略界面接触效应，模型算出的冷端温差往往偏大 10%～20%。工程上可先按经验值（如接触热阻 1×10⁻⁴ m²·K/W）设定，再通过验证环节反推修正。

第三步是求解与后处理。耦合求解得到温度场与电势场后，需要提取三类结果：冷端平均温度、制冷量 Q_c = ΠI − ½I²R − KΔT，以及 COP = Q_c/P_in。将这三项与需求指标逐项比对，即可判断模型是否可用。

第四步是验证与交付。验证的基准可以是厂商性能曲线，也可以是自己搭建的小型测试台实测数据。若冷端温度偏差超过 1 ℃ 或制冷量偏差超过 10%，应回到材料参数与接触效应两个环节排查。交付物通常包括：参数录入表、边界条件说明、关键结果曲线以及误差分析结论，使他人能够复现整个模型。