[논문 리뷰] An explicit P1 finite element scheme for Maxwell's equations with constant permittivity in a boundary neighborhood

Larisa Beilina, Vitoriano Ruas|arXiv (Cornell University)|2018. 08. 31.

Advanced Numerical Methods in Computational Mathematics참고 문헌 26인용 수 7

한 줄 요약

이 논문은 경계층 내에서 고정된 허용도를 가진 맥스웰 방정식을 해결하기 위해 표준 선형 요소와 명시적 시간 스텝을 사용하는 안정적인 명시적 P1 유한요소 스킴을 제안한다. 고전적인 CFL 조건과 발산 포착 안정화를 통해 전기장 근사에 대해 자연스러운 노름에서 최적의 O(h) 수렴이 증명된다.

ABSTRACT

This paper is devoted to the complete convergence study of the finite-element approximation of Maxwell's equations in the case where the magnetic permeability is constant. Standard linear finite elements for the space discretization are combined with a well-known explicit finite-difference scheme for the time discretization. The analysis applies to the particular case where the dielectric permittivity has a constant value outside a sub-domain, whose closure does not intersect the boundary of the problem-definition domain. Optimal convergence results are established in natural norms under reasonable assumptions, provided a classical CFL condition holds. A numerical validation of the theoretical results is provided.

연구 동기 및 목표

• 경계에 인접한 영역에서 고정된 허용도를 가진 맥스웰 방정식에 적용된 명시적 P1 유한요소법에 대한 완전한 수렴 분석을 수립하는 것.
• 특정 물성 조건 하에서 표준 연속 P1 유한요소가 명시적 시간 이산화와 조합될 경우 최적 수렴을 달성할 수 있음을 보여주는 것.
• 노드 기반 유한요소에서 발생하는 허위 해를 억제하기 위해 변분 형식에 발산 자유 조건을 통합하는 것.
• 허용도가 엄밀히 1이 아니지만 양수 상한 이하로 유계일 경우로의 수렴 분석을 확장하는 것.
• 허용도 함수에 대한 최소한의 정규성 가정 하에서도 스킴의 강건성을 검증하는 것, 특히 내부 영역에서 ε < 1 인 경우 포함.

제안 방법

• 비정렬 메esh에서 표준 연속 P1 유한요소를 통한 공간 이산화로 단순성과 효율성 확보.
• 명시적 중심 유한차분 스킴을 통한 시간 이산화로 매 단계에서 선형 시스템을 풀지 않아도 되는 명시적 시간 스텝 구현.
• 전기장의 발산 자유 조건을 강제하는 수정된 변분 형식을 통한 안정화로 허위 모드 방지.
• 계산 효율성 향상을 위해 납입 질량 행렬 사용으로 명시적 시간 적분 유지.
• 비다항식 허용도 함수의 수치적 적분에 삼각법칙 적용, 미약한 정규성 가정 하에서도 안정성 유지.
• ε < 1 인 경우를 다루기 위해 ε_min을 포함한 수정된 보조 문제를 도입하여 변분 형식을 적응시키며, 강건성과 수렴성 확보.

실험 결과

연구 질문

• RQ1 경계 근처에서 허용도가 일정할 경우, 명시적 P1 유한요소 스킴이 맥스웰 방정식에 대해 최적 수렴을 달성할 수 있는가?
• RQ2 시간 영역 맥스웰 해법에서 노드 기반 유한요소에서 발생하는 허위 해를 효과적으로 억제하는 방법은 무엇인가?
• RQ3 영역 내 ε < 1 인 경우에 안정성과 수렴성을 유지하기 위해 변분 형식에 어떤 수정이 필요한가?
• RQ4 이 명시적 P1 스킴에서 CFL 조건이 안정성에 충분한 조건이 되는 조건은 무엇인가?
• RQ5 허용도가 일정하지 않지만 양수 상한으로 유계일 경우로의 수렴 분석을 어느 정도까지 확장할 수 있는가?

주요 결과

• 허용도가 경계 근처에서 일정할 경우, 고전적인 CFL 조건 하에서 전기장에 대해 L2 노름에서 최적의 O(h) 수렴이 달성된다.
• ε_min > 2/3 일 경우, 내부에서 허용도가 1보다 작아도 스킴은 안정적이고 수렴성을 유지하며 보조 문제의 강건성 확보.
• 수정된 변분 형식을 통한 발산 자유 조건 강제가 효과적으로 이루어져, H(curl)-형식 요소를 사용하지 않아도 허위 해가 제거된다.
• 허용도가 일정하지 않지만 양수 상한으로 유계일 경우로의 수렴 결과는 확장 가능하며, 이 경우 수렴 속도는 보조 문제의 해의 정규성에 따라 달라진다.
• 비스무스 허용도 함수의 적분에 삼각법칙을 사용하더라도 안정성이 유지되며, 수렴 순서가 떨어지지 않는다.
• 분석은 두 차원으로 직접 일반화되며 본질적인 수정 없이 동일한 정성적 수렴 결과 유지.