Conforming/Non-conforming Virtual Elements and application to elasticity problems in curved three-dimensional domains

기술 요약: 곡선 3 차원 영역의 탄성 문제 적용을 위한 정합/비정합 가상 요소

문제 제기
가상 요소법 (VEM) 은 일반적인 다면체 메쉬에서 편미분 방정식 (PDE) 을 이산화하기 위한 유연한 프레임워크로 자리 잡았습니다. 그러나 곡선 영역이 포함된 3 차원 문제에 VEM 을 적용하는 데에는 여전히 상당한 과제가 남아 있습니다. 곡선 경계를 가진 2 차원 VEM 은 잘 정립되어 있지만, 3 차원 구현에서는 곡면에서의 공간 정의와 필요한 사영 (projection) 계산에 어려움이 있습니다. 기존 접근법은 종종 복잡한 마스터 - 슬레이브 구성, 약한 경계 조건 부과, 또는 불연속 변수를 가진 혼합 형식에 의존합니다. 고체 역학의 맥락에서 불연속 변위장을 생성하는 방법은 물리적으로 비현실적입니다. furthermore, 직선 모서리를 가진 표준 다면체 메쉬는 고차 스킴의 경우 해의 오차를 지배할 수 있는 기하학적 근사 오차를 도입합니다.

방법론
저자들은 메쉬 면을 정합 (conforming) 또는 비정합 (non-conforming) 으로 선택적으로 지정할 수 있는 새로운 정합/비정합 (C-NC) 가상 요소법을 제안합니다.

• 하이브리드 공간 구성: 이 방법은 다면체 $E$의 경계 면을 정합 집합 ($\mathcal{F}^c_E$) 과 비정합 집합 ($\mathcal{F}^n_E$) 으로 분할하여 다면체 $E$ 위의 국소 가상 공간 $W^k_h(E)$를 정의합니다.
  • 정합 면에서는 해가 연속 ($C^0$) 이며, 공간은 다항식 모멘트로 강화된 표준 노드 VEM 정의를 따릅니다.
  • 비정합 면에서는 해가 불연속일 수 있으며, 법선 방향 미분은 $k-1$ 차 다항식으로 제한됩니다.
  • 전역 공간은 정합 면을 가로질러만 연속성을 부과함으로써 구성됩니다. 이 "하이브리드" 요소는 정합 면에서는 디리클레, 비정합 면에서는 뉴만 경계 조건을 갖는 혼합된 국소 푸아송 문제를 통해 구축됩니다.
• 자유도 (DoFs): 공간은 정합 경계에서의 꼭짓점 및 모서리 점에 대한 점별 평가와 면에서의 함수 및 그 법선 방향 미분의 모멘트로 정의됩니다. 특히, 강성 행렬에 필요한 사영 (구체적으로 $L^2$ 및 에너지 사영) 은 혼합 경계 조건이 있더라도 이러한 DoFs 로만 계산 가능함이 입증되었습니다.
• 곡선 영역 적용: 이 방법은 곡선 영역의 3 차원 선형 탄성 문제에 적용됩니다. 전략은 다음과 같습니다:
  • 전역 변위 연속성을 보장하기 위해 영역 내부에는 평면 정합 면을 사용합니다.
  • 정확한 기하학을 표현하기 위해 경계 $\partial\Omega$ 에만 곡선 비정합 면을 사용합니다.
  • 매끄러운 매개변수화 $\gamma_F$를 통해 매핑된 다항식을 사용하여 곡면 위의 가상 공간을 정의합니다.
  • 사영 연산자에서의 경계 적분을 비정합 DoFs 를 사용하여 근사하며, 이는 사영의 법선 방향 미분이 매핑된 다항식 공간에 속하지 않을 때 필요한 근사임을 인정합니다.

주요 기여

