矢量有限元基函数的选择
在电磁场数值分析中，有限元法（FEM）是解决复杂几何结构下麦克斯韦方程组的核心工具。与标量场问题不同，电磁场通常由矢量势（如磁矢势 $\mathbf{A}$ 或电矢势 $\mathbf{F}$）描述。因此，选择合适的矢量有限元基函数对于保证解的收敛性、精度以及避免“虚假模”（Spurious Modes）至关重要。基函数的选择直接决定了离散系统的性质，是构建稳定数值求解器的基础。
标量基函数的局限性
在早期或简单的电磁模拟中，工程师常尝试使用标量拉格朗日（Lagrange）基函数来近似矢量场。然而，这种做法存在根本性的缺陷。标量基函数在单元内部是连续的，但在单元边界上，其法向分量可能不连续。根据电磁理论，磁感应强度 $\mathbf{B}$ 的法向分量在介质界面上必须连续，而磁矢势 $\mathbf{A}$ 的法向分量通常要求连续或满足特定的自然边界条件。
使用标量基函数会导致以下问题：
- 虚假模的产生：在低频或静态极限下，离散系统可能出现非物理的零频率解，即虚假模。这些模在物理上对应于零磁场，但在数值解中表现为非零的矢量场，严重干扰结果。
- 收敛性差：对于涉及旋度算子（$\nabla \times$）的问题，标量基函数无法正确逼近解的空间，导致误差收敛率低于理论预期。
Nédéc 边缘基函数
为了解决上述问题，Jean-Claude Nédéc 在 1980 年代提出了一组专门用于矢量场的基函数，现称为 Nédéc 基函数或边缘元（Edge Elements）。这类基函数是构造 $H(\text{curl})$ 有限元空间的标准选择，特别适用于求解涉及旋度的电磁问题，如时谐场和瞬态场分析。
Nédéc 基函数具有两个关键特性：
- 切向连续性：相邻单元共享边上的基函数值相等，确保了矢量场切向分量的连续性。这符合电磁场在边界上的物理约束。
- 法向不连续：基函数在单元法向方向上不强制连续，这允许解在界面处发生突变，符合物理实际。
一阶 Nédéc 基函数示例
在四面体单元中，一阶 Nédéc 基函数与单元的每条边相关联。每条边对应一个基函数，其值在该边的两个端点处分别为 1 和 0，而在其他边上为 0。这种构造方式直观地反映了矢量场沿边方向的变化。
对于六面体单元，一阶 Nédéc 基函数同样与边相关联，但形状函数更为复杂，通常涉及多项式插值。
高阶基函数与精度提升
当电磁场变化剧烈或需要更高精度时，一阶基函数可能不足以捕捉场的细节。此时，可采用高阶 Nédéc 基函数。高阶基函数通过在单元内部增加节点（如面节点或体节点）来提高多项式的阶数。
- 二阶及更高阶：高阶基函数能够更准确地表示曲率较大的场分布，减少网格细化带来的计算成本。
- 混合阶方法：在实际应用中，常采用混合阶方法，即在关键区域使用高阶基函数，在变化平缓区域使用低阶基函数，以平衡精度与计算效率。
基函数选择的工程建议
在实际电磁仿真软件或自研代码中，基函数的选择应遵循以下原则：
- 问题类型：对于涉及 $\nabla \times$ 的问题（如天线辐射、波导传输），必须使用 Nédéc 类基函数。对于静电或静磁问题，若仅求解标量势，可使用标量基函数；若求解矢量势，仍建议使用 Nédéc 基函数以避免虚假模。
- 网格类型：Nédéc 基函数支持四面体、六面体、楔形体等多种单元类型。选择时应确保基函数与网格拓扑结构兼容。
- 边界条件：正确施加边界条件依赖于基函数的性质。例如，在完美导体边界上，切向电场为零，这要求基函数在边界上的切向分量能够独立控制。
总结
矢量有限元基函数的选择是电磁场数值分析中的关键环节。Nédéc 基函数因其符合物理约束的连续性特性和良好的数值稳定性，成为解决旋度型电磁问题的标准工具。工程师应根据具体问题类型、精度要求和计算资源，合理选择基函数的阶数和类型，以确保仿真结果的准确性和可靠性。掌握基函数的数学性质及其物理意义，是深入理解和优化电磁仿真算法的基础。