[논문 리뷰] Convergence and Optimality of Adaptive Methods for Poisson's Equation in the FEEC Framework

Michael Holst, Yuwen Li외 2명|arXiv (Cornell University)|2013. 06. 08.

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

한 줄 요약

이 논문은 일반적인 위상과 임의의 공간 차원을 가진 도메인에서 피타고라스 방정식에 대한 적응형 혼합 유한요소법의 수렴성과 준최적성에 대해 유한요소 외부계산(FEEC) 프레임워크를 사용하여 확립한다. 신뢰할 수 있고 효율적인 오차 추정기와 새로운 준수직성 결과를 통해 저자들은 데이터 진동을 포함한 총 오차에 대해 적응 알고리즘의 수렴성을 증명하였으며, 진동에 대한 마킹이나 내부 노드 가정을 필요로 하지 않는다.

ABSTRACT

Finite Element Exterior Calculus (FEEC) was developed by Arnold, Falk, Winther and others over the last decade to exploit the observation that mixed variational problems can be posed on a Hilbert complex, and Galerkin-type mixed methods can then be obtained by solving finite-dimensional subcomplex problems. Chen, Holst, and Xu (Math. Comp. 78 (2009) 35-53) established convergence and optimality of an adaptive mixed finite element method using Raviart-Thomas or Brezzi-Douglas-Marini elements for Poisson's equation on contractible domains in two dimensions, which can be viewed as a boundary problem on the de Rham complex. Recently Demlow and Hirani (Found. Math. Comput. 14 (2014) 1337-1371) developed fundamental tools for a posteriori analysis on the de Rham complex. In this paper, we use tools in FEEC to construct convergence and complexity results on domains with general topology and spatial dimension. In particular, we construct a reliable and efficient error estimator and a sharper quasi-orthogonality result using a novel technique. Without marking for data oscillation, our adaptive method is a contraction with respect to a total error incorporating the error estimator and data oscillation.

연구 동기 및 목표

• R² 내의 수축 가능한 도메인을 넘어서 일반적인 위상과 임의의 공간 차원을 가진 도메인으로 적응형 혼합 유한요소법의 수렴성 및 최적성 결과를 확장하는 것.
• 데이터 진동을 고려한 FEEC 프레임워크 내에서 신뢰할 수 있고 효율적인 후행 오차 추정기의 수립.
• 내부 노드 성질과 같은 인위적인 가정을 피하기 위해 새로운 기법을 사용해 더 날카운 준수직성 결과를 도출하는 것.
• 이산 오차와 데이터 진동을 조합한 총 오차 노름에 대해 적응 알고리즘이 수렴함을 증명하여 진동에 대한 추가 마킹 단계 없이 수렴성을 확보하는 것.
• 이전 문헌의 APPROX 모듈이나 추가 가정에 의존하지 않고 근사 클래스의 관점에서 준최적성을 달성하는 것.

제안 방법

• Hilbert 복합체 위에서 혼합 변분 문제를 설정하기 위해 유한요소 외부계산(FEEC) 프레임워크를 사용하여 안정적인 혼합 방법의 체계적 구성이 가능하다.
• 국소 잔여항과 유량 재구성에 기반한 신뢰할 수 있고 효율적인 후행 오차 추정기를 구성하여 총 오차의 상한과 하한을 보장한다.
• 공형학, 경계 및 수직 부분공간과 관련된 성분으로 오차를 분해함으로써 FEEC 도구를 사용해 더 날카운 준수직성 추정을 도출하는 새로운 기법을 도입한다.
• 안정성 및 근사 상수에 따라 결정되는 임계값 θ* 이하의 마킹 매개수 θ < θ* 를 사용해 Dörfler 마킹을 적용하여 충분한 균형을 확보한다.
• 오차 추정기, 준수직성 및 마킹 전략을 조합하여 에너지 노름에서 총 오차(이산 오차 + 데이터 진동)의 수렴성을 증명한다.
• 메쉬 크기의 관점에서 수렴 속도를 특징짓기 위해 근사 클래스 A_s 와 A_s^o 를 사용하며, 표준 가정 하에 준최적성을 확보한다.

실험 결과

연구 질문

• RQ1 일반적인 위상과 임의의 공간 차원을 가진 도메인에서 FEEC 프레임워크를 사용하여 피타고라스 방정식에 대한 적응형 혼합 유한요소법의 수렴성과 최적성을 확립할 수 있는가?
• RQ2 FEEC 환경에서 데이터 진동을 고려한 신뢰할 수 있고 효율적인 후행 오차 추정기를 어떻게 구성할 수 있으며, 이를 위해 명시적인 진동 마킹이 필요로 하지 않는가?
• RQ3 FEEC 도구를 사용해 내부 노드 성질이나 APPROX 모듈과 같은 가정을 피할 수 있는 더 날카운 준수직성 결과를 유도할 수 있는가?
• RQ4 이산 오차와 데이터 진동을 조합한 총 오차 노름에 대해 적응 알고리즘이 수렴하는가? 이는 진동에 대한 추가 마킹 단계 없이도 보장되는가?
• RQ5 근사 클래스의 관점에서 이 적응 방법의 수렴 속도는 어떠한가? 그리고 인위적인 가정 없이 준최적성을 증명할 수 있는가?

주요 결과

• 적응형 혼합 유한요소법(AMFEM)은 진동에 대한 마킹이나 내부 노드 가정 없이도 총 오차 노름(이산 오차 + 데이터 진동)에 대해 수렴한다.
• 신뢰할 수 있고 효율적인 후행 오차 추정기를 구성하여 총 오차를 상한과 하한으로 둘러싸며 오차 추정의 강건성을 보장한다.
• FEEC 도구를 사용해 새로운 준수직성 결과를 도출하여 APPROX 모듈이나 인위적 가정에 의존하지 않고 오차의 수렴을 가능하게 한다.
• 마킹 매개수 θ 가 문제 상수에 따라 결정되는 임계값 θ* 이하로 선택될 경우, 총 오차에 대해 수렴성이 확보된다.
• 준최적성이 증명된다: 정확한 해와 데이터가 적절한 근사 클래스 A_s 와 A_s^o 에 속해 있을 경우 오차는 자유도 N 에 대해 최적 속도 O(N^{-s}) 로 감소한다.
• 수렴성 및 최적성 결과는 일반적인 공간 차원과 비트리비얼 위상을 가진 도메인에서도 성립하며, 이는 이전 결과가 R² 내의 수축 가능한 도메인에 국한되어 있던 것을 확장한다.