탄성파 고유모드 해석을 위한 고성능 격자 기반 솔버

초록

본 논문은 선형·등방성 재료의 2‑차원 단면을 갖는 음향 파이프라인(웨이브가이드)의 고유모드를 계산하기 위한 유한차분(FD) 스테거드 그리드 방식을 제안한다. 자유·고정·대칭·반대칭 경계조건과 완전 흡수층(PML)을 구현해 누설모드까지 효율적으로 해석한다. 전자기학에서 널리 쓰이는 Yee 스킴을 차용해 2차 정확도를 유지하면서도 희소 행렬 기반 고유값 해석으로 빠른 수치를 제공한다. 다양한 단순 기하와 상용 FEM 결과와의 비교를 통해 정확성을 검증하였다.

상세 분석

이 연구는 기존 3‑D 전자기·음향 전산해석이 갖는 높은 계산 비용을 회피하기 위해, 구조가 z‑축을 따라 불변인 경우를 2+1 D 문제로 축소한다는 점에서 의미가 크다. 저자들은 연속적인 탄성 파동 방정식을 ρ·∂²ₜ u = ∇·(C:∇ₛ u) 형태로 시작해, 질량 밀도 ρ와 강성 텐서 C를 각각 √ρ와 1/√ρ 로 정규화함으로써 Hermitian 연산자를 도출하고, ω²를 고유값으로 하는 표준 고유값 문제로 변환한다.

격자 설계는 전자기학의 Yee 스킴을 차용하되, 탄성 문제의 특성을 반영해 주응력(σ_xx, σ_yy, σ_zz)과 전단응력(σ_xy, σ_xz, σ_yz)을 서로 다른 격자 위치에 배치한다. 구체적으로, 셀 중심에 주응력을, 셀 면 중앙에 변위 성분(u_x, u_y, u_z)을, 그리고 면의 교차점에 전단응력을 두는 3‑D 스테거드 그리드를 2‑D 단면으로 압축한다. 이러한 배치는 미분 연산자를 중앙차분 형태로 구현하면서도 물성 파라미터가 격자 규모보다 급격히 변할 때도 2차 정확도를 유지한다.

경계조건 구현은 물리적 경계(고정, 자유)와 수학적 대칭(대칭, 반대칭) 조건을 격자점의 존재·삭제 혹은 값 고정으로 처리한다. 특히, PML은 복소 좌표 변환을 통해 z‑축 전파 상수 β에 대한 복소 지수 감쇠를 적용함으로써 방사형(누설) 모드의 복사 손실을 정확히 모델링한다.

수치 해법은 행렬 H를 희소 실수 대칭 행렬(손실이 없을 경우)로 구성하고, MATLAB의 ‘eigs’와 같은 Arnoldi 기반 전이역법을 이용해 원하는 β에 대한 몇 개의 저차 고유모드만을 효율적으로 추출한다. 이는 전통적인 FEM이 전체 모드 스펙트럼을 계산해야 하는 비효율성을 크게 개선한다.

검증 사례로는 직사각형, 원형, 그리고 복합 재료 코어‑클래드 구조의 단면을 선택해, 이론적 해석식(예: 라플라스·베셀 함수 해)과 상용 FEM(Comsol) 결과와 비교하였다. 격자 간격을 반으로 줄일 때 오차가 O(Δx²)로 수렴함을 보였으며, 10 GHz까지의 고주파에서도 수십 배 빠른 계산 속도와 1 % 이하의 상대 오차를 달성했다.

이 논문의 주요 기여는 (1) 탄성 파동 방정식에 대한 스테거드 격자 기반 고유모드 해석 프레임워크 제시, (2) 다양한 물리적·수학적 경계조건과 PML을 일관되게 구현, (3) 전자기학과의 격자 일치를 통해 다중물리 시뮬레이션(예: 광‑음향 상호작용)에서의 연동성을 확보한 점이다. 향후 비등방성 재료, 비선형 손실, 그리고 3‑D 변조 구조에 대한 확장 가능성도 논의된다.