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

Junping Wang, Xiu Ye|arXiv (Cornell University)|2011. 04. 14.

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

한 줄 요약

이 논문은 이차 타원 문제에 대해 불연속 함수 위에서 약하게 정의된 기울기 연산자를 사용함으로써 완전히 불연속적인 유한요소 공간을 사용할 수 있도록 하는 약한 갈레르킨 (WG) 유한요소 방법을 제안한다. 이 방법은 이산 $H^1$ 및 $L^2$ 노름에서 최적의 오차 추정을 달성하며, 이중 문제의 전체 $H^2$ 정규성 조건 하에서 초수렴이 관찰된다.

ABSTRACT

In this paper, authors shall introduce a finite element method by using a weakly defined gradient operator over discontinuous functions with heterogeneous properties. The use of weak gradients and their approximations results in a new concept called {\em discrete weak gradients} which is expected to play important roles in numerical methods for partial differential equations. This article intends to provide a general framework for operating differential operators on functions with heterogeneous properties. As a demonstrative example, the discrete weak gradient operator is employed as a building block to approximate the solution of a model second order elliptic problem, in which the classical gradient operator is replaced by the discrete weak gradient. The resulting numerical approximation is called a weak Galerkin (WG) finite element solution. It can be seen that the weak Galerkin method allows the use of totally discontinuous functions in the finite element procedure. For the second order elliptic problem, an optimal order error estimate in both a discrete $H^1$ and $L^2$ norms are established for the corresponding weak Galerkin finite element solutions. A superconvergence is also observed for the weak Galerkin approximation.

연구 동기 및 목표

• 이차 타원 문제에 대해 완전히 불연속적인 근사가 가능한 새로운 유한요소 방법을 개발하기 위해.
• 클래식 기울기의 대체로 디지털 약한 기울기 개념을 도입하기 위해.
• 이질적인 성질을 가진 함수에 미분 연산자를 적용하기 위한 일반적인 프레임워크를 수립하기 위해.
• 결과로 얻어진 약한 갈레르킨 유한요소 해에 대해 이산 $H^1$ 및 $L^2$ 노름에서 최적 수렴 속도를 증명하기 위해.
• 이중 문제의 전체 $H^2$ 정규성 조건 하에서 초수렴이 나타나는 조건을 규명하기 위해.

제안 방법

• 요소 내부와 경계에서 부분적 통합을 통해 정의된 약한 기울기 연산자를 사용하여, 유한요소 공간 내에서 불연속 함수를 허용한다.
• 다항식 공간 위에 국소 투영을 사용하여 디지털 약한 기울기를 구성함으로써, 불연속 함수 위에서 고전 기울기의 수치적 근사를 가능하게 한다.
• 약한 갈레르킨 형식은 고전 기울기를 디지털 약한 기울기로 대체함으로써 대칭적이고 정의된 양의 선형 시스템을 이끌어낸다.
• 각 요소가 별개의 내부 및 경계 자유도를 가지며, $v = \{v_0, v_b\}$로 표현되는 혼합 유한요소 유사한 구조를 사용한다.
• 오차 분석은 보간 추정, 쌍대성 원리 및 디지털 약한 기울기와 정확한 해의 투영 간 일致성 오차의 경계를 기반으로 한다.
• 이 프레임워크는 변수 계수와 혼합 경계 조건을 가진 일반적인 이차 타원 문제를 다룰 수 있도록 확장된다.

실험 결과

연구 질문

• RQ1 완전히 불연속적인 함수를 허용하면서도 최적 수렴을 유지할 수 있는 유한요소 방법을 구성할 수 있는가?
• RQ2 디지털 약한 기울기가 불연속 함수에 대해 고전 기울기를 안정화하고 근사하는 데 어떤 역할을 하는가?
• RQ3 이산 $H^1$ 및 $L^2$ 노름에서 약한 갈레르킨 유한요소 해에 대해 유도할 수 있는 오차 추정은 무엇인가?
• RQ4 약한 갈레르킨 방법이 초수렴을 보이는 정규성 조건은 무엇인가?
• RQ5 변수 계수와 비연속 해를 가진 문제에서 이 방법의 성능은 어떠한가?

주요 결과

• 약한 갈레르킨 유한요소 방법은 이산 $H^1$ 노름에서 최적 수렴 순서 $h^{m}$을 달성하며, $L^2$ 노름에서는 $h^{m+s}$의 수렴 순서를 보인다. 여기서 $m$은 다항식 차수이고 $s$는 이중 문제의 정규성 인덱스이다.
• 이중 문제의 전체 $H^2$ 정규성 조건이 성립하는 경우 ($s=1$)에, 이 방법은 $L^2$ 노름에서 수렴 순서 $h^{j+2}$로 초수렴을 보이며, 여기서 $j$는 각 요소에서의 다항식 차수이다.
• 이산 $H^1$ 노름에서의 오차 추정은 $C(h^m\|u\|_{m+1} + h^{1+s}\|f - Q_0f\|)$로 경계되며, 이는 부드러움 조건 하에서 최적 수렴을 보여준다.
• 이 방법은 대칭적이고 정의된 양의 계수 행렬과 변수 계수를 가진 일반적인 이차 타원 문제에 대해 안정적이고 수렴한다.
• 디지털 약한 기울기 연산자는 해가 요소 경계를 가로질러 불연속일 경우에도 고전 기울기의 일致성 있고 안정적인 근사를 가능하게 한다.
• 수치 결과는 이 방법이 최적 수렴 속도를 유지하고, 적절한 정규성 조건 하에서 초수렴을 보임을 확인한다.