[논문 리뷰] Optimal error estimation for H(curl)-conforming p-interpolation in two dimensions

Alex Bespalov, Norbert Heuer|ArXiv.org|2009. 03. 25.

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

한 줄 요약

이 논문은 삼각형 및 정사각형 기준 요소 위에서 H(curl)-적합한 p-보간에 대해 최적의 오차 추정을 수립한다. 주요 기여는 임의의 r > 0에 대해 H^r(curl, K)에 속하는 벡터장에 대해 최적의 p-근사 오차 경계를 증명하는 것으로, 정규화된 Poincaré 유형의 분해와 최적의 H^1-적합 추정을 기초로 한다.

ABSTRACT

In this paper we prove an optimal error estimate for the H(curl)-conforming projection based p-interpolation operator introduced in [L. Demkowicz and I. Babuska, p interpolation error estimates for edge finite elements of variable order in two dimensions, SIAM J. Numer. Anal., 41 (2003), pp. 1195-1208]. This result is proved on the reference element (either triangle or square) K for regular vector fields in H^r(curl,K) with arbitrary r>0. The formulation of the result in the H(div)-conforming setting, which is relevant for the analysis of high-order boundary element approximations for Maxwell's equations, is provided as well.

연구 동기 및 목표

• 기존의 H(curl)-적합한 p-보간 오차 추정에서의 격차를 메우기 위해, 이는 부분적으로 최적의 결과가 아니거나 정규성이 낮은 경우(r ∈ (0,1))에 국한되어 있었다.
• Demkowicz와 Babuška가 제안한 투영 기반 H(curl)-적합한 보간 연산자에 대해 삼각형 및 사각형 기준 요소에서 최적의 p-수렴 속도를 확립하기 위해.
• 이 보간 연산자의 적용 가능성을 정규성 r > 0의 모든 경우, 즉 저정규성 및 고정규성의 벡터장에까지 확장하기 위해.
• Maxwell 방정식을 위한 고차 유한요소 및 경계요소법의 이론적 기초를 제공하기 위해, H(curl)-적합 공간에서 최적의 수렴을 보장하기 위해.

제안 방법

• 벡터장 u ∈ H^r(curl, K)를 코일-フリー 성분과 정규성이 향상된 성분으로 분해하기 위해 정규화된 Poincaré 유형의 적분 연산자를 사용한다.
• H^1-적합한 투영 기반 p-보간 연산자에 대해 H^1-준노름에서 최적의 오차 추정을 먼저 유도하며, 이는 핵심 중간 결과가 된다.
• H(curl)-적합한 보간 연산자는 커muting de Rham 다이어그램 성질을 통해 분석되며, 코일 및 발산 연산자와의 일致성을 보장한다.
• 증명은 척도 조정 기법과 기준 요소(삼각형 또는 정사각형)의 형태 정규성 가정에 기반하며, Bramble-Hilbert 보조정리와 척도 조정을 통해 일반적인 준균일 메esh로 확장된다.
• H(div)-적합한 설정은 두 차원에서 코일과 발산 연산자 간의 동형성에 기반한 회전 대칭성을 이용해 유도된다.
• 분석은 기준 삼각형과 기준 정사각형 양쪽 모두에서 수행되며, 동일한 기법이 두 기하 구조에 모두 적용된다.

실험 결과

연구 질문

• RQ1 H(curl)-적합한 p-보간에 대해 삼각형 및 사각형 요소에서 최적의 p-수렴 속도를 달성할 수 있는가?
• RQ2 기존에 알려진 r ∈ (0,1) 범위를 초월하여 임의의 정규성 r > 0에 대해 오차 추정을 확장할 수 있는가?
• RQ3 투영 기반 보간 연산자가 최적의 수렴을 달성하면서도 de Rham 다이어그램 성질을 어떻게 유지하는가?
• RQ4 고차 유한요소법의 맥락에서 H^1-적합과 H(curl)-적합 보간 오차 추정 간의 관계는 무엇인가?
• RQ5 결과를 H(div)-적합한 보간으로 확장하여 경계요소법에 응용할 수 있는가?

주요 결과

• H^1-준노름에서 H^1-적합한 p-보간에 대해 최적의 오차 추정이 증명되었으며, u ∈ H^r(K), r > 0일 때 수렴 속도는 O(p^{-r})이다.
• H(curl)-적합한 p-보간 연산자에 대해 최적의 수렴 속도 ‖u − Π_p^{curl}u‖_{H(curl)} ≤ C p^{-r} ‖u‖_{H^r(curl)} 가 임의의 r > 0에 대해 성립한다.
• 결과는 기준 삼각형과 기준 정사각형 양쪽 모두에서 동일하게 성립하며, 동일한 기초 구조와 척도에 기반해 동일한 오차 경계를 가진다.
• 증명은 정규화된 Poincaré 유형의 연산자를 사용해 u ∈ H^r(curl, K)를 코일-フリー 성분과 정규성이 향상된 성분으로 분해하는 새로운 분해 기법에 기반한다.
• H(div)-적합한 p-보간 연산자는 회전을 통해 유도되며, 그 최적의 오차 추정은 r > 0에 대해 ‖u − Π_p^{div}u‖_{H(div)} ≤ C p^{-r} ‖u‖_{H^r(div)} 로 확립된다.
• 표준 척도 조정 및 Bramble-Hilbert 추론을 통해 결과는 준균일 메쉬에서의 hp-유한요소 및 경계요소법으로 확장된다.