[논문 리뷰] Convergence of adaptive FEM for some elliptic obstacle problem with inhomogeneous Dirichlet data

Michael Feischl, Marcus Page외 1명|arXiv (Cornell University)|2012. 07. 13.

Advanced Numerical Methods in Computational Mathematics참고 문헌 29인용 수 5

한 줄 요약

이 논문은 비균일 딜리클레 경계 조건과 일반적인 $ H^2 $ 장애물이 있는 타원형 장애물 문제에 대해 최소 순서 유한요소법의 적응형 방법의 수렴성을 증명한다. 경계 조건의 진동을 제어하기 위해 잔차 기반 오차 추정기를 확장하고, 수정된 에너지-추정기 결합을 도입함으로써, 이론적으로 이중 오차 양(에너지 오차, 추정기 오차, 경계 조건 진동 포함)의 수축 성질을 확립하여, 내재된 이산 공간이 필요 없이도 수렴성을 보장한다.

ABSTRACT

In this work, we show the convergence of adaptive lowest-order FEM (AFEM) for an elliptic obstacle problem with non-homogeneous Dirichlet data, where the obstacle χ is restricted only by χ in H^2(Ω). The adaptive loop is steered by some residual based error estimator introduced in Braess, Carstensen & Hoppe (2007) that is extended to control oscillations of the Dirichlet data, as well. In the spirit of Cascon et al. (2008), we show that a weighted sum of energy error, estimator, and Dirichlet oscillations satisfies a contraction property up to certain vanishing energy contributions. This result extends the analysis of Braess, Carstensen & Hoppe (2007) and Page & Praetorius (2010) to the case of non-homogeneous Dirichlet data as well as certain non-affine obstacles and introduces some energy estimates to overcome the lack of nestedness of the discrete spaces.

연구 동기 및 목표

• 비균일 딜리클레 경계 조건이 있는 장애물 문제에 대해 적응형 유한요소법의 수렴 분석을 확장한다.
• 비일관된 딜리클레 경계 조건 근사로 인해 이산 유한요소 공간이 중첩되지 않는 문제를 해결한다.
• 일반적인 $ H^2 $ 장애물, 비선형 경우를 포함하여 수렴성을 확립한다.
• 오차 추정기 및 분석 프레임워크에 딜리클레 경계 조건의 진동을 통합한다.
• 에너지 오차, 추정기 오차, 경계 조건 진동을 포함한 병합 오차 양의 수축 성질을 증명한다.

제안 방법

• 비균일 딜리클레 경계 조건의 진동을 제어할 수 있도록 확장된 잔차 기반 오차 추정기를 사용하여 메쉬를 적응적으로 세분화한다.
• 에너지 오차, 추정기 오차, 딜리클레 진동 항을 결합한 수정된 오차 양 $ \Delta_\ell $ 를 도입한다.
• 에너지 오차, 추정기 오차, 경계 조건 진동의 가중합을 사용하여, 소멸하는 외부 항을 고려한 수축 성질을 확립한다.
• 약한 신뢰성 추정을 적용하여 $ \varrho_\ell \to 0 $ 이면 $ \|u - U_\ell\|_{H^1} \to 0 $ 임을 보인다.
• 내부 노드 성질 대신 추정기 감소 성질을 사용하여, 메쉬 세분화 전략에 대해 강건한 결과를 얻는다.
• 이산 경계 조건 $ g_\ell $ 를 갖는 연속 해 $ u_\ell $ 를 도입하여, 이산 공간의 중첩되지 않는 문제를 회피한다.

실험 결과

연구 질문

• RQ1 비균일 딜리클레 경계 조건과 일반적인 $ H^2 $ 장애물이 있는 타원형 장애물 문제에 대해 적응형 유한요소법이 수렴할 수 있는가?
• RQ2 장애물 문제의 후행 오차 추정에서 딜리클레 경계 조건의 진동은 어떻게 제어할 수 있는가?
• RQ3 이산 공간이 중첩되지 않을 경우 병합 오차 양에 대해 수축 성질이 성립하는가?
• RQ4 수렴 증명을 애매한 장애물과 균일하지 않은 경계 조건의 경우로 확장할 수 있는가?
• RQ5 비선형 및 비균일 조건에서 적응형 유한요소법의 수렴 속도는 최적 성능을 보이는가?

주요 결과

• 비균일 딜리클레 경계 조건과 $ H^2 $ 장애물이 있는 타원형 장애물 문제에 대해, 내재된 이산 공간이 없더라도 적응형 유한요소법이 수렴한다.
• 병합 오차 양 $ \Delta_\ell $ 는 $ \Delta_{\ell+1} \leq \kappa \Delta_\ell + \alpha_\ell $ 를 만족하며, $ 0 < \kappa < 1 $ 이고 $ \alpha_\ell \to 0 $ 임을 보였다.
• 수치 실험을 통해 에너지 오차의 최적 수렴 속도 $ \mathcal{O}(N^{-1/2}) $ 와 딜리클레 진동의 $ \mathcal{O}(N^{-3/4}) $ 를 확인하였다.
• 추정기 $ \varrho\_\ell $ 는 신뢰성과 효율성 모두를 갖추었으며, $ \|u - U_\ell\|_{H^1} \simeq \varrho_\ell $ 임을 수치 결과의 평행 곡선을 통해 확인하였다.
• 적응형 메쉬 세분화는 균일 세분화보다 뛰어나며, 특히 모서리 특이점에서 최적 수렴 속도를 달성한다. 이는 적응성 매개변수 $ \theta \in \{0.4, 0.6, 0.8\} 에 관계없이 동일하게 성립한다.
• 메쉬 세분화는 해가 장애물에 닿지 않는 비활성 영역에 집중되며, 이는 예상된 결과이다.