[논문 리뷰] A posteriori error estimates for mixed virtual element methods

Andrea Cangiani, Mauricio Munar|arXiv (Cornell University)|2019. 04. 22.

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

한 줄 요약

이 논문은 혼합 경계 조건을 가진 2계 타원형 방정식에 적용된 혼합 가상원소법(mixed VEM)에 대해 잔차 기반의 잔여형 후행 오차 추정기를 제시한다. 이 추정기는 완전히 계산 가능하며 일반 다각형 메esh에서 안정적이며, VEM의 유량과 후처리된 유량에 대해 브로큰 $\mathrm{H}(\mathrm{div})$-노름에서 신뢰성 있고 효율적인 상한과 하한을 제공하여 효과적인 적응형 메쉬 구조 개선을 가능하게 한다.

ABSTRACT

We present an a posteriori error analysis for the mixed virtual element method (mixed VEM) applied to second order elliptic equations in divergence form with mixed boundary conditions. The resulting error estimator is of residual-type. It only depends on quantities directly available from the VEM solution and applies on very general polygonal meshes. The proof of the upper bound relies on a global inf-sup condition, a suitable Helmholtz decomposition, and the local approximation properties of a Clément-type interpolant. In turn, standard inverse inequalities and localization techniques based on bubble functions are the main tools yielding the lower bound. Via the inclusion of a fully local postprocessing of the mixed VEM solution, we also show that the estimator provides a reliable and efficient control on the broken $\mathrm{H}(\mathrm{div})$-norm error between the exact and the postprocessed flux. Numerical examples confirm the theoretical properties of our estimator, and show that it can be effectively used to drive an adaptive mesh refinement algorithm.

연구 동기 및 목표

• 혼합 경계 조건을 가진 2계 타원형 문제에 적용된 혼합 가상원소법(mixed VEM)에 대해 신뢰성 있고 효율적인 후행 오차 추정기를 개발하는 것.
• 정확한 해와 혼합 VEM 해 사이의 오차에 대해 브로큰 $\mathrm{H}(\mathrm{div})$-노름에서 상한과 하한을 확립하는 것, 이는 후처리된 유량을 포함한다.
• 오직 VEM 해로부터 유도된 양들에만 의존하는 완전히 局소적이고 계산 가능한 오차 추정기를 제공하여 적응형 메쉬 구조 개선에 적합하게 하는 것.
• 표준 VEM 유량에서 관찰된 부분 최적 수렴률을 초월하여, 후처리된 유량 근사에서 최적 수렴률을 증명하는 것.
• 이론적 결과를 수치 실험을 통해 검증하여 추정기가 적응형 메쉬 구조 개선을 효과적으로 이끄는지 확인하는 것.

제안 방법

• 오차 추정기는 잔차 유형이며, 오직 VEM 해로부터 계산 가능한 양들만을 사용하여 완전한 국소적 계산 가능성을 확보한다.
• 상한은 전역 inf-sup 조건, 오차의 헬름홀츠 분해, 클레멘트 유형의 내삽자 근사 성질을 이용해 유도된다.
• 하한은 역행 불등식과 버블 함수 기반 국소화 기법을 통해 확립된다.
• VEM 유량에 대해 후처리 기법을 적용하여 브로큰 $\mathrm{H}(\mathrm{div})$-노름에서 최적 수렴률을 달성하는 새로운 계산 가능한 근사값을 구성한다.
• 분석은 일반적인 다각형 메쉬, 포함하여 허니콤 노드가 있는 메쉬까지 적용 가능하며, VEM 프레임워크의 유연성을 활용한다.
• 오차 추정기를 기반으로 한 적응형 메쉬 구조 개선 알고리즘을 제안하였고, 특이한 경향이 있는 문제에 대한 수치적 검증을 수행하였다.

실험 결과

연구 질문

• RQ1 혼합 가상원소법에 대해 완전히 계산 가능한 후행 오차 추정기를 구성할 수 있는가, 이는 동시에 신뢰성 있고 효율적이어야 한다.
• RQ2 후처리된 유량에 대해 오차 추정기가 브로큰 $\mathrm{H}(\mathrm{div})$-노름에서 진짜 오차와 등가성을 가지는가?
• RQ3 제안된 추정기는 일반적인 다각형 메쉬에서 적응형 메쉬 구조 개선을 효과적으로 이끌 수 있는가?
• RQ4 후처리된 유량의 수렴 거동은 어떠한가? 브로큰 $\mathrm{H}(\mathrm{div})$-노름에서 최적 수렴률을 달성하는가?
• RQ5 특이성이 있는 문제, 예를 들어 L자형 도메인과 같은 경우에 이론적 오차 상한은 실질적으로 어떻게 작동하는가?

주요 결과

• 후행 오차 추정기는 브로큰 $\mathrm{H}(\mathrm{div})$-노름에서 오차에 대해 상한과 하한을 제공하며, 가상 불일치 항을 제외한 상태에서 성립한다.
• 후처리된 유량은 브로큰 $\mathrm{H}(\mathrm{div})$-노름에서 최적 수렴률을 달성하는 반면, 표준 VEM 유량은 부분 최적 수렴률에 머무른다.
• 추정기의 효과성 지수는 적응형 메쉬 구조 개선 단계 전반에 걸쳐 단일값에 매우 가까운 상태를 유지한다 (예: 0.87–0.95), 이는 신뢰할 수 있는 오차 제어를 의미한다.
• 수치 실험 결과, 추정기가 특이점 주변(예: L자형 도메인의 모서리)에서 메쉬 구조 개선을 효과적으로 이끌어내는 것으로 확인되었다.
• 2차 이상의 VEM($k=2$)를 사용한 적응 전략은 저차수 케이스보다 메쉬 구조 개선을 덜 강하게 이끌어내며, 이는 최적 수렴 거동을 보여준다.
• 왜곡된 다각형 메쉬에서도 추정기가 효과적으로 작동함을 확인하여 실용적인 적응형 환경에서의 안정성을 입증하였다.