超材料电磁分析

超材料是通过人工排列比波长小得多的金属图案（SRR或细线阵列）实现负介电常数或负透磁率的周期结构体。简单说，是由自然界的材料不可能实现的电磁响应。

完全的隐身斗篷还是理论上的梦想。Pendry（2006）的变换光学理论优美，但要实现广带、低损耗极其困难。不过已经实用化的应用有不少：

你抓住核心了。超材料分析的特点有3个：

是描述超材料电磁响应随频率变化的分散模型。有效介电常数可以写成Lorentz型：

类似地，SRR（分裂环共振器）这样的磁共鸣结构，有效透磁率也是Lorentz型：

完全正确。在$\omega_{0m} < \omega < \omega_{pm}$（磁等离子频率）范围内$\mu_{\text{eff}} < 0$。同时$\varepsilon_{\text{eff}} < 0$的范围重叠，折射率$n = \sqrt{\varepsilon_{\text{eff}} \mu_{\text{eff}}}$就变成负数——这就是Veselago（1968）预言的"左手超材料"。

形式本身不变。但结构关系（$\mathbf{D} = \varepsilon(\omega)\mathbf{E}$、$\mathbf{B} = \mu(\omega)\mathbf{H}$）具有频率分散性这是决定性差异。时间域FDTD中用ADE（附加微分方程）法或RC等效模型吸收分散性。

这里$k_0 = \omega/c$是自由空间的波数。$\varepsilon_r(\omega)$和$\mu_r(\omega)$遵循Drude-Lorentz模型变成复数的关键点。实部为负的区域，波要么不传播成为衰减波，要么作为左手的后向波传播。

聪明的疑问。这里就用Bloch-Floquet周期边界条件。只建模1个单元格，对相对的面加上相位差：

$\mathbf{a}$是格子矢量，$\mathbf{k}$是Bloch波数矢量。把$\mathbf{k}$沿着布里渊区的既约域（$\Gamma \to X \to M \to \Gamma$等）扫描，就能得到分散关系（能带图）。存在带隙的频率范围就是EBG（电磁带隙）。

完全是。光子晶体和超材料，数学上用的都是固体物理的带论框架。区别是，光子晶体的周期和波长可比，而超材料周期远小于波长（$a \ll \lambda$），所以能用等效连续介质描述为ε、μ。

用Nicolson-Ross-Weir（NRW）法的拓展。从超材料板的$S_{11}$（反射）和$S_{21}$（透射）反算阻抗$Z$和折射率$n$：

得到$n$和$Z$后，等效物性值就很简单：

聪明。最大的陷阱是分支选择的模糊性。$e^{jnk_0 d}$的反函数里有log函数，$n$的虚部出现$2m\pi/(k_0 d)$的任意整数$m$的不确定性。板厚$d$超过波长的一半就很难正确选择分支。对策在故障排查部分详述。

绝对不行。高频电磁场分析必须用边界元素（Nedelec/Whitney单元）否则会产生大量虚假模式（非物理的假解）。边界元素能自动保证电场的切向连续性，消除违反$\nabla \cdot \mathbf{D} = 0$的假解。

这样。边界元素中自由度放在单元的边上，直接补间矢量场的切向分量。弱形式是这样：

$\mathbf{N}_i$是边界元素的基函数。离散化后变成广义特征值问题：

$[S]$是curl-curl刚度矩阵，$[T]$相当于质量矩阵。代入周期边界条件，用Bloch波数$\mathbf{k}$作参数扫描就能得到能带图。

CST Studio Suite的时间域求解器就是FDTD。用ADE（附加微分方程）法处理分散性材料。Drude-Lorentz模型下，引入辅助变量$\mathbf{P}$（分极化电流）：

这个辅助微分方程在Yee网格上显式时间积分。FDTD的优点是广带S参数在1次计算中就能得到。入射高斯脉冲，时间域响应FFT就能一下子得到S11、S21的频特性。

步骤是这样的：

完全理解了。实务上不只算$\Gamma$点附近，要检查布里渊区全路径。有的方向只有部分带隙（方向带隙）和全方向都有隐隙（完全带隙），应用场景完全不同。

超材料特有的注意事项有几个：

超材料分析的标准工作流是这样：

10 GHz典型SRR设计参数是这样的：

共鸣频率的粗估计可以用LC回路模型：

完全同意。LC近似是初期设计用，最终还是要全波形仿真验证。实际上相互耦合、基板分散性、表面粗糙度损耗、制造偏差的影响，LC近似的10〜20%误差很常见。所以参数扫描很重要。

现场最常见的是下列3个：

由金属贴片+通孔+接地板组成的高阻抗表面（HIS）。反射相位为零的频率区间存在，该频带既不是PEC（完全电导体）也不是PMC（完全磁导体）而是特殊边界。天线能做成低轮廓，车载雷达和卫星通信天线已实用化。

超材料分析的检证有3个水平：

особенно注意的是，提取的$\varepsilon_{\text{eff}}$和$\mu_{\text{eff}}$的物理妥当性。下：

超材料分析对应的主要工具比较。各有擅长领域，用途选最合适的：

学术研究充分能用。MATLAB或Octave脚本建模风格，参数扫描自动化很方便。但没有GUI，能带图计算要自己实装。企业製品设计实际上HFSS或CST是标准。

那十成八九是虚假模式。原因和对策整理一下：

固有模式的電場分布可视化的最确実。物理的正電場共振结构体集中、虚假分析領域全体分布。$\nabla \cdot \mathbf{D}$処理计算、大虚假判定。

。NRW法屈折率$n$的実部$2m\pi/(k_0 d)$的不定性。以下的対策有効：

最安全的「板厚薄 + 条件」的組合。厚结构的S参数取場合、薄1層正$n$確認、初期値厚结构的$n$連続的追跡方法。Smith et al.（2005）的実装広使。

超材料分析周期境界条件、分散性材料模型、S参数反演3的特殊要素絡合、一確実理解積上的大事。分聞。