[논문 리뷰] A New Primal-Dual Weak Galerkin Finite Element Method for Ill-posed Elliptic Cauchy Problems

Chunmei Wang|arXiv (Cornell University)|2018. 09. 12.

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

한 줄 요약

이 논문은 약한 기울기와 대칭적이며 잘 정의된 이산 체계를 활용하여 존재할 경우 정확한 해를 유지하는, 악조건의 타원형 캐시 문제를 해결하기 위한 새로운 원시-쌍대 약한 갈레르킨 유한요소법을 제안한다. 이 방법은 스케일링된 잔차 노름에서 최적 수렴 속도를 달성하며, 문제적일 수 있는 경우에도 안정적이고 수렴하는 수치적 근사를 보여준다.

ABSTRACT

A new numerical method is devised and analyzed for a type of ill-posed elliptic Cauchy problems by using the primal-dual weak Galerkin finite element method. This new primal-dual weak Galerkin algorithm is robust and efficient in the sense that the system arising from the scheme is symmetric, well-posed, and is satisfied by the exact solution (if it exists). An error estimate of optimal order is established for the corresponding numerical solutions in a scaled residual norm. In addition, a mathematical convergence is established in a weak $L^2$ topology for the new numerical method. Numerical results are reported to demonstrate the efficiency of the primal-dual weak Galerkin method as well as the accuracy of the numerical approximations.

연구 동기 및 목표

• 원인이 불안정하고 데이터에 대해 연속적인 의존성이 없는 악조건의 타원형 캐시 문제의 수치적 해를 해결한다.
• 기초가 되는 PDE 문제의 악조건성에도 불구하고 안정적이고 수렴하는 유한요소법을 개발한다.
• 원시-쌍대 사다리꼴 구조를 가진 약한 갈레르킨 공식화를 사용하여 대칭적이며 잘 정의된 이산 체계를 구성한다.
• 메쉬에 의존하는 스케일링된 잔차 노름에서 최적의 오차 추정치를 확립하고, 수치적 해에 대한 약한 $L^2$ 수렴을 보장한다.
• 정확한 해가 존재하는 기준 문제에 대한 수치 실험을 통해 방법의 강건성과 정확성을 입증한다.

제안 방법

• 약한 기울기와 일부 경계 조건이 있는 테스트 함수를 포함하는 약한 갈레르킨 혼합 변분 형식을 사용하여 타원형 캐시 문제를 수식화한다.
• 다각형/다면체 메쉬에서 원시 변수와 그 약한 기울기 모두에 대해 다항식 근사 근사를 사용하여 약한 형식을 이산화한다.
• 일반화된 inf-sup 조건을 통해 안정성을 확보하기 위해 이산 약한 형식과 쌍대 방정식을 결합하여 원시-쌍대 약한 갈레르킨 스킴을 구성한다.
• 쌍대 변수를 변분 형식에 통합하여 유도된 선형 체계의 대칭성과 잘 정의됨을 보장한다.
• 오차를 측정하기 위해 스케일링된 잔차 노름을 사용하며, 이는 악조건 문제의 맥락에서 수렴 분석에 필수적이다.
• 일반화된 inf-sup 조건과 새로운 수학적 도구를 적용하여 이산 스킴의 안정성과 수렴성을 증명한다.

실험 결과

연구 질문

• RQ1 악조건의 타원형 캐시 문제를 안정화하고 정확하게 해결할 수 있는 원시-쌍대 약한 갈레르킨 유한요소법을 설계할 수 있는가?
• RQ2 제안된 방법이 정확한 해가 존재할 경우 이를 만족하는 대칭적이며 잘 정의된 이산 체계를 제공하는가?
• RQ3 스케일링된 잔차 노름에서 수치적 해의 수렴 속도는 무엇이며, 최적 순서를 달성하는가?
• RQ4 해당 방법은 $L^2$ 및 $H^1$ 노름에서 어떻게 성능을 보이며, 일부 경우에서 수렴 속도가 최적을 하회하는 이유는 무엇인가?
• RQ5 원래의 캐시 문제에 해가 없거나 악조건일 경우에도 이 방법이 유일하고 수렴하는 수치적 해를 제공할 수 있는가?

주요 결과

• 제안된 원시-쌍대 약한 갈레르킨 유한요소법은 정확한 해가 존재할 경우 이를 만족하는 대칭적이며 잘 정의된 이산 체계를 생성한다.
• 스케일링된 잔차 노름에서 $\mathcal{O}(h)$ 순서의 최적 오차 추정치가 확립되어 이론적 수렴 속도를 확인한다.
• 수치 결과는 기울기의 $L^2$ 노름에서 약 $\mathcal{O}(h)$, 해의 경우 약 $\mathcal{O}(h^2)$의 수렴 속도를 보이며, 이는 이론적 기대와 일치한다.
• 스케일링된 잔차 노름 $|||e_h|||_{h,\Gamma_d}$ 에서는 여러 테스트 케이스, 즉 $u_1 = \sin(x)\sin(y)$, $u_2 = \cos(x)\cos(y)$, $u_3 = \cos(x)\sin(y)$ 에서도 $\mathcal{O}(h)$의 수렴 속도가 관찰되었으며, 메쉬가 세밀해질수록 순서가 약 1.0에 수렴한다.
• 이 방법은 강건성을 보이며, 원래 문제에 해가 없거나 악조건일 경우에도 유일한 수치적 해를 생성한다.
• 약한 $L^2$ 수렴은 수학적으로 입증되었지만, 이산 체계의 조건 수 문제나 악조건성으로 인해 $L^2$ 및 $H^1$ 노름에서는 다소 감소된 수렴 속도가 나타날 수 있다.