[논문 리뷰] $L^\infty$- and $W^{1,\infty}$-error estimates of linear finite element method for Neumann boundary value problems in a smooth domain

Takahito Kashiwabara, Tomoya Kemmochi|arXiv (Cornell University)|2018. 04. 02.

Advanced Numerical Methods in Computational Mathematics참고 문헌 18인용 수 3

한 줄 요약

이 논문은 매끄러운 도메인에서 두 번째 차수 타원형 뉴먼 문제에 선형 유한요소법을 적용할 때, 다각형 근사에 의한 도메인 변형을 고려하여 점별 오차 추정을 최적화한다. 정규화된 그린 함수와 이중 원판에서의 국소적 $H^1$- 및 $L^2$-추정을 조합하여, 비일관성(경계 근사로 인한)에도 불구하고 $L^\infty$-노름에서 $O(h^2|\text{log}~h|)$, $W^{1,\infty}$-노름에서 $O(h)$ 수렴 속도를 증명한다.

ABSTRACT

Pointwise error analysis of the linear finite element approximation for $-Δu + u = f$ in $Ω$, $\partial_n u = τ$ on $\partialΩ$, where $Ω$ is a bounded smooth domain in $\mathbb R^N$, is presented. We establish $O(h^2|\log h|)$ and $O(h)$ error bounds in the $L^\infty$- and $W^{1,\infty}$-norms respectively, by adopting the technique of regularized Green's functions combined with local $H^1$- and $L^2$-estimates in dyadic annuli. Since the computational domain $Ω_h$ is only polyhedral, one has to take into account non-conformity of the approximation caused by the discrepancy $Ω_h eq Ω$. In particular, the so-called Galerkin orthogonality relation, utilized three times in the proof, does not exactly hold and involves domain perturbation terms (or boundary-skin terms), which need to be addressed carefully. A numerical example is provided to confirm the theoretical result.

연구 동기 및 목표

• 매끄러운 도메인에서 뉴먼 경계 조건을 가진 선형 유한요소법에 대한 점별 오차 추정을 제공한다.
• 다각형 도메인 근사($\Omega_h \neq \Omega$)와 경계 불일치($\Gamma_h \neq \Gamma$)로 인한 비일관성 문제를 다룬다.
• 표준 $L^\infty$ 및 $W^{1,\infty}$ 오차 분석을 도메인 변형 하에서 뉴먼 문제에까지 확장한다.
• 도메인 불일치로 인한 갈레르킨 수직성의 붕괴에도 불구하고 최적 수렴 속도를 확립한다.
• 이론적 결과를 수치 예제로 검증하여 유도된 오차 한계를 확인한다.

제안 방법

• 정규화된 그린 함수를 활용하여 $L^\infty$ 및 $W^{1,\infty}$ 오차 분석을 그린 함수 근사의 $W^{1,1}$-오차 추정으로 환원한다.
• 관심 점을 중심으로 한 이중 원판 $A_j$로 도메인을 분해하여 국소적 오차 제어를 가능하게 한다.
• 각 원판 내에서 국소적 $H^1$- 및 $L^2$-추정을 적용하여 정규화된 그린 함수 근사의 오차를 제한한다.
• $\Omega_h \neq \Omega$ 및 $\Gamma_h \neq \Gamma$로 인한 도메인 변형 항(경계 스킨 효과)을 가중 노름에서 정확히 추정하여 처리한다.
• 비일관성을 고려하기 위해 수정항을 포함한 갈레르킨 수직성 관계를 사용하여 안정성과 수렴성을 확보한다.
• 스케일링 힌트와 이중 분해를 적용하여 그린 함수의 오차 기울기의 $L^1$-노름에 대한 날카운 경계를 도출한다.

실험 결과

연구 질문

• RQ1 계산 도메인이 다각형일 때 매끄러운 도메인에서 뉴먼 문제에 대해 선형 유한요소법의 최적 점별 수렴 속도는 무엇인가?
• RQ2 경계 근사로 인한 도메인 변형은 유한요소법에서 $L^\infty$ 및 $W^{1,\infty}$ 오차 추정에 어떤 영향을 미치는가?
• RQ3 정규화된 그린 함수와 이중 원판을 효과적으로 조합하여 비일관성 하에서 날카운 $L^\infty$ 및 $W^{1,\infty}$ 오차 경계를 유도할 수 있는가?
• RQ4 $\Omega_h \neq \Omega$ 및 $\Gamma_h \neq \Gamma$일 때 경계 스킨 항은 오차 분석에서 어떤 역할을 하는가?
• RQ5 유도된 오차 추정은 조건부 경우의 표준 결과와 어떻게 비교되며, $L^\infty$-오차에서 로그 인자 $|\text{log}~h|$의 영향은 무엇인가?

주요 결과

• 논문은 뉴먼 문제의 선형 유한요소 해에 대해 $L^\infty$-노름에서 $O(h^2|\text{log}~h|)$ 오차 경계를 확립한다.
• $W^{1,\infty}$-노름에서 $O(h)$ 오차 경계가 증명되었으며, 이는 매끄러운 해에 대해 최적 속도와 일치한다.
• 도메인 근사로 인한 비일관성에도 불구하고, $\Omega_h \neq \Omega$에서 기인하는 경계 스킨 항을 철저히 다루어 수렴 속도를 달성한다.
• 특이점 근처의 점별 오차를 제어하기 위해 정규화된 그린 함수와 이중 원판 내 국소적 추정에 기반한 분석이다.
• 보조정리들은 경계 스킨층 기여가 $\delta = C_0 E h^2$로 제어되며, 기울기의 $L^1$-노름이 적절히 감소함을 확인한다.
• 수치 예제를 통해 이론적 수렴 속도가 검증되었으며, 유도된 오차 추정이 타당함을 확인한다.