[논문 리뷰] Superconvergence of high order finite difference schemes based on variational formulation for elliptic equations

Hao Li, Xiangxiong Zhang|arXiv (Cornell University)|2019. 04. 02.

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

한 줄 요약

이 논문은 $C^0$-$Q^k$ 유한요소법에서 유도된 고차 유한차법에 $(k+1)\times(k+1)$ 가우스-로바토 quadrature를 적용하여 초수렴(superconvergence)을 확립한다. 변분형식 내의 적분을 이 quadrature로 대체함으로써, 딜리클레 경계조건을 가진 타원형 방정식에 대해 이산 $L^2$-노름에서 $(k+2)$-차 정확도를 달성하며, 이는 가우스-로바토 점에서의 함수 값이 초최적의 속도로 수렴한다는 것을 보여준다.

ABSTRACT

The classical continuous finite element method with Lagrangian $Q^k$ basis reduces to a finite difference scheme when all the integrals are replaced by the $(k+1) imes (k+1)$ Gauss-Lobatto quadrature. We prove that this finite difference scheme is $(k+2)$-th order accurate in the discrete 2-norm for an elliptic equation with Dirichlet boundary conditions, which is a superconvergence result of function values.

연구 동기 및 목표

• 고차 유한차법이 $C^0$-$Q^k$ 유한요소법에서 유도될 때, 가우스-로바토 점에서의 함수 값의 초수렴을 확립하는 것.
• 변분형식에 $(k+1)\times(k+1)$ 가우스-로바토 quadrature를 사용하면, 딜리클레 경계조건을 가진 타원형 문제에서 이산 $L^2$-노름에서 $(k+2)$-차 정확도를 갖는 유한차법이 도출된다는 것을 보여주는 것.
• 초수렴을 유지하면서 정확도나 대칭성을 희생시키지 않는 실용적이고 효율적인 유한차 구현을 제공하는 것.
• 변수 계수 문제, 대류 항, 그리고 3차원 문제로 분석를 확장하여 동일한 수렴 속도를 확인하는 것.

제안 방법

• 모든 적분을 $C^0$-$Q^k$ 유한요소 형식에서 $(k+1)\times(k+1)$ 가우스-로바토 quadrature로 대체함으로써 유한차법을 유도하며, 이는 $Q^k$ 기저의 자유도와 정확히 일치한다.
• 이산 문제를 $A_h(u_h, v_h) = \braket{f, v_h}_h$ ($\forall v_h \to V_0^h$)로 설정하며, $A_h$와 $\braket{f, v_h}_h$는 이차형식과 선형형식에 대해 가우스-로바토 quadrature를 사용한다.
• 해와 계수에 대한 매끄러움 조건을 가정할 경우, 유도된 방법이 이산 $L^2$-노름에서 $(k+2)$-차 정확도를 갖는다는 것을 증명한다.
• 변분 분석 및 이산 그린 함수 기법을 사용하여 가우스-로바토 점에서의 오차를 유 bounds하며, quadrature 노드의 초수렴 성질을 활용한다.
• 2차원 및 3차원에서 포아송 방정식, 변수 계수 문제, 대류-확산 방정식에 대해 딜리클레 및 노이만 조건을 적용한 수치적 실험을 통해 이론적 결과를 검증한다.
• 섹션 7.4에서 보여준 바와 같이, 텐서곱 quadrature를 통한 행렬 조립을 활용하여 효율적인 구현을 수행하고, 반복적 또는 직접적 해법기를 사용해 선형 시스템을 해결한다.

실험 결과

연구 질문

• RQ1 가우스-로바토 quadrature를 통해 $C^0$-$Q^k$ 유한요소법에서 유도된 고차 유한차법이 quadrature 점에서의 함수 값에 대해 초수렴을 달성할 수 있는가?
• RQ2 $(k+1)\times(k+1)$ 가우스-로바토 quadrature를 사용하면, 딜리클레 경계조건을 가진 타원형 문제에서 이산 $L^2$-노름에서 $(k+2)$-차 정확도가 유지되는가?
• RQ3 변수 계수 문제, 대류 항, 그리고 3차원 문제에서 이 방법은 어떻게 작동하며, 초수렴 속도는 유지되는가?
• RQ4 유도된 방법이 고차 정확도와 대칭성을 유지하면서 표준 유한차 스텐실로 효율적으로 구현될 수 있는가?
• RQ5 가우스-로바토 점에서의 최대 노름 초수렴 행동은 어떠한가? 그리고 이는 $L^2$-노름 결과와 어떻게 비교되는가?

주요 결과

• $(k+1)\times(k+1)$ 가우스-로바토 quadrature를 사용하여 $C^0$-$Q^k$ FEM에서 유도된 유한차법은 딜리클레 경계조건을 가진 타원형 방정식에 대해 이산 $L^2$-노름에서 $(k+2)$-차 정확도를 달성한다.
• 특히 $k=2$일 경우, 2차원 및 3차원 포아송 문제에 대한 수치 실험을 통해 $L^2$ 및 $L^\frown$ 노름 모두에서 네 번째 차 정확도를 확보한다.
• 이 방법은 강성 행렬의 대칭성을 유지하며, 효율적인 행렬 조립을 가능하게 하여 고차 시뮬레이션에 실용적이다.
• 수치 결과는 $k=2$일 때 이산 $L^2$-노름에서 약 3.95–3.96, $L^\frown$-노름에서 3.92–3.94의 수렴 차수를 보이며, 이는 이론적 $(k+2)$-차 초수렴을 확인한다.
• 초수렴 결과는 3차원 문제 및 변수 계수 문제로까지 확장되며, $k=2$일 때 이산 $L^2$-노름에서 약 4.00–4.04, $L^\frown$-노름에서 정확히 4.00의 수렴 차수를 관측한다.
• 비압축성 속도장을 가진 대류-확산 문제에서도 이 방법은 $L^2$ 및 $L^\frown$ 노름 모두에서 $O(h^{k+2})$ 수렴을 달성하며 효과적이다.