[논문 리뷰] Virtual Element Methods for general second order elliptic problems on polygonal meshes

L. Beirão da Veiga, F. Brezzi외 2명|arXiv (Cornell University)|2014. 12. 08.

Advanced Numerical Methods in Computational Mathematics인용 수 7

한 줄 요약

이 논문은 다각형 메esh에서 변수 계수를 가진 일반적인 제2차 타원형 문제를 위한 새로운 가상 요소 방법(Virtual Element Method, VEM)을 제안한다. 정확한 $L^2$-투영 연산자를 활용하여 최적 수렴성을 유지한다. 기존 VEM들과 달리 라플라스 연산자에 대해 고전적 방법으로 축소되지 않으며, 다양한 다각형 메esh, 특히 극도로 기울인 요소를 포함하여 $L^2$ 및 $H^1$ 노름에서 강건하고 최적의 $h$-수렴성을 보인다.

ABSTRACT

We consider the discretization of a boundary value problem for a general linear second-order elliptic operator with smooth coefficients using the Virtual Element approach. As in [59] the problem is supposed to have a unique solution, but the associated bilinear form is not supposed to be coercive. Contrary to what was previously done for Virtual Element Methods (as for instance in [9]), we use here, in a systematic way, the L^2-projection operators as designed in [1]. In particular, the present method does not reduce to the original Virtual Element Method of [9] for simpler problems as the classical Laplace operator (apart from the lowest order cases). Numerical experiments show the accuracy and the robustness of the method, and they show as well that a simple-minded extension of the method in [9] to the case of variable coefficients produces, in general, sub-optimal results.

연구 동기 및 목표

• 기존 접근 방식이 최적 수렴성을 상실할 수 있는 변수 계수를 가진 일반적인 제2차 타원형 문제로 VEM을 확장하는 데 도전한다.
• 이차형식이 강제적이지 않은 경우에도 최적 수렴 속도를 유지하는 체계적인 프레임워크를 개발한다.
• 임의의 다각형 메esh, 특히 극도로 기울이고 볼록하지 않은 요소를 포함하여 강건성과 정확성을 확보한다.
• 기본 VEM을 변수 계수로 직접 확장하는 것이 고차수에서 최적 수렴성을 유지하지 못함을 보여주며, 특히 고차수에서 문제가 있음을 밝힌다.
• VEM 해를 다항식 공간에 $L^2$-투영하여 내부 점의 해 값 계산을 신뢰할 수 있는 방법으로 제공한다.

제안 방법

• 부드럽고 변수 계수를 가진 일반 선형 제2차 타원형 연산자의 변분 형식을 수립한다.
• 이차형식 계산의 핵심 구성 요소로 다항식 공간에 대한 $L^2$-투영 연산자를 사용하는 안정화된 VEM 형식을 도입한다.
• 정확한 다항식 차수의 $L^2$-투영과 안정화 항을 조합하여 이산 이차형식 $B_h(\cdot, \cdot)$를 정의한다.
• VEM 해의 $L^2$-투영 $\Pi^0_k p_h$를 사용하여 내부 점의 점별 오차, 특히 정확한 해의 최대값을 평가한다.
• 보간 및 투영 오차 한계를 사용하여 사전 오차 추정을 수립하여 $\|p - p_h\|_0 \leq C h^{k+1}(\|p\|_{k+1} + |f|_k)$를 도출한다.
• L^2-투영 프레임워크를 통해 일致성 및 안정성 항을 제어함으로써 방법의 일致성과 안정성을 확보한다.

실험 결과

연구 질문

• RQ1 변수 계수를 가진 일반적인 제2차 타원형 문제로 VEM을 체계적으로 확장할 수 있으며, 최적 수렴성을 유지할 수 있는가?
• RQ2 기존의 상수 계수 전용 VEM에 대해 $L^2$-투영 연산자를 도입함으로써, 변수 계수로 확장할 경우 최적 수렴성이 손상되지 않는가?
• RQ3 극도로 기울이거나 비볼록인 다각형 메쉬에서 이 방법은 표준 FEM과 비교해 어떻게 성능을 발휘하는가?
• RQ4 해가 요소 내에서 명시적으로 알려져 있지 않을 때, 내부 점에서 VEM 해의 정확도는 어떠한가?
• RQ5 기존 VEM(예: [9])을 변수 계수로 단순히 확장하면 고차수 방법에서 최적 수렴성을 유지할 수 있는가?

주요 결과

• 제안된 VEM은 수치 실험으로 확인된 lin $L^2$ 및 $H^1$ 노름에서 최적의 $h$-수렴성을 보이며, 다항식 차수 $k$에 대해 수렴 차수 $k+1$를 확보한다.
• 다음과 같은 메쉬 유형에서 $k=1$일 때 최적 수렴성을 보였다: 무작위 바르노이, Lloyd-100(CVT), 정사각형, 오목 메쉬.
• $k=4$일 때도 $L^2$ 및 $H^1$ 노름에서 최적 수렴성을 유지하지만, 단순한 확장(4.18)은 수렴 차수에 심각한 손실을 입는다.
• $\mathbb{P}_k$에 대한 VEM 해의 $L^2$-투영은 내부 점에서의 해 근사치, 특히 정확한 해의 최대값 점에서 정확한 근사를 제공한다.
• 수치 결과는 강건성을 확인한다: 극도로 기울인 비볼록 다각형 요소에서도 우수한 성능을 보이며, 메쉬 형태에 대한 독립성을 보여준다.
• 이 방법은 고차수에서 [9]의 원래 VEM 라플라스 연산자에 대해 축소되지 않으며, 더 일반적이고 독립적인 공식화임을 나타낸다.