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

L. Beirão da Veiga, Franco Brezzi외 2명|arXiv (Cornell University)|2015. 06. 24.

Advanced Numerical Methods in Computational Mathematics참고 문헌 48인용 수 23

한 줄 요약

이 논문은 다각형 메esh에서 변동 계수를 가진 일반적인 2계 타원 문제를 위한 혼합 가상요소법(Mixed Virtual Element Method, VEM)을 제안한다. 비상수 계수 처리를 위해 체계적인 L²-투영 연산자를 사용하며, 최적 수렴 속도를 달성하고 왜곡된 메쉬, 다각형 메쉬(바르노이 및 볼록형이 아닌 요소 포함)에서도 뛰어난 성능을 보인다.

ABSTRACT

In the present paper we introduce a Virtual Element Method (VEM) for the approximate solution of general linear second order elliptic problems in mixed form, allowing for variable coefficients. We derive a theoretical convergence analysis of the method and develop a set of numerical tests on a benchmark problem with known solution.

연구 동기 및 목표

• 기존의 삼각형 메쉬보다 더 유연한 다각형 메쉬를 사용하여 일반적인 2계 타원 문제를 해결하는 데 도전하는 것.
• 비상수 계수가 존재하는 상황에서도 최적 수렴 속도를 유지하는 가상요소법의 혼합형식을 개발하는 것.
• 비볼록형 및 바르노이 유형의 메쉬를 포함한 극도로 왜곡된 다각형 요소에서도 안정성과 정확성을 확보하는 것.
• 이론적 수렴 분석과 알려진 해가 존재하는 기준 문제에 대한 수치적 검증을 제공하는 것.
• 다양한 메쉬 유형(비볼록형, 중심점 기반 바르노이 테세레이션 포함)에 대한 테스트를 통해 방법의 강건성을 입증하는 것.

제안 방법

• 유한요소의 다각형 영역에서 정의된 이차형식을 사용하여 유량과 압력 변수를 포함하는 혼합 문제의 약한 형태를 수립한다.
• 비상수 계수 처리를 위해 다항식 공간 위에 L²-투영 연산자를 도입하여 방법을 안정화시킨다. 이는 가상요소 공간 내에서 명시적인 형상 함수가 없을 경우에 유용하다.
• 안정화 항을 내부 자유도를 제외하고 경계 자유도만을 사용하여 정의함으로써 계산 비용을 감소시킨다.
• 지역 투영과 요소 경계 및 모멘트에 관련된 자유도를 사용하여 이산 변분형식을 구성한다.
• 정확한 유량이 요소 내에서 계산될 수 없기 때문에, 수치적 유량을 $\mathbb{P}_k^2$에 투영한 결과를 오차 평가에 사용한다.
• 변동 계수 확산 연산자를 투영 기반 안정화를 통해 체계적으로 다루는 방법을 사용하여 다각형 메쉬에서 방법을 구현한다.

실험 결과

연구 질문

• RQ1 일반적인 다각형 메쉬에서 변동 계수를 가진 2계 타원 문제에 대해 혼합 가상요소법이 최적 수렴 속도를 달성할 수 있는가?
• RQ2 비볼록형 또는 비볼록형 다각형 요소, 바르노이 및 중심점 기반 바르노이 테세레이션을 포함한 극도로 왜곡된 메쉬에서 이 방법은 어떻게 성능을 보이는가?
• RQ3 공간적으로 변화하는 계수 상황에서 L²-투영 연산자가 최적 수렴을 유지하는 데 어떤 역할을 하는가?
• RQ4 이산 해와 그 L²-투영 간의 초수렴 현상이 존재하는가? 이는 총 오차 행동에 어떤 영향을 미치는가?
• RQ5 구조적, 무작위, 비볼록형 메쉬와 같은 다양한 메쉬 유형이 혼합 VEM의 수렴 속도에 어떤 영향을 미치는가?

주요 결과

• 이론적으로 예측된 바와 같이, 압력에 대해 $L^2$-노름에서 수렴 속도가 $k+1$이며, 유량에 대해서도 $L^2$-노름에서 수렴 속도가 $k+1$임을 확인하였다.
• k=1인 경우, 거친 메쉬에서는 $p_h$와 그 $L^2$-투영 $p_I$ 사이의 초수렴으로 인해 압력의 $L^2$ 오차 기울기가 $k+2=3$을 보이며, 메쉬 크기가 감소함에 따라 기울기가 $k+1=2$로 전이된다.
• 비볼록 다각형 및 무작위 바르노이 테세레이션을 포함한 극도로 왜곡된 메쉬에서도 최적 수렴을 유지함을 입증하였다.
• 안정화 항은 경계 자유도에만 의존하며, 내부 자유도는 명시적인 안정화가 필요 없어 구현을 단순화시킨다.
• 수치 결과는 모든 테스트된 메쉬 유형(Lloyd-0, Lloyd-100, 정사각형, 비볼록형)에서 일관된 수렴 행동을 보이며 방법의 강건성을 확인한다.
• k=4인 경우 오차 분해에서 초수렴 효과가 뚜렷하게 관찰되며, $\|p_I - p_h\|_0 \leq Ch^{k+2} = Ch^6$로 이론적 예측을 확인한다.