3. 热电耦合控制方程推导
热电耦合问题的数学本质，是电荷输运与热量输运在同一介质中相互影响、彼此制约。要建立可求解的物理模型，必须从两个基本守恒定律出发：电荷守恒给出电流连续性，能量守恒给出热传导与焦耳热、汤姆孙热、珀尔帖热的耦合关系。本节将逐步推导稳态与瞬态下的热电耦合偏微分方程组，并明确各物理量的量纲与符号约定。
一、符号约定与量纲体系
在推导之前，先统一符号与单位制。本文采用国际单位制（SI），温度取绝对温标，电流正方向与电场正方向一致。关键物理量如下表所示。
| 符号 | 物理量 | 量纲（SI） | 说明 |
|---|---|---|---|
| φ | 电势 | V（伏特） | 标量场，φ=φ(x,t) |
| T | 绝对温度 | K（开尔文） | 标量场，T=T(x,t) |
| J | 电流密度矢量 | A/m² | 单位面积电流 |
| q | 热流密度矢量 | W/m² | 单位面积热流 |
| σ | 电导率 | S/m（A/(V·m)） | σ=1/ρ，ρ为电阻率 |
| κ | 热导率 | W/(m·K) | 傅里叶导热系数 |
| S | 塞贝克系数 | V/K | 又称热电动势率，S=dφ/dT（开路） |
| ρe | 电荷密度 | C/m³ | 瞬态电荷积累项 |
| cp | 定压比热容 | J/(kg·K) | 瞬态热惯性 |
| ρm | 质量密度 | kg/m³ | 与cp共同构成热容 |
符号约定需特别注意：塞贝克系数S的正负号取决于材料类型——P型半导体S>0，N型半导体S<0。电场E与电势的关系为E=−∇φ，电流正方向沿E方向。全文及后续图示均保持此约定。
二、电荷守恒方程
电荷守恒的积分形式为：单位时间内流出闭合面的电荷量等于该体积内电荷的减少率。写成微分形式：
∂ρ_e/∂t + ∇·J = 0在热电材料中，载流子浓度极高，电荷松弛时间约为10⁻¹⁴ s量级，远小于热过程的时间尺度（秒至分钟）。因此在绝大多数工程问题中，可忽略电荷积累项，取准稳态近似：
∇·J = 0 （稳态/准稳态电荷守恒）这一近似是热电耦合方程推导的关键简化：电场分布由瞬时边界条件决定，而温度场则按热扩散时间演化。两者时间尺度分离，正是分步求解或弱耦合迭代的依据。
三、本构关系：电流密度与热流密度
热电材料中，电流不仅由电场驱动（欧姆定律），还由温度梯度驱动（塞贝克效应）。热流不仅由温度梯度驱动（傅里叶定律），还由电流携带（珀尔帖效应）。完整本构关系为：
J = σ(−∇φ − S∇T)
q = −κ∇T + S·T·J第一式中，−∇φ为电场驱动力，−S∇T为温度梯度产生的热电动势。第二式中，−κ∇T为傅里叶热流，S·T·J为载流子携带的可逆热流（珀尔帖热流）。注意S·T·J项中T为当地绝对温度，不可用温差代替。
四、能量守恒方程
能量守恒的微分形式为：单位体积内能变化率等于热流散度的负值加上源项。对热电介质，源项包括焦耳热和汤姆孙热：
ρ_m·c_p·∂T/∂t = −∇·q + J·J/σ − τ·J·∇T其中τ=T·dS/dT为汤姆孙系数。将本构关系q=−κ∇T+S·T·J代入并展开，可得完整瞬态能量方程：
ρ_m·c_p·∂T/∂t = ∇·(κ∇T) − T·(dS/dT)·J·∇T + J·J/σ右端三项分别对应：傅里叶导热（可逆）、汤姆孙热（可逆，与电流和温度梯度点积相关）、焦耳热（不可逆，恒为正）。稳态时左端为零，得到稳态能量方程：
∇·(κ∇T) − T·(dS/dT)·J·∇T + J·J/σ = 0五、耦合方程组的完整形式与求解逻辑
将电荷守恒与能量守恒联立，得到瞬态热电耦合偏微分方程组：
∇·[σ(−∇φ − S∇T)] = 0
ρ_m·c_p·∂T/∂t = ∇·(κ∇T) − T·(dS/dT)·J·∇T + J·J/σ稳态形式则去掉时间导数项。该方程组的特点是非线性：σ、κ、S均随温度变化，焦耳热项J·J/σ与电势梯度平方相关，汤姆孙项含T与∇T的乘积。工程上通常采用有限元法离散，以电势φ和温度T为节点自由度，通过牛顿-拉夫逊迭代处理非线性。
六、边界条件与量纲一致性检查
方程组需配合边界条件才能定解。常见边界包括：电学侧给定电势（Dirichlet）或电流密度（Neumann）；热学侧给定温度或热流密度，对流边界取−κ∂T/∂n=h(T−T∞)。以典型Bi₂Te₃热电腿为例，取σ=1×10⁵ S/m、κ=1.5 W/(m·K)、S=2×10⁻⁴ V/K、腿长3 mm、截面积1 mm²，在冷端300 K、热端400 K、电流0.5 A条件下，焦耳热密度约为J²/σ≈(5×10⁵)²/10⁵≈2.5×10⁶ W/m³，与傅里叶导热项κ∇²T≈1.5×(100/0.003²)≈1.7×10⁷ W/m³相比不可忽略，说明耦合项必须保留。上述数值为教学示例假设，非实测数据。
量纲检查是验证推导正确性的重要手段：能量方程左端ρmcp∂T/∂t的量纲为(kg/m³)(J/(kg·K))(K/s)=W/m³，右端∇·(κ∇T)为(W/(m·K))(K/m²)=W/m³，J·J/σ为(A²/m⁴)/(A/(V·m))=W/m³，三项一致。若推导中出现量纲不匹配，通常意味着遗漏了T因子或S的符号处理有误。
上图将两条守恒链的推导路径与耦合关系并列呈现：左侧电荷守恒经准稳态近似得到∇·J=0，右侧能量守恒代入本构关系后展开为导热、汤姆孙热与焦耳热三项。中间的绿色双向箭头标明两个物理场的相互作用方向——电流密度J通过焦耳热进入能量方程，温度梯度∇T通过塞贝克效应进入电荷方程。底部给出完整方程组与“先电后热”的迭代求解逻辑，便于在建模软件中对照实现。