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 中电磁热耦合通过“多物理场”节点下的 电磁热 接口实现。其核心作用是把电磁接口计算出的损耗密度,自动映射为传热接口的热源项。建立过程如下:
- 在“多物理场”下添加 电磁热 节点,选择电磁场接口(如
ec、mf、emw)与传热接口(ht)。 - 在电磁热节点设置中,指定“电磁损耗”类型:电阻损耗、磁损耗、介电损耗,或全部。
- 在传热接口中,热源项会自动引用电磁损耗变量,例如
ec.Qrh + mf.Qml + ec.Qdl。 - 若需要温度反馈影响材料属性,还需在材料属性中设置温度依赖,并启用双向耦合。
耦合节点建立后,电磁场求解得到的损耗密度会作为体积热源进入传热方程:
ρ Cp ∂T/∂t + ρ Cp u · ∇T + ∇ · (-k∇T) = Q_em
其中 Q_em = Q_rh + Q_ml + Q_dl
若为稳态问题,左侧时间项为零,方程退化为 ∇ · (-k∇T) = Q_em。此时温度分布完全由损耗密度与边界散热条件决定。
三、双向耦合策略
电磁热耦合存在单向与双向两种策略,选择依据是温度是否显著改变电磁参数。
单向耦合:先求解电磁场,得到损耗密度,再求解传热场。适用于温度变化对电导率、磁导率、介电常数影响较小的场景,如低温升的母线、小功率电感。计算量小,收敛稳定。
双向耦合:电磁场与传热场交替迭代,温度更新后重新计算材料属性,再回代电磁场。适用于高温升场景,如大功率变压器、感应加热、微波消融。此时电导率随温度上升而下降(金属),或随温度上升而上升(半导体),磁导率与介电常数也会随温度漂移。
在 COMSOL 中实现双向耦合,通常采用 全耦合 或 分离式 求解器。全耦合将电磁与传热自由度放在同一方程组中求解,收敛快但内存占用高;分离式则交替求解两个物理场,适合强非线性问题。建议先以单向耦合获得初始温度场,再切换为双向耦合,可显著提高收敛性。
四、耦合关系可视化
下图展示了电磁热耦合中三类损耗的输入、耦合节点映射与温度反馈路径,帮助理解单向与双向耦合的差异。
图中左侧为电磁场接口,中间为三类损耗密度,右侧为多物理场耦合节点,下方为传热接口。绿色实线表示损耗作为热源输入传热场,红色虚线表示温度反馈回电磁场更新材料属性。单向耦合只保留绿色路径,双向耦合则同时激活红色反馈路径。
五、实操要点与判断依据
建立耦合后,建议按以下顺序检查:
- 确认电磁接口输出的损耗变量名与耦合节点中引用的变量一致,避免热源为零。
- 检查材料属性是否已设置为温度的函数,否则双向耦合退化为单向耦合。
- 观察求解日志中的残差曲线,若温度场残差下降缓慢,可尝试降低松弛因子或改用分离式求解器。
- 后处理时,分别绘制
ec.Qrh、mf.Qml、ec.Qdl的体图,判断哪类损耗占主导,据此决定是否需要简化模型。
以高频变压器为例,若铁芯磁损耗密度约为 5×10⁴ W/m³,绕组电阻损耗密度约为 2×10⁵ W/m³,则绕组温升会明显高于铁芯,散热设计应优先考虑绕组。若温度升高后铜电导率下降约 20%,电阻损耗随之下降,双向耦合结果会比单向耦合低约 10%~15%。这一差异在高温升场景中不可忽略。
掌握损耗表达与耦合节点建立后,下一步即可进入具体电磁接口的建模设置,把本节的热源映射关系落实到实际模型中。