A Staggered Discontinuous Galerkin Method for linear elasticity problem on Polytopal Meshes

기술 요약: 다각형 메쉬 상의 선형 탄성학을 위한 스태거드 불연속 갤러킨 방법 (A Staggered Discontinuous Galerkin Method for Linear Elasticity on Polytotal Meshes)

문제 정의
본 논문은 헬링거-라이스너 혼합 변분 원리(Hellinger-Reissner mixed variational principle)를 사용하여 선형 탄성 문제를 수치적으로 근사하는 문제를 다룬다. 목표는 대칭 응력 텐서 $\sigma \in H(\text{div}, \Omega; S)$와 변위장 $u \in L^2(\Omega; \mathbb{R}^d)$를 동시에 근사하는 것이다. 이 혼합 정식화(mixed formulation)를 이산화할 때 발생하는 주요 과제는 세 가지이다: (i) 안정성을 보장하기 위해 인프-섭(inf-sup) 조건을 만족시키는 것, (ii) 응력 텐서의 대칭성과 요소 경계에서의 법선 연속성(normal continuity)을 모두 강제하는 것, (iii) 최적의 수렴률을 달est하는 것. 기존 방법들은 종종 정점 자유도(vertex degrees of freedom, DoFs)에 의존하여 하이브리드화를 방해하거나, 약한 대칭성(weak symmetry)을 채택하여 국소적인 각운동량 균형을 엄격하게 만족시키지 못할 수 있다. 또한, 많은 방법들이 비압축 한계($\lambda \to \infty$)에서 "락킹(locking)" 현상을 겪는 데 어려움이 있다.

방법론
저자들은 임의의 다각형 메쉬(polytopal meshes) 상에서 정의된 새로운 스태거드 불연속 갤러킨(Staggered Discontinuous Galerkin, SDG) 방법을 제안한다. 이 접근 방식은 프라이멀-듀얼(primal-dual) 메쉬 구조를 활용한다:

1. 메쉬 구성: 프라이말 메쉬 $K_h$는 일반적인 다각형 요소들로 구성된다. 듀얼 메쉬 $K^*_h$는 각 프라이말 요소의 무게중심을 요소의 정점들과 연결함으로써 정교화된 삼각 분할 $T_h$를 생성하여 구축된다.
2. 스태거드 연속성 (Staggered Continuity): 이 방법은 이전의 SDG 탄성학 연구들과는 반대로 연속성 할당을 역전시킨다.
   • 응력 ($\sigma_h$): 정교화된 삼각 분할 $T_h$ 상의 조각별 다항식으로 정의된다. 결정적으로, 응력은 **강한 대칭성(strongly symmetric)**을 가지며 프라이말 요소($K \in K_h$)의 경계에서 법선 연속성을 갖는다.
   • 변위 ($u_h$): 프라이말 요소 $K_h$ 상의 조각별 다항식으로 근사된다. 변위는 듀얼 요소($\omega_F \in K^*_h$)의 경계에 대해 연속적이다.
3. 유한 요소 공간:
   • 응력 공간 $\Sigma_h$는 강한 대칭성을 가진 $H(\text{div})$-적합(conforming) 공간으로 구성된다. 국소적 DoF는 정점 기반의 미지수를 피하기 위해 프라이말 면(face) 상의 법선 플럭스와 내부 모멘트를 통해 정의된다.
   • 변위 공간 $U_h$는 벡터 값 다항식을 갖는 프라이말 셀 상의 함수들로 구성된다.
   • 최저 차수($k=0$)의 경우, 잘 정의된 문제(well-posedness)를 보장하기 위해 응력 공간은 강체 운동(rigid motion) 공간으로 확장된다.
