PML — CAE术语解说

好问题。PML是"Perfectly Matched Layer"的缩写，是在计算领域的边界处放置的人工层，用来消除"非物理反射"和"计算领域内物理波动"之间的不匹配，即将反射降至接近零。例如，从天线辐射的电磁波在分析区域的边缘突然遇到墙壁时，会产生现实中不存在的反射波，污染结果。PML的作用就像海绵一样吸收这种反射。

完全正确，从根本上就不同。普通的衰减材料会因为阻抗不匹配而产生反射。PML的核心是对波动方程进行"复数坐标变换"。通过将实数坐标扩展到复数，可以得到波在PML中传播时呈指数衰减的解。这种变换产生了"匹配"。例如，在一维波动方程中，坐标x可以变换为

在实现层面，复数坐标变换的效果被建模为具有"复介电常数和复磁导率"或"复折射率"的媒质。也就是说，PML层被定义为具有"损耗"和"色散"特性的虚拟材料。在Ansys HFSS中，作为"PML边界条件"，用户设置层的厚度（例如中心波长的10%至20%）和衰减廓线（多项式或几何级数）。后台会自动计算这个等效媒质的参数。

必须的。PML作为"层"添加到计算域外侧，那里当然也需要网格。在FEM中，PML区域的单元被赋予前面所说的复杂材料特性。重要的是保持PML内网格质量。特别是当使用急剧变化的衰减廓线时，粗网格无法正确表现衰减，反射低减性能会下降。根据经验，PML层的厚度方向至少需要3～5个单元。

是有好的起点的。厚度根据所处理的波长λ确定。电磁波分析中通常为λ/4到λ。衰减系数σ的最大值设置成使PML出口处波动振幅足够小。经验法则是，PML内传播一个波长时的衰减量应≥60dB（振幅为1/1000）。COMSOL Multiphysics的"吸收边界条件"中，用户可以调整多项式次数和缩放系数，但多数情况下默认设置（例如多项式次数3、理论反射率-60dB）就足够了。

敏锐的指出。完全正确，PML区域的单元具有复数材料常数，所以整体系数矩阵也是复数。与实数矩阵相比，存储容量约增加2倍，计算时间也会增加。但是，替代方案"吸收边界条件(ABC)"的反射低减性能较差（特别是斜入射和低频）。这是个折衷。在大规模问题中，只在需要的方向应用PML（例如辐射问题只在辐射方向）可以降低成本。

首先是确认PML是否真的必要。对于封闭的腔体或波导，边界是完全导体且物理上有反射，就不需要PML。PML必要的情况是"开放边界"模拟，即天线辐射、散射截面、声学辐射等问题。其次是放置位置。至少应离光源或散射体λ/2以上。过近的话，PML无法正确处理近场的衰落波，结果会不准确。

有两个现实的方法。第一是"与参考解比较"。用非常大的计算域（反射可以忽略）代替PML求解，把它作为正解，与用PML的小域解比较。第二是"监测残余反射"。在Ansys HFSS中可以用"Field Override"定义PML前的面，计算该面上后向的电功率流（坡印廷矢量），看相对于入射功率的比值。实用上目标是-50dB以下。

主要有三个。1)共振频率偏移：虽然不是闭合系统，但特定频率出现尖锐共振峰。这是边界反射产生驻波的证据。2)辐射方向图失真：天线指向性相比理论或测量变得左右不对称。3)收敛不好：频率扫描中S参数呈锯齿状，或者即使网格加细结果仍不稳定。这些都表明PML的反射低减性能不足，波动"被困"在计算域内。

概念相同，但界面和应用范围不同。Ansys HFSS中PML作为3D对象定义（"Create"→"Boundary"→"PML"）。COMSOL中则作为"吸收边界条件"的一种选择"PML"，并指定应用的"域"。CST Studio Suite中在"边界条件"设置中选择"PML"，可以细致地控制层数和廓线。另外，CST基于FIT（有限积分法），所以PML的实现在内部上与基于FEM的软件不同，需要注意。

也用。是在声场（压力音场）和弹性波（固体中振动）的辐射和散射问题中。COMSOL的"声学模块"和"结构力学模块"有专用的"PML"功能。但因为控制方程不同，定义的等效复杂材料特性与电磁的情况不同。例如声学中定义复数密度和复数体积弹性率。Ansys Mechanical中也用"无限单元"达到类似目的，但数学基础有时不同于PML。

相比商业软件难度高。需要用户自己定义PML域，并在该域内实现基于适当复杂坐标变换的弱形式。例如用FEniCS进行电磁波模拟时，需要在弱形式中明确加入PML复数坐标变换的雅可比值。OpenFOAM的声学求解器"acousticFoam"虽有PML功能，但需要正确编写设置文件（`fvOptions`）。对教育目的或需要定制化的研究很强大，但实务应用的进入门槛较高。

首先怀疑的是"数值不稳定性"。如果把PML的衰减系数σ_max设置得太大，PML层内波动急剧衰减。这种急峻变化网格跟不上，离散化误差增大，求解器可能发散。对策是三点：降低σ_max（从默认值开始），增加PML厚度（用更平缓的梯度达到同样衰减），PML层内网格加细。检查Ansys日志中是否有"matrix is ill-conditioned"这样的警告。

完全正确。PML按照波长设计厚度，所以低频时相对PML变薄，吸收不足。具体对策两个。第一是"PML厚度不用波长比，用绝对长度确保"。例如，最低频率1MHz（空气中波长300m）也要确保0.5m厚度。第二是"用频率相关的衰减系数廓线"。COMSOL等中可设置衰减系数σ与频率成正比（

需要谨慎。经典PML理论基于各向同性线性媒质。面向各向异性材料的PML理论上无法完全匹配，会有多少反射产生。商业软件（Ansys、COMSOL）内部试图应对，但性能可能不如各向同性情况。在非线性媒质是主角的问题（如高强度光的非线性光学效应）中，PML应用非常困难，通常不用。此时要么物理上扩大计算域，要么考虑其他开放边界方法。使用前务必查阅软件手册中的支持情况。