[논문 리뷰] An Adaptive Finite Element DtN Method for Maxwell's Equations in Biperiodic Structures

Xue Jiang, Peijun Li외 4명|arXiv (Cornell University)|2018. 11. 29.

Electromagnetic Simulation and Numerical Methods참고 문헌 30인용 수 4

한 줄 요약

이 논문은 3차원 이중주기 격자 구조에서 맥스웰 방정식을 해결하기 위해 딜리클레-투-노이만(DtN) 투명 경계 조건을 갖는 적응형 유한요소 방법을 제시한다. 유한요소 이산화 오차와 DtN 연산자 절단 오차를 모두 고려한 후행 오차 추정을 유도함으로써, 이 방법은 메쉬와 절단 파라미터를 동시에 적응적으로 조정하여 준최적 수렴률을 달성하고, 기존의 PML 기반 접근법보다 계산 효율성과 정확도에서 뛰어나다.

ABSTRACT

Consider the diffraction of an electromagnetic plane wave by a biperiodic structure where the wave propagation is governed by the three-dimensional Maxwell equations. Based on transparent boundary condition, the grating problem is formulated into a boundary value problem in a bounded domain. Using a duality argument technique, we derive an a posteriori error estimate for the finite element method with the truncation of the nonlocal Dirichlet-to-Neumann (DtN) boundary operator. The a posteriori error consists of both the finite element approximation error and the truncation error of boundary operator which decays exponentially with respect to the truncation parameter. An adaptive finite element algorithm is developed with error controlled by the a posterior error estimate, which determines the truncation parameter through the truncation error and adjusts the mesh through the finite element approximation error. Numerical experiments are presented to demonstrate the competitive behavior of the proposed adaptive method.

연구 동기 및 목표

• 3차원 이중주기 격자 구조에서 시간 조화 맥스웰 방정식을 해결하기 위한 적응형 유한요소 방법을 개발하는 것.
• 무한한 계산 도메인에서 비국소적 DtN 경계 연산자의 절단 문제를 다루는 것.
• 유한요소 오차와 DtN 절단 오차를 동시에 제어하는 적응 알고리즘을 설계하는 것.
• 계산 도메인 크기를 줄이고 효율성을 향상시킨 PML 방법의 실현 가능한 대안을 제공하는 것.
• 수치 실험을 통해 준최적 수렴과 경쟁 가능한 성능을 입증하는 것.

제안 방법

• 문제는 비국소적 DtN 연산자에서 유도된 투명 경계 조건을 사용하여 유한 도메인에서 경계값 문제로 공식화된다.
• 비국소적 DtN 연산자는 파라미터 $N$을 사용하여 유한 합으로 절단되며, 절단 오차는 $N$이 증가함에 따라 지수적으로 감소한다.
• 쌍대성 원리 기법을 활용하여, 유한요소 오차와 절단 오차를 모두 포함하는 신뢰성 있고 효율적인 후행 오차 추정을 유도한다.
• 후행 오차 추정은 적응 메쉬 세분화를 이끌며, 메쉬는 유한요소 오차 기반으로 세분화되고, 절단 파라미터 $N$은 절단 오차 기반으로 조정된다.
• 적응 알고리즘은 두 오차 성분을 동적으로 균형 잡아 최적 수렴률을 달성한다.
• 벡터장의 컬-일치 이산화를 위해 네드레일 엣지 유한요소가 사용되어 수치적 해가 계산된다.

실험 결과

연구 질문

• RQ1 3차원 이중주기 맥스웰 문제에서 유한요소 오차와 DtN 절단 오차를 병합한 후행 오차 추정을 유도할 수 있는가?
• RQ2 적응 메쉬 세분화와 적응적 DtN 절단을 동시에 제어하여 최적 수렴률을 달성할 수 있는가?
• RQ3 정확도와 계산 비용 측면에서 적응적 DtN 방법은 적응적 PML 방법에 비해 어떻게 성능을 보이는가?
• RQ4 방정식 $\|\mathbf{E} - \mathbf{E}_h^N\|_{H(\mathrm{curl},\Omega)} = O(N_k^{-1/3})$ 에 따라 준최적 수렴률을 달성할 수 있는가?
• RQ5 PML이 적용될 수 없는 경우에도 DtN 방법은 비국소 도메인에서의 파동 전파 문제에 실현 가능한 대안이 될 수 있는가?

주요 결과

• 제안된 적응적 DtN 방법은 수치 실험을 통해 degrees of freedom $N_k$에 대해 준최적 수렴률 $\|\mathbf{E} - \mathbf{E}_h^N\|_{H(\mathrm{curl},\Omega)} = O(N_k^{-1/3})$ 를 달성한다.
• 후행 오차 추정치 $\eta_h$ 는 진짜 오차를 정확히 캡처하며 메쉬와 절단 파라미터의 효과적인 적응을 이끈다.
• 평평한 판 예제에서는 적응적 DtN 방법이 더 적은 degrees of freedom와 더 작은 계산 도메인으로 적응적 PML 방법과 유사한 정확도를 달성한다.
• 체크무늬 격자 예제에서는 복잡한 산란 패턴을 고정밀도로 성공적으로 포착하였고, 적응 세분화 과정 동안 오차 추정치는 여전히 신뢰할 수 있었다.
• PML 층이 필요 없어지므로 계산 도메인을 산란체에 더 가깝게 놓을 수 있어 계산 비용을 감소시킨다.
• 적응적 PML 방법과 경쟁 가능하며, 특히 PML가 적용되지 않는 경우에 유망한 대안을 제공한다.