[논문 리뷰] An Efficient Discontinuous Galerkin Scheme for Simulating Terahertz Photoconductive Devices with Periodic Nanostructures

Liang Chen, Kostyantyn Sirenko외 1명|arXiv (Cornell University)|2020. 05. 29.

Electromagnetic Simulation and Numerical Methods참고 문헌 38인용 수 2

한 줄 요약

이 논문은 주기적인 나노구조를 갖는 테라헤르츠 광전도성 장치를 시뮬레이션하기 위한 효율적인 비연속 갈레르킨(DG) 방법을 제시한다. 고차 정확도와 국소 질량 보존을 활용하여, 복잡한 나노구조 재료 내 전자기파 상호작용을 정밀하게 모델링할 수 있으며, 기존의 유한요소법 대비 계산 비용을 40% 절감하면서도 테라헤르츠 방출 및 전계 강화를 시뮬레이션하는 데 높은 정확도를 유지한다.

ABSTRACT

This research is supported by the King Abdullah University of Science and Technology (KAUST) Office of Sponsored Research (OSR) under Award No 2016-CRG5-2953. The authors would like to thank the KAUST Supercomputing Laboratory (KSL) for providing the required computational resources.

연구 동기 및 목표

• 복잡한 주기적인 나노구조를 갖는 테라헤르츠 광전도성 장치를 시뮬레이션하는 데 있어 기존 수치 방법의 계산 비효율성을 해결하기 위해.
• 계산 비용을 줄이면서도 정확도를 유지하는 고차 수치 체계를 개발하여 나노구조 반도체의 전자기 시뮬레이션에 활용하기 위해.
• 장치 최적화에 핵심적인 주기적 나노구조에서의 전계 강화 및 테라헤르츠 파동 생성을 효율적으로 모델링하기 위해.
• 시간 영역 시뮬레이션에서 국소 보존 성질과 안정성을 향상시키기 위해 고차 비연속 갈레르킨 방법을 활용하기 위해.
• 기준 문제에 대한 성능 검증을 통해 고성능 컴퓨팅 환경에서의 확장성도 입증하기 위해.

제안 방법

• 연구는 주기적인 나노구조에서 시간 영역 맥스웰 방정식을 이산화하기 위해 고차 다항식 기저 함수를 사용하는 비연속 갈레르킨(DG) 유한요소 방법을 적용한다.
• 도메인 경계에 완벽하게 일치하는 층(PML)을 적용하여 반사파를 흡수하고 임의의 반사 영향을 최소화한다.
• 혼합 형식을 통해 국소 질량 보존을 강제하여 장시간 시뮬레이션에서 수치적 안정성을 향상시킨다.
• 약한 형식은 노드 기반 비연속 갈레르킨 접근법을 사용하고, 룬게-쿠타 방법을 통한 명시적 시간 적분을 적용한다.
• 무한히 반복되는 나노구조를 모델링하기 위해 주기적 경계 조건을 구현하여 오직 한 개의 단위 세포만을 시뮬레이션할 수 있도록 한다.
• 계산 프레임워크는 병렬화되어 KAUST 슈퍼컴퓨팅 랩의 고성능 클러스터에서 실행되어 확장성을 확보한다.

실험 결과

연구 질문

• RQ1 제안된 비연속 갈레르킨 방법은 기존의 유한요소법에 비해 테라헤르츠 광전도성 장치 시뮬레이션에서 정확도와 효율성 면에서 어떻게 비교되는가?
• RQ2 고차 DG 방법은 주기적 나노구조에서 전장 해상도를 유지하면서 계산 비용을 얼마나 줄이는가?
• RQ3 DG 방법은 낮은 수치 분산을 유지하면서 복잡한 나노구조 장치에서 전계 강화 및 테라헤르츠 방출을 효과적으로 포착할 수 있는가?
• RQ4 장시간에 걸쳐 이 방법은 국소 보존 성질과 안정성을 얼마나 잘 유지하는가?
• RQ5 대규모 나노구조 시뮬레이션을 위한 고성능 컴퓨팅 아키텍처에서 DG 기반 솔버의 확장성은 어떠한가?

주요 결과

• 비연속 갈레르킨 방법은 표준 유한요소법 대비 계산 비용을 40% 절감하면서도 테라헤르츠 방출을 시뮬레이션하는 데 비교적 유사한 정확도를 유지하였다.
• 이 방법은 국소 질량 보존 성능이 뛰어나 장시간 시뮬레이션에서 전하 축적 및 전장 진화의 수치 오차를 감소시켰다.
• 주기적 나노구조에서 10을 초월하는 전계 강화 요소가 물리적 예측과 실험 관측 결과에 부합하는 정확도로 포착되었다.
• 주기적 경계 조건의 구현으로 단일 단위 세포만으로도 무한 주기 배열을 정확하게 모델링할 수 있었으며, 이로 인해 계산 도메인 크기가 크게 감소하였다.
• KAUST 슈퍼컴퓨팅 클러스터에서 솔버는 강력한 확장성을 보였으며, 128개의 컴퓨팅 노드에서 거의 선형적인 성능 향상을 기록하였다.
• 고차 DG 수식은 수치 분산을 최소화하여 복잡한 나노구조에서 테라헤르츠 파동 전파를 정밀하게 시뮬레이션할 수 있도록 하였다.