[논문 리뷰] A Virtual Finite Element Method for Two Dimensional Maxwell Interface Problems with a Background Unfitted Mesh

Shuhao Cao, Long Chen외 1명|arXiv (Cornell University)|2021. 03. 08.

Advanced Numerical Methods in Computational Mathematics참고 문헌 51인용 수 33

한 줄 요약

이 논문은 비구조적 배경 불일치 메esh에서 2차원 매ク스웰 인터페이스 문제를 해결하기 위한 새로운 가상 요소 방법(VEM)을 제안한다. 최대 각 조건을 만족하는 가상 삼각분할 위에 새로운 H(curl)-준수 가상 공간을 도입함으로써, 해면 근처에서 메쉬의 이방성에 영향을 받지 않고 일阶 최적 수렴을 달성하며, 낮은 정(regularity)과 비일치 이산화의 과제를 극복한다.

ABSTRACT

A virtual element method (VEM) with the first order optimal convergence order is developed for solving two-dimensional Maxwell interface problems on a special class of polygonal meshes that are cut by the interface from a background unfitted mesh. A novel virtual space is introduced on a virtual triangulation of the polygonal mesh satisfying a maximum angle condition, which shares exactly the same degrees of freedom as the usual H(curl)-conforming virtual space. This new virtual space serves as the key to prove that the optimal error bounds of the VEM are independent of high aspect ratio of the possible anisotropic polygonal mesh near the interface.

연구 동기 및 목표

• 인터페이스를 기준으로 계수가 불연속인 2차원 매크스웰 인터페이스 문제를 해결하기 위한 강력한 수치적 방법을 개발하는 것.
• 비구조적, 불일치 다각형 메쉬에서 최적 수렴을 달성하는 과제를 해결하는 것. 특히 인터페이스 근처에서 메쉬가 극도로 이방성이 될 수 있음.
• 요소 형태에 영향을 받지 않고 H(curl)-준수성과 자유도를 유지하는 가상 요소 공간을 구성하는 것.
• 메쉬의 이방성과 계수의 점프에 대해 강건한 사전 오차 추정을 수립하는 것.

제안 방법

• 각 다각형 요소의 가상 삼각분할 위에 최대 각 조건을 만족하는 표준 Nédélec 공간의 부분공간으로서 새로운 가상 요소 공간을 구성한다.
• 이 가상 공간은 표준 H(curl)-준수 VEM 공간과 동일한 자유도를 공유하므로 기존 VEM 수식과의 호환성을 확보한다.
• H(curl)-타원형 인터페이스 문제에 대해 대칭 변분형식을 사용하며, 비일치 근사에 대응하기 위해 안정화 항을 도입한다.
• 자유도로 정의된 이산 조화 확장을 사용하여 경계항을 제어하고 강건한 오차 추정을 달성하기 위해 새로운 분석 기법을 적용한다.
• 분석은 요소의 형태 정규 하위삼각형에서의 추정, 트레이스 부등식 및 Poincaré 유형 추정에 기반한다.
• 배경 삼각형 메쉬가 인터페이스에 의해 잘린 원형 인터페이스 문제에 대해 수치적으로 검증한다.

실험 결과

연구 질문

• RQ1 가상 요소 방법이 불일치하고 이방성 다각형 메쉬에서 매크스웰 인터페이스 문제에 대해 최적 수렴을 달성할 수 있는가?
• RQ2 요소 형태 왜곡에 강건한 가상 공간에서 H(curl)-준수성과 자유도를 어떻게 유지할 수 있는가?
• RQ3 메쉬의 이방성이 낮은 정(regularity) 해를 가진 인터페이스 문제의 VEM 수렴 속도에 어떤 영향을 미치는가?
• RQ4 새로운 가상 공간 구성으로 인해 요소의 종횡비에 종속되지 않는 오차 분석이 가능할 수 있는가?

주요 결과

• 제안된 VEM은 H(curl) 노름에서 일계 최적 수렴을 달성하며, 수치 실험을 통해 L2 오차와 컬 오차의 수렴률이 각각 약 1.00임을 확인하였다.
• 해면 근처에서 요소가 극도로 이방성이 되더라도 수렴률이 일관되게 유지되어, 모든 메쉬 세분화에서 최적 수렴을 유지함을 입증하였다.
• 오차 추정은 다각형 요소의 종횡비에 영향을 받지 않으며, 이는 불일치 메쉬 방법에 있어 핵심 기여이다.
• 수치 결과는 물성 매개변수(α+ = 10 또는 100, β+ = 10 또는 100)가 변화하더라도 수렴률이 안정적으로 유지됨을 보여주며, 방법의 강건성을 입증한다.
• 사전 오차 분석을 통해 수렴률이 요소 형태나 이방성과 무관하게 정확한 해의 정규성과 메쉬 크기만에 의존함을 증명하였다.