[논문 리뷰] A Weak Galerkin Mixed Finite Element Method for Second-Order Elliptic Problems

Junping Wang, Xiu Ye|arXiv (Cornell University)|2012. 02. 16.

Advanced Numerical Methods in Computational Mathematics참고 문헌 14인용 수 6

한 줄 요약

이 논문은 임의의 다각형/다면체 메시에서 불연속 다항식을 사용함으로써 제2차 타원형 문제에 대한 새로운 약한 갈레르킨 혼합 유한요소법(WG-MFEM)을 제안한다. 이 방법은 약한으로 정의된 이산 발산 연산자를 사용하여 이산 $H^1$ 및 $L^2$ 노름에서 최적의 오차 추정을 가능하게 하며, 메시 설계와 다항식 근사에서 높은 유연성을 확보하면서 국소 질량 보존을 보장한다.

ABSTRACT

A new weak Galerkin (WG) method is introduced and analyzed for the second order elliptic equation formulated as a system of two first order linear equations. This method, called WG-MFEM, is designed by using discontinuous piecewise polynomials on finite element partitions with arbitrary shape of polygons/polyhedra. The WG-MFEM is capable of providing very accurate numerical approximations for both the primary and flux variables. Allowing the use of discontinuous approximating functions on arbitrary shape of polygons/polyhedra makes the method highly flexible in practical computation. Optimal order error estimates in both discrete $H^1$ and $L^2$ norms are established for the corresponding weak Galerkin mixed finite element solutions.

연구 동기 및 목표

• 제2차 타원형 문제에 대한 탄력적인 유한요소법을 개발하여 임의의 다각형/다면체 메시와 불연속 근사를 허용한다.
• 약한 갈레르킨 프레임워크를 이산 약한 발산 연산자를 도입하여 혼합 유한요소 형식으로 확장한다.
• 이산 $H^1$ 및 $L^2$ 노름에서 유량 및 주요 변수 근사에 대해 최적의 오차 추정을 수립한다.
• 다공질 유체 흐름 및 이동 시뮬레이션과 같은 응용 분야에서 중요한 국소 질량 보존 성질을 유지한다.
• 기존 혼합 유한요소 방법의 연속성 제약을 분포 도함수 기반 약한 공식화를 통해 극복한다.

제안 방법

• 유량 $\mathbf{q}$와 스칼라 변수 $u$를 포함하는 1차 체계로 제2차 타원형 문제를 공식화하며, $\alpha\mathbf{q} + \nabla u = 0$ 및 $\nabla \cdot \mathbf{q} = f$를 만족시킨다.
• 각 요소에서 국소적으로 정의된 이산 약한 발산 연산자 $\nabla_w \cdot$를 도입하여 불연속 벡터값 함수를 다룬다.
• 일般적인 다각형/다면체 분할에서 $\mathbf{q}_h$ 및 $u_h$에 대해 불연속 다항식을 사용하여 약한 갈레르킨 혼합 유한요소 형식을 구성한다.
• 변분 형식을 유도한다: 모든 테스트 함수에 대해 $\mathbf{q}_h \in \mathcal{V}_h$, $u_h \in \mathcal{W}_h$ 를 찾기 위한 조건으로 $(\alpha \mathbf{q}_h, \mathbf{v}) - (\nabla_w \cdot \mathbf{v}, u_h) = \langle g, \mathbf{v} \cdot \mathbf{n} \rangle_{\partial\Omega}$ 및 $(\nabla_w \cdot \mathbf{q}_h, w) = (f, w)$를 만족시킨다.
• 지역 노름을 유계로 하고 오차 한계를 유도하기 위해 형태 정규 단순체 기반의 역부등식과 도메인 역부등식 추정을 사용한다.
• 이산 에너지 노름에서 안정성, 일致성 및 역부등식 추론을 결합하여 최적 수렴률을 확립한다.

실험 결과

연구 질문

• RQ1 제2차 타원형 문제에 대해 임의의 다각형/다면체 메시에서 불연속 근사를 允許하는 약한 갈레르킨 혼합 유한요소법을 구성할 수 있는가?
• RQ2 국소 질량 보존과 최적 수렴을 가능하게 하는 이산 약한 발산 연산자의 구조와 안정성은 어떠한가?
• RQ3 WG-MFEM 프레임워크에서 $H^1$ 및 $L^2$ 노름의 최적 오차 추정은 메시의 형상과 다항식 차수에 어떻게 의존하는가?
• RQ4 메시 설계의 고도의 유연성과 함께 높은 정밀도를 유지하면서 국소 질량 보존을 유지할 수 있는가?
• RQ5 약한 기울기와 약한 발산 기반 혼합 유한요소 계획의 최적 수렴을 증명하기 위해 필요한 이론적 도구는 무엇인가?

주요 결과

• WG-MFEM는 유량 $\mathbf{q}$ 및 주요 변수 $u$에 대해 이산 $H^1$ 및 $L^2$ 노름에서 최적 수렴률을 달성한다.
• 일般적인 유한요소 분할 조건, 즉 임의의 다각형 및 다면체 요소를 포함하여 최적의 오차 추정이 수립된다.
• 약한 발산과 약한 갈레르킨 공식화를 사용함으로써, 불연속 근사가 있더라도 국소 질량 보존이 보장된다.
• 이산 약한 발산 연산자 $\nabla_w \cdot$는 불연속 다항식 공간에서 잘 정의되고 안정적이며, 이는 방법의 탄력성에 기여한다.
• 도메인 역부등식이 증명된다: 형태 정규 단순체 $K$ 상의 차수 $n$ 다항식 $v$에 대해, 반지름이 $h_K$ 비례하는 구 $S \subset K$ 에 대해 $\|v\|_K^2 \leq C(\varsigma_*, n) \|v\|_S^2$ 가 성립한다.
• 각 요소 $T$ 상의 차수 $n$ 다항식 $\varphi$에 대해 $\|\nabla \varphi\|_T \leq C(n) h_T^{-1} \|\varphi\|_T$ 라는 역부등식이 성립하며, 이는 오차 분석에 필수적이다.