[논문 리뷰] Low Regularity Primal-Dual Weak Galerkin Finite Element Methods for Ill-Posed Elliptic Cauchy Problems

Chunmei Wang|arXiv (Cornell University)|2020. 11. 24.

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

한 줄 요약

이 논문은 정확한 해에 대한 최소한의 정규성 가정 하에 악조건의 타원형 캐시 문제를 위한 새로운 원시-쌍대 약한 갈레르킨 유한요소 방법(PDWG)을 소개한다. 원시 방정식과 그의 수반(쌍대) 방정식을 변분 형식을 통해 결합함으로써, 낮은 정규성 조건에서도 최적 수렴 속도를 달성하며, 이는 이산 $L^2$ 노름에서 $\mathcal{O}(h)$를 초월하는 수렴 차수를 보여주는 수치 실험을 통해 검증된다.

ABSTRACT

A new primal-dual weak Galerkin (PDWG) finite element method is introduced and analyzed for the ill-posed elliptic Cauchy problems with ultra-low regularity assumptions on the exact solution. The Euler-Lagrange formulation resulting from the PDWG scheme yields a system of equations involving both the primal equation and the adjoint (dual) equation. The optimal order error estimate for the primal variable in a low regularity assumption is established. A series of numerical experiments are illustrated to validate effectiveness of the developed theory.

연구 동기 및 목표

• 정확한 해가 연속적으로 데이터에 의존하지 않는 악조건의 타원형 캐시 문제를 다루며, 특히 정규성이 낮거나 노이즈가 있는 경계 조건에서의 문제를 해결한다.
• 정확한 해가 높은 정규성을 갖지 않더라도 효과적으로 유지되는 안정적이고 수렴하는 유한요소 방법을 개발한다.
• 약한 갈레르킨 형식을 통해 원시 방정식과 쌍대 방정식을 결합한 수학적 프레임워크를 구축하여 안정성과 정확도를 향상시킨다.
• 최소한의 정규성 가정 하에 최적의 오차 추정을 확립하여 실제 노이즈 또는 불확실한 데이터가 있는 문제에까지 적용 가능성을 넓힌다.

제안 방법

• 원시-쌍대 약한 갈레르킨(PDWG) 유한요소 방법을 사용하여 타원형 캐시 문제를 변분 형식을 통해 원래 편미분방정식과 그 수반 방정식과 결합한다.
• 약한 갈레르킨 기법을 활용하여 도메인 내에서 도함수를 약하게 정의한 불연속 유한요소 공간을 사용함으로써 정규성이 낮은 해를 다룰 수 있다.
• 원시 방정식과 쌍대 방정식을 조합하여 대칭적이고 잘 정의된 시스템을 구성함으로써, 정확한 해가 존재할 경우 일致성을 확보한다.
• Nitsche 유형의 마orting을 사용하여 경계 조건을 약하게 적용함으로써 일반적인 다각형/다면체 메esh와 경계 처리의 유연성을 확보한다.
• 원시 변수와 쌍대(수반) 변수를 포함하는 사다리꼴 시스템을 유도하기 위해 이산 변분 형식을 구성한다.
• 원시 해와 그의 쌍대 해를 동시에 풀 수 있는 혼합 유한요소 프레임워크를 사용하여 안정성과 수렴성을 향상시킨다.

실험 결과

연구 질문

• RQ1 정확한 해에 대한 최소한의 정규성 가정 하에 원시-쌍대 약한 갈레르킨 유한요소 방법이 타원형 캐시 문제에 대해 최적의 수렴 속도를 달성할 수 있는가?
• RQ2 원시 방정식과 쌍대 방정식의 결합이 노이즈가 있거나 정규성이 낮은 데이터 상황에서 안정성과 정확도를 어떻게 향상시키는가?
• RQ3 정확한 해가 제한된 소볼레프 정규성($H^1$ 유사 또는 그 이하 정규성 영역)을 가질 경우 PDWG 방법의 수렴 행동은 어떠한가?
• RQ4 계수 텐서가 일정하지 않거나 경계 데이터가 혼합형(Dirichlet 및 Neumann)일 경우에도 방법이 최적의 수렴 차수를 유지할 수 있는가?
• RQ5 해가 급격한 기울기를 보이거나 부드럽지 않을 경우, 기존의 유한요소 방법과 비교해 본 방법의 수치적 성능은 어떠한가?

주요 결과

• 정확한 해에 대한 최소한의 정규성 가정 하에 이산 $L^2$ 노름에서 원시 변수의 최적 수렴 차수를 달성하며, 수치 실험에서 $\mathcal{O}(h)$를 초월하는 수렴 속도를 보인다.
• 정확한 해 $u = xy(1-x)(1-y)$에 대해 오차 $\|e_h\|_0$는 약 $\mathcal{O}(h^{1.2})$에서 $\mathcal{O}(h^{1.6})$의 수렴 속도를 보이며, 예상되는 $\mathcal{O}(h)$ 속도를 초월한다.
• $u = e^{xy}$의 경우 $\|e_h\|_0$의 수렴 차수는 $\mathcal{O}(h^{1.1})$에서 $\mathcal{O}(h^{1.3})$ 사이로 나타나 실질적으로 초선형 수렴을 보인다.
• 쌍대 변수 오차 $\|\lambda_0\|_0$, $\|\lambda_b\|_0$, $\|\lambda_n\|_0$는 $\mathcal{O}(h^{1.8})$에서 $\mathcal{O}(h^{2.5})$ 사이의 수렴 차수를 보이며, 수반 해의 높은 정확도를 나타낸다.
• 수치 실험을 통해 계수 텐서가 일정하지 않거나 혼합 경계 조건이 존재하는 경우에도 방법이 안정적이고 수렴함을 확인하였다.
• 노이즈가 있거나 불확실한 경계 측정이 있는 경우에도 정규성이 낮은 데이터에 대해 강건함을 입증하여, 노이즈나 불확실성이 있는 역문제에의 적용 가능성을 입증한다.