4. 网格划分策略

在 CST 的时域求解器(Time Domain Solver)中,网格并不是“画完模型后随手点一下”的附属步骤,而是直接决定仿真精度、内存占用与计算时间的核心参数。CST 时域求解器基于有限积分技术(Finite Integration Technique, FIT),其网格是六面体(hexahedral)的,与频域有限元常用的四面体网格有本质区别。理解 FIT 网格的构造逻辑,是制定合理划分策略的前提。

一、时域有限积分网格的基本构造

FIT 把麦克斯韦方程组在空间上离散到一对正交的网格上:电场分量定义在网格棱边(edge)上,磁场分量定义在网格面元(face)上,二者通过网格的拓扑关系相互耦合。因此 CST 的网格是“结构化六面体网格”,每个网格单元都是长方体,网格线沿 x、y、z 三个方向分布。这种结构带来两个直接后果:一是网格生成快、内存效率高;二是对曲面和斜面的贴合能力弱于四面体网格,需要用阶梯近似(staircase)或薄片技术(PBA,Perfect Boundary Approximation)来逼近曲面。

网格密度由“每波长网格数”(cells per wavelength)控制。CST 默认在最高频率处保证约 10~15 个网格单元每波长,这是保证数值色散误差可接受的经验下限。对于天线问题,最高频率通常取工作频段上限,若关注带外特性或宽带匹配,则需按实际关心的最高频率设定。网格太粗会导致谐振频率偏移、方向图畸变;网格太细则内存和时间急剧上升,且收益递减。

FIT 网格的“结构化”意味着网格线是全局贯通的,局部加密会沿整条网格线传播。因此局部加密要谨慎,避免在无关区域引入大量冗余网格。

二、自适应网格加密的工作机制

CST 提供自适应网格加密(Adaptive Mesh Refinement, AMR),其基本流程是:先用较粗的初始网格求解一次,然后根据场分布或 S 参数的变化估计各区域的误差,在误差大的区域自动加密网格,再求解、再评估,直到满足收敛判据或达到最大迭代次数。自适应加密通常以 S 参数的收敛为主要判据,例如要求相邻两次迭代的 S11 幅度差小于设定阈值(如 0.01~0.02),或谐振频率偏移小于设定值。

自适应加密的优点是减少人工试错,尤其适合结构复杂、场分布事先难以预判的天线。但它也有代价:每次迭代都要重新求解,总时间可能是单次求解的数倍;而且自适应判据基于全局 S 参数,可能对某些局部场(如近场、SAR、表面电流)不够敏感。因此对于需要精确近场或电流分布的问题,不能完全依赖自适应,仍需结合人工判断。

三、网格收敛性检查与优化

网格收敛性检查的核心思想是:逐步加密网格,观察关键结果是否趋于稳定。若加密后结果变化很小,说明网格已足够;若变化仍明显,说明网格不足。实际操作中,可以固定其他设置,只改变网格密度(如从 10 cells/λ 逐步提高到 15、20 cells/λ),记录谐振频率、S11 最小值、增益、方向图等指标,绘制“结果—网格密度”曲线,找到曲线趋于平坦的拐点。

优化网格的策略包括:

  • 局部加密优先于全局加密:在馈电点、缝隙、细导线、介质边缘等场变化剧烈处设置局部网格加密(Local Mesh Refinement),而不是整体提高密度。
  • 利用对称性:若天线结构具有对称面,可设置对称边界(Electric/Magnetic Symmetry),把计算域减半甚至减到四分之一,等效于提高网格利用率。
  • 合理设置频率范围:网格密度按最高关心频率设定,不要盲目按远高于工作频段的频率划分。
  • 检查网格统计:求解后查看 Mesh Statistics,关注网格总数、最小/最大网格尺寸、长宽比(aspect ratio)。长宽比过大(如超过 10:1)会降低精度,应尽量控制在合理范围。
  • 使用 PBA 与薄片技术:对曲面和薄结构启用 PBA,可在不显著增加网格的情况下改善几何逼近精度。

下面这张图概括了从初始网格到收敛判断的完整决策关系,以及各环节的关键参数与观察对象。

FIT 网格划分与收敛性检查决策关系 初始网格 10~15 cells/λ @ fmax 时域求解 记录 S11 / 谐振频率 / 增益 收敛判断 ΔS11 < 0.01 或 Δf < 0.5% 未收敛 自适应加密 / 局部加密 馈电点、缝隙、介质边缘、细导线 检查长宽比 < 10:1,启用 PBA 重新求解 收敛 输出结果 S 参数、方向图、增益、效率 网格优化要点 • 局部加密优先于全局加密 • 利用对称面减半计算域 • 按最高关心频率定密度 • 查看 Mesh Statistics 统计 • 控制长宽比,启用 PBA • 记录多组网格密度对比 收敛趋势示意 网格密度 (cells/λ) f₀ 拐点:网格已足够

四、实操示例:2.4 GHz 微带贴片天线的网格设置

以一个 2.4 GHz 矩形微带贴片天线为例,介质基板为 FR4(εr≈4.3,厚度 1.6 mm),贴片尺寸约 38 mm × 29 mm。设定频率范围为 2.0~2.8 GHz,最高频率 2.8 GHz,对应自由空间波长约 107 mm。若按 15 cells/λ 设定,则最大网格尺寸约 7 mm;但介质中波长缩短为约 51.6 mm,介质内网格应更密,约 3.4 mm。CST 的 Mesh Expert 会自动根据材料分配网格,但建议手动检查介质层厚度方向是否有足够网格:1.6 mm 厚度若只有 1 层网格,会严重低估介质效应,通常应保证厚度方向至少 3~5 层,即网格尺寸约 0.3~0.5 mm。

馈电点附近电流集中,可设置局部加密,把馈电探针周围 2~3 mm 区域的网格尺寸降到 0.2 mm。贴片边缘和缝隙处也建议局部加密。完成设置后,先跑一次粗网格,记录 S11 最小值频率 f₀ 和增益 G;再把全局网格密度提高到 20 cells/λ,重跑,比较 f₀ 和 G 的变化。若 f₀ 偏移小于 0.5%、G 变化小于 0.2 dB,可认为收敛。若偏移仍大,继续加密或检查是否存在网格长宽比过大、曲面逼近不足等问题。

自适应加密虽然方便,但不要把它当作“一键解决”。若初始网格过粗,自适应可能收敛到错误的局部极值;若模型存在细小结构(如 0.1 mm 的缝隙),初始网格必须能分辨它,否则自适应也无从加密。

五、常见问题与判断依据

实际使用中,网格问题常表现为:谐振频率偏高(网格过粗,等效介电常数偏低)、S11 曲线毛刺多(网格不均匀或长宽比过大)、方向图不对称(对称面网格不一致)、计算时间异常长(局部加密过度或全局密度过高)。遇到这些现象,应先查看 Mesh Statistics 中的网格总数和尺寸分布,再决定是加密、减密还是调整局部区域。网格收敛性检查不是一次性的,而应作为仿真流程中的固定环节,尤其在结果用于设计决策时,必须给出网格收敛的证据。

总结来说,FIT 网格是结构化的六面体网格,密度由每波长网格数控制;自适应加密能减少人工试错,但不能替代对关键区域的判断;收敛性检查通过逐步加密观察结果稳定性,是保证仿真可信度的必要步骤。掌握这三者的关系,才能在精度与效率之间找到平衡。