[논문 리뷰] Gradient recovery for elliptic interface problem: I. body-fitted mesh

Hailong Guo, Xu Yang|arXiv (Cornell University)|2016. 07. 20.

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

한 줄 요약

이 논문은 체류형 메esh에서 타원형 인터페이스 문제를 위한 새로운 임마ージ드 기울기 복구 방법을 제안하며, 인터페이스에서 해의 불규칙성으로 인해 기존 복구 기법이 실패하는 문제를 해결한다. 인터페이스에 의해 분할된 각 부분영역에서 국소적으로 다항식 보존 복구(PPR)를 적용함으로써, 초수렴을 달성하고 이론적 분석 및 경미하게 비정규화된 메쉬와 적응형 메쉬에서의 수치 예제를 통해 검증된 점차 정확한 후행 오차 추정기를 제공한다.

ABSTRACT

In this paper, we propose a novel gradient recovery method for elliptic interface problem using body-fitted mesh in two dimension. Due to the lack of regularity of solution at interface, standard gradient recovery methods fail to give superconvergent results, and thus will lead to overrefinement when served as a posteriori error estimator. This drawback is overcome by designing an immersed gradient recovery operator in our method. We prove the superconvergence of this method for both mildly unstructured mesh and adaptive mesh, and present several numerical examples to verify the superconvergence and its robustness as a posteriori error estimator.

연구 동기 및 목표

• 해가 인터페이스에서 국소적으로 부드럽지 않을 경우 기존 기울기 복구 기법이 초수렴을 달성하지 못하는 문제를 해결하기 위해.
• 해가 전역적으로 매끄럽지 않더라도 여전히 초수렴을 유지하는 기울기 복구 기법을 개발하기 위해.
• 적응형 유한요소법을 위한 강인하고 점차 정확한 후행 오차 추정기를 구축하기 위해.
• 경미하게 비정규화된 메쉬 및 적응적으로 메쉬가 세분화된 체류형 메쉬로 초수렴 이론을 확장하기 위해.

제안 방법

• 인터페이스로 정의된 각 부분영역에서 별도로 다항식 보존 복구(PPR)를 적용하는 임마ージ드 기울기 복구 연산자를 설계하기 위해.
• 인터페이스에서 떨어진 노드에 대해서는 표준 PPR을 사용하고, 인터페이스에 인접한 노드에 대해서는 각 부분영역 내의 샘플링 점들만을 사용하여 최소제곱법으로 2차 다항식을 피팅하기 위해.
• 각 인터페이스 노드에 대해 두 개의 별개의 기울기 근사치를 구성하여 해가 인터페이스를 넘어가면서 미분 가능하지 않음을 반영하기 위해.
• 이전 연구에서 확보된 초근접성 결과를 활용하여, 경미하게 비정규화된 메쉬 조건 하에서도 초수렴을 입증하기 위해.
• 실용적인 가정과 [41]에서 제시된 초근접성 추정치를 활용하여, 적응적으로 세분화된 메쉬에서의 초수렴을 증명하기 위해.
• 회복된 기울기를 기반으로 한 복구 유형의 후행 오차 추정기를 유도하여 점차 정확성을 보장하기 위해.

실험 결과

연구 질문

• RQ1 해가 전역적으로 규칙적이지 않을 경우 기울기 복구 기법이 타원형 인터페이스 문제에서 초수렴을 달성할 수 있는가?
• RQ2 부분영역 국소 PPR 접근법이 기존 복구 기법의 한계를 극복하여 인터페이스 문제에서 성능을 향상시키는가?
• RQ3 제안된 기울기 복구 연산자가 경미하게 비정규화된 체류형 메쉬에서 초수렴을 달성하는가?
• RQ4 적응적으로 세분화된 메쉬에서 메쉬가 초수렴을 유지하고 점차 정확한 후행 오차 추정기를 제공할 수 있는가?
• RQ5 특이성이 있는 인터페이스 문제에서 복구 기반 오차 추정기는 전통적인 추정기보다 강인성과 정확도 면에서 뛰어나게 작동하는가?

주요 결과

• 제안된 임마ージ드 기울기 복구 방법은 경미하게 비정규화된 메쉬와 적응적으로 세분화된 체류형 메쉬 모두에서 초수렴을 달성한다.
• 균일 메쉬에서 β⁻ = 1000 및 β⁻ = 10000 조건 하에서 에너지 오차의 수치적 결과는 O(N^0.95)의 초수렴 속도를 확인한다.
• 이론적으로 예측된 바와 같이 특이성이 존재하는 경우 에너지 오차에 대해 최적의 수렴 속도 O(N^{-0.5})를 달성한다.
• 수치 실험에서 후행 오차 추정기의 유효 지수는 균일 메쉬 및 적응형 메쉬에서 항상 1에 가까운 것으로 확인되어 점차 정확성을 확인한다.
• 복구 기반 오차 추정기는 과다 메쉬 세분화 없이 특이성을 정확히 포착하며, 적응형 메쉬 세분화에서 기존의 SPR 및 PPR 추정기보다 뛰어난 성능을 보인다.
• 해가 전역적으로 부드럽지 않은 경우(예: u ∈ H^{1+ε})에도 불구하고, 추정기는 점차 정확하지 않지만 강인성을 유지하며 실용적인 신뢰성을 입증한다.