[论文解读] Elastic Wave Eigenmode Solver for Acoustic Waveguides

Nathan Dostart, Yangyang Liu+1 more|arXiv (Cornell University)|Jul 19, 2017

Mechanical and Optical Resonators参考文献 47被引用 3

一句话总结

本文提出了一种用于非均匀截面声学波导中弹性波的有限差分特征模态求解器，采用类似于电磁学中Yee方法的交错网格格式。该方法可精确模拟导波、漏泄模态及辐射模态，结合完美匹配层（PMLs），经与解析解和三维FEM求解器对比验证，且针对光机械器件设计中与电磁求解器的集成进行了优化。

ABSTRACT

A numerical solver for the elastic wave eigenmodes in acoustic waveguides of inhomogeneous cross-section is presented. Operating under the assumptions of linear, isotropic materials, it utilizes a finite-difference method on a staggered grid to solve for the acoustic eigenmodes of the vector-field elastic wave equation. Free, fixed, symmetry, and anti-symmetry boundary conditions are implemented, enabling efficient simulation of acoustic structures with geometrical symmetries and terminations. Perfectly matched layers are also implemented, allowing for the simulation of radiative (leaky) modes. The method is analogous to eigenmode solvers ubiquitously employed in electromagnetics to find waveguide modes, and enables design of acoustic waveguides as well as seamless integration with electromagnetic solvers for optomechanical device design. The accuracy of the solver is demonstrated by calculating eigenfrequencies and mode shapes for common acoustic modes in several simple geometries and comparing the results to analytical solutions where available or to numerical solvers based on more computationally expensive methods.

研究动机与目标

• 解决在光机械和射频MEMS应用中设计芯片级声学波导对高效、精确数值工具的需求。
• 通过开发一种2+1D特征模态公式，克服全三维求解器的局限性，降低计算成本，同时保持精度。
• 利用完美匹配层（PMLs）实现漏泄模态和辐射损耗的模拟，以实现对真实波导端面和能量泄漏的建模。
• 通过采用类似Yee的交错网格，实现与电磁求解器的无缝集成，支持声光与光机械器件的协同设计。
• 为各向同性、非均匀介质中的矢量弹性波方程提供一种鲁棒的二阶精度有限差分方法。

提出的方法

• 将线性、各向同性弹性波方程表述为广义特征值问题，以传播常数β为输入，模态频率ω为输出。
• 采用交错网格有限差分方法对矢量波方程进行离散化，保持所有场分量的二阶精度。
• 实现自由、固定、对称与反对称边界条件，以利用几何对称性，减小仿真域尺寸。
• 在域边界引入完美匹配层（PMLs），以吸收出射波，实现无反射的漏泄模态模拟。
• 采用2+1D公式将三维波导问题简化为二维特征值问题，聚焦于给定β下的截面模态。
• 应用类似Yee网格的交错排列，与电磁求解器对齐，支持光机械系统中多物理场协同仿真。

实验结果

研究问题

• RQ1 能否有效将有限差分特征模态求解器适配于具有任意截面的非均匀、各向同性波导中的弹性波方程？
• RQ2 与解析解和三维FEM解相比，该求解器在计算基模态态频率和模态形状方面的精度如何？
• RQ3 该求解器中PMLs在模拟声学波导中辐射损耗和漏泄模态方面的准确性如何？
• RQ4 该求解器在对称波导几何结构中能否高效处理对称与反对称边界条件，以降低计算成本？
• RQ5 该求解器在使用Yee网格时，与电磁求解器的集成效果如何，适用于光机械器件的协同设计？

主要发现

• 求解器在空间离散化中实现了二阶精度，通过收敛性研究与简单几何结构的解析解对比得到验证。
• 在二氧化硅-硅波导中，准-Love模态与类Rayleigh模态的特征频率和模态形状计算精度极高，与解析预测的误差在0.1%以内。
• PML的实现成功吸收了出射波，实现了对漏泄模态及辐射损耗的精确模拟，经与三维FEM求解器（COMSOL）对比验证。
• 在硅障碍层宽度可变的漏泄波导中，模式求解器与三维FEM求解器在传播损耗计算上表现出极佳的一致性，证实了PML的准确性。
• 对漏泄模态辐射场的傅里叶分析揭示了两个主要辐射分量，分别对应于二氧化硅板中的准-Love波与准-Rayleigh波。
• 求解器成功模拟了声能通过硅障碍层向辐射模态的隧穿过程，波前与位移分布与理论预期一致。