Certified anisotropic meshing of thin, high-aspect-ratio regions in device modeling

技术摘要：器件建模中薄、高长宽比区域的认证各向异性网格划分

问题陈述
器件仿真经常遇到嵌入在较大主体中的薄层（例如半导体薄膜、隔离片），这些区域具有极高的长宽比（$A = W/e_1$），范围从 $10^3$ 到 $10^5$。标准的器件网格生成器通常生成近各向同性的、符合 Delaunay 准则的单元，以满足经典的误差估计，但这对于扁平或针状单元会导致性能退化。为了通过 $m$ 个各向同性单元来解析厚度为 $e_1$ 的薄层，网格步长 $h$ 在所有方向上都必须很小。因此，未知量数量的缩放关系为 $N_{iso} \sim m^3 A^2$。对于一个典型的参考结构（10 mm 宽，30 µm 厚薄膜），各向同性网格需要大约 $1.6 \times 10^7$ 个节点，这使得计算成本因网格各向同性而非物理必要性而变得难以承受。

方法论
作者提出了一种结合各向异性结构化网格与严谨、无参考（reference-free）误差认证的两部分策略：

1. 各向异性结构化单元： 该方法并非采用各向同性细化，而是使用由粗糙横向网格（$n_x \times n_y$）和分级纵向网格（层内 $n_{z1}$ 个单元，主体内 $n_{z2}$ 个单元）组成的张量积进行域离散化。六面体通过 Kuhn–Freudenthal 细分法被分解为四面体。这种方法允许在层内使用高长宽比（10:1 至 1000:1）的扁平单元。至关重要的是，总节点数 $N_{aniso}$ 与长宽比 $A$ 无关，将该参考结构的需求降低至 567 个节点。

2. 基于互补能量的认证包络： 为了克服此类拉伸单元缺乏可计算误差界限的问题，该方法采用了基于 Prager–Synge 超圆（对偶能量法）的“认证”机制。在同一网格上计算两个解：

   • 一个符合 $P_1$ 连续性的电势 $V_h$（原问题解），它提供电导的上界（$G_{up}$）。
   • 一个平衡的、无散的 Raviart–Thomas 通量 $j_h$（对偶问题解，$RT_0$），它提供电导的下界（$G_{lo}$）。

   这两个界限包络了精确电阻 $R$，使得 $R_{lo} \le R \le R_{up}$。这两个界限的相对宽度 $\eta = (R_{up} - R_{lo})/R_{up}$ 作为一种保证的离散误差度量，无需参考解即可使用。该宽度取决于场的平面性，而非单元形状。

关键结果
该方法通过闭式解进行了验证，并应用于复杂的器件结构：

• 层平行流： 对于平行于薄层的流动，该方法无论单元如何拉伸都是精确的。对于高达 1000:1 的长宽比，认证宽度 $\eta$ 仍保持在舍入误差水平（$10^{-11}$ 至 $10^{-13}$）。各向异性网格仅用 567 个节点就实现了精确的面电阻，而各向同性网格则需要 $1.6 \times 10^7$ 个节点（未知量减少了 $2.8 \times 10^4$ 倍）。
• 二维场（点源）： 对于位于基底上 30 nm 薄膜上的点源，在考虑接地边界后，各向异性解与闭式镜像级数解的匹配度在 0.14% 以内。拉伸单元（高达 49:1）产生的误差仅约为 0.1%。
• 接触电阻认证： 对于薄膜-基底结构上的接触点，认证区间是嵌套的，并在细化过程中从 23% 缩小到 6.9%。收敛速率是由接触边缘奇异性驱动的，而非单元拉伸。
• 三层堆叠应用： 应用于一个包含 30 nm 隔离片的堆叠结构（占总厚度的 0.5%），该方法解析了掺杂对电势降的影响。高掺杂层承担了 0.13% 的压降，而低掺杂层承担了 13.4%。这两个情况的接触电阻认证包络在最粗糙的网格上也是不相交的，这证明该方法能够区分那些如果不经过校准，紧凑模型可能会忽略的物理效应。

意义与主张
论文声称，薄且宽的层是一个“数值瓶颈”而非物理瓶颈。其主要贡献在于证明了各向异性结构化单元可以实现计算成本与层长宽比的解耦，同时保持严谨的误差控制。

不同于现有的依赖于网格-解对齐常数的各向异性后验估计器（残差型或恢复型），该方法提供了一个保证的包络，且不含未知常数。该包络对于任何长宽比下的层平行流都是精确的，并且仅在场发生重新分布的地方（如接触点附近）报告残差误差。作者断言，这种方法允许在器件建模中进行电阻和面电阻的认证提取，而无需使用各向同性网格那令人望而生畏的节点数，也无需使用紧凑界面模型那无界的误差。该方法使用 Python 实现，并配有确定性脚本，所有结果均可通过提供的代码进行复现。