[논문 리뷰] Weak Galerkin Finite Element Methods for Second-Order Elliptic Problems on Polytopal Meshes

Lin Mu, Junping Wang외 1명|arXiv (Cornell University)|2012. 04. 16.

Numerical methods in engineering인용 수 3

한 줄 요약

이 논문은 이차 타원 문제에 대해 임의의 다면체 메esh에서 이산 약한 기울기 연산자를 비연속 조각다항식에 적용하는 약한 갈레르킨 유한요소법(WG-FEM)을 제안한다. 이 방법은 이산 $H^1$ 및 $L^2$ 노름에서 최적의 오차 추정치를 달성하여 일반적인 다면체 기반 이산화에 대해 강건성, 신뢰성 및 정확성을 입증한다.

ABSTRACT

This paper introduces a new weak Galerkin (WG) finite element method for second order elliptic equations on polytopal meshes. This method, called WG-FEM, is designed by using a discrete weak gradient operator applied to discontinuous piecewise polynomials on finite element partitions of arbitrary polytopes with certain shape regularity. The paper explains how the numerical schemes are designed and why they provide reliable numerical approximations for the underlying partial differential equations. In particular, optimal order error estimates are established for the corresponding WG-FEM approximations in both a discrete $H^1$ norm and the standard $L^2$ norm. Numerical results are presented to demonstrate the robustness, reliability, and accuracy of the WG-FEM. All the results are derived for finite element partitions with polytopes. Allowing the use of discontinuous approximating functions on arbitrary polytopal elements is a highly demanded feature for numerical algorithms in scientific computing.

연구 동기 및 목표

• 이차 타원 문제에 대해 임의의 다면체 요소에서 비연속 근사를 지원하는 유한요소법을 개발하는 것.
• 형상 정규성 이외의 형상 제약 조건 없이 일반적인 다면체 메쉬에서 수치적 신뢰성과 정확성을 보장하는 것.
• 제안된 방법에 대해 이산 $H^1$ 및 $L^2$ 노름에서 최적 수렴 속도를 확립하는 것.
• 임의의 다면체 요소와 비연속 형상 함수를 허용함으로써 과학 계산의 유연성을 향상시키는 프레임워크를 제공하는 것.

제안 방법

• 이 방법은 다면체 요소에서 비연속 조각다항식 근사를 다루기 위해 이산 약한 기울기 연산자를 사용한다.
• 비연속 및 비일관 형상 함수의 사용을 가능하게 하는 약한 도함수 기반 변분 형식을 사용한다.
• 지역 요소별 투영과 안정화 항을 사용하여 일致성과 안정성을 확보하는 이산 약한 기울기를 구성한다.
• 요소 경계를 넘는 약한 연속성을 강제하는 혼합 변분 형식으로 방법을 설정한다.
• 요소 기하학에 대한 최소한의 가정을 유지하면서 최적 수렴 속도를 유지하도록 설계되었으며, 형상 정규성만 요구한다.
• 이론적 오차 추정치를 검증하고 강건한 성능을 입증하기 위해 수치 실험을 실시한다.

실험 결과

연구 질문

• RQ1 약한 갈레르킨 유한요소법은 비연속 근사를 포함한 임의의 다면체 메쉬에서 최적 수렴을 달성할 수 있는가?
• RQ2 형상 일치 제약 조건이 없는 상황에서 이산 약한 기울기 연산자는 어떻게 안정성과 정확성을 보장하는가?
• RQ3 이산 $H^1$ 및 $L^2$ 노름에서 WG-FEM에 대해 엄밀하게 확립할 수 있는 오차 추정치는 무엇인가?
• RQ4 이 방법은 다양한 다면체 메쉬 유형과 다항식 차수에서 수치적으로 어떻게 성능을 발휘하는가?

주요 결과

• WG-FEM는 다면체 메쉬에서 이차 타원 문제에 대해 이산 $H^1$ 노름에서 최적 순서의 오차 추정치를 달성한다.
• 표준 $L^2$ 노름에서도 최적 수렴 속도가 확립되어 방법의 강건성과 정확성을 확인한다.
• 임의의 다면체 요소에서 비연속 다항식 근사를 사용하더라도 방법은 안정적이고 정확하다.
• 수치 실험을 통해 이론적 오차 추정치가 확인되었으며, 다양한 메쉬 구성에서 방법의 신뢰성이 입증되었다.