2. HFSS自适应网格收敛原理
在HFSS中,有限元求解的精度并不取决于初始网格有多密,而取决于网格是否“足够贴合”当前模型的电磁场分布。自适应网格细分(Adaptive Mesh Refinement)正是HFSS用来逐步逼近这一目标的核心机制:它先以较粗的初始网格求解,再根据场解误差估计自动加密网格,反复迭代,直到结果收敛或达到迭代上限。理解这一机制,是判断仿真结果是否可信的第一步。
一、自适应迭代的完整闭环
HFSS的每一次自适应迭代(Adaptive Pass)都包含四个连续动作:求解、误差估计、网格加密、收敛判断。求解阶段在现有网格上计算S参数与场分布;误差估计阶段并不依赖用户判断,而是由求解器基于相邻单元之间的场解连续性来估计数值误差,误差大的区域意味着该处场变化剧烈而网格分辨率不足;网格加密阶段只针对误差较大的单元进行细分,通常将四面体单元一分为多,而不是全局均匀加密;收敛判断阶段则比较本次与上一次迭代的Delta S,决定是否继续。
这种“局部加密”策略是自适应方法效率远高于全局均匀加密的根本原因。对于微带线边缘、过孔、缝隙、介质交界面等场集中区域,网格会自动变密;而对于场变化平缓的均匀区域,网格保持稀疏,从而在精度与计算量之间取得平衡。
二、收敛判据Delta S的含义
Delta S是HFSS自适应收敛的核心判据,其定义为:本次迭代与上一次迭代之间,所有端口S参数幅度的最大变化量。用公式表达即:
Delta S = max | |S(pass N)| - |S(pass N-1)| |
注意这里取的是绝对值之差的最大值,遍历所有端口的S参数。例如某二端口网络,Pass 3得到S11=-12.35 dB、S21=-0.82 dB,Pass 4得到S11=-12.41 dB、S21=-0.83 dB,则两次的幅度差分别为0.06 dB和0.01 dB,Delta S取0.06 dB。若用户设定的收敛判据为0.02,则本次尚未收敛,需要继续迭代。
Delta S的物理意义是:当网格加密后S参数几乎不再变化,说明当前网格已经能够稳定地描述电磁响应,继续加密对结果影响可忽略。因此Delta S越小,结果越可信,但代价是迭代次数和计算时间增加。工程上常用的收敛判据范围是0.01~0.02,对于高Q谐振结构、窄带滤波器等对精度敏感的问题,可收紧到0.005甚至更小;对于宽带粗略评估,0.02~0.05也可接受。
三、最大迭代次数与收敛失败
HFSS默认最大迭代次数为6,很多用户会直接沿用。但默认值并不总是足够:当模型包含高对比度介质、细小缝隙、强谐振结构时,Delta S可能在6次迭代后仍未降到判据以下,此时求解器会停止并给出“未收敛”提示。这时得到的S参数并非完全错误,但精度无法保证,不能直接作为最终结论。
遇到未收敛时,正确的处理顺序是:先检查Delta S随迭代次数的下降趋势,若仍在稳定下降,可适当增加最大迭代次数(如10~15);若Delta S在某次迭代后反而上升或剧烈波动,则往往说明网格加密策略与模型特征不匹配,需要检查是否存在过细的几何特征、材料参数突变或端口设置问题,而不是盲目增加迭代次数。
四、收敛过程的关键观察量
在HFSS的求解数据中,与自适应收敛直接相关的观察量包括:每次迭代的Delta S、网格单元总数、求解矩阵规模、单次求解时间以及累计求解时间。随着迭代推进,网格单元数通常逐次增加,矩阵规模随之增大,单次求解时间也会上升。如果发现某次迭代网格数暴增而Delta S改善很小,说明加密效率低,可能需要调整初始网格或收敛设置。
下表给出一个典型微带贴片天线在收敛判据0.02下的自适应迭代示例(数据为教学示例,非实测):
| 迭代次数 | Delta S | 网格单元数 | 是否收敛 |
|---|---|---|---|
| 1 | 0.185 | 8,420 | 否 |
| 2 | 0.072 | 14,760 | 否 |
| 3 | 0.031 | 23,180 | 否 |
| 4 | 0.014 | 35,640 | 是 |
可以看到,Delta S从0.185逐步降到0.014,网格单元数从8千增长到3.5万,第4次迭代满足0.02判据后停止。这一过程说明:自适应收敛不是“一次算准”,而是通过多轮“求解—估计—加密”逐步逼近稳定解。理解这一点,才能在结果异常时判断问题出在网格、设置还是模型本身。