[논문 리뷰] Virtual elements for Maxwell's equations

L. Beirão da Veiga, Franco Dassi외 2명|arXiv (Cornell University)|2021. 02. 01.

Advanced Numerical Methods in Computational Mathematics참고 문헌 41인용 수 8

한 줄 요약

이 논문은 일반 다면체 메esh를 사용하는 시간에 의존하는 맥스웰 방정식을 위한 저차수 가상요소법(VEM)을 제시한다. 이 방법은 노드, 변, 면 가상요소 공간의 de Rham 수열을 활용한다. 반연속 및 전연속 스킴 모두 최적 수렴 속도를 달성하며, 왜곡되고 다면체 메쉬에서의 수치 실험을 통해 자석장의 수렴성과 기계 정밀도 수준의 정확한 수렴성 유지가 검증되었다.

ABSTRACT

We present a low order virtual element discretization for time dependent Maxwell's equations, which allow for the use of general polyhedral meshes. Both the semi- and fully-discrete schemes are considered. We derive optimal a priori estimates and validate them on a set of numerical experiments. As pivot results, we discuss some novel inequalities for de Rahm sequences of nodal, edge, and face virtual element spaces.

연구 동기 및 목표

• 일반 다면체 메쉬를 지원하는 시간에 의존하는 맥스웰 방정식을 위한 저차수 가상요소 이산화 방법을 개발한다.
• 3차원에서 H(curl) 및 H(div) 적합 근사에 대해 가상요소 공간에 de Rham 수열 프레임워크를 확장한다.
• 반연속 및 전연속 수식에 대해 최적의 사전 오차 추정치를 수립한다.
• 왜곡되고 다면체 메쉬에서의 수치 실험을 통해 방법의 타당성을 검증하며, 수렴 특성과 자석장의 무차원성 성질 유지 여부를 평가한다.
• 특이점이나 자기유체역학(MHD) 문제에 대한 고차수 확장 및 응용을 위한 기초를 마련한다.

제안 방법

• 이 방법은 3차원 다면체 요소에 대해 노드(H1), 변(H(curl)), 면(H(div)) 가상요소 공간의 de Rham 복합체를 기반으로 한 가상요소 수식을 사용한다.
• 특히 컬 및 발산 연산자에 대해 일관성과 안정성을 확보하기 위해 안정화 항을 구성한다.
• 반연속 스킴은 공간적 VEM 근사를 사용하고, 전연속화를 위해 후진 오일러 방법을 사용한 시간 적분을 적용한다.
• 방법은 맥스웰 방정식의 혼합 변분 수식을 기반으로 하며, 이는 이산 수준에서 연속 문제의 구조를 유지한다.
• 가상요소 공간은 다항식 부분공간 위로의 국소 사영을 통해 구성되며, 자유도는 모멘트와 경계 적분을 통해 정의된다.
• 보간 성질과 가상요소 맥락에서의 de Rham 수열에 대한 새로운 부등식을 사용하여 사전 오차 추정치를 유도한다.

실험 결과

연구 질문

• RQ1 일반 다면체 메쉬에서 시간에 의존하는 맥스웰 방정식에 대해 최적 수렴성을 갖는 저차수 가상요소 방법을 구성할 수 있는가?
• RQ2 3차원에서 H(curl) 및 H(div) 적합 구조를 유지하기 위해 가상요소 공간의 de Rham 수열을 어떻게 설계할 수 있는가?
• RQ3 이산 변분 수식의 안정성과 일관성을 보장하기 위해 어떤 안정화 기법을 사용할 수 있는가?
• RQ4 이 방법은 자석기구장 B의 무차원성 성질을 고정밀도로 유지하는가?
• RQ5 공간과 시간을 동시에 균일하게 미세화할 경우, 이 방법의 수렴 특성은 어떠한가?

주요 결과

• 반연속 및 전연속 가상요소 스킴이 전기장과 자석장 모두에 대해 L2 노름에서 최적 수렴 속도를 달성하며, 수치 실험을 통해 확인되었다.
• 전기장 E의 경우, 공간 수렴 속도는 약 1.0이며 시간 수렴 속도 역시 약 1.0이며, h 및 τ가 미세해질수록 오차가 감소한다.
• 자석장 B의 경우, 공간 수렴 속도는 약 1.0이며 시간 수렴 속도 역시 약 1.0이며, 모든 메쉬 유형에서 오차가 일관되게 감소한다.
• 이산 자석장 B_h는 발산이 약 10^-13에서 10^-15 수준을 유지하여 기계 정밀도 수준에서 거의 정확한 무차원성 성질을 보존함을 시사한다.
• 형태가 규칙적인 voro 메쉬의 경우, 무작위 메쉬보다 시스템 행렬의 조건수(조건수)가 낮아지며, 이는 더 작은 대수적 오차와 더 나은 발산 보존을 초래한다.
• 자석장 오차는 공간 이산화에 의해 주로 결정되며, 실험에서 시간 해상도가 비교적 낮기 때문에 시간 이산화 오차는 덜 두드러진다.