4. 电磁热耦合机制与多物理场接口

电磁热耦合的本质,是把电磁场求解得到的损耗密度,作为热源项代入传热方程。理解这条“能量转换链”,是正确建立 COMSOL 多物理场耦合节点的前提。本节先梳理三类损耗的数学表达,再说明耦合节点的建立方式与双向耦合策略。

一、三类电磁损耗的数学表达

在频域电磁场中,总电磁损耗密度通常由三部分构成,它们对应不同的物理机制,也对应 COMSOL 中不同的输出变量。

1. 电阻损耗(传导电流损耗)

由传导电流密度 J 与电场 E 共同决定,单位体积的时均损耗密度为:

Q_rh = 1/2 * real(J · conj(E))     [W/m^3]

其中 J = σE + jωD + J_ext
若忽略位移电流与外部电流源:
Q_rh = 1/2 * σ * |E|^2 = 1/2 * |J|^2 / σ

在 COMSOL 中,该量对应 ec.Qrh(电磁热源,电阻性)。对于母线、线圈、导体等以传导电流为主的部件,这一项是主要热源。

2. 磁损耗(磁滞与涡流损耗)

磁损耗包含磁滞损耗与涡流损耗。在频域磁场接口中,涡流损耗密度为:

Q_ml = 1/2 * real(jωB · conj(H))     [W/m^3]

对线性各向同性材料:
Q_ml = 1/2 * ω * μ'' * |H|^2
其中 μ = μ' - jμ'',μ'' 表征磁损耗

磁滞损耗通常通过材料的 B-H 曲线面积或损耗角正切 tanδ_m 引入。在 COMSOL 中,磁损耗对应 mf.Qrh(磁场接口的电阻损耗)与 mf.Qml(磁损耗)。对于变压器铁芯、电感磁芯,磁损耗往往不可忽略。

3. 介电损耗

由介质极化弛豫引起,损耗密度为:

Q_dl = 1/2 * ω * ε'' * |E|^2     [W/m^3]

其中 ε = ε' - jε'',ε'' 为介电损耗因子
也可写为 Q_dl = 1/2 * ω * ε' * tanδ_e * |E|^2

在 COMSOL 中,介电损耗对应 ec.Qdl 或 emw.Qh(电磁波接口)。微波加热、射频消融、电容器介质发热等场景以介电损耗为主。

三类损耗并非总是同时显著。工频线圈以电阻损耗为主,高频变压器需同时考虑电阻损耗与磁损耗,微波腔体则以介电损耗为主。建模时应根据工作频率与材料属性判断主导项,避免无谓地开启全部损耗计算而增加求解负担。

二、多物理场耦合节点的建立

COMSOL 中电磁热耦合通过“多物理场”节点下的 电磁热 接口实现。其核心作用是把电磁接口计算出的损耗密度,自动映射为传热接口的热源项。建立过程如下:

  1. 在“多物理场”下添加 电磁热 节点,选择电磁场接口(如 ec、mf、emw)与传热接口(ht)。
  2. 在电磁热节点设置中,指定“电磁损耗”类型:电阻损耗、磁损耗、介电损耗,或全部。
  3. 在传热接口中,热源项会自动引用电磁损耗变量,例如 ec.Qrh + mf.Qml + ec.Qdl。
  4. 若需要温度反馈影响材料属性,还需在材料属性中设置温度依赖,并启用双向耦合。

耦合节点建立后,电磁场求解得到的损耗密度会作为体积热源进入传热方程:

ρ Cp ∂T/∂t + ρ Cp u · ∇T + ∇ · (-k∇T) = Q_em

其中 Q_em = Q_rh + Q_ml + Q_dl

若为稳态问题,左侧时间项为零,方程退化为 ∇ · (-k∇T) = Q_em。此时温度分布完全由损耗密度与边界散热条件决定。

三、双向耦合策略

电磁热耦合存在单向与双向两种策略,选择依据是温度是否显著改变电磁参数。

