为衰减地质雷达数值模拟时因截断区域厚度或者参数影响产生的倏逝波,提出了一种利用复频移完美匹配层(complex frequency shifted perfectly matched layer,CFS-PML)作为吸收边界,并结合频域高阶有限元算法(higher order finite element...为衰减地质雷达数值模拟时因截断区域厚度或者参数影响产生的倏逝波,提出了一种利用复频移完美匹配层(complex frequency shifted perfectly matched layer,CFS-PML)作为吸收边界,并结合频域高阶有限元算法(higher order finite element method,HO-FEM)求解电磁总场分量分布情况的新方法.该方法改善了截断域介质本构张量矩阵的特性,利用矩阵频移参量α吸收倏逝波,提高掠射角处吸收性能.数值实验结果表明,与传统PML相比,结合CFS-PML的HO-FEM方法精度得到明显提高,反射误差额外降低10~30 dB,节省了17%~30%的计算时间.此外,该方法能应用于复杂地质结构模型的电磁总场计算,为地质雷达数值模拟提供了一种新的方法.展开更多
文中将CPML引入3维弱无条件稳定算法HIE-FDTD中,详细推到了CPML在3维弱无条件稳定HIE-FDTD中的差分公式。为了验证CPML在3维HIE-FDTD中的吸波性能,建立了数值计算模型,并将CPML的吸波性能同其它几种常用的吸收边界条件进行了比较。结果...文中将CPML引入3维弱无条件稳定算法HIE-FDTD中,详细推到了CPML在3维弱无条件稳定HIE-FDTD中的差分公式。为了验证CPML在3维HIE-FDTD中的吸波性能,建立了数值计算模型,并将CPML的吸波性能同其它几种常用的吸收边界条件进行了比较。结果显示,当将CPML层数设置为8时,其最大反射误差为-72 d B,远低于传统FDTD方法的反射误差。另外,当匹配层参数设置为α=0.05,可以在一个较大范围内选取κmax和σmax来实现最佳误差,从而使得在选值时易于预测反射情况。展开更多
文摘为衰减地质雷达数值模拟时因截断区域厚度或者参数影响产生的倏逝波,提出了一种利用复频移完美匹配层(complex frequency shifted perfectly matched layer,CFS-PML)作为吸收边界,并结合频域高阶有限元算法(higher order finite element method,HO-FEM)求解电磁总场分量分布情况的新方法.该方法改善了截断域介质本构张量矩阵的特性,利用矩阵频移参量α吸收倏逝波,提高掠射角处吸收性能.数值实验结果表明,与传统PML相比,结合CFS-PML的HO-FEM方法精度得到明显提高,反射误差额外降低10~30 dB,节省了17%~30%的计算时间.此外,该方法能应用于复杂地质结构模型的电磁总场计算,为地质雷达数值模拟提供了一种新的方法.
文摘文中将CPML引入3维弱无条件稳定算法HIE-FDTD中,详细推到了CPML在3维弱无条件稳定HIE-FDTD中的差分公式。为了验证CPML在3维HIE-FDTD中的吸波性能,建立了数值计算模型,并将CPML的吸波性能同其它几种常用的吸收边界条件进行了比较。结果显示,当将CPML层数设置为8时,其最大反射误差为-72 d B,远低于传统FDTD方法的反射误差。另外,当匹配层参数设置为α=0.05,可以在一个较大范围内选取κmax和σmax来实现最佳误差,从而使得在选值时易于预测反射情况。
