这篇开始介绍理想匹配层数学模型
理想匹配层,也可以叫做理想吸收边界, 按照字面意思理解:就是在这个边界上无反射,完美被吸收(在任何 角度和任意频率的电磁波),并且电磁波会在里面做衰减,不会透射出去.
先回顾下一些基本知识:
斯涅尔定理:
菲涅尔定理:
如果介质1和介质2都是LHI材料, 如果要满足
因为是LHI材料的
由此可见,理想匹配层不可能是LHI材料,
如果材料不是isotropic的话,我们就需要改写Maxwell equation,因为在各个方向上的介电常数/磁导率是不一样的,我们需要引入tensor的概念来描述(传输矩阵那篇有做过介绍),再写下:
在晶格里面,我们一般只看x,y,z三个方向介电常数,(介电常数只出现在对角线上) 定义有三种介质:
在anisotropic材料中,菲涅尔定理和斯涅尔定理改写成如下:
如果我们选择的材料
这样可以得到和入射角度无关的项;
如果a=b的话,
所以理想匹配层的材料是uniaxial材料, 并且满足a=b=1/c
其中a,b,c是复数
上面的图只是在z方向做的计算,我们一般需要在x,y,z三个方向上都满足理想匹配层条件,我们先定义一个参数s:
整合后得到:
3D方向的PML层可以这么看:)
接下来需要把PML整合进Maxwell方程组:
没有整合理想匹配层Maxwell方程可以表示如下:
把PML边界加进去后,在含PML边界的Maxwell方程可以表示为:
然后我们用计算机计算[s]
,简单解释下上面这些方程的意思:
在x<0时,
x>0时,会出现随x增大而增大的衰减项;
下一篇会先简单介绍下FDM(Finite difference method)有限差分法基本原理 主要是泰勒展开式和矩阵算符.这个是仿真软件CST用的最基本算法.
参考资料: Dr.Raymond Rumpf在youtube上的CEM课程
文章知识点与官方知识档案匹配,可进一步学习相关知识算法技能树首页概览34187 人正在系统学习中
声明:本站部分文章及图片源自用户投稿,如本站任何资料有侵权请您尽早请联系jinwei@zod.com.cn进行处理,非常感谢!