[논문 리뷰] Numerical Approximation in Riemannian Manifolds by Karcher Means

Stefan W. von Deylen|arXiv (Cornell University)|2015. 05. 14.

Geometric Analysis and Curvature Flows참고 문헌 7인용 수 6

한 줄 요약

이 논문은 비정상적인 곡률을 가진 리만 다양체를 근사하기 위한 수치적 방법을 제안하며, 카르커 평균을 사용해 중심단순체를 정의함으로써 결과적인 사상이 형태 정규성 조건 하에 거의 등거리사상이 되도록 보장한다. 주요 기여는 중심단순체 사상의 미분과 그 기울기가 유계임을 증명한 것으로, 이는 곡률 다양체 위에서 갈레르킨 유한요소 방법을 적용하고 고전적인 유클리드 방법과 유사한 오차 추정치를 확보할 수 있게 한다.

ABSTRACT

(1) For a compact Riemannian manifold without boundary $(M,g)$ containing $n+1$ points $p_i$ and the $n$-dimensional standard simplex $Δ$, the miniser of \[ E: M imes Δ o {\mathbf R}, (a,λ) \mapsto λ^0 d^2(a,p_0) + \dots + λ^n d^2(a,p_n) \] is considered as point with "barycentric coordinates" $λ_i$ within the so-called Karcher simplex (or Riemannian simplex or geodesic finite element) defined by vertices $p_i$. In the small, existence and uniqueness is well-known. Now suppose $Δ$ carries a flat Riemannian metric $g^e$ induced by edge lengths $d(p_i,p_j)$, where $d$ is the geodesic distance in $M$. If all edge lengths are small than $h$ and $vol(Δ,g^e) \geq αh^n$ for some $α> 0$, then we can show that \begin{equation} |(x^*g - g^e)(v,w)| \leq c h^2 |v| |w|, \qquad |( abla^{x^*g} - abla^{g^e})_v w| \leq c h |v| |w| \end{equation} with some constant $c$ depending only on the curvature tensor $R$ of $(M,g)$ and $α$. From this we derive several estimates for Finite Element calculations in which $(M,g)$ is replaced by a piecewise flat realised simplicial complex. (2) Let $M$ be the geometric realisation of a simplicial complex $K$. The simplicial cohomology $(C^k(K), \partial^*)$ has been interpreted as "discrete outer calculus" (DEC) in the literature. We define spaces $P^{-1}Ω^k \subset L^\inftyΩ^k$ and outer differentials and give an isometric cochain map $C^k o P^{-1}Ω^k$. This reduces the computation of variational problems in discrete outer calculus to variational problems in a trial space of non-conforming differential forms. We investigate the approximation properties of $P^{-1}Ω^k$ in $H^1Ω^k$ and compare the solutions to variational problems in both spaces.

연구 동기 및 목표

• 비정상적인 곡률을 가진 리만 다양체에서 계산적으로 실현 가능하고 기하학적으로 의미 있는 단순체의 개념을 정의하기 위해.
• 카르커 중심단순체 사상이 단사적이며 타당한 중심좌표를 유도할 수 있는 조건을 설정하기 위해.
• 중심단순체 사상의 왜곡에 대한 정량적 추정치를 제공하여, 이 사상이 곡률에 따라 결정되는 상한 내에서 등거리사상처럼 행동함을 보장하기 위해.
• 조각별 평탄한 단순체 근사화된 다양체 위에서 갈레르킨 유한요소 방법을 적용하고 고전적인 유클리드 유한요소 이론과 유사한 오차 제어를 가능하게 하기 위해.
• 내재된 기하학적 구성에 기반해 이산 기하 방법과 곡률 다양체 위의 변분 PDE 사이의 격차를 메우기 위해.

제안 방법

