[논문 리뷰] Weak Galerkin Finite Element Methods for Optimal Control Problems Governed by Second Order Elliptic Partial Differential Equations

Chunmei Wang, Junping Wang외 1명|arXiv (Cornell University)|2023. 10. 02.

Advanced Numerical Methods in Computational MathematicsEngineering인용 수 3

한 줄 요약

이 논문은 이차 타원 PDE에 의해 지배되는 최적 제어 문제를 위한 약한 갈레르킨(WG) 유한요소법을 개발한다. 이는 불연속 유한요소를 사용하여 유연성과 정확도를 향상시킨다. 에너지 노름과 $L^2$ 노름에서 최적의 오차 추정을 수립하였으며, 다양한 정규화 파라미터를 가진 다각형 메쉬에서의 수치적 수렴 연구를 통해 검증하였다.

ABSTRACT

This paper is concerned with the development of weak Galerkin (WG) finite element method for optimal control problems governed by second order elliptic partial differential equations (PDEs). It is advantageous to use discontinuous finite elements over the traditional $C^1$ finite elements here. Optimal order error estimates are established and confirmed by some numerical tests.

연구 동기 및 목표

• 이차 타원 PDE에 의해 지배되는 최적 제어 문제를 위한 약한 갈레르킨 유한요소법을 개발한다.
• 더 높은 유연성을 확보하기 위해 전통적인 $C^1$ 유한요소의 한계를 극복하기 위해 불연속 근사법을 도입한다.
• 수치적 해에 대해 에너지 노름과 $L^2$ 노름에서 최적의 수렴 속도를 확립한다.
• 구조적 및 비구조적 다각형 메쉬에서 광범위한 수치 실험을 통해 이론적 결과를 검증한다.
• 이전에 결과가 없었던 최적 제어 설정에 약한 갈레르킨 방법론을 확장한다.

제안 방법

• 상태 방정식과 노이만 경계 제어를 포함한 비용 기능을 가진 제약 최소화 문제로 최적 제어 문제를 공식화한다.
• 다각형 또는 다면체 요소 위에서 정의된 약한 라플라시안 연산자를 사용한 약한 갈레르킨 변분 형식을 도입한다.
• 상태 변수와 그 약한 기울기의 경우 불연속 유한요소 공간을 사용하여 일반적인 메쉬 유형과 고차 정확도를 가능하게 한다.
• 이산 WG 스킴을 유도하고, inf-sup 안정성 분석을 통해 이산 해의 존재성과 유일성을 증명한다.
• 투영 연산자를 사용하여 정확한 해와 이산 해를 비교함으로써 오차 방정식을 유도한다.
• WG 투영의 안정성 및 근사 성질을 활용하여 에너지 노름과 $L^2$ 노름에서 사전 오차 추정을 적용한다.

실험 결과

연구 질문

• RQ1 약한 갈레르킨 유한요소법은 이차 타원 PDE에 의해 지배되는 최적 제어 문제에 효과적으로 확장될 수 있는가?
• RQ2 이러한 제어 문제에 대해 WG 방법의 에너지 노름과 $L^2$ 노름에서의 수렴 행동은 어떠한가?
• RQ3 정규화 파라미터 $\alpha$는 수치적 해의 조건수와 정확도에 어떤 영향을 미치는가?
• RQ4 일반적인 다각형 또는 다면체 메쉬에서 WG 방법이 최적 순서의 수렴을 달성할 수 있는가?
• RQ5 WG 유한요소 공간에서의 다양한 다항식 차수는 정확도와 수렴 속도에 어떤 영향을 미치는가?

주요 결과

• $P_k$-$P_{k-2}$ WG 유한요소 쌍에 대해 에너지 노름에서 $h^{k}$, $L^2$ 노름에서 $h^{k+1}$의 최적 수렴 속도를 달성한다.
• $P_4$-$P_2$ WG 방법의 경우, 에너지 노름 오차는 메쉬 정밀도와 $\alpha$에 따라 2.8–3.0의 순서로 수렴하고, $L^2$ 오차는 2.7–4.9의 순서로 수렴한다.
• $P_5$-$P_3$ WG 방법의 경우, 에너지 노름 오차는 5.9–6.0의 순서로 수렴하고, $L^2$ 오차는 6.0–6.2의 순서로 수렴하여 최적 수렴을 확인한다.
• 수치 결과는 $\alpha \to 0$ 또는 $\alpha \to \infty$일 때 해가 해석적 극한에 수렴하며, 메쉬 정밀도가 향상될수록 오차가 일관되게 감소함을 확인한다.
• 모든 $\alpha$ 값, 특히 극단적인 값($10^{-8}$에서 $10^9$)에서도 방법은 안정적이고 정확하여 강건성을 입증한다.
• 수치 해의 시각화 결과는 $\alpha$가 증가함에 따라 평탄하거나 구조적인 프로파일로 수렴하는 평탄한 수렴을 보이며, 이는 이론적 기대와 일치한다.