4. 변분 정식화: 이 스킴은 1차 선형 탄성계로부터 유도된다. 여기에는 응력의 발산(divergence)과 정교화된 삼각 분할의 내부 면에서의 점프(jump)를 포함하는 안정화 항이 포함된다. 이러한 항들은 연속성을 강제하기보다는 락킹을 방지하기 위해 설계되었다.
5. 하이브리드화 (Hybridization): 이 방법은 하이브리드화 가능한 구조를 갖는다. 요소 경계에 라그랑주 승수(Lagrange multipliers)를 도입함으로써 응력 변수를 정적 축소(statically condensed)할 수 있으며, 결과적으로 변위의 경계 DoF만을 포함하는 전역 시스템을 얻게 된다.

핵심 기여

• 다각형 상의 강한 대칭성: 본 논문은 정점 DoF에 의依赖하지 않고 임의의 다각형 메쉬 상에서 법선 연속성을 갖는 대칭 응력 공간을 구축한다. 이는 응력 공간의 특정한 기하학적 분해와 DoF의 신중한 선택을 통해 달성된다.
• 역전된 연속성: 약한 대칭성을 강제하는 기존의 SDG 탄성학 연구들과 달리, 본 연구는 프라이말 및 듀얼 메쉬에 대한 응력과 변위의 연속성 패턴을 교체함으로써 강한 대칭성을 강제한다.
• 락킹 프리(Locking-Free) 특성: 본 방법은 락킹이 없는 것으로 증명되었다. 안정성과 수렴성은 응력의 트레이스(trace)를 라메 상수(Lamé parameters)와 독립적으로 제어하는 안정화 항 덕분에 비압축 한계($\lambda \to \infty$)에서도 견고하게 유지된다.
• 국소적 보존: 이 정식화는 선형 운동량의 국소적 균형(발산 방정식을 통해)과 각운동량의 균형(응력 텐서의 강한 대칭성을 통해)을 보존한다.
• 하이브리드 구조: 정점 DoF의 부재는 하이브리드화 절차를 가능하게 하여, 전역 선형 시스템의 크기를 크게 줄이고 구현을 단순화하며 계산 비용을 절감한다.

결과

• 이론적 분석: 저자들은 이산적 인프-섭 조건과 코어시비티(coercivity)를 통해 스킴의 적절성(well-posedness)을 증명하는 포괄적인 이론적 분석을 제공한다. 최적의 사전적(a priori) 오차 추정치가 도출되었으며, 이는 에너지 노름에서의 응력과 이산 $H^1$-유사 노름에서의 변위에 대해 $O(h^{k+1})$의 수렴율을, $L^2$ 노름에서의 변위에 대해 $O(h^{k+2})$의 수렴율을 보여준다.
• 수치 실험:
  • 락킹 프리 검증: 다양한 라메 파라미터($\lambda = 1, 10^2, 10^4, 10^6$)를 가진 삼각형 메쉬에서의 수치 테스트는 오차가 $\lambda$에 독립적임을 보여주며, 이를 통해 방법의 락킹 프리 특성을 확인한다.
  • 수렴율: 삼각형 및 일반 다각형 메쉬 모두에서의 실험은 $k=0$ 및 $k=1$에 대한 이론적 수렴 차수를 확인시켜 준다.
  • 다각형 유연성: 본 방법은 일반적인 다각형 메쉬를 성공적으로 처리하여 표준 심플리셜(simplicial) 또는 사각형 그리드 이상의 적용 가능성을 입증한다.

의의
본 논문은 포아송 방정식[16]에 적용되었던 SDG 프레임워크를 더 복잡한 선형 탄성학 분야로 확장했다고 주장한다. 본 연구의 의의는 강한 대칭성, 법선 연속성, 그리고 하이브리드화 가능성을 결과적으로 정점 DoF 없이 일반 다각형 메쉬 상에서 결합한 방법을 제공한다는 데 있다. 이는 탄성학의 혼합 유한 요소법에서 직면해 온 오랜 난제, 즉 대칭성과 연속성을 동시에 만족시키면서도 정점 DoF를 요구하는 복잡한 요소나 약한 대칭성을 피하는 문제를 해결한다. 제안된 방법은 이론적으로 견고하면서도 하이브리드 구조 덕분에 계산적으로 효율적인 락킹 프리 대안을 제시한다.