[논문 리뷰] Virtual Elements for the Navier-Stokes problem on polygonal meshes

L. Beirão da Veiga, Carlo Lovadina외 1명|arXiv (Cornell University)|2017. 03. 01.

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

한 줄 요약

이 논문은 일반 다각형 메esh에서 2차원 라우아-스토크스 방정식을 위한 가상 요소 방법(VEM)의 가족을 제안한다. 이 방법은 수렴 속도가 최적이고, 비연결 다각형 메쉬, 특히 복합적인 메쉬 구성에서 강건한 성능을 보이며, 점별로 수렴 자유인 속도장과 대각 행렬 형태의 강성 행렬을 통해 실현된다. 비선형 대류 항을 다루기 위해 수렴 자유인 이산 속도 공간과 향상된 다항식 자유도를 사용한다.

ABSTRACT

A family of Virtual Element Methods for the 2D Navier-Stokes equations is proposed and analysed. The schemes provide a discrete velocity field which is point-wise divergence-free. A rigorous error analysis is developed, showing that the methods are stable and optimally convergent. Several numerical tests are presented, confirming the theoretical predictions. A comparison with some mixed finite elements is also performed.

연구 동기 및 목표

• 일반 다각형 메쉬에서 2차원 라우아-스토크스 방정식을 위한 적합하고 수렴 자유인 가상 요소 방법을 개발한다.
• 비선형 대류 항의 과제를 해결하기 위해 적절한 프로젝터와 향상된 이산 속도 공간을 도입한다.
• 가상 요소 방법 프레임워크에 맞춘 엄밀한 오차 분석을 통해 안정성과 최적 수렴을 확보한다.
• 복잡한 메쉬 구성, 특히 비연결 및 다면체 메쉬를 포함한 다양한 메쉬에서의 강건성과 유연성을 입증한다.
• 특히 독립된 메쉬를 사용하는 연계된 스토크스-다르시 흐름 문제에서 수치적으로 접근법을 검증한다.

제안 방법

• 속도에 대해 수렴 자유인 이산 해를 보장하는 자유도를 갖는 가상 요소 공간을 사용한다.
• 대류 삼선형 형식은 안정성을 유지하기 위해 비대칭형 변형을 사용하여 근사화한다. 이는 프로젝터로 인해 비대칭성이 손실되더라도 유지된다.
• 비선형 항을 처리하기 위해 프로젝터 기반 접근법을 사용하여, 형상 함수의 점별 지식 없이도 강성 행렬을 계산할 수 있다.
• 압력은 요소 간 연속성 요구 없이 국소 다항식으로 근사한다.
• 압력-속도 결합 행렬에 대각형 행렬 구조를 활용하여 정적 축소 없이도 효율적인 해법 전략을 가능하게 한다.
• 혼합 유한 요소 방법을 사용한 스토크스 문제와의 원활한 결합이 가능하며, 독립된 메쉬를 사용하는 다르시-스토크스 경계 문제에 적용 가능하다.

실험 결과

연구 질문

• RQ1 가상 요소 방법은 다각형 메쉬에서 비선형 라우아-스토크스 방정식으로 확장될 수 있는가? 이때 수렴 자유 속도 근사가 유지되는가?
• RQ2 삼선형 형식의 선택(기본형 대비 비대칭형)은 라우아-스토크스 방정식의 VEM에서 안정성과 수렴에 어떤 영향을 미치는가?
• RQ3 수렴 자유 성질은 낮은 점성도 또는 연계된 흐름 문제에서 강건성을 얼마나 향상시키는가?
• RQ4 헤어진 노드와 서로 다른 하위영역에서 독립적인 분할을 가진 비연결 다각형 메쉬를 VEM 프레임워크가 다룰 수 있는가?
• RQ5 복잡한 기하 구조와 경계 문제에서 이 방법은 최적 수렴 속도와 정확도를 유지하는가?

주요 결과

• 제안된 VEM은 구조적 및 비구조적 다각형 메쉬에서의 수치 실험을 통해 속도 및 압력 필드에서 최적 수렴 속도를 달성한다.
• 점별 수렴 자유 성질로 인해 이산 압력과 속도 오차가 거의 분리되어, 낮은 점성도 영역에서 강건성이 향상된다.
• 독립적인 메쉬 분할(예: 다르시 영역에서는 삼각형, 스토크스 영역에서는 정사각형)이 있는 메쉬에서도 최적 수렴을 유지하여, 다물리스 커플링에서의 강건성을 입증한다.
• 경계를 따라 계산된 유량은 기준 해에 대해 2차 수렴 정확도를 보이며, 조밀한 메쉬에서도 상대 오차가 1% 미만으로 유지된다.
• 압력-속도 강성 행렬의 대각형 행렬 구조 덕분에 내부 자유도를 정적 축소 없이 제거할 수 있어 효율적인 해법 전략이 가능하다.
• 특히 수직류 및 다르시-스토크스 문제에서 표준 혼합 유한 요소 방법보다 속도 정확도와 강건성이 뛰어나 성능이 뛰어나다.