为衰减地质雷达数值模拟时因截断区域厚度或者参数影响产生的倏逝波,提出了一种利用复频移完美匹配层(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%的计算时间.此外,该方法能应用于复杂地质结构模型的电磁总场计算,为地质雷达数值模拟提供了一种新的方法.展开更多
为了快速求解电磁散射问题中具有震荡性、奇异性、慢收敛性的索末菲积分,提出了一种利用分部外推算法加速索末菲尾部积分计算,并结合稳定双共轭快速傅里叶变换(stabilized biconjugate gradient fast Fourier transform,BCGS-FFT)算法...为了快速求解电磁散射问题中具有震荡性、奇异性、慢收敛性的索末菲积分,提出了一种利用分部外推算法加速索末菲尾部积分计算,并结合稳定双共轭快速傅里叶变换(stabilized biconjugate gradient fast Fourier transform,BCGS-FFT)算法求解电磁散射问题场分布情况的新方法.首先给出电场积分方程(electric field integral equation, EFIE)的表达形式,且在求解过程的索末菲积分中应用一种便捷的椭圆积分路径来最小化索末菲积分的震荡性与奇异性,在索末菲尾部积分使用Levin分部外推法来提高积分收敛速度,以此来快速填充并矢格林函数矩阵.然后对新方法进行了多种数值实验,验证算法的精确度,并对比了新方法与传统BCGS-FFT方法的计算效率,发现在保持相同计算精度的条件下,新方法可节省20%~37%的计算时间.该方法能应用于复杂散射体嵌入多层空间的电磁散射计算,为快速求解目标区域的电磁散射场提供了一种新的方法.展开更多
文摘为衰减地质雷达数值模拟时因截断区域厚度或者参数影响产生的倏逝波,提出了一种利用复频移完美匹配层(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%的计算时间.此外,该方法能应用于复杂地质结构模型的电磁总场计算,为地质雷达数值模拟提供了一种新的方法.
文摘为了快速求解电磁散射问题中具有震荡性、奇异性、慢收敛性的索末菲积分,提出了一种利用分部外推算法加速索末菲尾部积分计算,并结合稳定双共轭快速傅里叶变换(stabilized biconjugate gradient fast Fourier transform,BCGS-FFT)算法求解电磁散射问题场分布情况的新方法.首先给出电场积分方程(electric field integral equation, EFIE)的表达形式,且在求解过程的索末菲积分中应用一种便捷的椭圆积分路径来最小化索末菲积分的震荡性与奇异性,在索末菲尾部积分使用Levin分部外推法来提高积分收敛速度,以此来快速填充并矢格林函数矩阵.然后对新方法进行了多种数值实验,验证算法的精确度,并对比了新方法与传统BCGS-FFT方法的计算效率,发现在保持相同计算精度的条件下,新方法可节省20%~37%的计算时间.该方法能应用于复杂散射体嵌入多层空间的电磁散射计算,为快速求解目标区域的电磁散射场提供了一种新的方法.