[논문 리뷰] A stabilized cut finite element method for partial differential equations on surfaces: The Laplace-Beltrami operator

Erik Burman, Peter Hansbo외 1명|arXiv (Cornell University)|2013. 12. 04.

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

한 줄 요약

이 논문은 삼각형 메시에 임bed된 표면에서 라플라스-베르트라미 방정식을 해결하기 위해 안정화된 컷 유한요소 방법을 제안한다. 연속적인 조각별 선형 함수를 표면에 제한하여 사용하며, 정상적인 기울기 점프를 제어하는 일관된 안정화 항을 도입함으로써 표면의 위치와 관계없이 최적의 조건수 범위를 확보하고 에너지 및 $L^2$ 노름에서 최적의 $O(h^2)$ 수렴을 가능하게 한다.

ABSTRACT

We consider solving the Laplace-Beltrami problem on a smooth two dimensional surface embedded into a three dimensional space meshed with tetrahedra. The mesh does not respect the surface and thus the surface cuts through the elements. We consider a Galerkin method based on using the restrictions of continuous piecewise linears defined on the tetrahedra to the surface as trial and test functions. The resulting discrete method may be severely ill-conditioned, and the main purpose of this paper is to suggest a remedy for this problem based on adding a consistent stabilization term to the original bilinear form. We show optimal estimates for the condition number of the stabilized method independent of the location of the surface. We also prove optimal a priori error estimates for the stabilized method.

연구 동기 및 목표

• 배경 메시에 대한 표면의 위치에 민감한 표준 컷 유한요소 방법의 심각한 악조건화 문제를 해결한다.
• 배경 메시에 대한 표면의 위치에 따라 조건수가 변하는 민감성을 해소한다.
• 최적의 수렴 속도를 유지하면서 수치적 안정성을 향상시키는 일관된 안정화 항을 개발한다.
• 표면 기하학에 영향을 받지 않는 조건수 및 오차 추정의 이론적 보장을 수립한다.
• 복잡한 표면 PDE 문제를 위한 강력한 프레임워크를 제공하며, 향후 헬름홀츠 방정식 및 결합된 밀도-표면 문제로의 확장 가능성을 고려한다.

제안 방법

• 배경의 삼각형 요소에 정의된 연속적인 조각별 선형 유한요소를 표면에 제한한다.
• 표면과 교차하는 요소 면에서 정규 기울기 점프를 제어하는 안정화 항을 도입한다.
• 조건수 개선을 위해 안정화된 이차형식을 사용하는 일관된 갈레르킨 방법을 설정한다.
• 이산 피카르에 유사한 추정을 활용하여 표면의 위치와 관계없이 조건수를 제한한다.
• 오차 분석을 위해 표면의 이산 함수를 연속 표면으로 옮기는 리프팅 연산자를 적용한다.
• 접선 미적분학과 표면 미분 연산자(라플라스-베르트라미)를 사용하여 표면에서의 약한 설정을 정의한다.

실험 결과

연구 질문

• RQ1 표면에서 불일치하는 유한요소 방법의 악조건화를 제거할 수 있는 일관된 안정화 항을 설계할 수 있는가?
• RQ2 안정화된 방법이 라플라스-베르트라미 문제에 대해 에너지 및 $L^2$ 노름에서 최적의 수렴 속도를 달성하는가?
• RQ3 이산 시스템의 조건수가 배경 메시에 대한 표면의 위치와 관계없이 유계로 유지될 수 있는가?
• RQ4 실제로 안정화된 방법은 안정화되지 않은 방법에 비해 조건수 및 수렴 행동에 어떤 영향을 미치는가?
• RQ5 이론적 결과를 헬름홀츠 방정식과 같은 더 복잡한 표면 PDE로 확장할 수 있는가?

주요 결과

• 안정화된 방법은 배경 메시 내 표면의 위치와 관계없이 최적의 조건수 상한 $O(h^{-2})$ 를 확보한다.
• 이산 해에 대해 에너지 및 $L^2$ 노름에서 최적의 $O(h^2)$ 수렴 속도가 증명된다.
• 수치 실험은 안정화된 방법이 최적의 수렴 속도를 유지하는 반면, 안정화되지 않은 방법은 악조건화 및 허위 영 고유값을 보임을 확인한다.
• 안정화 없이 조건수는 $O(h^{-2})$ 로 증가하지만, 안정화로 인해 조건수가 크게 향상되며, 특히 대각 스케일링과 조합될 경우 더욱 뚜렷하다.
• 안정화 항은 이산 함수의 법선 도함수를 효과적으로 제어하며, 표면의 임의의 위치에 의한 악조건화를 방지한다.
• 표면이 요소를 가로질러도 방법은 강건하며, 안정화 항은 약한 설정과 일관된다.