[논문 리뷰] A discrete Weber inequality on three-dimensional hybrid spaces with application to the HHO approximation of magnetostatics

Florent Chave, Daniele A. Di Pietro외 1명|arXiv (Cornell University)|2020. 07. 07.

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

한 줄 요약

이 논문은 다면체 메esh에 대해 정의된 3차원 하이브리드 다항식 공간에서 이산 웨버 부등식을 수립하여 자석계 문제에 대한 두 가지 하이브리드 고차수 방법(HHO)의 수렴 분석을 가능하게 한다. 이 방법들은 자석장 및 벡터 포텐셜 수식 모두에서 에너지 오차 수렴률이 최적의 $k+1$을 달성하며, 입방체 및 삼각형 메쉬에서 이론적 수렴률을 확인하는 수치적 검증이 수행되었다.

ABSTRACT

We prove a discrete version of the first Weber inequality on three-dimensional hybrid spaces spanned by vectors of polynomials attached to the elements and faces of a polyhedral mesh. We then introduce two Hybrid High-Order methods for the approximation of the magnetostatics model, in both its (first-order) field and (second-order) vector potential formulations. These methods are applicable on general polyhedral meshes, and allow for arbitrary orders of approximation. Leveraging the previously established discrete Weber inequality, we perform a comprehensive analysis of the two methods. We finally validate them on a set of test-cases.

연구 동기 및 목표

• 일반적인 다면체 메쉬에 대해 3차원 하이브리드 다항식 공간에서 첫 번째 웨버 부등식의 이산형을 수립하는 것.
• 자석계 문제에 대한 두 가지 하이브리드 고차수 방법(HHO)을 개발하고 분석하는 것: 하나는 일阶장 수식이고, 다른 하나는 이阶벡터 포텐셜 수식이다.
• 에너지 노름과 $L^2$-유사 노름에서 두 방법 모두 최적의 $k+1$ 수렴률을 증명하는 것. 이는 임의의 다항식 차수 $k \geq 0$에 대해 유효하다.
• 입방체 및 정규 삼각형 메쉬에서 수치적으로 방법을 검증하여 이론적 수렴률을 확인하는 것.

제안 방법

• 각 요소별 다항식 벡터 함수의 분해를 이용하여 메쉬 요소와 면에서 생성되는 다항식 벡터로 이루어진 하이브리드 공간에서 이산 웨버 부등식을 유도한다 (보조정리 2).
• 일阶장 수식과 자석계 문제의 일반화된 벡터 포텐셜 수식을 위한 두 개의 HHO 방법을 구성한다.
• 국소 재구성과 하이브리드화를 통해 슈어 여유법을 이용하여 요소 기반 자유도를 제거함으로써 압축된 전역 시스템의 효율적 해법을 가능하게 한다.
• 최근에 증명된 이산 웨버 부등식과 이산 인프-서프 조건을 사용하여 방법의 안정성과 수렴성을 확립한다.
• 라그랑주 승수를 사용한 혼합 유한요소 유사 프레임워크를 도입하여 이산 설정에서 일致성과 호환성을 확보한다.
• 입방체 메쉬에서 입방체 및 삼각형 요소를 사용하여 수치 실험을 수행하고, 압축된 시스템을 SparseLU로 풀며 에너지, $L^2$, $L^2$-유사 노름에서 수렴성을 측정한다.

실험 결과

연구 질문

• RQ1 일반적인 다면체 메쉬에 대해 3차원 하이브리드 다항식 공간에서 연속 사례와 유사한 이산 웨버 부등식을 수립할 수 있는가?
• RQ2 제안된 HHO 방법이 자석계 문제에서 일阶장 수식과 이阶벡터 포텐셜 수식 모두에서 에너지 노름에서 최적의 수렴률 $k+1$을 달성하는가?
• RQ3 이산 이차형식 $d_h$를 제거하면 삼각형 메쉬에서 HHO 방법의 정확도와 조건수에 어떤 영향을 미치는가?
• RQ4 이산 정확한 복합체나 전역 일치성에 의존하지 않고도 이산 웨버 부등식을 사용하여 안정성과 수렴성을 증명할 수 있는가?
• RQ5 다양한 메쉬 유형과 다항식 차수에서 HHO 방법의 실용적 성능은 수렴률과 해법 시간 측면에서 어떻게 나타나는가?

주요 결과

• 다면체 메쉬에서 하이브리드 다항식 공간에 대한 이산 웨버 부등식이 엄밀히 증명되었으며, 이는 수렴 분석의 기초가 되었다.
• 일阶장 수식을 위한 HHO 방법은 $k \geq 0$에 대해 최적의 에너지 오차 수렴률 $k+1$을 달성한다.
• 벡터 포텐셜 수식을 위한 HHO 방법 역시 최적의 $k+1$ 에너지 오차 수렴률을 달성하며, 자석장에 대한 $L^2$ 오차 수렴률은 $k+2$이다.
• 일치하는 삼각형 메쉬에서 이차형식 $d_h$를 제거하면 자석장의 근사가 略적으로 더 정확해지며, $k+1$ 수렴률은 유지된다.
• 수치 결과는 입방체 및 삼각형 메쉬에서 예측된 수렴률을 확인하였으며, 자유도의 수에 따라 해법 시간이 적절히 스케일링됨을 보였다.
• 이 방법은 다항식 차수 $k = 0, 1, 2$ 전반에서 강건하며, 오차 측정치는 잘 조절되어 있고 이론적 기대와 일致한다.