1. 새로운 하이브리드 요소: 이 논문은 면의 정합/비정합 특성이 사용자 정의 선택 사항인 VEM 형식을 도입하여, 고전적 정합 VEM, 비정합 VEM, 또는 둘의 혼합을 복원할 수 있는 "하이브리드" 요소를 가능하게 합니다.
2. 이론적 분석: 저자들은 C-NC 공간에 대한 엄밀한 안정성 및 수렴 분석을 제공합니다. $H^1$-반노름과 DoF 기반 안정화 사이의 동등성을 증명하고, 보간 추정을 수립하며, 선형 탄성 문제에 대한 최적 오차 상한을 유도합니다.
3. 곡선 3 차원 탄성: 이 논문은 정확한 기하학 표현을 가진 곡선 영역에서 3 차원 탄성 문제를 해결하는 실용적 방법론을 보여줍니다. 불연속 변위를 초래하거나 복잡한 구성을 요구하는 이전 3 차원 곡선 VEM 접근법과 달리, 이 방법은 정확한 곡선 경계를 포착하는 동시에 벌크 (bulk) 에서 연속적인 변위장을 보장합니다.
4. 수치적 검증: 이 방법은 세 가지 시나리오에서 테스트되었습니다:
   • 표준 직선 기하학 (수렴 속도 검증).
   • 정확한 메쉬 근사를 가진 곡선 기하학 (원통).
   • 선형 및 2 차 기하학적 근사를 사용하여 임의의 비평면 다각형으로 테셀레이션된 곡선 기하학 (구).

결과
수치 실험은 이론적 예측을 확인합니다:

• 수렴 속도: 평면 및 곡선 영역 모두에서 이 방법은 다항식 차수 $k$에 대해 $L^2$-노름에서 $O(h^{k+1})$, $H^1$-반노름에서 $O(h^k)$의 최적 수렴 속도를 달성합니다.
• 곡선 경계 성능: 정확한 곡선 경계 표현을 가진 원통의 경우, 이 방법은 기하학적 근사 오차 ($O(h^2)$ 및 $O(h^{3/2})$) 로 고통받는 사면체 메쉬 위의 표준 2 차 유한 요소 스킴보다 우수한 최적 속도를 유지합니다.
• 기하학적 근사 한계: 구 예제의 경우, 이 방법은 임의의 다각형 테셀레이션을 성공적으로 처리합니다. 그러나 저자들은 $k=2$의 경우 기하학적 근사가 충분히 정확할 때만 최적 수렴 속도가 회복됨을 관찰합니다. 2 차 기하학적 근사 (Quad) 를 사용할 때, $L^2$ 오차는 처음에는 더 높은 차수의 수렴을 보이지만 결국 정체됩니다. 저자들은 이는 $k=2$ 요소에 대한 최적 $L^2$ 수렴에 필요한 $O(h^3)$ 기하학적 오차를 보장하지 않는 세렌디피티 (serendipity) 기저의 사용으로 인해 기하학적 근사 오차 (특히 표면적 및 부피 오차) 가 지배적이 되기 때문이라고 설명합니다.

의의 및 주장
이 논문은 제안된 C-NC VEM 이 곡선 영역의 3 차원 탄성 문제를 위한 유연하고 견고한 프레임워크를 제공한다고 주장합니다. 그 주요 의의는 연속적인 변위 해를 제공 (고체 역학에 필수적) 하면서 비정합 면을 사용하여 곡선 경계를 정확하게 표현할 수 있다는 점에 있으며, 이는 기존 표준 VEM 또는 FEM 접근법으로는 달성하기 어려웠던 조합입니다. 저자들은 이 방법이 곡면 위의 경계 적분 사영에 근사를 도입하지만, 기하학적 근사가 적절하다면 이는 방법의 안정성이나 최적 수렴 속도를 훼손하지 않는다고 강조합니다. 이 작업은 비물리적인 불연속 해나 지나치게 복잡한 마스터 - 슬레이브 구성에 의존하지 않고 복잡한 3 차원 기하학을 처리할 수 있음을 보여주는 개념 증명입니다.