[논문 리뷰] Reliable and efficient a posteriori error estimates of DG methods for a frictional contact problem

Fei Wang, Weimin Han|arXiv (Cornell University)|2014. 01. 20.

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

한 줄 요약

이 논문은 제2종 타입 타원형 변분부등식으로 서술된 마찰 접촉 문제에 적용된 불연속 갈레르킨(DG) 방법을 위한 신뢰성 있고 효율적인 사후 오차 추정기를 제시한다. 변분부등식의 오차를 관련된 선형 문제의 오차와 연결함으로써, 저자들은 잔차 유형의 추정기를 유도하고, 이들이 신뢰성과 효율성 모두를 만족함을 증명하며, 국소 잔차, 해의 점프, 자료 진동성에 의존하는 경계를 제시한다.

ABSTRACT

A posteriori error estimators are studied for discontinuous Galerkin methods for solving a frictional contact problem, which is a representative elliptic variational inequality of the second kind. The estimators are derived by relating the error of the variational inequality to that of a linear problem. Reliability and efficiency of the estimators are shown.

연구 동기 및 목표

• 제2종 타입 타원형 변분부등식으로 제어되는 마찰 접촉 문제에 적용된 불연속 갈레르킨(DG) 방법에 대한 사후 오차 추정기를 개발하는 것.
• 제안된 오차 추정기의 신뢰성과 효율성을 확립하여, 적응형 유한요소법에서 진짜 오차를 정확히 반영함을 보장하는 것.
• 이중성과 잔차 기법을 활용하여 선형 문제의 사후 오차 분석 프레임워크를 변분부등식으로 확장하는 것.
• 적응 알고리즘에서 메쉬 정밀화를 이끌어내는 계산 가능한 오차 지표를 제공하는 것, 특히 해가 매끄럽지 않거나 미분 가능하지 않은 문제에 대해.
• DG 방법의 유연성 덕분에 메쉬의 비정규성, 특히 매달린 노드가 있는 경우에도 오차 추정기가 강인함을 보장하는 것.

제안 방법

• 변분부등식의 오차를 관련된 선형 타원형 문제의 오차와 연결함으로써 사후 오차 추정기를 유도하는 것.
• 요소 잔차, 요소 간면에서의 점프 항, 자료 진동성을 기반으로 한 잔차 유형 오차 지표를 구성하는 것.
• 잔차와 수치적 유량의 국소 투영을 이용한 복구 기법을 통해 진짜 오차를 근사하는 것.
• 코시-슈바르츠 부등식과 역추정식을 적용하여 오차를 국소 잔차와 해의 점프로 경계하는 것.
• 이중 문제 접근법과 안정성 추정식을 활용하여 오차 지표의 일致성과 신뢰성을 제어하는 것.
• 진짜 오차, 자료 진동성, 이중 해의 점프를 이용하여 국소 오차 지표를 아래에서 경계함으로써 효율성을 증명하는 것.

실험 결과

연구 질문

• RQ1 불연속 갈레르킨 방법을 마찰 접촉 문제에 적용할 때, 신뢰성 있고 효율적인 사후 오차 추정기를 구성할 수 있는가?
• RQ2 제2종 타입 변분부등식의 오차를 관련된 선형 문제의 오차와 어떻게 연결하여 오차 추정을 가능하게 할 수 있는가?
• RQ3 DG 방법을 사용할 때 국소 오차 지표에 기여하는 핵심 구성요소(잔차, 점프, 자료 진동성)는 무엇인가?
• RQ4 제안된 추정기가 적응형 메쉬 정밀화에서 진짜 오차를 어느 정도 정확히 반영하는가?
• RQ5 기본적인 유한요소 가정 하에 추정기의 효율성이 이론적으로 증명될 수 있는가?

주요 결과

• 제안된 사후 오차 추정기는 신뢰성 있다. 즉, 전역 오차가 국소 오차 지표의 합에 대한 상수배로 위에서 경계된다.
• 추정기는 효율성도 있다. 국소 오차 지표가 진짜 오차에 대한 상수배로 아래에서 경계되므로 과도한 메쉬 정밀화가 방지된다.
• 효율성 경계에는 오차의 에너지 노름, 이중 해의 점프, 자료 진동성, 국소 잔차 항이 포함된다.
• 오차 추정기는 다음 경계를 만족한다: $ \eta_K \leq C\left(|u - u_h|_{\omega_K} + \sum_{e \in \mathcal{E}(K) \cap \mathcal{E}_2} |\lambda - \lambda_h|_{*,e} + h_K \|f - f_h\|_{\omega_K} + \sum_{e \in \mathcal{E}(K) \cap \mathcal{E}_2} h_e \|\lambda_h - \overline{\lambda}_h\|_e^2 \right) $, 여기서 $ C $ 는 메쉬 크기와 무관하다.
• 분석 결과, 오차 추정기가 메쉬 정밀화에 대해 강인하며, DG 방법의 유연성 덕분에 매달린 노드가 있는 일반적인 메쉬에도 적용 가능함을 확인하였다.