[논문 리뷰] FEM-analysis on graded meshes for turning point problems exhibiting an interior layer

Simon Becher|arXiv (Cornell University)|2016. 03. 15.

Differential Equations and Numerical Methods참고 문헌 4인용 수 3

한 줄 요약

이 논문은 내부 전환점과 코초형 내부 계면을 가진 특성화된 두점 경계값 문제에 대해 층에 적합한 격자에서 유한요소해석을 제시한다. Liseikin의 격자 생성 함수를 사용한 고차 유한요소를 적용하여, 에너지 노름에서 $\varepsilon$-균일 오차 추정과 선형 요소에 대해 최적의 $O(N^{-2})$ 수렴 속도를 증명함으로써, 강한 경계층과 내부 계면을 가진 문제에 대해 수치적 안정성과 정확성이 뛰어나다는 것을 입증한다.

ABSTRACT

We consider singularly perturbed boundary value problems with a simple interior turning point whose solutions exhibit an interior layer. These problems are discretised using higher order finite elements on layer-adapted graded meshes proposed by Liseikin. We prove $ε$-uniform error estimates in the energy norm. Furthermore, for linear elements we are able to prove optimal order $ε$-uniform convergence in the $L^2$-norm on these graded meshes.

연구 동기 및 목표

• 단순한 내부 전환점과 코초형 내부 계면을 가진 특성화된 두점 경계값 문제를 해결하는 데 있어 수치적 과제를 해결하기 위해.
• 퍼티urbation 매개변수 $\varepsilon$에 관계없이 수렴이 보장되는 층에 적합한 격자 위의 유한요소법을 개발하고 분석하기 위해.
• 기존 결과를 선형 유한요소를 넘어서, Liseikin의 격자에서 고차 유한요소에 대해 $\varepsilon$-균일 수렴을 증명함으로써 확장하기 위해.
• 이러한 격자 위에서 선형 유한요소에 대해 최적의 $L^2$-노름 오차 추정을 확립하여 이론적 수렴 속도를 확인하기 위해.

제안 방법

• 연구는 Liseikin의 격자 생성 함수를 사용한 고차 라그랑주 요소를 적용한 유한요소법(FEM)을 사용하며, 이는 내부 계면을 해석하기 위해 설계되어 있다.
• 격자는 전환점 $x_0 = 0$ 근처에서 조밀하게 집중되도록 보장하는 변환 $\varphi(\xi,\varepsilon)$를 사용하여 구성되며, 계면 두께와 일치시키기 위해 그라디에이션 매개변수 $\alpha \leq \lambda$를 사용한다.
• 오차 측정과 $\varepsilon$-균일 수렴을 증명하기 위해 가중 에너지 노름 $|||\cdot|||_{\varepsilon} = \left(\varepsilon|\cdot|_1^2 + \|\cdot\|^2\right)^{1/2}$ 를 사용한다.
• 선형 요소의 경우 초밀착성 결과(superconvergence)를 도출하고, 후처리 및 보간 추정을 통해 $L^2$-노름에서 최적의 $O(N^{-2})$ 수렴을 증명한다.
• 핵심 기술 도구로는 격자 의존 역행렬 부등식, 보간 오차 추정, 격자 비율과 격자 함수의 이阶도 도함수에 대한 경계를 포함한다.
• 이론적 분석은 등분포 원칙과 주요량 기반의 계면 감쇠 변환을 활용하여 해의 도함수와 격자 그라디에이션을 제어한다.

실험 결과

연구 질문

• RQ1 Liseikin의 격자에서 고차 유한요소가 내부 전환점과 코초형 계면을 가진 특성화된 문제에 대해 $\varepsilon$-균일 수렴을 달성할 수 있는가?
• RQ2 이러한 격자에서 선형 유한요소의 $L^2$-노름에서 최적의 수렴 속도는 무엇인가?
• RQ3 격자 그라디에이션 매개변수 $\alpha$와 계면 매개변수 $\lambda$는 수렴 행동에 어떤 영향을 미치는가?
• RQ4 초밀착성 성질을 활용하여 선형 요소에 대해 최적의 $L^2$-노름 오차 추정을 유도할 수 있는가?

주요 결과

• Liseikin의 격자에서 고차 요소를 사용한 유한요소법은 요소 차수 $k$에 대해 에너지 노름에서 $O(N^{-k})$의 순서로 $\varepsilon$-균일 수렴을 달성한다.
• 선형 요소($k=1$)의 경우, $L^2$-노름에서 최적의 $\varepsilon$-균일 수렴 속도 $O(N^{-2})$를 달성한다.
• 오차 추정은 퍼티urbation 매개변수 $\varepsilon$에 독립적이며, 작은 $\varepsilon$에 대해 강건함을 확인한다.
• 분석 결과, 최적의 수렴을 유지하기 위해 격자 그라디에이션 매개변수 $\alpha$는 $\alpha \leq \lambda$를 만족해야 한다는 것이 확인된다.
• 수치 실험은 이론적 수렴 속도를 검증하여 예측된 순서와 일치하는 결과를 보였다.
• 초밀착성 결과는 이산 해가 정확한 해의 보간과 관련되어 있음을 이용하여 최적의 $L^2$-노름 오차 추정을 도출하는 데 기여한다.