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

Chunmei Wang, Junping Wang|arXiv (Cornell University)|2018. 06. 05.

Numerical methods in engineering참고 문헌 31인용 수 5

한 줄 요약

이 논문은 제약 최소화와 약한 갈레르킨 유한요소 설정을 결합하여 악조건의 타원형 코시 문제를 해결하기 위한 원시-이중 약한 갈레르킨 유한요소 방법을 제안한다. 이 방법은 정확한 해에 의해 만족되는 대칭적이고 잘 정의된 시스템을 도출하며, 이는 이산 $H^2$ 노름에서 $O(h)$의 최적 수렴 속도와 $L^2$, $L^1$, $L^∞$ 노름에서 $O(h^2)$의 수렴 속도를 보이며, 모델 문제에 대한 수치 실험을 통해 검증된다.

ABSTRACT

The authors propose and analyze a well-posed numerical scheme for a type of ill-posed elliptic Cauchy problem by using a constrained minimization approach combined with the weak Galerkin finite element method. The resulting Euler-Lagrange formulation yields a system of equations involving the original equation for the primal variable and its adjoint for the dual variable, and is thus an example of the primal-dual weak Galerkin finite element method. This new primal-dual weak Galerkin algorithm is consistent in the sense that the system is symmetric, well-posed, and is satisfied by the exact solution. A certain stability and error estimates were derived in discrete Sobolev norms, including one in a weak $L^2$ topology. Some numerical results are reported to illustrate and validate the theory developed in the paper.

연구 동기 및 목표

• 연속적인 데이터에 대한 의존성이 없는 악조건의 타원형 코시 문제를 위한 안정적이고 잘 정의된 수치적 스킴을 개발한다.
• 결측된 경계 조건을 가진 타원형 편미분방정식을 제약 최소화 접근법을 통해 해결하는 데 도전한다.
• 원시 및 이중 변수에 대해 약한 갈레르킨 프레임워크를 통해 대칭적이고 일致적이며 안정적인 유한요소 방법을 구축한다.
• 이산 소볼레프 노름, 특히 약한 $L^2$ 위상에서의 엄밀한 오차 추정을 유도하여 방법의 수렴성을 검증한다.
• 정확한 및 훼손된 경계 데이터를 가진 수치 실험을 통해 방법의 성능을 시연한다.

제안 방법

• 이 방법은 원시 변수(해)와 이중 변수(수반)를 연결하기 위해 제약 최소화 설정을 사용하여 대칭적 시스템을 도출한다.
• 약한 갈레르킨 유한요소 방법을 적용하여 원시 및 이중 방정식의 약한 형태를 이산화하며, 약한 기울기와 약한 함수를 사용한다.
• 제약 최소화에서 유도된 오일러-라그랑주 방정식은 원래 PDE와 그 수반을 포함하는 혼합 시스템을 형성하며, 이는 원시-이중 약한 갈레르킨 설정이 된다.
• 수치적 플럭스와 라그랑주 승수를 통해 경계 조건을 약하게 강제하여 안정성과 유연성을 향상시킨다.
• 이산 소볼레프 노름을 사용하여 오차 추정을 도출하며, 이는 약한 $L^2$ 위상에서의 추정을 포함하여 강력한 수렴 분석을 보장한다.
• 삼각형 또는 사각형 메esh에서 다항식의 조각별 근사로 수치 스킴을 구현하며, 메쉬 정밀도 조절을 통해 수렴 속도를 검증한다.

실험 결과

연구 질문

• RQ1 제약 최소화를 통해 악조건의 타원형 코시 문제를 안정화할 수 있는 원시-이중 약한 갈레르킨 유한요소 방법을 구성할 수 있는가?
• RQ2 유도된 시스템이 대칭성, 잘 정의됨, 그리고 정확한 해와의 일致성을 유지하는가?
• RQ3 이 방법의 수렴 속도는 $H^2$, $L^2$, $L^1$, $L^∞$ 노름을 포함한 다양한 이산 노름에서 최적인가?
• RQ4 노이즈가 있는 경계 자료에서 방법의 성능은 어떠하며, 자료의 훼손에 얼마나 민감한가?
• RQ5 수치 결과에 따르면 실제 적용에서 이론적 예측을 초월하는 높은 수렴 속도를 달성할 수 있는가?

주요 결과

• 원시-이중 약한 갈레르킨 유한요소 방법은 정확한 해에 의해 만족되는 대칭적이고 잘 정의된 시스템을 생성하며, 일致성을 보장한다.
• 해 $u$에 대해 이산 $H^2$ 노름에서 $O(h)$의 최적 수렴이 관찰된다.
• $L^2$, $L^1$, $L^∞$ 노름에서 근사값 $u_0$는 $O(h^2)$로 수렴하며, 모든 테스트 케이스에서 수치적 수렴 차수는 약 2에 가깝다.
• 이중 변수 $u_n$은 이산 $H^1$ 및 $W^{1,1}$ 노름에서 $O(h^2)$로 수렴하며, 이는 이론적 기대와 일치한다.
• 수치 실험 결과, 이론적 예측을 초월하여 방법이 더 잘 작동하는 것으로 나타났으며, 여러 노름에서 수렴 차수는 $O(h^2)$에 매우 가깝게 유지된다.
• 경계 자료의 소규모 무작위 훼손(예: $0.005*(0.5-Rand)$)에 매우 민감하게 반응하여, 노이즈가 있는 환경에서는 추가 정규화가 필요함을 시사한다.