대규모 PMUT 어레이 시뮬레이션을 위한 혼합 차원 축소와 고정밀 DG 스펙트럴 엘리먼트 프레임워크

초록

본 논문은 PMUT(피에조 전기 마이크로머신 초음파 트랜스듀서) 어레이의 전기·기계·음향 복합 거동을 효율적으로 모델링하기 위해, 각 소자를 고유 진동 모드 기반의 차원 축소 모델로 표현하고 이를 비정형 메쉬와 고차원 DG 스펙트럴 엘리먼트 방법으로 구현한 통합 시뮬레이션 프레임워크를 제안한다. SPEED 오픈소스 플랫폼에 구현된 이 방법은 METIS 기반 도메인 분할, MPI 병렬 처리 등을 활용해 대규모 어레이에서도 높은 정확도와 확장성을 확보한다.

상세 분석

이 연구는 PMUT 어레이 시뮬레이션에서 가장 큰 병목인 전·기·음향 결합 해석을 두 단계로 나누어 해결한다. 첫 번째 단계는 개별 PMUT의 3차원 구조-전기-기계 고유값 문제를 풀어 주요 진동 모드와 고유 주파수를 추출하고, 이를 선형 결합 형태의 저차원 모델(모드 좌표 qₘ)로 축소한다. 이렇게 얻은 차원 축소 모델은 전압 구동 시 전기‑기계 결합 계수 ηₘ과 전기 용량 C를 통해 전압 입력 ϕ(t)와 직접 연결되며, 전송(TX) 단계에서는 전압에 의해 압력 파동을 발생시키고, 수신(RX) 단계에서는 압력에 의해 전압이 유도되는 양방향 동작을 동일한 방정식 체계로 기술한다. 두 번째 단계는 이러한 저차원 PMUT 모델을 2차원 멤브레인 경계 조건으로 변환하여, 대규모 3차원 음향 영역에 적용한다. 여기서 DG-SEM(Discontinuous Galerkin Spectral Element Method)을 사용해 비정형 메쉬와 비정합 인터페이스를 자유롭게 구성할 수 있다. 특히, 내부 음향 영역(Ω_in)과 외부 음향 영역(Ω_out)을 DG 인터페이스(Γ_I)로 연결함으로써 서로 다른 폴리노미얼 차수와 메쉬 해상도를 동시에 적용할 수 있다. 이는 고주파 파동 전파에 필요한 고차원 근사와, 멀리 떨어진 경계에서는 저차원 근사만으로도 충분히 정확한 해를 얻어 계산 비용을 크게 절감한다. 또한 PML(Perfectly Matched Layer) 구현을 위한 추가 변수(ψ, Φ)와 감쇠 프로파일 ζ_i(x_i)를 도입해 인공 반사 없이 무한 공간을 효과적으로 모사한다. 수치 해석은 명시적 Newmark 시간 적분 스킴을 사용해 고성능 병렬 환경에서 확장성을 확보한다. 메쉬 분할은 METIS 기반 그래프 파티셔닝으로 수행되며, MPI 통신 최적화를 통해 수천 개 코어에서도 거의 선형적인 스케일링을 달성한다. 실험 결과는 (1) 차원 축소 모델이 전체 3D FEM 결과와 비교해 0.5 % 이하의 오차를 보이며, (2) DG-SEM 기반 음향 해석이 전파 속도와 직접ivity 패턴을 정확히 재현하고, (3) 코어 수를 128에서 4096까지 증가시켜도 효율적인 병렬 가속을 유지함을 입증한다. 이러한 접근은 기존 고해상도 FEM이 감당하기 어려운 수천 개 이상의 PMUT을 포함한 대형 어레이 설계와 최적화에 필수적인 도구가 될 것으로 기대된다.