Yee 网格结构及其坐标定义
在电磁场数值分析与计算方法中，有限差分时间域法（FDTD）因其直接、高效且易于实现的特点，成为求解麦克斯韦方程组的主流方法之一。而构建 FDTD 模拟的基础，正是 Yee 网格结构。该结构由美国物理学家 John Yee 于 1966 年提出，其独特的交错网格设计不仅满足了数值稳定性要求，更在物理本质上实现了电场与磁场分量的自然耦合。
Yee 网格的核心特征在于将空间划分为均匀的立方体单元，并在每个单元的中心分别放置电场和磁场分量。与传统的笛卡尔网格不同，电场分量（$E$）和磁场分量（$H$）并非位于同一节点，而是相互交错排列。具体而言，电场矢量被放置在立方体的面心位置，而磁场矢量则位于立方体的边心位置。这种布局使得任意一个电场分量在空间上都被四个磁场分量所包围，反之亦然。
为了准确描述这种几何关系，必须建立明确的三维坐标定义体系。假设模拟区域在 $x, y, z$ 三个方向上的网格步长均为 $\Delta$，时间步长为 $\Delta t$。我们将空间原点 $(0,0,0)$ 设为参考点，定义第 $i, j, k$ 个网格单元的中心坐标分别为 $(i\Delta, j\Delta, k\Delta)$。
在 Yee 网格中，电场分量通常沿坐标轴方向分布。以 $x$ 方向的电场分量 $E_x$ 为例，它位于 $y-z$ 平面的中心。若考虑第 $i$ 个单元，其 $E_x$ 分量的空间坐标可表示为：
$$ (x, y, z) = (i\Delta, (j+0.5)\Delta, (k+0.5)\Delta) $$
同理，$y$ 方向的电场分量 $E_y$ 位于 $x-z$ 平面的中心，其坐标为：
$$ (x, y, z) = ((i+0.5)\Delta, j\Delta, (k+0.5)\Delta) $$
$z$ 方向的电场分量 $E_z$ 位于 $x-y$ 平面的中心，其坐标为：
$$ (x, y, z) = ((i+0.5)\Delta, (j+0.5)\Delta, k\Delta) $$
磁场分量的分布则完全交错于电场之间。以 $x$ 方向的磁场分量 $H_x$ 为例，它位于 $y-z$ 棱的中点，其空间坐标为：
$$ (x, y, z) = (i\Delta, (j+0.5)\Delta, k\Delta) $$
这种交错关系可以概括为：若电场分量位于坐标的“半整数”位置（即 $0.5$），则相邻的磁场分量位于“整数”位置，或者反之。例如，$E_x$ 的 $y$ 和 $z$ 坐标包含 $0.5$，而与其耦合计算 $H_y$ 和 $H_z$ 的 $y$ 和 $z$ 坐标则取整数。
这种特殊的坐标定义直接决定了数值计算的更新顺序。由于电场和磁场在空间上交错，它们在时间上也需要进行半拍（half-step）的更新。标准的 Yee 算法采用“偶 - 奇”时间推进法：
- 在 $n$ 时刻，先更新所有电场分量 $E^n$。
- 在 $n+0.5$ 时刻，利用更新后的电场计算磁场分量 $H^{n+0.5}$。
- 在 $n+1$ 时刻，利用更新后的磁场计算下一时刻的电场 $E^{n+1}$。
这种时间上的交错与空间上的交错相结合，使得麦克斯韦方程组中的旋度项能够自然地通过中心差分近似得到体现。例如，计算 $H_y$ 的更新时，需要用到 $E_x$ 和 $E_z$ 在空间上的差分，由于 $E_x$ 和 $E_z$ 在 $y$ 方向上交错排列，其差分恰好对应于 $H_y$ 所在位置的磁场变化率，从而在离散形式下保持了物理方程的对称性和守恒律。
理解 Yee 网格的坐标定义对于正确编写 FDTD 代码至关重要。在实际编程中，必须严格区分空间索引的奇偶性，以避免数值误差累积导致的不稳定。例如，在计算 $E_x$ 时，必须确保读取的是 $y$ 和 $z$ 方向上相邻的 $H$ 分量，这些 $H$ 分量的空间位置必须与 $E_x$ 的半整数坐标精确匹配。任何坐标定义的偏差，都可能导致模拟结果出现非物理的振荡或能量发散。
综上所述，Yee 网格结构通过巧妙的空间交错和半拍时间更新，为电磁波传播问题的数值求解提供了坚实的理论基础。掌握其坐标定义与更新机制，是从事电磁仿真、天线设计等领域数值计算工作的必备技能。