Trefftz Discontinuous Galerkin methods for scattering by periodic structures

기술 요약: 주기적 구조에 의한 산란을 위한 트레프츠 불연속 갤러킨 방법 (Trefftz Discontinuous Galerkin Methods)

문제 정식화
본 논문은 주기적 회절 격자(diffraction gratings)에 의한 평면파 산란의 수치 근사 문제를 다룬다. 물리적 모델은 횡전기(TE) 또는 횡자기(TM) 편광 가정을 따르는 시간 조화 맥스웰 방정식으로부터 유도된 2차원 헬름홀츠 방정식이다. 영역은 하나의 방향으로 불변하는 주기적 구조로 구성되며, 투과 가능한 영역(조각별로 일정한, 복소수 값을 가질 수 있는 유전 매질)과 투과 불가능한 영역(완전 전기 도체, PEC)을 모두 포함한다.

문제는 절단된 단위 셀 $\Omega = (0, L) \times (-H, H) \setminus D$ 상의 경계값 문제(BVP)로 정식화된다. 경계 조건은 다음과 같다:

• 구조와 입사파의 주기성을 고려하기 위한 수직 측면($\Gamma_{left}, \Gamma_{right}$)에서의 준주기성(Quasi-periodicity).
• 투과 불가능한 장애물 $D$ 상의 디리클레(Dirichlet) 조건.
• 송출되는 파동이 산란체로부터 멀어지도록 보장하며, 소머펠트 방사 조건을 정확하게 대체하기 위한 수평 절단 경계($\Gamma_{\pm H}$) 상의 디리클레-투-네만(Dirichlet-to-Neumann, DtN) 조건.

방법론
저자들은 트레프츠 불연속 갤러킨(Trefftz Discontinuous Galerkin, TDG) 방법, 구체적으로는 DtN-TDG 스킴을 제안한다. 주요 방법론적 특징은 다음과 같다:

1. 이산 공간(Discrete Space): 근사 공간은 헬름홀츠 방정식의 조각별 해(트레프츠 함수)로 구성된다. 구체적으로, 이 방법은 요소별로 정의된 평면파 기저 함수(복소 지수 함수)를 채택한다.
2. 수치 플럭스(Numerical Fluxes): 약한 연속성(weak continuity)을 강제하기 위해 요소 인터페이스 상에서 수치 플럭스를 활용한다.
   • 내부 면(interior faces)에서는 불연속적인 재료 계수에 적응된 표준 플럭스가 사용된다.
   • 인공 경계 $\Gamma_{\pm H}$ 상에서는 절단된 DtN 연산자($T^{\pm}_M$)를 포함하는 새로운 플럭스가 유도된다.
   • PEC 경계 $\Gamma_D$ 상에서는 표준 디리클레 플럭스가 적용된다.
3. 해석적 조립(Analytical Assembly): 이 방법의 중요한 특징은 모든 선형 시스템 항목(행렬 및 로드 벡터)을 해석적으로 계산할 수 있다는 점이다. 기저 함수가 복소 지수 함수이고 DtN 연산자가 푸리에 영역에서 대각 성분을 가지기 때문에, 메쉬 면에 대한 적분과 DtN 맵을 위한 푸르리에 계수 계산에 수치 적분(numerical quadrature)이 필요하지 않다.
4. 절단(Truncation): DtN 연산자를 정의하는 무한 급수는 차수 $M$에서 절단된다. 논문은 이 절단이 연속 연산자의 특성을 보존하는 조건을 확립한다.

주요 기여 및 이론적 결과

