[논문 리뷰] A Cut Finite Element Method with Boundary Value Correction

Erik Burman, Peter Hansbo외 1명|arXiv (Cornell University)|2015. 07. 11.

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

한 줄 요약

이 논문은 경계 값 보정을 통한 컷 유한요소법을 제안하며, 기하학적 복잡성을 극복하기 위해 단지 조각별 선형 경계 근사만을 사용함으로써 최적 수렴성을 달성한다. 법선 방향의 테일러 전개를 도입하여 경계 조건을 보정함으로써, 높은 정확도의 경계 표현을 요구하지 않고도 에너지 노름과 $L^2$ 노름에서 최적의 오차 추정을 달성한다.

ABSTRACT

In this contribution we develop a cut finite element method with boundary value correction of the type originally proposed by Bramble, Dupont, and Thomee. The cut finite element method is a fictitious domain method with Nitsche type enforcement of Dirichlet conditions together with stabilization of the elements at the boundary which is stable and enjoy optimal order approximation properties. A computational difficulty is, however, the geometric computations related to quadrature on the cut elements which must be accurate enough to achieve higher order approximation. With boundary value correction we may use only a piecewise linear approximation of the boundary, which is very convenient in a cut finite element method, and still obtain optimal order convergence. The boundary value correction is a modified Nitsche formulation involving a Taylor expansion in the normal direction compensating for the approximation of the boundary. Key to the analysis is a consistent stabilization term which enables us to prove stability of the method and a priori error estimates with explicit dependence on the meshsize and distance between the exact and approximate boundary.

연구 동기 및 목표

• 컷 유한요소법(CutFEM)에서 컷 요소에 대한 정확한 적분 계산의 계산적 과제를 해결하기 위해.
• 경계를 단지 조각별 선형 근사만 사용할 때에도 최적 수렴 속도를 달성할 수 있도록 하기 위해.
• 경계 근사 오차가 존재하는 상황에서도 안정성과 최적 수렴성을 유지하는 안정화된 니츠에 형식을 개발하기 위해.
• 메쉬 크기와 경계 근사 오차에 명시적인 의존성을 갖는 사전 오차 추정을 증명하기 위해.
• 표준 CutFEM에서 요구하는 $O(h^{p+1})$ 수준의 경계 정확도를 회피하면서도 최적 수렴이 가능함을 보여주기 위해.

제안 방법

• 정확한 경계에서의 해를 근사 경계에서의 값과 도함수를 사용하여 법선 방향의 테일러 전개를 통해 근사하기 위해 경계 값 보정을 도입한다.
• 근사 경계에서 딜레르레 조건을 니츠에 유형의 약한 강제 방식으로 구현하며, 기하 오차를 보정하기 위해 테일러 전개를 수정한다.
• 경계 근처에서 함수의 변동을 제어하는 일관된 안정화 항을 적용하여 강제성과 최적 조건 수치 $O(h^{-2})$를 보장한다.
• 가장 가까운 점 사상법을 사용하여 정확한 경계를 정의하고 근사 도메인에서의 일관된 적분을 가능하게 한다.
• 조각별 선형 거리 함수를 사용하여 이산 도메인을 정의하여 기하 계산을 단순화한다.
• 니츠에 방법과 경계 보정 및 안정화를 결합한 대칭 변분형식을 유도하고 분석한다.

실험 결과

연구 질문

• RQ1 CutFEM에서 경계를 단지 조각별 선형 근사로만 사용할 경우 최적 수렴을 달성할 수 있는가? 고차 기하 재구성은 필요하지 않은가?
• RQ2 법선 방향의 테일러 전개가 컷 유한요소법에서 딜레르레 경계 조건의 정확도를 어떻게 향상시키는가?
• RQ3 경계가 요소를 임의로 가로지르는 상황에서 안정성과 최적 조건 수치를 보장하는 안정화 메커니즘은 무엇인가?
• RQ4 경계 근사 오차가 존재하는 상황에서도 제안된 방법이 에너지 노름과 $L^2$ 노름 모두에서 최적 오차 추정을 유지하는가?
• RQ5 최소한의 테일러 전개 항으로도 고차 다항식 근사(예: $P^2$, $P^3$)에 대해 최적 수렴을 달성할 수 있는가?

주요 결과

• 경계가 단지 조각별 선형 근사로만 이루어져 있어도 에너지 노름과 $L^2$ 노름 모두에서 최적 순서의 수렴을 달성한다.
• $P^2$ 및 $P^3$ 요소의 경우, 테일러 전개의 첫 두 항만으로도 최적 수렴 속도를 달성할 수 있다.
• 안정화 항은 선형 시스템의 조건 수가 배경 메쉬에 대한 경계의 위치와 무관하게 $O(h^{-2})$로 유지됨을 보장한다.
• 수치 결과는 $P^2$ 요소에 대해 수렴 속도 2, $P^3$ 요소에 대해 수렴 속도 4를 확인하였으며, 이는 이론적 예측과 일치한다.
• 요소가 부분적으로 컷되어 있어도 활성 도메인을 초월한 이산 해의 연장이 필요 없이 안정성과 정확성을 유지한다.
• 오차 추정은 메쉬 크기 $h$와 정확한 경계와 근사 경계 사이의 거리에 대해 명시적으로 표현되어 정밀한 수렴 분석이 가능하다.