求解波动电磁场问题时，您可能会希望模拟一个带有开放边界的域，也就是说，计算域的边界支持电磁波以无反射的方式通过。针对这一问题，COMSOL 提供了几种解决方案。今天，我们将分析如何使用散射边界条件和完美匹配层来截断域，并讨论它们各自的适用范围。

为什么要截断域？

通常我们会对自由空间中的辐射对象产生模拟兴趣，比如天线。我们可以开发一个模型来模拟太空深处卫星上的天线，或者更常见的安装在吸波测试暗室中的天线。

无限自由空间中的天线，我们仅希望模拟环绕天线的一小块区域。

这类模型可通过 RF 模块或波动光学模块的‘电磁波，频域’接口开发。这些模块提供了相似的接口，可用于通过有限元方法求解频域形式的麦克斯韦方程组。（有关这些模块的主要差别，请阅读下方的文章。）

本文我们将只考察二维问题，其中电磁波沿 x-y 平面传播，电场沿 z 方向极化。我们还将假定模拟域为完全真空，因此频域麦克斯韦方程组可简化为：

其中， 

是电场、真空中的相对磁导率和介电常数为 

， 

是波数。

利用有限元方法求解上述方程需要一个有限大小的模拟域和一组边界条件。外部对所有辐射透明，因此我们希望沿外部应用一些边界条件，从而把此域截断为对自由空间的合理近似。我们还希望将该域截断得尽量小，因为对模型大小的控制将能帮我们降低计算成本。

现在来看一下 COMSOL Multiphysics 仿真环境中用于截断模拟域的两个选项：散射边界条件和完美匹配层。

散射边界条件

其中 r 是辐射轴。

当模拟域的边界位于源的无限远处时，这一边界条件就可以提供无反射传输，只是我们无法模拟一个无限大的域。因此，虽然我们无法精确应用 Sommerfeld 辐射条件，但可以应用它的一个合理近似边界。

现在让我们看看这一边界条件：

您可以清楚看到这一条件与 Sommerfeld 条件极其相似。本边界更正式的叫法应该是一阶散射边界条件 (SBC)，可以轻松在 COMSOL Multiphysics 中执行。实际上，它就是一个带复值系数的 Robin 边界条件。

不过，该条件有一个非常大的局限：它只有在辐射精确沿法向入射到边界上时才会无反射。所有非法向入射到 SBC 上的波都会发生部分反射。下图绘制了不同入射角平面波在一阶 SBC 上的反射系数。

一阶 SBC 上不同入射角平面波的反射。

从上图中我们可以观察到，随着入射波接近切向入射，波几乎会被完全反射掉。在 60° 的入射角下，反射约为 10%，显然，我们希望能有更好的边界条件。

COMSOL Multiphysics 还提供了二阶 SBC：

方程增加了二阶项，即电场沿边界的二阶切向导数，这在 COMSOL 软件架构中也非常容易执行。

我们来对比一阶和二阶 SBC 的反射系数：

不同入射角平面波在一阶和二阶 SBC 上的反射情况。

我们可以看到二阶 SBC 的一致性更好。我们现在可以让入射角约为 75° 时才使反射达到 10%。情况有所改善，但我们还可以做得更好。现在让我们将注意力从边界条件转到完美匹配层。

完美匹配层

回想一下，我们正尝试模拟类似于电磁波吸波暗室中的天线这一情景，这是一个壁面包含锥体楔形辐射吸收材料的房间，能够最小化所有的反射信号。我们可以对完美匹配层 (PML) 进行这样一个物理类比：相比边界条件，它其实更像我们沿模型外部增加的一块域，能够吸收所有的出射波。

从数学上讲，PML 只是一个带有各向异性和复值介电常数及磁导率的域。虽然从理论上讲 PML 无反射，但由于数值离散，即网格的关系，它们还是会表现出一些反射。为了最小化这一反射，我们希望在 PML 中使用一个与材料属性中的各向异性对齐的网格。下方显示了适用于二维圆形和三维球形域的 PML 网格。产品文档中还讨论了迪卡尔坐标、球形 PML，以及它们的正确使用。

适合二维和三维球形 PML 的网格。

现在让我们看看与 SBC 相比，PML 在不同入射角下的反射情况：

不同入射角下一阶和二阶 SBC 及 PML 的反射情况。

我们可以看到 PML 能在最广泛的范围内反射最少波。当波几乎平行于边界传播时，仍存在反射，不过幸运的是，实际中类似的场景并不多见。PML 的另一项特征是它不仅能吸收辐射波，还能够吸收倏逝波，不过，本文不对此进行过多介绍。因此，从物理视角来看，PML 的确可以作为一种拥有完美吸收特征的材料。

那么，您该选择哪种？

显然，根据文中的介绍，PML 是最佳方案。但与 SBC 相比，PML 会的内存使用更多。

因此，如果您还处于早期模拟阶段，希望开发一个计算强度略低的模型，二阶 SBC 是不错的选择。当您有理由相信任何发生在 SBC 上的反射并不会严重影响您关心的结果时，您就可以使用 SBC。

出于和软件之前版本兼容性的考虑，目前的缺省选项是一阶 SBC，但在 COMSOL Multiphysics 4.4 及以上版本中，则使用了二阶 SBC。我们这里只介绍了平面波形式的 SBC，其实还存在柱面波和球面波（3D）形式的一阶和二阶 SBC。虽然它们使用的内存较低，但与 PML 相比，都表现出更多的反射。

如果您事先对边界处的物理场了解不多，那 SBC 和 PML 将是很适合的条件。另一方面，如果您希望模拟的边界上的物理场具有某种模式时，比如表征了波导的边界，‘端口’和‘集总端口’会更加适合。在接下来的博客中，我们将会讨论这些条件。

声明：本文来源于互联网，麦涛网本着传播知识、相互交流与学习的目的进行转载，并为网友免费提供。文章已尽力标明转载出处，如有著作权人或出版方提出异议，或有疑问，请随时联系（tel：0731-89579179，email：office@maitaonet.com）并告知我们，以便及时删除和纠正，谢谢。