• 명시적 안정성 추정(Explicit Stability Estimate): 저자들은 헬름홀츠 해에 대한 새로운 완전 명시적 안정성 추정치(정리 3.5)를 유도한다. 이 유계(bound)는 재료 점프에 대해 강건하며, 파수(wavenumber), 기하학적 구조, 그리고 가장 가까운 레일리-우드 이상점(Rayleigh–Wood anomaly, DtN 연산자의 단사성을 상실하는 지점)까지의 거리 $\delta$에 명시적으로 의존한다.
• Rellich 항등식: 이 DtN 경계를 갖는 준주기적 설정에 특화된 새로운 Rellich 항등식(23)이 유도된다. 이 항등식은 문제의 적절성(well-posedness)과 안정성 추정을 증명하는 초석이 된다. 이는 레일리-우드 이상점과 복소 재료 매개변수를 처리하기 위해 기존의 결과를 확장한 것이다.
• 적절성(Well-Posedness): 비포획(non-trapping) 가정(장애물에 대한 기하학적 조건 및 유전율의 단조성)과 파수가 레일리-우드 이상점이 아니라는 가정 하에, 연속 문제와 절단된 이산 문제가 모두 적절함을 증명한다.
• 강성(Coercivity) 및 준최적성(Quasi-Optimality): 이산 방법이 특정 메쉬 의존 세미 노름(seminorm)에서 강성을 가짐을 보이며, 이를 통해 오차에 대한 준최적성 유계(proposition 4.2)를 도출한다.
• 수렴 분석: 평면파 공간의 근사 특성을 논하며, 이들이 매끄러운 해에 대해 조각별 다항식 공간보다 점근적으로 더 빠른 수렴율을 제공함을 언급한다. 또한, 특이점(모서리)에서 발생하는 수치적 불안정성을 다루고, 이에 대한 해결책으로 감쇠 평면파(evanescent plane waves)의 사용을 제안한다.

수치 결과
이 방법은 MATLAB으로 구현되었으며 다양한 구성에 대해 테스트되었다:

• 매끄러운 해(Smooth Solutions): 균질한 재료 간의 평평한 인터페이스(무손실 및 손실 매질)와 다층 구조 문제를 대상으로 했을 때, 평면파 방향의 수($p$)에 대해 **지수적 수렴(exponential convergence)**을 보여준다.
• 특이 해(Singular Solutions): 모서리 특이점이 포함된 문제(예: 모서리가 있는 인터페이스, 투과 불가능한 장애물)의 경우, 수렴 속도가 느려지기는 하지만 방법론이 여전히 효과적임을 보여준다. 오차는 요소 경계와 기하학적 특이점 근처에 집중된다.
• 비유일성 시나리오(Non-Unique Scenarios): BVP가 무수히 많은 해(유도 모드)를 갖는 경우를 조사한다. 수치 결과는 이 방법이 해의 아핀 공간(affine space) 내의 하나의 해로 수렴함을 나타내며, 특정 해는 이산화 방식에 따라 달라지며 임계값 근처의 입사각에 민감하다.
• 절단 오차(Truncation Error): DtN 절단($M$)에 의한 오차가 충분히 큰 $M$에 대해 이산화 오차에 의해 지배되며 빠르게 감소함을 실험을 통해 확인한다.

의의 및 주장
본 논문은 DtN-TDG 방법을 주기적 산란 문제에 대한 강건하고 효율적인 도구로 포지셔닝한다. 이 방법의 주요 의의는 다음과 같다:

1. 쿼드러처 프리 조립(Quadrature-Free Assembly): 모든 시스템 항목을 해석적으로 계산할 수 있는 능력은 수치 적분 오차를 제거하고, 특히 고주파 문제에서 계산 비용을 줄여준다.
2. 명시적 안정성: 재료 점프와 이상점 근접성을 고려한 명시적 안정성 유계(32)의 도출은 실제 구현에서 종종 결여되는 이론적 보장을 제공한다.
3. 복잡한 기하학적 구조 처리: 이 방법은 불연속적인 재료 매개변수를 가진 투과 가능 및 투과 불가능 주기적 구조를 성공적으로 처리한다.

저자들은 이 방법이 매끄러운 해에 대해 다항식 방법보다 빠른 수렴율을 제공하지만, 특이점이 존재하는 경우 수치적 불안정성에 주의해야 하며, 이를 위해 감쇠 평면파가 잠재적인 개선책으로 제안된다고 언급한다. 재현성을 높이기 위해 코드는 무료로 공개되어 있다.