• 표준 단순체에서 다양체로의 중심단순체를 카르커 평균을 통해 정의된 중심단순체 사상 $ x: riangle \to M $ 의 이미지로 구성한다. 이 때 주어진 점 $ p_i \notin M $ 의 리만 중심을 사용한다.
• 리만 중심(카르커 평균)을 사용하여 표준 단순체에서 다양체로의 유일하고 연속적인 중심단순체 사상을 정의한다.
• 다양체와 점의 구성에 대한 형태 정규성 조건을 활용하여 미분 $ dx $ 와 그 기울기 $ \nabla dx $ 의 유계성을 확보함으로써 거의 등거리사상 성질을 보장한다.
• 조스와 카르커(1982)의 $ \nabla dx $ 추정 이론을 응용하고 갈레르킨 근사 이론과 융합하여 오차 추정치를 유도한다.
• 유한요소 공간의 $ L^2 $-최적 근사 성질과 파oincaré 유형 부등식을 활용하여 갈레르킨 해의 오차를 제어한다.
• 특히 $ \nabla W $ 가 $ Q = dx - Pdy $ 의 $ L^2 $-최적 근사임을 이용하여, 보간과 기울기 제어를 통해 오차 추정치를 도출한다.

실험 결과

연구 질문

• RQ1 카르커 중심단순체 사상 $ x $ 가 그 이미지 위에서 미분동형사상이 되는 조건은 무엇인가? 즉, 단순체 위에 타당한 중심좌표를 유도할 수 있는가?
• RQ2 중심단순체 사상 $ x $ 의 왜곡은 일반 좌표계와 비교해 어떻게 되며, 곡률과 메쉬 크기의 함수로 유계일 수 있는가?
• RQ3 조각별 평탄한 단순체 근사화된 리만 다양체 위의 변분 문제에 대해 갈레르킨 유한요소 방법을 엄밀하게 적용할 수 있는가?
• RQ4 카르커 기반의 단순체 메esh를 사용할 때, 다양체 위의 PDE 근사에 대한 갈레르킨 근사의 오차 추정치는 무엇인가?
• RQ5 특히 정확한 지오데식 거리 계산의 비용이 높은 점을 감안할 때, 중심단순체의 암묵적 정의를 어떻게 계산적으로 실현 가능하게 할 수 있는가?

주요 결과

• 카르커 중심단순체 사상의 미분 $ dx $ 는 $ |dx| \triangleq \text{ess sup} |\nabla dx| \triangleq C_{0,1}^{\text{primes}} h $ 를 만족하며, 여기서 $ h $ 는 메쉬 크기이고 $ C_{0,1}^{\text{primes}} $ 는 곡률과 형태 정규성에 의존한다.
• 미분의 기울기 $ \nabla dx $ 는 $ L^\nu $-노름에서 유계이며, $ \nabla dx \triangleq \text{ess sup} |\nabla dx| \triangleq C_{0,1}^{\text{primes}} h $ 로 표현되어 사상이 거의 등거리사상임을 보장한다.
• 카르커 단순체 위의 갈레르킨 해는 최적 오차 추정치를 달성한다: $ \rho_0(x,y) + \rho_1(x,y) \triangleq \rho_{0,\text{bdry}}(x,y) + h \nabla dx + C_{0,1}^{\text{primes}} h $, 상수들은 다양체 기하학과 차원에 따라 달라진다.
• 이 방법은 중심단순체 사상이 변분 공식화에 필요한 필수 기하학적 구조를 유지함으로써, 곡률 다양체 위에서 표준 유한요소 기법을 적용할 수 있도록 한다.
• 형태 정규성과 곡률 조건이 만족되면, 갈레르킨 근사의 오차 추정치는 고전적인 유클리드 유한요소 이론과 동일한 순서의 오차를 보인다.
• 제곱 거리 함수 $ \textbf{d}^2 $ 의 근사화에 대해 이 방법이 강건함을 시사하며, 크랭크 등(2013)의 빠른 지오데식 거리 해법을 통해 실용적인 구현이 가능하다.