计算电磁学中的能量动量模拟

在计算电磁学中，能量与动量模拟的核心是麦克斯韦方程组的守恒形式。对于无耗、线性、各向同性介质，电磁能量密度可写为：

在计算电磁学中，能量与动量模拟的核心是麦克斯韦方程组的守恒形式。对于无耗、线性、各向同性介质，电磁能量密度可写为：

[
u = \frac{1}{2}(\mathbf{E}\cdot\mathbf{D} + \mathbf{B}\cdot\mathbf{H})
]

能量流动由坡印廷矢量描述：

[
\mathbf{S} = \mathbf{E}\times\mathbf{H}
]

其守恒方程为：

[
\frac{\partial u}{\partial t} + \nabla\cdot\mathbf{S} = -\mathbf{J}\cdot\mathbf{E}
]

右端项表示焦耳损耗或外部源做功。动量方面，真空中的电磁动量密度为：

[
\mathbf{g} = \varepsilon_0 \mathbf{E}\times\mathbf{B} = \frac{\mathbf{S}}{c^2}
]

在介质中，Minkowski 动量密度 (\mathbf{D}\times\mathbf{B}) 与 Abraham 动量密度 (\mathbf{E}\times\mathbf{H}/c^2) 存在差异，数值模拟中需根据所研究的力类型选择一致的定义。

动量守恒通过麦克斯韦应力张量表达。真空中的应力张量为：

[
T_{ij} = \varepsilon_0\left(E_iE_j - \frac{1}{2}\delta_{ij}E^2\right)

  • \frac{1}{\mu_0}\left(B_iB_j - \frac{1}{2}\delta_{ij}B^2\right)
    ]

对闭合曲面 (A) 积分可得电磁力：

[
\mathbf{F} = \oint_A \mathbf{T}\cdot\mathbf{n},dA
]

该式是光镊、光操控和纳米光子学中计算光学力的基础。

数值离散中的守恒性挑战

计算电磁学常用时域有限差分（FDTD）、有限元法（FEM）、有限积分技术（FIT）和间断伽辽金时域法（DGTD）。不同方法对能量与动量的离散保持能力不同。

  • FDTD：Yee 网格天然交错，电场位于棱边，磁场位于面心。蛙跳时间积分在无耗、无源区域可近似保持能量守恒，但 PML 吸收边界、色散材料和非线性材料会引入人工损耗或能量交换。
  • FEM：需采用稳定混合元或矢量元，避免伪解。能量范数可用于自适应网格加密，动量通量通常通过后处理应力张量获得。
  • FIT/DGTD：基于离散外微分可精确保持离散斯托克斯定理，适合复杂几何与高阶精度模拟。

离散误差主要来自空间插值、时间积分和材料模型。例如，坡印廷矢量 (\mathbf{S}=\mathbf{E}\times\mathbf{H}) 中 (\mathbf{E}) 与 (\mathbf{H}) 在 Yee 网格上不同位置，必须插值到同一节点或面心，否则会引入相位误差。麦克斯韦应力张量涉及场分量乘积，插值误差会直接影响力计算精度。

典型模拟流程

  1. 几何与材料建模：定义计算域、边界条件、材料色散模型。
  2. 选择离散格式：根据频带、几何复杂度和精度要求选择 FDTD、FEM 或 DGTD。
  3. 场更新：按麦克斯韦旋度方程推进电场与磁场。
  4. 能量与动量后处理：计算 (u)、(\mathbf{S})、(\mathbf{g}) 和 (\mathbf{T})，并进行体积分或面积分。
  5. 验证与收敛性分析：检查能量守恒误差、动量守恒误差和网格收敛阶。

示例：FDTD 计算光镊力

考虑一个介质微粒置于高斯光束中。模拟步骤如下：

  • 使用 FDTD 更新三维电磁场，记录稳态周期内的场分布。
  • 在微粒外选取一个闭合立方曲面，确保曲面不穿过材料突变区。
  • 将 Yee 网格上的 (\mathbf{E}) 和 (\mathbf{H}) 插值到曲面面心。
  • 计算每个面元上的麦克斯韦应力张量 (\mathbf{T})。
  • 对曲面进行数值积分：

[
\mathbf{F} \approx \sum_k \mathbf{T}_k \cdot \mathbf{n}_k \Delta A_k
]

  • 与偶极子近似或米氏理论结果对比，验证力的大小与方向。

伪代码示意如下：

for n in range(Nt):
    update_H(E, H, dt, mu)
    update_E(E, H, dt, eps)
    if n % save_step == 0:
        E_s = interpolate_to_surface(E)
        H_s = interpolate_to_surface(H)
        T = maxwell_stress(E_s, H_s, eps, mu)
        F = integrate_surface(T, normals, areas)

跨学科应用与前沿

能量动量模拟已广泛用于多个领域：

  • 纳米光子学：计算光学力、光镊、光马达和微流控操控。
  • 等离子体物理：粒子模拟（PIC）中电磁场与带电粒子动量交换。
  • 超材料与拓扑光子学：分析负折射、动量带隙和自旋轨道耦合。
  • 无线输能与能量收集：评估近场能量传输效率和辐射动量。
  • 生物电磁学：研究细胞膜电场力与机械效应。

前沿方向包括：量子真空中的卡西米尔力计算、时空调制介质中的动量交换、以及基于离散外微分的保结构算法。这些方法强调在离散层面精确保持能量动量守恒，从而提升长期模拟稳定性与物理可信度。

验证与误差控制

为保证模拟可靠，建议：

  • 监测总能量随时间变化，无耗系统相对误差应低于 (10^{-4})。
  • 对动量通量进行曲面无关性测试，改变积分面位置结果应稳定。
  • 进行网格收敛性分析，力误差应随网格加密按预期阶数下降。
  • 谨慎处理 PML 与色散材料，必要时将材料能量纳入总能量核算。

能量动量模拟不仅是数值后处理，更是连接麦克斯韦理论与工程应用的关键桥梁。正确离散、一致插值和严格验证，是获得可信物理结果的前提。