2. 电磁仿真基础

在CST电磁仿真中，网格划分并非孤立的前处理步骤，而是与电磁场求解的数学本质紧密耦合。理解麦克斯韦方程组如何被离散化，以及有限元法（FEM）、有限积分法（FIT）和矩量法（MOM）各自对网格的“偏好”，是正确设置网格参数、判断收敛性的前提。本章从控制方程出发，剖析三种主流数值方法在CST中的实现差异，并揭示网格密度、网格类型与求解精度、资源消耗之间的内在联系。

2.1 麦克斯韦方程组：连续域到离散域的桥梁

电磁仿真的本质是在计算机上求解麦克斯韦方程组。该方程组在连续空间中具有无穷多个自由度，而计算机只能处理有限维度的代数方程。因此，所有数值方法的核心思想都是将连续的偏微分方程或积分方程转化为离散的矩阵方程。这一转化过程直接决定了网格的形态与分布。

在CST中，不同求解器对麦克斯韦方程组的处理方式截然不同：FEM（频域求解器）基于微分形式的波动方程，将计算域划分为有限个四面体单元，在每个单元内用基函数近似场分布；FIT（时域求解器）则直接对积分形式的麦克斯韦方程组进行空间离散，在交错网格上定义电场和磁场分量；MOM（积分方程求解器）将边界条件转化为积分方程，仅对导体表面或介质分界面进行离散。这三种数学路径决定了网格的维度（体网格或面网格）、单元形状（四面体或六面体）以及加密策略。

数值方法	控制方程形式	离散对象	典型网格	CST求解器
FEM	微分形式（波动方程）	整个计算域（体）	四面体（自适应）	频域求解器
FIT	积分形式（Maxwell网格方程）	整个计算域（体）	六面体（正交网格）	时域求解器
MOM	积分方程（边界积分）	导体表面/介质分界面	三角形面网格	积分方程求解器

2.2 FEM与四面体网格：精度与自由度的平衡

有限元法（FEM）在CST频域求解器中占据主导地位。其核心思想是将计算域划分为有限个互不重叠的四面体单元，在每个单元内使用矢量基函数（如Whitney基函数）来近似电场分布。四面体网格的最大优势在于其几何适应性——能够精确贴合任意复杂的曲面结构，如弯曲的微带线、共形天线罩或螺旋电感。

然而，这种灵活性也带来了代价：四面体网格的单元数量通常远大于同等精度下的六面体网格，导致矩阵规模庞大。CST通过自适应网格加密策略来缓解这一问题——初始网格较粗，求解后根据场梯度分布自动在电场变化剧烈的区域（如导体边缘、缝隙处）加密网格，迭代数次直至S参数收敛。

2.3 FIT与六面体网格：时域效率的基石

有限积分法（FIT）是CST时域求解器的数学基础。它将麦克斯韦方程组离散到一对相互正交的网格系统上（主网格与对偶网格），电场定义在主网格的棱边上，磁场定义在对偶网格的棱边上。这种交错布局天然满足旋度方程的离散形式，保证了电荷守恒和能量守恒的数值稳定性。

FIT使用六面体（正交）网格，其优势在于：① 网格生成简单，内存占用低；② 采用显式时间步进（如蛙跳格式），无需求解大型矩阵方程；③ 一次时域仿真即可获得宽频响应。但六面体网格对曲面结构存在“阶梯近似”误差，需要局部网格加密来缓解。CST提供局部网格细化和共形网格技术，在保持六面体结构的同时提高曲面拟合精度。

// 时域求解器网格设置示例（CST宏语言）
With Mesh
    .Type = "Hexahedral"
    .StepsPerWavelength = 15
    .EdgeLengthRatio = 1.5
    .UseLocalRefinement = True
    .RefineAtEdges = True
    .RefineAtSurfaces = True
End With

2.4 MOM与面网格：开放问题的自然选择

矩量法（MOM）在CST中主要用于电大尺寸的辐射与散射问题，如天线阵列、雷达散射截面（RCS）分析。与FEM和FIT不同，MOM仅需对导体表面或介质分界面进行离散，将三维体积问题降维为二维面问题，极大减少了未知量数量。其网格单元为三角形面片，能够精确描述任意曲面。

