[논문 리뷰] A C^0-Weak Galerkin Finite Element Method for the Biharmonic Equation

Lin Mu, Junping Wang외 2명|arXiv (Cornell University)|2012. 12. 02.

Advanced Numerical Methods in Computational Mathematics참고 문헌 22인용 수 34

한 줄 요약

이 논문은 2차원 및 3차원에서 빈도방정식(biharmonic equation)을 위한 C⁰-약한 갈레르킨 유한요소법을 제안하며, C⁰ 연속 유한요소 함수를 위한 새로운 약한 라플라시안 공식을 사용한다. 이 방법은 대칭적이며, 양의 정부호이고, 매개변수 없이 작동하며, 이산 H² 및 L² 노름에서 최적 수렴률을 달성한다. 오차 추정치는 체적 및 표면 질량을 차수 (k+1−d) 및 (k+2−d)까지 유지하는 보완된 Scott-Zhang 보간 연산자를 통해 검증된다.

ABSTRACT

A C^0-weak Galerkin (WG) method is introduced and analyzed for solving the biharmonic equation in 2D and 3D. A weak Laplacian is defined for C^0 functions in the new weak formulation. This WG finite element formulation is symmetric, positive definite and parameter free. Optimal order error estimates are established in both a discrete H^2 norm and the L^2 norm, for the weak Galerkin finite element solution. Numerical results are presented to confirm the theory. As a technical tool, a refined Scott-Zhang interpolation operator is constructed to assist the corresponding error estimate. This refined interpolation preserves the volume mass of order (k+1-d) and the surface mass of order (k+2-d) for the P_{k+2} finite element functions in d-dimensional space.

연구 동기 및 목표

• C¹ 연속 유한요소가 필요하지 않은 빈도방정식을 위한 C⁰-약한 갈레르킨 유한요소법을 개발하기 위해.
• C⁰ 유한요소 함수를 위한 새로운 약한 라플라시안 연산자를 정의하여 대칭적이고 양의 정부호인 공식화를 가능하게 하기 위해.
• 결과로 얻어진 약한 갈레르킨 방법에 대해 이산 H² 및 L² 노름에서 최적의 오차 추정치를 확립하기 위해.
• d차원 공간에서 체적 질량을 차수 (k+1−d)까지, 표면 질량을 차수 (k+2−d)까지 유지하는 보완된 Scott-Zhang 보간 연산자를 구성하기 위해.
• C⁰ 유한요소 공간에 대해 최적의 근사 성질을 보장하는 새로운 보간 연산자를 기반으로 한 수렴 분석을 제공하기 위해.

제안 방법

• 요소 내 약한 도함수를 포함하는 변분 공식을 통해 정의된, C⁰ 유한요소 함수를 위한 새로운 약한 라플라시안 연산자 Δ_w를 도입한다.
• 새로운 약한 라플라시안을 사용하여 대칭적이며, 양의 정부호이고, 매개변수 없는 빈도방정식에 대한 약한 갈레르킨 유한요소 공식을 구성한다.
• 요소 내부 및 면에서 모멘트 조건을 강제하는 보완된 Scott-Zhang 보간 연산자를 사용하여 다항식 모멘트를 차수 (k+1−d) 및 (k+2−d)까지 유지한다.
• 별개의 요소 스타일 스텐실에서 다항식 차수 k+2를 정확히 재현함을 증명하여 보간 연산자의 최적 근사 성질을 확립한다.
• 약한 갈레르킨 공식의 안정성과 일致성과 보간 오차 한계를 조합하여 이산 H² 노름과 L² 노름에서 오차 추정치를 유도한다.
• 보완된 보간 연산자에 표준 스케일링 추론과 소볼레프 부등식을 적용하여 H¹ 및 H² 노름에서 순서 h^{k+3}의 최적 수렴률을 증명한다.

실험 결과

연구 질문

• RQ1 C¹-연속 요소의 복잡성을 피하면서도 빈도방정식에 대해 C⁰-약한 갈레르킨 유한요소법을 구성할 수 있는가?
• RQ2 제안된 C⁰ 함수를 위한 약한 라플라시안 공식이 대칭적이며, 양의 정부호이고, 매개변수 없이 작동하는가?
• RQ3 d차원 공간에서 체적 질량을 차수 (k+1−d)까지, 표면 질량을 차수 (k+2−d)까지 유지하는 보완된 Scott-Zhang 보간 연산자를 설계할 수 있는가?
• RQ4 C⁰-약한 갈레르킨 방법이 이산 H² 및 L² 노름에서 최적 수렴률을 갖는가?
• RQ5 보완된 보간 연산자가 오차 분석에 어떻게 기여하며, 최적의 근사 성질을 보장하는가?

주요 결과

• 제안된 C⁰-약한 갈레르킨 유한요소법은 대칭적이며, 양의 정부호이며, 안정화 매개변수 없이 작동한다.
• 이산 H² 노름과 L² 노름에서 최적 순서의 오차 추정치가 확립되었으며, H² 노름에서 해의 수렴률은 O(h^{k+3})이며, L² 노름에서도 O(h^{k+3})이다.
• 보완된 Scott-Zhang 보간 연산자는 체적 질량을 차수 (k+1−d)까지, 표면 질량을 차수 (k+2−d)까지 유지하며, 최소 요구 조건보다 한 단계 높은 수준이다.
• 보간 연산자는 모든 p ∈ P_{k+1−d}(T)에 대해 ∫_T (v−Q₀v)p dx = 0 이고, 모든 p ∈ P_{k+2−d}(E)에 대해 ∫_E (v−Q₀v)p dE = 0 를 만족하여 고차 모멘트 보존을 보장한다.
• 요소 스텐실에서 다항식 P_{k+2}를 국소적으로 재현함으로써 최적의 근사 성질을 달성하여 날카운 오차 한계를 가능하게 한다.
• 수치적 결과는 이론적 수렴률을 확인하며, 삼각형 및 테트라헤드론 메esh에서 제안된 방법의 강건성과 정확성을 검증한다.