[논문 리뷰] Stable broken H(curl) polynomial extensions and p-robust a posteriori error estimates by broken patchwise equilibration for the curl-curl problem

Théophile Chaumont-Frelet, Alexandre Ern외 1명|arXiv (Cornell University)|2020. 05. 28.

Advanced Numerical Methods in Computational Mathematics참고 문헌 21인용 수 4

한 줄 요약

이 논문은 삼각형 메esh의 모서리 패치 위에서 조각별 다항식 데이터를 불연속 H(curl) 공간에 대해 안정적이고 다항식 차수에 영향을 받지 않는 확장 기법을 수립하여, 새로운 불연속 국소 균형화 기법을 통해 컬-컬 문제에 대한 p-로버스트 후행 오차 추정을 가능하게 한다. 주요 기여는 다항식 차수에 의존하지 않는 상수를 가진 안정성 결과로, 전역적으로 H(curl)-일치하는 유량이 필요 없이 효율적이고 신뢰할 수 있으며 局소적으로 효율적인 오차 추정기를 제공한다.

ABSTRACT

We study extensions of piecewise polynomial data prescribed in a patch of tetrahedra sharing an edge. We show stability in the sense that the minimizers over piecewise polynomial spaces with prescribed tangential component jumps across faces and prescribed piecewise curl in elements are subordinate in the broken energy norm to the minimizers over the broken H(curl) space with the same prescriptions. Our proofs are constructive and yield constants independent of the polynomial degree. We then detail the application of this result to the a posteriori error analysis of the curl-curl problem discretized with N\\'ed\\'elec finite elements of arbitrary order. The resulting estimators are reliable, locally efficient, polynomial-degree-robust, and inexpensive. They are constructed by a broken patchwise equilibration which, in particular, does not produce a globally H(curl)-conforming flux. The equilibration is only related to edge patches and can be realized without solutions of patch problems by a sweep through tetrahedra around every mesh edge. The error estimates become guaranteed when the regularity pick-up constant is explicitly known. Numerical experiments illustrate the theoretical findings.

연구 동기 및 목표

• 공통 모서리를 공유하는 삼각형 패치 위에서 불연속 H(curl) 공간 내 조각별 다항식 데이터에 대해 안정적이고 p-로버스트 확장을 수립하기.
• 컬-컬 문제의 고차 Nédélec 유한요소 해법에 대한 p-로버스트 후행 오차 추정의 이론적 기반을 확립하기.
• 전역적으로 H(curl)-일치하는 유량이 필요 없이 효율적이고 신뢰할 수 있으며 局소적으로 효율적인 오차 추정기를 구성하기.
• 전역 패치 문제를 해결하지 않고도 실용적인 오차 추정을 가능하게 하기 위한 국소적이고 스weep 기반 알고리즘을 제공하기.

제안 방법

• 모서리 패치에서 국소 최소화를 통해 불연속 Nédélec 공간 내 국소 유량 재구성을 정의함으로써 불연속 국소 균형화를 구성하기.
• 다항식 차수 p에 의존하지 않는 상수를 가진 불연속 H(curl) 공간 내 이산 최적 근사 오차의 안정성을 증명하기 (연속 최적 근사 오차에 상대적으로).
• 지정된 컬과 면을 가로질러서 발생하는 탄성 도약을 만족시키는 이산 잠재함수를 구성하는 구축적 접근법을 사용하여 에너지 노름이 유계임을 보장하기.
• 균형화 기법을 적용하여, 정규화 회복 상수가 알려져 있을 경우 보장되는 후행 오차 추정기를 유도하기.
• 각 메쉬 모서리 주변의 단순체를 국소적으로 스캔하는 방식으로 방법을 구현하여 전역 해법을 피하기.
• 모서리 패치의 구조와 Raviart–Thomas 장을 활용하여 오차 추정기의 국소 효율성과 신뢰성을 확보하기.

실험 결과

연구 질문

• RQ1 삼각형 메쉬의 모서리 패치 위에서 불연속 H(curl) 공간 내 조각별 다항식 데이터에 대해 안정적이고 p-로버스트 확장을 구성할 수 있는가?
• RQ2 이러한 확장은 고차 Nédélec 해법에 대해 신뢰할 수 있고 국소적으로 효율적인 후행 오차 추정기를 도출하는 데 사용될 수 있는가?
• RQ3 전역적으로 H(curl)-일치하는 유량이 필요 없거나 전역 패치 문제를 해결할 필요 없이 p-로버스트 오차 추정기를 구성할 수 있는가?
• RQ4 정규화 회복 상수가 후행 분석에서 보장된 오차 범위를 확보하는 데 어떤 역할을 하는가?
• RQ5 균형화 과정은 국소 연산만을 사용하여 실용적으로 어떻게 효율적으로 구현할 수 있는가?

주요 결과

• 모서리 패치 위에서 불연속 H(curl) 공간 내 이산 최적 근사 오차는 다항식 차수 p에 관계없이 연속 최적 근사 오차로 유계이며, 이 상수는 p에 영향을 받지 않는다.
• 제안된 오차 추정기는 신뢰성 있고 국소적으로 효율적이며 p-로버스트하며, 정규화 회복 상수가 명시적으로 알려져 있을 경우 보장된 오차 범위를 가진다.
• 균형화 과정은 국소적이며 각 메쉬 모서리 주변의 단순체를 스캔하는 방식으로 구현 가능하여 전역 패치 해법을 피할 수 있다.
• 이 방법은 이산 해의 컬과 면을 가로질러 발생하는 탄성 도약을 만족시키는 유량을 구성하여 최적의 오차 제어를 보장한다.
• 수치 실험은 이론적 결과를 확인하며, 다양한 다항식 차수에서 제안된 추정기의 효과성과 안정성을 입증한다.
• 불연속 공간 내 최소화를 통한 이산 잠재함수의 구성은 안정성과 p-로버스트성을 보장하며, 핵심 부등식 (1.1)은 형상 정규성에 따라 결정되지만 다항식 차수에 무관한 상수로 성립한다.