[논문 리뷰] Dual Raviart-Thomas mixed finite elements

François Dubois|arXiv (Cornell University)|2010. 12. 08.

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

한 줄 요약

이 논문은 2차원 푸아송 방정식을 위한 파라미터화된 이중 라비아르트-로마노프스키 혼합 유한요소를 사용한 페트로프-갈레르킨 유한체적 방법을 제안한다. 여기서 스칼라 미지수는 조각별로 일정하고, 운동량은 최저 차수의 라비아르트-로마노프스키 요소로 근사된다. 이 방법은 두 개의 자유 매개변수를 가진 이중 기저 함수의 가족을 통해 안정성과 최적 수렴성을 보장하며, 각 변에 대해 여섯 개의 삼각형 스텐실을 사용하여 국소 보존성과 일관된 이산 기울기를 확보한다.

ABSTRACT

For an elliptic problem with two space dimensions, we propose to formulate the finite volume method with the help of Petrov-Galerkin mixed finite elementsthat are based on the building of a dual Raviart-Thomas basis.

연구 동기 및 목표

• 혼합 유한요소 공식화와 페트로프-갈레르킨 변분 원리를 융합한 2차원 푸아송 방정식을 위한 유한체적 방법을 개발한다.
• 운동량 미지수를 위한 이중 라비아르트-로마노프스키 기저를 구성하여 국소 보존성과 일관된 이산 기울기를 확보한다.
• 유한체적 스킴의 안정성과 최적 수렴성을 보장하기 위한 충분 조건을 수립한다.
• 유한체적 방법의 수학적 공식화와 혼합 유한요소 이론 간의 다리를 놓으며, 특히 박스 스킴과 인프-스터프 조건과의 관계를 다룬다.
• 전체 평면의 경우에 대해 명시적인 보간 함수를 구성할 수 있는 프레임워크를 제공하여 향후 안정성 분석을 가능하게 한다.

제안 방법

• 이 방법은 스칼라 해 $ u $ 가 $ L^2 $ 에 속하고, 운동량 $ p = \nabla u $ 가 $ H({\rm div}, \Omega) $ 에 속하는 제품 공간 $ V = L^2(Ω) \times H({\rm div}, \Omega) $ 를 사용한 혼합 변분 공식화를 채택한다.
• 이산 공간은 각 삼각형에서 일정한 $ L_{\cal T}^2(\Omega) $ 와 최저 차수의 라비아르트-로마노프스키 요소를 사용한 $ H_{\cal T}({\rm div}, \Omega) $ 로 정의된다.
• 시험 함수가 이중 라비아르트-로마노프스키 기저에서 유도된 페트로프-갈레르킨 변분 공식화를 도입하여, 스킴에 대해 두 개의 자유 매개변수를 가진 기저 함수의 가족을 도출한다.
• 이산 기울기 연산자는 각 변을 둘러싼 여섯 개의 삼각형 스텐실을 기반으로 구성되어 국소 보존성과 기울기의 일관된 재구성 가능성을 보장한다.
• 이중 기저에 관련된 충분 조건 하에서 안정성이 증명되며, 스킴이 잘 정의된 해를 보장하기 위해 필요한 인프-스터프 조건을 만족함을 보였다.
• 각 삼각형에 대해 보존 법칙을 통합하여, 라비아르트-로마노프스키 기저 함수를 통해 유도된 통량을 통해 균형을 강제한다.

실험 결과

연구 질문

• RQ1 운동량에 대해 이중 라비아르트-로마노프스키 기저를 기반으로 한 페트로프-갈레르킨 혼합 유한요소를 사용한 유한체적 방법은 어떻게 공식화할 수 있는가?
• RQ2 이러한 혼합 유한체적 스킴의 안정성과 최적 수렴성을 보장하기 위한 조건은 무엇인가?
• RQ3 삼각형 메esh에서 각 변에 대해 여섯 개의 삼각형 스텐실을 사용할 때 이산 기울기 연산자는 어떻게 일관되게 재구성할 수 있는가?
• RQ4 두 개의 자유 매개변수를 가진 이중 라비아르트-로마노프스키 기저 함수의 가족은 안정적이고 수렴성이 보장되는 스킴을 구성하는 데 어떤 역할을 하는가?
• RQ5 이 공식화는 박스 스킴과 기존의 혼합 유한요소 이론과 어떻게 관련되어 있으며, 수학적 일관성과 보존성 측면에서 어떤 의미를 갖는가?

주요 결과

• 제안된 페트로프-갈레르킨 유한체적 스킴은 라비아르트-로마노프스키 기저 함수를 사용하여 각 삼각형 요소에 대해 운동량 방정식을 통합함으로써 통량의 국소 보존성을 보장한다.
• 이중 라비아르트-로마노프스키 기저에 관련된 충분한 안정성 가정 하에서 스킴은 최적 수렴률을 달성하며, 이 기저는 두 개의 자유 계수로 매개변수화된다.
• 이산 기울기 연산자는 각 변에 대해 여섯 개의 삼각형 스텐실을 기반으로 구성되어 기울기 장의 일관되고 보존적인 재구성을 가능하게 한다.
• 이 방법은 혼합 공식화의 잘 정의된 해를 보장하기 위해 필요한 인프-스터프 조건을 만족한다. 이는 이산 공간 내에서 유일한 해가 존재함을 보장한다.
• 안정성에 대한 충분 조건이 유도되었으며, 이는 이중 라비아르트-로마노프스키 기저의 선택에 의존한다. 또한, 이 방법은 연속 변분 문제와 일관됨을 보였다.
• 이 프레임워크는 $ \Omega = \mathbb{R}^2 $ 의 경우에 대해 명시적인 보간 함수를 구성할 수 있도록 하여 향후 안정성 분석과 구현에 기반을 마련한다.