摘要
为保证瞬变电磁场晚期信号计算正确,在时域有限差分正演中通常采用低频近似,即将磁场散度方程显式的包含在电磁场迭代方程中,因而H_z的求解方式与H_x、H_y不同,导致常规CPML吸收边界无法施加.本文重新推导了CPML介质中的迭代方程,并将其与常规介质中的迭代方程在形式上进行了统一,提出了一种适用于瞬变电磁低频近似下三维FDTD的CPML边界条件及施加方法.首先采用均匀半空间模型来验证本方法的有效性,发现模拟结果与解析解吻合较好,且模拟用时为采用Dirichlet边界FDTD用时的一半;而且,反射误差实验验证了本文的CPML边界吸收效率较高.之后,模拟了典型的三层模型,发现本方法的模拟结果与线性数字滤波解吻合程度较好.最后采用本方法计算了航空瞬变电磁带地形的复杂模型,实验证明了本文提出的CPML具有广泛的适用性.
- 单位