单向耦合:先求解电磁场,得到损耗密度,再求解传热场。适用于温度变化对电导率、磁导率、介电常数影响较小的场景,如低温升的母线、小功率电感。计算量小,收敛稳定。

双向耦合:电磁场与传热场交替迭代,温度更新后重新计算材料属性,再回代电磁场。适用于高温升场景,如大功率变压器、感应加热、微波消融。此时电导率随温度上升而下降(金属),或随温度上升而上升(半导体),磁导率与介电常数也会随温度漂移。

在 COMSOL 中实现双向耦合,通常采用 全耦合 或 分离式 求解器。全耦合将电磁与传热自由度放在同一方程组中求解,收敛快但内存占用高;分离式则交替求解两个物理场,适合强非线性问题。建议先以单向耦合获得初始温度场,再切换为双向耦合,可显著提高收敛性。

双向耦合中,材料属性若设置为温度的函数,必须保证函数连续可导。若使用实验数据插值,建议先做平滑处理,否则容易导致迭代振荡或求解器不收敛。

四、耦合关系可视化

下图展示了电磁热耦合中三类损耗的输入、耦合节点映射与温度反馈路径,帮助理解单向与双向耦合的差异。

电磁热耦合:损耗 → 热源 → 温度反馈 电磁场接口 ec / mf / emw 求解 E、H、J、B 输出损耗密度 Q_rh、Q_ml、Q_dl 三类损耗密度 电阻损耗 Q_rh = ½σ|E|² 磁损耗 Q_ml = ½ωμ''|H|² 介电损耗 Q_dl = ½ωε''|E|² 单位:W/m³ 多物理场耦合节点 电磁热 映射损耗 → 热源 Q_em = Q_rh+Q_ml+Q_dl 自动变量引用 传热接口 ht ρCp∂T/∂t + ∇·(-k∇T) = Q_em 求解温度场 T(x,y,z,t) 边界:对流、辐射、固定温度 热源 Q_em 温度反馈 σ(T)、μ(T)、ε(T) 单向耦合: 电磁场 → 损耗 → 传热场(温度不反馈) 双向耦合: 电磁场 ⇄ 损耗 ⇄ 传热场(温度反馈更新材料属性)

图中左侧为电磁场接口,中间为三类损耗密度,右侧为多物理场耦合节点,下方为传热接口。绿色实线表示损耗作为热源输入传热场,红色虚线表示温度反馈回电磁场更新材料属性。单向耦合只保留绿色路径,双向耦合则同时激活红色反馈路径。

五、实操要点与判断依据

建立耦合后,建议按以下顺序检查:

  • 确认电磁接口输出的损耗变量名与耦合节点中引用的变量一致,避免热源为零。
  • 检查材料属性是否已设置为温度的函数,否则双向耦合退化为单向耦合。
  • 观察求解日志中的残差曲线,若温度场残差下降缓慢,可尝试降低松弛因子或改用分离式求解器。
  • 后处理时,分别绘制 ec.Qrh、mf.Qml、ec.Qdl 的体图,判断哪类损耗占主导,据此决定是否需要简化模型。

以高频变压器为例,若铁芯磁损耗密度约为 5×10⁴ W/m³,绕组电阻损耗密度约为 2×10⁵ W/m³,则绕组温升会明显高于铁芯,散热设计应优先考虑绕组。若温度升高后铜电导率下降约 20%,电阻损耗随之下降,双向耦合结果会比单向耦合低约 10%~15%。这一差异在高温升场景中不可忽略。

判断是否启用双向耦合的实用准则:若预估温升引起的电导率变化超过 5%,或磁导率/介电常数变化超过 3%,建议启用双向耦合;否则单向耦合即可满足工程精度。

掌握损耗表达与耦合节点建立后,下一步即可进入具体电磁接口的建模设置,把本节的热源映射关系落实到实际模型中。