4. 边界条件与界面条件设定

热电耦合模型的可靠性，很大程度上取决于边界条件是否与真实工况一致。偏微分方程描述的是器件内部的普遍规律，而边界条件决定了这个规律在具体器件、具体环境中的特解。本节围绕电学边界、热学边界和热电偶臂界面连续性条件三部分展开，重点说明它们的数学表达形式以及在建模中如何落地。

一、电学边界：电流、电压与接触电阻

电学边界通常施加在热电偶臂的冷端和热端金属化层外表面。常见形式有三类：第一类为给定电势（Dirichlet 条件），如 φ = V₀，适合描述与理想电压源直接相连的端面；第二类为给定电流密度（Neumann 条件），如 J·n = J₀，其中 n 为外法向单位矢量，适合恒流源驱动场景；第三类为混合条件，通过接触电阻把端面电势与外部电路联系起来。

接触电阻是实际建模中最容易被忽略、却对结果影响显著的因素。在金属化层与热电材料之间，接触电阻 R_c 可等效为接触电阻率 ρ_c（单位 Ω·m²）除以接触面积。其数学表达为端面电势的跃变：

φ_metal − φ_te = ρ_c · (J·n)

也就是说，电流穿过界面时会产生一个与电流密度成正比的电势差。若 ρ_c 取 1×10⁻⁹ Ω·m²、工作电流密度为 1×10⁵ A/m²，则界面电势跃变约为 1×10⁻⁴ V，即 0.1 mV。单看数值不大，但当器件串联数十对偶臂、且追求高转换效率时，这部分焦耳损耗会直接压低输出电压。建模时应在金属化层与热电臂之间插入一层薄界面单元，或直接在端面施加上述跃变关系，而不是把金属与热电材料简单设为理想欧姆接触。

此外，电学边界还需注意参考电位。整个模型必须指定一个零电势参考点，否则电势场不唯一。通常把冷端金属化层外表面设为 φ = 0，热端施加电流或电压，这样得到的端电压即为器件输出电压。

二、热学边界：温度、热流、对流与辐射

热学边界比电学边界更复杂，因为实际器件同时存在传导、对流和辐射三种换热方式。常见表达如下：

• 给定温度：T = T₀，适用于与恒温冷板、恒温炉膛紧密贴合的端面，是最简单也最常用的第一类边界。
• 给定热流：−k ∂T/∂n = q₀，适用于已知输入热功率的加热面，如薄膜加热器或激光加热。
• 对流换热：−k ∂T/∂n = h(T − T_∞)，其中 h 为对流换热系数，T_∞ 为环境温度。自然对流下 h 约为 5～10 W/(m²·K)，强制风冷可达 20～100 W/(m²·K)。
• 辐射换热：−k ∂T/∂n = εσ(T⁴ − T_∞⁴)，其中 ε 为发射率，σ 为斯特藩–玻尔兹曼常数。高温工况下辐射不可忽略，例如热端 500 K、环境 300 K、ε = 0.8 时，辐射热流密度约为 3.6 kW/m²，已与中等对流水平相当。

实际建模中，端面往往是对流与辐射并存的混合边界，可写成 −k ∂T/∂n = h(T − T_∞) + εσ(T⁴ − T_∞⁴)。需要注意，辐射项对温度是四次非线性，求解器需要开启非线性迭代；若把 T⁴ 线性化为 4T_m³(T − T_m) 的等效辐射换热系数，可加快收敛，但要在结果中核对线性化误差是否可接受。

三、热电偶臂界面的连续性条件

热电偶臂由 P 型与 N 型两种材料组成，二者通过端部金属化层串联。在材料交界面上，需要同时满足电学和热学的连续性条件。对电势而言，若忽略接触电阻，则界面两侧电势连续：φ₁ = φ₂；电流密度法向分量连续：J₁·n = J₂·n。对温度而言，温度连续：T₁ = T₂；热流法向分量连续：−k₁ ∂T₁/∂n = −k₂ ∂T₂/∂n。

但真实界面存在接触热阻 R_t（单位 m²·K/W），此时温度不再连续，而是产生跃变：

T₁ − T₂ = R_t · (−k ∂T/∂n)

接触热阻与接触电阻往往同时存在，且都集中在金属化层附近几微米范围内。建模时若网格在该区域过粗，会把界面跃变“抹平”，导致算出的端电压偏低、热端温度偏高。建议在界面附近至少布置 3～5 层渐变网格，并在后处理中单独提取界面两侧的温度差和电势差，与 R_t、ρ_c 的设定值对照，验证界面条件是否被正确施加。

四、边界设定的检查要点

完成边界与界面条件设定后，建议做三项自检。第一，检查边界类型与物理工况是否匹配：恒温冷板对应第一类温度边界，恒流源对应第二类电流边界，风冷或自然对流对应第三类换热边界，不能混用。第二，检查界面参数是否与材料数据一致：ρ_c 和 R_t 应来自工艺实测或可靠文献，且单位统一为 Ω·m² 与 m²·K/W。第三，检查参考点与方向约定：零电势点、外法向方向、电流正方向应在整个模型中保持一致，避免出现同一端面既设正电流又设负电流的矛盾。

以一个典型示例说明：某热电模块热端温度设定为 450 K，冷端为 300 K，热端输入电流 2 A，偶臂截面积 4 mm²，接触电阻率取 5×10⁻⁹ Ω·m²，接触热阻取 2×10⁻⁵ m²·K/W。建模时热端施加电流边界、冷端设零电势，热端与冷端金属化层外表面分别施加恒温边界，侧面施加 h = 8 W/(m²·K) 与 ε = 0.8 的混合换热。求解后应重点观察界面两侧的电势差是否接近 ρ_c·J、温度差是否接近 R_t·q，若偏差明显，说明界面单元或网格需要调整。这些数值属于教学示例假设，实际建模应以器件工艺数据为准。

边界条件不是模型的附属设置，而是把物理方程“锚定”到真实器件上的关键环节。电学边界决定电流如何进出，热学边界决定热量如何交换，界面条件决定两种材料之间如何耦合。三者设定一致、量纲统一、方向明确，后续的耦合求解才有意义。