[논문 리뷰] Estimate of the Convergence Rate of Finite Element Solutions to Elliptic Equations of Second Order with Discontinuous Coefficients

Jinchao Xu|arXiv (Cornell University)|2013. 11. 17.

Advanced Mathematical Modeling in Engineering참고 문헌 1인용 수 29

한 줄 요약

이 논문은 2차원에서 계수가 불연속인 제2형 타원형 문제에 선형 유한요소법을 적용할 때 최적의 오차 추정을 수립한다. 변수 기반 테일러 전개와 가중치 적분 추정에 기반한 새로운 분석 기법을 사용하여 수렴 속도 ‖u−uₕ‖₁,Ω ≤ Ch|ln h|¹ᐟ²‖u‖₂,Ω₁₊Ω₂ 를 증명함으로써, 계수의 불연속성이 부드러운 경우에 비해 수렴을 약간 떨어뜨릴 뿐임을 보여준다.

ABSTRACT

In this paper, we consider elliptic boundary value problems with discontinuous coefficients and obtain the asymptotic optimal error estimate $\|u-u_k\|_{1,Ω}\leqslant Ch|\ln h|^{1/2}\|u\|_{2,Ω_1+Ω_2}$ for triangle linear elements.

연구 동기 및 목표

• 계수가 불연속인 제2형 타원형 편미분방정식에 대한 선형 유한요소법의 수렴 행동을 분석하는 것.
• 2차원에서 계수의 불연속성이 유한요소 오차 추정에 미치는 영향을 정량화하는 것.
• 조각별로 스무스인 경계 S가 하위영역 Ω₁과 Ω₂를 분리하는 상황에서, 점차적으로 최적인 오차 추정을 유도하는 것.
• 계수의 점프로 인해 전역 H² 정규성의 결여가 발생하더라도 여전히 거의 최적인 수렴 속도를 확립하는 것.

제안 방법

• 계수 B와 1차 항 σ가 불연속인 모델 문제의 변분형식을 이원형 형식 a(u,v)를 통해 정의한다.
• 경계 S를 고려한 정규 삼각분할 위의 연속적이며 조각별 선형 함수로 이루어진 유한요소 공간 Sₕ 를 적용한다.
• 경계 S와 교차하는 비정상적인 요소를 다루기 위해 가변 기저점 기반의 테일러 전개에 기반한 새로운 오차 분석 기법을 활용한다.
• 기울기의 경계 근처 행동을 제어하기 위해 변수 변환과 가중치 적분 추정을 도입한다.
• 두 가지 핵심 보조정리를 활용한다: 하나는 비대칭 커널을 가진 특이 적분의 L² 추정을 위한 것이고, 另一个是 측도 의존성과 함께 소부분집합에서의 L² 노름을 위한 소볼레프 유형의 삽입 정리이다.
• 오차 추정을 최적화하기 위해 ε = 1/(2|ln h|) 를 선택하여 h와 |ln h|에 대한 의존도를 최소화함으로써 최종 수렴 속도를 도출한다.

실험 결과

연구 질문

• RQ1 계수의 불연속성이 2차원에서 선형 유한요소법의 수렴 속도에 어떤 영향을 미치는가?
• RQ2 정확한 해가 계수 점프로 인해 전역 H² 정규성이 결여된 경우, H¹-일치 유한요소법에 대해 최적의 오차 추정을 도출할 수 있는가?
• RQ3 경계 S가 조각별 스무스이고 계수 B가 S를 기준으로 불연속일 때, 오차가 메쉬 크기 h에 대해 어떻게 정밀하게 의존하는가?
• RQ4 계수가 불연속일 경우 기존의 유한요소 오차 추정이 얼마나 깨지거나 악화되는가?

주요 결과

• 유한요소 해 uₕ 는 H¹-노름에서 정확한 해 u 로 수렴하며, 수렴 속도는 ‖u−uₕ‖₁,Ω ≤ Ch|ln h|¹ᐟ²‖u‖₂,Ω₁₊Ω₂ 이다.
• 로그 인자 |ln h|¹ᐟ² 는 경계 S 근처의 비정상적인 요소 분석에서 유래되며, 주어진 가정 하에서 최적임이 입증된다.
• 계수 B의 불연속성에도 불구하고 오차 추정은 점차적으로 최적을 유지함을 보여주며, 부드러운 계수 경우에 비해 수렴 속도가 약간 떨어지는 것으로 나타난다.
• 분석을 통해 u가 H²(Ω)에 속하지 않더라도 H¹-오차는 H²-준노름 ‖u‖₂,Ω₁₊Ω₂ 를 상수배로 곱한 h|ln h|¹ᐟ² 로 유계임을 증명한다.
• 오차를 정상적인 요소와 비정상적인 요소로 분해하고 각각에 맞는 특화된 추정을 적용함으로써, 전역 H² 정규성의 부재를 성공적으로 다룬다.