2. 电磁场基础与COMSOL电磁接口

在COMSOL中做电磁热耦合，第一步不是急着点“电磁波，频域”或“磁场”，而是先判断你面对的物理问题属于哪一类：是辐射波在波长尺度上传播，还是低频电流在导体中扩散，抑或是两者在同一模型中并存。这个判断直接决定了你该选RF模块、AC/DC模块还是波动光学模块，也决定了后续网格、边界条件和耦合方式的合理性。本节从麦克斯韦方程组出发，梳理本构关系与边界条件，再落到COMSOL三个电磁接口的适用场景与选择依据。

麦克斯韦方程组：四个方程，四种物理图像

微分形式的麦克斯韦方程组是电磁仿真的出发点：

∇×E = -∂B/∂t          （法拉第定律）
∇×H = J + ∂D/∂t        （安培-麦克斯韦定律）
∇·D = ρ                （高斯定律）
∇·B = 0                （磁通连续性）

这四个方程本身并不足以求解，还需要本构关系把场量联系起来。在线性、各向同性介质中：

D = εE = ε₀εᵣE
B = μH = μ₀μᵣH
J = σE

其中εᵣ、μᵣ、σ分别是相对介电常数、相对磁导率和电导率。这三个参数决定了材料对电磁场的响应方式，也决定了COMSOL中材料节点的填写内容。例如，铜的σ≈5.998×10⁷ S/m，在低频下表现为良导体；而石英玻璃的σ≈10⁻¹⁷ S/m，几乎可视为理想介质。当σ很大时，位移电流∂D/∂t相对传导电流J可以忽略，方程退化为磁准静态形式；当频率升高到微波或光学波段，位移电流重新变得重要，必须保留完整波动方程。

边界条件：把无限空间截断成可计算问题

数值仿真只能计算有限区域，边界条件就是把无限空间“截断”的数学手段。常见的电磁边界条件包括：

• 理想电导体（PEC）：n×E=0，电场切向分量为零，用于模拟金属壁或接地导体。
• 理想磁导体（PMC）：n×H=0，磁场切向分量为零，常用于对称面或周期边界。
• 阻抗边界条件：用表面阻抗Zs近似薄导体层，避免对趋肤深度网格的过度细化。
• 散射边界条件（SBC）：一阶吸收边界，允许波几乎无反射地离开计算域。
• 完美匹配层（PML）：在计算域外围加一层各向异性吸收介质，把出射波完全吸收，是波动问题中最常用的截断方式。
• 端口边界：用于馈入或提取电磁能量，可指定功率、电压或模式场。

边界条件选错，往往表现为结果“看起来有场”，但能量不守恒或反射异常。例如，在微波加热模型中，如果波导端口用PEC封闭而不用端口边界，入射波无法进入，整个模型就没有能量输入。

三个电磁接口的适用场景

COMSOL的电磁产品线按频率和尺度划分，核心接口如下：

选择依据可以归纳为三条：看尺度——当结构尺寸与波长可比时，必须用全波接口；当尺寸远小于波长时，准静态接口足够。看材料——高电导率、低频率优先考虑AC/DC；低损耗介质、高频率优先考虑RF或波动光学。看输出——需要S参数、辐射方向图、模式分析时用RF；需要磁场分布、电感、涡流损耗时用AC/DC。

从电磁到热的耦合路径

电磁热耦合的核心是电磁损耗作为热源。在频域电磁中，体积损耗密度为：

Q = ½ Re(J·E*) + ½ Re(iωB·H*)

第一项是电阻损耗（焦耳热），第二项是磁损耗。在COMSOL中，这个量通常以“电磁热源”多物理场节点自动传递到“传热”接口的“热源”项。对于微波加热，Q主要来自介电损耗；对于感应加热，Q主要来自涡流和磁滞损耗。理解这一点，才能在耦合时正确选择“电磁热”多物理场节点，而不是手动输入一个常数热源。

实操中的选择示例

假设要仿真一个微波炉腔体加热一块土豆（介电材料），频率2.45 GHz，腔体尺寸约30 cm，波长约12.2 cm。此时结构尺寸与波长可比，必须用“电磁波，频域”接口，边界用PEC模拟金属壁，用端口馈入微波功率。土豆的介电损耗会产生Q，通过“电磁热”多物理场耦合到“传热”接口，求解温度分布。如果换成感应加热一块钢坯，频率50 Hz，钢坯尺寸10 cm，波长约6000 km，远大于结构尺寸，此时用“磁场”接口计算涡流，再用“电磁热”耦合到传热。两个案例的接口选择截然不同，但耦合路径一致。

最后强调一点：COMSOL中的“电磁热”多物理场节点默认是单向耦合（电磁→热），即温度不影响电磁参数。如果材料电导率或介电常数随温度显著变化，需要手动开启双向耦合，在材料属性中把σ、εᵣ写成温度T的函数，并确保求解器能处理这种非线性。这一步在感应加热和微波烧结中尤为关键，也是电磁热耦合从“能算”走向“算得准”的分水岭。