정규분해 기반 연속 라그랑주 유한요소를 이용한 맥스웰 고유값 문제 해결

초록

본 논문은 정규분해 기법을 활용해 (H_0^s(\mathrm{curl};\Omega)) 공간을 (H^{s+1}(\Omega)) 벡터 포텐셜과 (H^{s+1}(\Omega)) 스칼라 포텐셜의 합으로 분해한다. 이를 기반으로 연속 (C^0) 라그랑주 고차 유한요소를 사용한 새로운 맥스웰 고유값 해법을 제안하고, 스펙트럴 수렴과 최적 오류 추정(특히 (H(\mathrm{curl})) 노름에서의 차수 (k) 및 (k+1) 요소에 대한 수렴률)을 엄밀히 증명한다. 수치 실험을 통해 제안 방법이 스푸리어스 모드 없이 정확히 수렴함을 확인한다.

상세 분석

이 연구는 기존의 Nédélec 계열 고차 엣지 요소가 구현 복잡성과 자유도와의 이중성 문제를 안고 있다는 점에 주목한다. 저자들은 (H_0^s(\mathrm{curl};\Omega)) 공간을 (H^{s+1}(\Omega)) 벡터 포텐셜 (u)와 (H^{s+1}(\Omega)) 스칼라 포텐셜 (p) 의 그래디언트로 분해하는 새로운 정규분해 정리를 증명한다. 핵심은 경계조건을 (u\times n=0) (접선 제로)로 완화하고, (s>1/2) 인 경우에만 (u,p\in H^{s+1}(\Omega)) 임을 보이는 것이다. 이 정규분해는 기존 (H^1_0) 분해와 달리 포텐셜 함수가 더 높은 정규성을 가지므로, 연속 라그랑주 요소로 직접 근사화가 가능해진다.
수치 스키마는 차수가 (k) 인 라그랑주 요소로 (u)를, 차수가 (k+1) 인 요소로 (p)를 각각 근사한다. 이렇게 하면 자유도는 전통적인 엣지 요소와 비슷하거나 더 적게 유지하면서도, 스펙트럴 이론에 따라 연산자 (T) (역연산자와 정규분해 연산자의 조합)가 콤팩트함을 이용해 고유값 근사에 대한 무스푸리어스 보장을 얻는다. 오류 분석에서는 (|E-E_h|_{H(\mathrm{curl})}\le C h^{\min(k,s)}) 와 같은 최적 차수 수렴을 도출하고, 고유값 오차는 (O(h^{2\min(k,s)})) 임을 증명한다.
또한, 저자들은 복합 메쉬나 가중 발산 안정화와 같은 기존의 복잡한 기법을 전혀 사용하지 않고, 표준 FEM 패키지에 바로 적용 가능한 형태를 제시한다는 점에서 실용성이 높다. 수치 실험에서는 구형, L형, 그리고 복합 곡률을 가진 다면체에 대해 (k=1,2,3) 차수의 요소를 적용했으며, 이론적 수렴률과 거의 일치하는 로그‑로그 수렴 그래프를 확인한다. 특히, 스푸리어스 0 모드가 전혀 나타나지 않아 스펙트럴 정확도와 물리적 신뢰성을 동시에 확보한다는 결론을 얻는다.