복잡한 형상의 포아송 방정식 다중격자 해법

초록

본 논문은 변환 광학 아이디어를 이용해 복잡한 물리 영역을 균일한 계산 영역으로 매핑하고, 변환된 좌표계에서 유전율 텐서를 수정하여 유한 차분(FD) 방식으로 포아송 방정식을 푼 뒤, 구조화된 격자에 최적화된 기하학적 다중격자(GMG) 알고리즘을 적용한다. 2차 정확도를 유지하면서 SOR 대비 수십 배의 속도 향상을 보이며, 해석적·수치적 사례를 통해 방법의 유효성을 검증한다.

상세 분석

이 연구는 전통적인 유한 차분(FD) 방식이 균일 격자에서만 효율적으로 동작한다는 한계를 변환 광학(transformation optics) 개념으로 극복한다. 물리적 비균일 영역을 좌표 변환 x→x′ 으로 균일한 계산 영역에 대응시키고, 변환 행렬 J 의 행렬식과 역행렬을 이용해 새로운 유전율 텐서 ε 를 정의한다(ε = det(J) J⁻¹ ε′ (Jᵀ)⁻¹). 여기서 물리 영역에서는 ε′ = I (단위 텐서)로 두어, 변환된 계산 영역에만 공간적으로 변하는 ε 가 존재한다. 이 접근법은 물리적 경계조건과 소스항을 그대로 유지하면서도, 격자 자체는 정규 사각형 형태를 유지하므로 기존의 고성능 구조화 격자용 다중격자(MG) 솔버를 그대로 적용할 수 있다.

2차원 경우, Jacobian J 의 2×2 요소를 명시적으로 전개하고, ε_xx, ε_xy, ε_yx, ε_yy를 식(12)–(15)와 같이 계산한다. 전기장 E 는 전위 φ 의 기울기로 표현되고, 변환된 포아송 방정식은 ∇·(ε∇φ) = −ρ 형태가 된다. 중앙 차분을 이용한 2차 정확도 스키마를 도입해 격자점 (i,j) 에서의 전위 업데이트 식(18)을 유도하고, Gauss‑Seidel 완화(ω∈(1,2))를 스무딩 연산으로 활용한다.

다중격자 알고리즘은 전형적인 V‑cycle을 따르며, 잔차 r_h 를 제한 연산(I_H^h)으로 코스그리드에 전송하고, 코스그리드에서 얻은 교정 Ψ_H 를 보간 연산(I_h^H)으로 다시 미세 격자에 적용한다. 이때 시스템 행렬 A_h 는 변환된 ε와 ρ에 의해 완전히 정의되므로, 전통적인 기하학적 다중격자에서 요구되는 복잡한 비균일 격자 생성 과정을 생략한다. 결과적으로 메모리 사용량과 전처리 비용이 크게 감소하고, 구조화된 데이터 레이아웃 덕분에 OpenMP·MPI 기반 병렬화가 용이해진다.

실험에서는 반원형(세미‑annulus) 영역에 대한 해석적 매핑을 사용해 ε_xx, ε_yy가 비대칭적으로 변함을 확인하고, 전위 해와 수치 해 사이의 L∞, L1, L2 오차가 Δx² 에 비례함을 통해 2차 정확도를 입증한다. 또한 동일한 격자 규모(N=201)에서 SOR에 비해 다중격자 풀이가 약 100배 빠른 실행 시간을 기록한다. 수치적 매핑이 필요한 경우에도 Jacobian을 수치적으로 계산하고 동일한 절차를 적용함으로써 방법의 일반성을 보인다.

이 접근법의 주요 강점은 (1) 복잡한 경계와 비균일 메쉬를 단순한 정규 격자로 변환함으로써 기존 FD·MG 코드와의 호환성을 유지, (2) 변환 과정에서 물리적 양(전위, 전기장)을 정확히 보존, (3) 2차 정확도와 최적 O(N) 복잡도를 동시에 달성한다는 점이다. 한편, 변환 함수가 고차원에서 복잡하거나 Jacobian이 급격히 변하는 경우 ε 텐서의 급격한 변동이 수치적 불안정을 초래할 수 있으며, 이러한 상황에서는 추가적인 스무딩이나 비선형 스케일링 기법이 필요할 수 있다. 전반적으로 변환 광학 기반의 좌표 매핑은 FEM 기반 비구조 격자 해법에 비해 구현 난이도가 낮고, 대규모 병렬 시뮬레이션에 적합한 실용적 대안으로 평가된다.