[논문 리뷰] Elastic Wave Eigenmode Solver for Acoustic Waveguides

Nathan Dostart, Yangyang Liu외 1명|arXiv (Cornell University)|2017. 07. 19.

Mechanical and Optical Resonators참고 문헌 47인용 수 3

한 줄 요약

이 논문은 유한차분 고유모드 해석기법을 제안하며, 비균일한 횡단면을 가진 탄성파 가로파이프에서 Yee 방법에 유사한 격자배치를 사용하여, 유한차분 기반의 해석기법을 통해 유도 모드, 누출 모드 및 복사 모드를 정확하게 시뮬레이션할 수 있다. 이는 완벽하게 매칭된 계층(PML)을 활용하여 검증되었으며, 해석적 해석 및 3차원 유한요소법(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.

연구 동기 및 목표

• 옵티오메카닉스 및 RF MEMS 응용 분야에서 칩스케일의 음향 가로파이프를 설계하기 위한 효율적이고 정확한 수치 도구의 필요성을 해결한다.
• 전체 3차원 해석기기의 한계를 극복하기 위해 계산 비용을 줄이면서도 정확도를 유지하는 2+1D 고유모드 수식을 개발한다.
• 완벽하게 매칭된 계층(PML)을 사용하여 누출 모드와 복사 손실을 시뮬레이션함으로써 실제 가로파이프 종단 및 에너지 누출 모델링을 가능하게 한다.
• Yee 유사 격자배치를 사용함으로써 전자기 해석기기와의 원활한 통합을 가능하게 하여, 음향-광학 및 옵티오메카닉스 장치의 공동 설계를 지원한다.
• 이sovotropic 비균일 매질에서의 벡터형 탄성파 방정식을 위한 강력하고 두 번째 차수 정확도를 갖춘 유한차분 방법을 제공한다.

제안 방법

• 선형적이고 등방성인 탄성파 방정식을 고유값 문제로 재구성하며, 전파 상수 β를 입력으로, 고유모드 주파수 ω를 출력으로 한다.
• 모든 장성분에 대해 두 번째 차수 정확도를 유지하는 격자배치를 사용한 스태그너드-그리드 유한차분 방법을 적용하여 벡터형 파동 방정식을 이산화한다.
• 자기대칭, 반대칭 경계 조건을 포함한 자유, 고정, 대칭 및 반대칭 경계 조건을 도입하여 기하학적 대칭성을 활용하고 시뮬레이션 영역 크기를 줄인다.
• 영역 경계에 완벽하게 매칭된 계층(PML)을 도입하여 반사 없이 외부로 향하는 파동을 흡수하고, 반사 없이 누출 모드를 시뮬레이션할 수 있도록 한다.
• 3차원 가로파이프 문제를 2차원 고유값 문제로 줄이기 위해 2+1D 수식을 적용하여 특정 β에서의 횡단면 모드에 집중한다.
• 전자기 해석기기와의 통합을 위해 Yee 격자 유사 스태그너드 배열을 적용하여 옵티오메카닉스 시스템에서의 다물리학 공동 시뮬레이션을 가능하게 한다.

실험 결과

연구 질문

• RQ1 비균일하고 등방성인 매질에서 임의의 횡단면을 가진 가로파이프에 대해, 고유모드 해석기기의 유한차분 기법을 탄성파 방정식에 효과적으로 적용할 수 있는가?
• RQ2 해석적 해석 및 3차원 FEM 해석기기와 비교했을 때, 이 해석기기가 기본 음향 모드의 고유주파수와 모드 형상에 대해 얼마나 정확하게 계산할 수 있는가?
• RQ3 이 해석기기에서 PML가 음향 가로파이프의 복사 손실과 누출 모드를 얼마나 정확하게 시뮬레이션할 수 있는가?
• RQ4 이 해석기기가 대칭 및 반대칭 경계 조건을 효율적으로 처리하여 대칭적인 가로파이프 기하구조에서 계산 비용을 줄일 수 있는가?
• RQ5 Yee 격자 구조를 사용함으로써 전자기 해석기기와의 통합이 얼마나 잘 이루어지는가? 이는 옵티오메카닉스 장치의 공동 설계에 얼마나 유리한가?

주요 결과

• 해석기기는 공간 이산화에서 두 번째 차수 정확도를 확보하였으며, 단순 기하구조에서의 수렴 연구와 해석적 해석과의 비교를 통해 검증되었다.
• 실리카-실리콘 가로파이프에서의 준-로브 및 레일리 유사 모드에 대해 고유주파수와 모드 형상을 높은 정확도로 계산하였으며, 해석적 예측과 0.1% 이내의 오차로 일치하였다.
• PML의 구현이 외부로 향하는 파동을 효과적으로 흡수하여, 복사 손실이 있는 누출 모드의 정확한 시뮬레이션을 가능하게 하였으며, 3차원 FEM 해석기기(COMSOL)와의 비교를 통해 검증되었다.
• 변동하는 실리콘 장벽 폭을 가진 누출 가로파이프에서의 전파 손실은 고유모드 해석기기와 3차원 FEM 해석기기 간에 뛰어난 일치를 보였으며, PML 정확도를 확인하였다.
• 누출 모드의 복사장에 대한 푸리에 분석 결과, 실리카 필름에서의 준-로브 및 준-레일리 파에 해당하는 두 개의 주요 복사 성분이 확인되었다.
• 해석기기는 실리콘 장벽을 통해 음향 에너지가 복사 모드로 터널링되는 것을 성공적으로 시뮬레이션하였으며, 이론적 기대와 일치하는 波front 및 변위 패턴을 보였다.