[发明专利]基于时域有限差分的双曲超材料的完美匹配层方法

说明书

基于时域有限差分的双曲超材料的完美匹配层方法，属于计算电磁仿真领域。包括以下步骤：1)设置FDTD算法的仿真参数；2)确定仿真精度和离散网格数；3)在双曲介质中加入正弦波源；4)更新电场和磁场分量；5)绘制磁场分布，并分析稳定性；6)数据后处理，若步骤5)中得到的磁场场量的结果不收敛，或误差较大，则说明PML不稳定，需要对PML进行修正。在双曲材料中传播的电磁波在传统PML和COMSOL频域PML仿真中存在数值发散的问题，本发明解决时域传统PML和COMSOL中频域PML不能吸收双曲材料中的电磁波的问题，并表现出很强的稳定性，数值结果进一步验证了PML改进技术的有效性。

技术领域

本发明属于计算电磁仿真领域，尤其是涉及一种基于时域有限差分的双曲超材料的完美匹配层方法。

背景技术

近年来，随着微纳加工技术的不断进步，在亚波长尺度下物质与电磁波的相互作用引起了人们的广泛关注。电磁超构材料是一种由特征尺寸远小于工作波长的结构单元所构成，具有自然界材料不存在的电磁特性的人工材料。作为电磁超构材料的一个重要分支，双曲超材料由于其独特的近场电磁波操控特性成为人们研究的焦点。通过改变组成双曲超材料内部的表面等离子体激元激发强度和方向的调控，从而使双曲超材料呈现出独特的色散特性。双曲超材料已经在很多领域得到应用，包括亚波长成像，广场局域和增强自发辐射等。另一方面，双曲超表面作为一种特殊新型的平面超材料，具有双曲色散特性，在理论和应用上也与双曲超材料有着诸多共性。与块体结构的双曲超材料相比，由于其纵向维度尺寸的大幅度减小，双曲超表面可以将电磁波限制在二维平面上，易于片上器件的集成，因而在近些年受到研究人员的极大关注。

在电磁波传播的无界空间中，完美匹配层(PML)被开发为一种有效的技术来截断计算域，PML可以吸收传入的电磁波，PML在计算区域边界引入了虚拟各向异性有损介质，在一定条件下，计算区域空间与虚拟有损介质完全匹配，而不会在吸收边界和常规材料的交界处引起反射，计算空间中的外部电磁波可以进入PML层而无需反射，介质中的电磁波逐渐衰减，从而有效地吸收向外传播的电磁波。经典PML在不同类型的介质中表现出了超高的吸收效率，适用于线性色散材料，理论上也适用于负折射率材料的模拟，但是通过分析，PML在很多具有负折射率的材料中是不稳定的。比如在双曲超材料中，由于波的群速和相速的矢量积在PML的方向为负，所以在双曲超材料中传播的电磁波进入完美匹配层时会表现出不稳定性。

发明内容

本发明的目的在于针对现有技术存在的上述问题，提供基于Drude模型的双曲材料PML吸收算法，可提高双曲材料中电磁波吸收的稳定和效率的基于时域有限差分(Finitedifference time domain，FDTD)的双曲超材料的完美匹配层方法。

本发明包括以下步骤：

1)设置FDTD算法的仿真参数；

2)确定仿真精度和离散网格数；

3)在双曲介质中加入正弦波源；

4)更新电场和磁场分量；

5)绘制磁场分布，并分析稳定性；

6)数据后处理，若步骤5)中得到的磁场场量的结果不收敛，或误差较大，则说明PML不稳定，需要对PML进行修正。

在步骤1)中，所述设置FDTD算法的仿真参数，设置仿真时间步数为1500，PML层厚度为9网格数，仿真频率为1.5×10¹⁴Hz。

在步骤2)中，所述确定仿真精度和离散网格数的具体步骤可为：每最小波长采样点数为N_lambda为10，x和y方向上的离散网格数分别为：3λ×6λ。

在步骤3)中，所述在双曲介质中加入正弦波源，并更新磁场的具体步骤和原理为：