[논문 리뷰] Virtual Element Methods for hyperbolic problems on polygonal meshes

Giuseppe Vacca|arXiv (Cornell University)|2016. 02. 18.

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

한 줄 요약

이 논문은 다각형 메esh에서 2계 초구형 문제, 특히 선형 파동 방정식을 해결하기 위한 가상 요소 방법(Virtual Element Method, VEM)을 개발한다. $H^1$- 및 $L^2$-준수 이산 이차형식을 사용하여 안정적인 준연속 및 완전 연속 스킴을 구성함으로써, 이 방법은 $H^1$-반노름과 $L^2$-노름 모두에서 최적 수렴 속도를 달성하며, 수치적 검증을 통해 정확성과 강건성을 입증한다. 특히 시간 적분에 비트 방법(Bathe method)을 뉴마크(Newmark)보다 선호할 수 있다.

ABSTRACT

In the present paper we develop the Virtual Element Method for hyperbolic problems on polygonal meshes, considering the linear wave equations as our model problem. After presenting the semi-discrete scheme, we derive the convergence estimates in H^1 semi-norm and L^2 norm. Moreover we develop a theoretical analysis on the stability for the fully discrete problem by comparing the Newmark method and the Bathe method. Finally we show the practical behaviour of the proposed method through a large array of numerical tests.

연구 동기 및 목표

• 일반적인 다각형 메쉬에서 초구형 문제, 특히 선형 파동 방정식에 대해 가상 요소 방법(Virtual Element Method, VEM)을 확장하는 것.
• 다항식 투영을 이용한 계산 가능한 이산 이차형식을 사용한 준연속 VEM 형식을 개발하는 것.
• 준연속 문제에 대해 $H^1$-반노름과 $L^2$-노름에서의 수렴성을 분석하는 것.
• 뉴마크(Newmark) 및 비트(Bathe) 시간 적분 방법을 사용한 완전 연속 스킴의 안정성과 수렴성 조사.
• 바이어노비 메쉬(Voronoi meshes)를 활용한 다각형 다항식 차수와 시간 간격을 변화시켜 광범위한 수치 실험을 통해 이론적 결과를 검증하는 것.

제안 방법

• 이 방법은 다항식 공간으로의 투영과 자유도를 통해 정의된 가상 요소 공간을 활용하며, 비다항 기저 함수의 명시적 지식 없이도 이산 이차형식을 계산할 수 있도록 한다.
• 이산 $H^1$-이차형식은 기울기 투영을 통해 구성되며, $L^2$-이차형식은 [1]의 기법을 활용하여 최적 정확도를 확보한다.
• 준연속 형식은 변분적 의미에서 파동 방정식을 약하게 강제함으로써 유도되며, 시간에 대한 상미분 방정식의 연립계로 이어진다.
• 완전 연속화를 위해 뉴마크와 비트 방법이 시간 적분기로 적용되며, 양자 모두에 대해 안정성과 수렴성 분석을 수행한다.
• 이 방법은 패치 테스트를 충족하며, 볼록이 아닌 다각형 요소에서도 최적 수렴을 유지한다.
• 수치 실험은 바이어노비 메쉬를 사용하며, $k=1$ 및 $k=2$ 다항식 차수에 대해 오차 노름 $E^{1}_{h, au}$ 및 $E^{0}_{h, au}$ 를 계산한다.

실험 결과

연구 질문

• RQ1 가상 요소 방법(Virtual Element Method)은 다각형 메쉬에서 선형 파동 방정식과 같은 초구형 문제로 성공적으로 확장될 수 있는가?
• RQ2 다각형 요소에 대해 VEM 준연속 스킴의 수렴 특성은 $H^1$-반노름과 $L^2$-노름에서 어떻게 되는가?
• RQ3 완전 연속 VEM 형식에서 뉴마크와 비트 시간 적분 스킴은 안정성과 정확성 측면에서 어떻게 비교되는가?
• RQ4 비안정화된 질량 행렬을 사용할 경우 VEM은 여전히 최적 수렴을 유지하는가?
• RQ5 특히 명시적 시간 적분에서 임의의 고주파 진동을 효과적으로 다룰 수 있는가?

주요 결과

• VEM 준연속 스킴은 $H^1$-반노름과 $L^2$-노름 모두에서 최적 수렴 속도를 달성하며, 다항식 차수 $k=1,2$ 에서 오차가 $h^k$ 비율로 감소함으로써 이론적 예측을 확인한다.
• 다항식 차수 $k=1$ 에서는 $\tau=1/40$ 에서 $H^1$-오차 $E^{1}_{h, au}$ 가 약 $3.38 \times 10^{-2}$ 수준에서 안정화되며, $k=2$ 에서는 $7.93 \times 10^{-4}$ 로 감소하여 최적 수렴을 나타낸다.
• $k=2$ 에서 $L^2$-오차 $E^{0}_{h, au}$ 는 $\tau=1/5$ 에서 $8.14 \times 10^{-3}$ 에서부터 $\tau=1/40$ 에서는 $5.77 \times 10^{-6}$ 로 감소하여 최적 수렴 순서를 확인한다.
• 비트 방법은 작은 $\mathrm{CFL}$ 값에서도 해에서 임의의 고주파 진동을 제거함으로써 뉴마크의 삼각형 법칙보다 뛰어난 성능을 보인다.
• 비안정화된 이차형식은 안정화된 버전과 매우 유사한 오차를 낼 수 있으나, 뉴마크 스킴에서 $\beta<1/4$ 일 경우 질량 행렬이 잠재적으로 특이해져 사용에 제한이 있다.
• 수치 실험 결과 VEM은 다각형 메쉬, 특히 볼록이 아닌 요소에서도 정확성과 안정성을 유지하며, 수렴 순서를 잃지 않고 임의의 다각형 요소를 지원한다.