[논문 리뷰] Obtaining higher-order Galerkin accuracy when the boundary is polygonally approximated

Todd F. Dupont, Johnny Guzmán외 1명|arXiv (Cornell University)|2020. 01. 09.

Advanced Numerical Methods in Computational Mathematics참고 문헌 9인용 수 4

한 줄 요약

이 논문은 다각형으로 근사된 영역에서 포아송 방정식을 해결하기 위한 대칭적 고차원 갈레르킨 방법을 제안한다. 표준 유한요소 방법에서의 수렴 속도 저하를 유발하는 기하 오차를 보정함으로써, 최적의 수렴 속도를 달성한다. 진짜 경계와의 부호 거리 함수를 기반으로 한 경계층 보정을 변분 공식에 적용하여, 대칭성을 유지하면서도 최적의 수렴 속도(예: 2차 요소에 대해 $O(h^3)$)를 확보한다. 이는 이전의 비대칭적 접근 방식(예: Bramble-Dupont-Thomée 방법)의 한계를 극복한다.

ABSTRACT

We study two techniques for correcting the geometrical error associated with domain approximation by a polygon. The first was introduced some time ago \cite{bramble1972projection} and leads to a nonsymmetric formulation for Poisson's equation. We introduce a new technique that yields a symmetric formulation and has similar performance. We compare both methods on a simple test problem.

연구 동기 및 목표

• 매끄러운 영역을 다각형으로 근사한 경우 표준 갈레르킨 유한요소 방법의 수렴 속도가 최적치에 못 미치는 문제를 해결하기 위해.
• 고차 다항식 근사(예: 2차 및 3차 요소)에 대해 최적의 수렴 속도를 복원하는 대칭적 변분 공식을 개발하기 위해.
• 기존의 고차 보정 방법에서 나타나는 비대칭 이차형식 문제를 해결하기 위해.
• 영역 근사로 인한 기하 오차가 부호 거리 함수와 법선 미분을 포함하는 경계 적분 항을 통해 보정될 수 있음을 보여주기 위해.

제안 방법

• 부호 거리 $\delta$와 법선 미분을 포함하는 경계항을 추가함으로써 표준 갈레르킨 공식을 수정하는 대칭 이차형식 $M_h(u,v)$를 도입한다.
• 일致성과 대칭성을 확보하기 위해 $\int_{\partial\Omega_h} \left(\gamma \frac{\delta}{h} - 1\right) \left( \delta \frac{\partial u}{\partial n} \frac{\partial v}{\partial n} + \frac{\partial u}{\partial n} v + \frac{\partial v}{\partial n} u \right) ds$ 와 같은 보정 항을 포함한다.
• 경계 근처에서 해의 타일러 전개를 이용하여 다각형 경계와 진짜 경계 간의 격차를 반영하는 보정 항을 유도한다.
• 체적 및 경계 적분을 포함하는 변분 공식을 사용하여 딜레르흐 경계 조건을 가진 포아송 방정식에 적용한다.
• 안정성과 최적 수렴을 확보하기 위해 안정화 파rameter $\gamma$를 선택한 로빈 유형 경계 조건을 변분 형식에 적용한다.
• 곡선 경계를 다각형으로 근사한 매끄러운 영역에서 수치적으로 방법을 검증하며, 표준 방법 및 비대칭 보정 방법과 비교한다.

실험 결과

연구 질문

• RQ1 영역이 다각형으로 근사된 경우에도 최적의 수렴 속도를 복원할 수 있는 대칭적 변분 공식을 구성할 수 있는가?
• RQ2 대칭성을 유지하면서 고차 정확도를 확보할 수 있는 최적의 경계 보정 항의 구조는 무엇인가?
• RQ3 대칭적 고차 보정 방법을 사용할 경우, 표준 방법 또는 비대칭 방법과 비교해 오차 분포는 어떻게 변화하는가?
• RQ4 부호 거리 함수 $\delta$를 사용할 경우, 표준 Nitsche 유형 방법에 비해 수렴 속도 향상 정도는 어느 정도인가?

주요 결과

• 제안된 대칭 방법 $M_h$는 다각형 경계 근사에도 불구하고(piecewise quadratic finite elements에 대해) $H^1$-노름에서 최적의 $O(h^3)$ 수렴 속도를 달성한다.
• 수치 실험 결과, 대칭 방법은 비대칭 Bramble-Dupont-Thomée 방법과 동일한 수렴 속도를 보이며, $k=2$일 경우 $L^2$ 오차는 $O(h^3)$, $k=3$일 경우 $O(h^4)$로 감소한다.
• 대칭 방법의 오차 분포는 경계 부근에 집중되어 있으며, 내부 오차는 제로로 감소함을 보여, 오염 효과의 효과적인 억제를 나타낸다.
• 환형 영역 테스트 케이스에서 $k=2$일 때 $h_{\text{max}} = 0.036$에서 $L^2$ 오차가 $1.54 \times 10^{-5}$로 측정되어 최적 수렴이 확인된다.
• 이 방법은 $\delta^\ell \partial^n u / \partial n^n$ 항을 통해 고차 보정을 통합하면서도 안정성과 대칭성을 유지하며, 고차 정확도를 달성한다.
• 비대칭 시스템의 계산적 단점 없이도 대칭 형식을 유지함으로써, 반복 해법 및 대규모 시뮬레이션에 더 적합하다.