MOM的网格密度要求与工作频率直接相关：通常每个波长需要8-10个三角形单元。对于线状结构（如偶极子天线），可使用线网格进一步降低自由度。但MOM生成的矩阵是稠密的，求解复杂度为O(N²)至O(N³)，因此对内存需求较高。CST采用多层快速多极子算法（MLFMA）加速矩阵-向量乘积，使电大尺寸问题的求解成为可能。

2.5 网格与收敛判据的联动机制

网格划分不是一次性的操作，而是与收敛判据形成闭环迭代。CST的典型流程为：初始网格 → 求解 → 评估误差 → 局部加密 → 重新求解 → 检查收敛。收敛判据通常包含两类指标：目标量稳定性（如S参数随网格加密的变化量）和场分布一致性（如关键位置的电场幅值变化）。

在实际工程中，建议采用“两步验证法”：第一步，使用默认网格求解，观察S参数曲线趋势；第二步，将全局网格密度提高50%，对比两次结果。若S参数变化小于0.05dB且谐振频率偏移小于0.1%，则认为网格收敛。对于宽带结构，还需检查多个频点处的场分布是否一致。

收敛指标	阈值建议	检查方法	不满足时的措施
S参数幅度变化	< 0.01 dB	对比两次加密结果	加密场强区域网格
谐振频率偏移	< 0.1%	提取S11最小值频率	细化敏感结构
远场方向图差异	主瓣 < 0.1 dB	对比方向图曲线	增加辐射边界网格
能量残差	< 1e-6	查看求解器日志	增加迭代次数

2.6 案例拆解：微带贴片天线的网格收敛性分析

以一个2.45 GHz矩形微带贴片天线为例，基板为Rogers RT5880（εr=2.2，厚度1.575 mm），贴片尺寸38mm×28.5mm。使用CST频域求解器（FEM）进行仿真，目标为S11 < -10dB的阻抗带宽。

初始设置：全局网格密度设为每波长12个四面体单元，开启自适应加密，最大迭代次数为10次。第一次求解后，S11在2.45GHz处为-18.3dB，但谐振频率偏移至2.52GHz（偏移2.9%），明显不收敛。

分析判断：谐振频率偏高说明贴片电尺寸偏小，可能是基板厚度方向网格不足导致有效介电常数计算偏低。检查网格剖面发现，基板厚度方向仅有1层四面体，无法准确模拟垂直场分布。

操作步骤：① 在基板厚度方向设置最小网格数为3层；② 对贴片边缘（电流集中区域）设置局部加密，网格尺寸从2mm细化至0.5mm；③ 重新求解并对比S11曲线。经过3次自适应迭代后，谐振频率稳定在2.45GHz±0.01GHz，S11幅度变化小于0.005dB。

结果对比：最终网格数量从初始的12万增加至35万，求解时间从4分钟延长至15分钟，但仿真结果与实测数据高度吻合（谐振频率偏差0.3%）。复盘经验：对于含薄介质层的结构，必须保证厚度方向至少3层网格；自适应加密应结合人工判断，不能完全依赖自动流程。

2.7 常见误区与总结

误区一：网格越密越精确。实际上，当网格密度超过每波长20个单元后，精度提升趋于饱和，而内存和时间的开销呈指数增长。过度加密还可能引入数值色散误差。正确做法是采用局部加密，仅在场梯度大的区域细化。

误区二：所有求解器使用相同的网格策略。FEM适合复杂曲面和介质体，FIT适合宽频和电大尺寸，MOM适合开放辐射问题。网格类型必须与求解器匹配，例如在时域求解器中使用四面体网格会导致时间步长过小而无法实际计算。

误区三：忽略网格质量指标。除了密度，网格的畸变度（Skewness）、长宽比（Aspect Ratio）和雅可比行列式也直接影响求解精度。CST的网格检查工具会标记质量较差的单元，应优先修复这些区域。

总结：网格划分是电磁仿真中连接物理模型与数值算法的关键纽带。理解麦克斯韦方程组如何被FEM、FIT、MOM离散化，能够帮助工程师从数学本质出发，合理选择网格类型、密度和加密策略，从而在有限的计算资源下获得高可信度的仿真结果。收敛判据不是简单的“跑完即止”，而是通过网格-求解-验证的闭环迭代，确保仿真结果逼近真实物理行为。