[논문 리뷰] H(div) and H(curl)-conforming VEM

L. Beirão da Veiga, Franco Brezzi외 2명|arXiv (Cornell University)|2014. 07. 25.

Advanced Numerical Methods in Computational Mathematics참고 문헌 47인용 수 45

한 줄 요약

이 논문은 일반 다각형 및 다면체 요소에 대해 H(div) 및 H(curl)-형 가상 요소 방법(Virtual Element Method, VEM)을 제안하며, 기존의 Raviart-Thomas 및 BDM 유한요소를 삼각형/체적 기반 메esh에서 벗어나 일반 다각형/다면체 메쉬로 확장한다. 안정적인 다항식 보존 사영 연산자들을 구성하고, 정확한 VEM 공간의 순서를 확립함으로써, 요소별 PDE를 직접 풀지 않고도 비정규 메쉬에서 고차 혼합 형식을 구현할 수 있도록 한다.

ABSTRACT

In the present paper we construct Virtual Element Spaces that are $H({ m div})$-conforming and $H({ m \bf curl})$-conforming on general polygonal and polyhedral elements; these spaces can be interpreted as a generalization of well known Finite Elements. We moreover present the basic tools needed to make use of these spaces in the approximation of partial differential equations. Finally, we discuss the construction of exact sequences of VEM spaces.

연구 동기 및 목표

• 기존의 단순체 기반 방법의 한계를 극복하고, 일반 다각형 및 다면체 요소로 H(div) 및 H(curl)-형 유한요소 방법을 확장하는 것.
• 다항식 일致성을 유지하고 고차 근사화를 지원하는 안정적이고 계산 가능한 국소 VEM 공간을 개발하는 것.
• 요소별 PDE를 풀지 않고도 다항식 부분공간 위로의 L²-사영 연산자를 구성함으로써 효율적인 행렬 조립을 가능하게 하는 것.
• 혼합 형식에 대한 호환성을 보장하기 위해 2D 및 3D에서 VEM 공간의 정확한 순서를 확립하는 것.
• 면, 모서리 및 요소에서의 다항식 차수를 독립적으로 제어할 수 있도록 하여 특정 PDE에 맞게 조정 가능한 다양한 VEM 공간 구조를 제공하는 융통성 있는 프레임워크를 마련하는 것.

제안 방법

• 면, 모서리 및 요소의 모멘트를 기반으로 한 자유도를 사용하여 다각형 및 다면체 요소에서 H(div) 및 H(curl)-형 VEM 공간을 정의한다.
• VEM 공간에서 다항식 부분공간 위로의 L²-사영 연산자를 구성함으로써, 기저 함수의 명시적 지식 없이도 국소 강성 행렬을 계산할 수 있도록 한다.
• 최적 수렴성을 보장하기 위해 관련 다항식 부분공간(예: Raviart-Thomas 또는 BDM 유사 다항식)이 VEM 공간에 포함되도록 보장한다.
• 변분 범죄 및 일致성 논증을 통해 국소 사영의 존재성과 안정성을 증명함으로써 최적의 근사 성질을 확보한다.
• 2D 및 3D에서 이산 de Rham 복합체의 호환성을 검증함으로써 VEM 공간의 정확한 순서를 확립한다. 이는 발산 및 컬 연산자를 포함한다.
• 면에서의 다항식 차수(k_b), 발산에 대한 다항식 차수(k_d), 컬에 대한 다항식 차수(k_r)를 독립적으로 조정할 수 있도록 일반화된 프레임워크를 제안함으로써 특정 PDE에 맞는 최적화가 가능하도록 한다.

실험 결과

연구 질문

• RQ1 일반 다각형 및 다면체 요소에서 H(div) 및 H(curl)-형 가상 요소 공간을 어떻게 구성할 수 있을까? 이때 최적 수렴성을 유지할 수 있는가?
• RQ2 L²-사영 연산자가 요소별 기저 함수 평가 없이 국소 강성 행렬을 조립하는 데 어떤 역할을 하는가?
• RQ3 2D 및 3D 모두에서 VEM 공간의 정확한 순서를 구성할 수 있는가? 이는 de Rham 복합체의 구조를 유지하는 데 필수적이다.
• RQ4 이 프레임워크는 면, 모서리 및 요소에서의 다항식 차수를 독립적으로 제어할 수 있도록 일반화되어 있어 특정 응용에 맞게 최적화할 수 있는가?
• RQ5 이러한 프레임워크는 복잡한 기하구조에서 혼합 유한요소 형식에 어떤 영향을 미치는가?

주요 결과

• 제안된 H(div) 및 H(curl)-형 VEM 공간은 표준 형태 조건을 만족할 경우 일반 다각형 및 다면체 요소에서 안정적이고 최적 수렴성을 보인다.
• 다항식 부분공간 위로의 L²-사영 연산자의 구성은 요소별 PDE를 풀지 않고도 효율적이고 정확한 국소 행렬 조립을 가능하게 한다.
• 2D 및 3D 모두에서 VEM 공간의 정확한 순서가 확립되어 이산 de Rham 복합체가 유지됨을 보장하며, 이는 혼합 형식에 매우 중요하다.
• 면(k_b), 발산(k_d), 컬(k_r)에 대한 다항식 차수를 독립적으로 조정할 수 있도록 프레임워크가 허용되어 다양한 변형을 생성할 수 있다.
• 기존의 Raviart-Thomas 및 BDM 요소들이 임의의 다각형 및 다면체 메쉬로 일반화되며, 다항식 일치성 및 최적 근사성과 같은 핵심 성질을 유지한다.
• 이론적 결과는 이산 해가 혼합 형식에서 수렴을 보장하기 위해 필수적인 조건(예: inf-sup 안정성)을 만족함을 확인한다.