[논문 리뷰] On the Design of Locking Free Ghost Penalty Stabilization and the Relation to CutFEM with Discrete Extension

Erik Burman, Peter Hansbo외 1명|arXiv (Cornell University)|2022. 05. 03.

Numerical methods for differential equations인용 수 4

한 줄 요약

이 논문은 컷 유한요소법(CutFEM)에 대한 새로운 유령 보정 안정화 기법을 제안하며, 이는 이전의 접근 방식을 일반화하여 임의의 함수형식(예: 개별 자유도)의 안정화를 가능하게 하고, 면 이웃 외의 요소 간 연결도 허용한다. 주요 기여는 절점 기반 안정화가 강건성 조건을 만족함을 증명함으로써, 안정화 계수를 무한대로 갈수록 락킹(locking) 없이도 가능하게 하여, 컷유한요소법에서 이산 확장 연산자와의 직접적인 등가성을 확립한다.

ABSTRACT

In this note, we develop a new stabilization mechanism for cut finite element methods that generalizes previous approaches of ghost penalty type in two ways: (1) The quantity that is stabilized and (2) The choice of elements that are connected in the stabilization. In particular, we can stabilize functionals of the discrete function such as finite element degrees of freedom. We subsequently show that the kernel of our ghost penalty operator defines a finite element space based on discrete extensions in the spirit of those introduced in Burman, E.; Hansbo, P. and Larson, M. G., CutFEM Based on Extended Finite Element Spaces, arXiv2101.10052, 2021.

연구 동기 및 목표

• 커트 유한요소법에서 표준 요소 기반 안정화를 넘어서는 일반화된 유령 보정 안정화 프레임워크를 개발하는 것.
• 해의 이산화된 함수에 대한 임의의 함수형식(예: 개별 자유도)의 안정화를 가능하게 하여 불안정한 모드를 더 정밀하게 제어하는 것.
• 일반화된 안정화가 안정화 계수를 무한대로 갈수록 락킹을 피함으로써 강건성을 확보하는지 보여주는 것.
• 유령 보정 안정화와 컷유한요소법에서 이산 확장 연산자 간의 직접적인 이론적 및 알고리즘적 연결 고리를 확립하는 것.
• 두 번째 순서 타원 문제에 대해 최적 순서 수렴성과 조건수 추정치를 지원하는 통합된 프레임워크를 제공하는 것.

제안 방법

• 이산 함수의 기능적 형태(예: 절점 자유도 포함)를 안정화 대상으로 하는 일반화된 유령 보정 형태를 제안한다.
• 메쉬 크기 비례 거리 이내의 요소를 통해 경계에 의해 잘린 요소를 안정화할 수 있는 융통성 있는 결합 전략을 도입한다. 이는 이웃 면 요소 외의 요소 간 안정화를 허용한다.
• 자유도에 직접 작용하는 안정화 형태를 정의함으로써, 안정화 계수 파arameter가 무한대에 갈수록 제약 조건을 강하게 구현할 수 있도록 한다.
• 이산 확장 연산자 프레임워크를 사용하여 유령 보정이 확장 기반 안정화와 관련되며, 안정화 계수 파arameter가 매우 클 경우 등가성을 보여준다.
• 표준 확장에 대한 역부등식과 안정성 추정치를 적용하여 안정화 형태의 경계를 유도한다.
• 보정 항의 국소화 및 강도에 충분한 조건이 만족될 경우, 이산 양선형 형식의 강제성(coercivity)을 확립한다.

실험 결과

연구 질문

• RQ1 유령 보정 안정화는 표준 요소 기반 양상 외에, 이산 해의 임의의 함수형식(예: 개별 자유도)을 안정화할 수 있는가?
• RQ2 안정화 계수 파arameter가 무한대로 갈수록 일반화된 유령 보정 형식이 여전히 강건하고 락킹이 없는가?
• RQ3 제안된 절점 기반 안정화는 컷유한요소법에서 이산 확장 연산자 프레임워크와 어떻게 관련이 있는가?
• RQ4 요소 간 연결 조건이 컷유한요소법에서 안정성과 최적 수렴성을 보장하기 위해 어떤 조건을 가져야 하는가?
• RQ5 일반화된 유령 보정 형태가 안정성과 최적 오차 추정치를 위한 추상적 프레임워크를 만족하는지 입증할 수 있는가?

주요 결과

• 제안된 절점 기반 안정화는 안정화 계수 파arameter가 무한대로 갈수록 강건성과 최적 수렴 조건을 만족한다.
• 유령 보정 연산자의 핵(kernel)은 이산 확장으로 정의된 유한요소 공간과 정확히 일치하며, 이는 안정화와 확장 접근법 간의 직접적인 등가성을 수립한다.
• 일반화된 유령 보정 형태는 두 번째 순서 타원 문제에 대해 최적 순서 사전 오차 추정치와 조건수 경계를 보장한다.
• 필요한 강제성 및 일致성 조건을 만족함으로써, 안정화 계수 파arameter가 무한대에 갈수록 락킹 없이 안정성과 최적 수렴성을 달성한다.
• 이론적 프레임워크는 비교적 동반 논문들과의 비교를 통해 고차수 타원 문제로 자연스럽게 확장 가능하다.
• 분석 결과, 결합의 국소화가 충분하고 보정 파arameter가 적절히 선택될 경우 안정화 형태가 강제성을 만족함을 확인하였다.