[논문 리뷰] Adaptive and anisotropic finite element approximation : Theory and algorithms

Jean‐Marie Mirebeau|arXiv (Cornell University)|2011. 01. 07.

Computational Geometry and Mesh Generation참고 문헌 49인용 수 5

한 줄 요약

이 박사학위 논문은 고정된 원소 수를 가진 삼각형 메쉬 위에서 근사 오차를 최소화하는 방식으로 적응형 이방성 유한요소 근사의 이론적 및 알고리즘적 기반을 수립한다. 부드러운 함수에 대해 渐近적 오차 감쇠율을 증명하고, 최적의 메쉬 적응과 관련된 함수 클래스를 규명하며, 이산 최적화 문제의 연속적 일반화를 도출하고, 근사적으로 최적의 메쉬를 생성하는 계층적 정밀화 전략을 개발한다.

ABSTRACT

Mesh adaption procedures for finite element approximation allows one to adapt the resolution, by local refinement in the regions of strong variation of the function of interest. This procedure plays a key role in numerous applications of scientific computing. The use of anisotropic triangles allows to improve the efficiency of the procedure by introducing long and thin triangles that fit in particular the directions of the possible curves of discontinuity. Given a norm X of interest and a function f to be approximated, we formulate the problem of optimal mesh adaptation, as minimizing the approximation error over all (possibly anisotropic) triangulations of prescribed cardinality. We address the four following questions related to this problem: I. How does the approximation error behave in the asymptotic regime when the number of triangles N tends to infinity, when f is a smooth function ? II. Which classes of functions govern the rate of decay of the approximation error as N grows, and are in that sense naturally tied to the problem of optimal mesh adaptation? III. Could this optimization problem, which is posed on triangulations of a given cardinality N, be replaced by an equivalent more tractable problem posed on a continuous object? IV. Is it possible to produce a near-optimal sequence of triangulations using a hierarchical refinement procedure?

연구 동기 및 목표

• 이방성 삼각형 메쉬를 사용한 최적의 메쉬 적응을 위한 엄밀한 수학적 프레임워크를 개발하기 위해.
• 원소 수가 증가함에 따라 근사 오차가 어떻게 渐近적으로 감소하는지 이해하기 위해.
• 이방성 메쉬 적응과 자연스럽게 부합하는 함수의 정(regularity)이 무엇인지 규명하기 위해.
• 이산 메쉬 최적화 문제를 더 다룰 수 있는 연속 문제로 재구성하기 위해.
• 근사적으로 최적의 삼각형 메쉬를 생성하는 계층적 정밀화 알고리즘을 설계하기 위해.

제안 방법

• 주어진 원소 수 N에 대해 최적의 메쉬 적응 문제를 모든 삼각형 메쉬의 집합에서 근사 오차를 최소화하는 문제로 수식화한다.
• N → ∞ 일 때 부드러운 함수에 대한 근사 오차의 渐近적 행동을 분석한다.
• Sobolev 공간과 가중치가 있는 Sobolev 공간을 근사 오차 감쇠율을 결정짓는 함수 클래스로 규명한다.
• 이방성 특성을 표현하기 위해 리만 메트릭을 사용하여 이산 메쉬 최적화 문제의 연속적 일반화를 도입한다.
• 연속적 일반화를 바탕으로 후행 오차 추정량을 유도하여 계층적 정밀화를 안내한다.
• 연속적 수식에서 유도된 오차 지표에 기반해 요소를 적응적으로 정밀화하는 계층적 정밀화 알고리즘을 제안한다.

실험 결과

연구 질문

• RQ1 부드러운 함수에 대해 원소 수 N이 무한대에 가까워질 때 근사 오차는 어떻게 행동하는가?
• RQ2 근사 오차 감쇠율을 결정짓는 함수 클래스는 무엇이며, 이는 왜 이방성 메쉬 적응과 자연스럽게 부합하는가?
• RQ3 이산 메쉬 최적화 문제를 더 다룰 수 있는 연속 문제로 대체할 수 있는가?
• RQ4 계층적 정밀화 전략을 사용하여 근사적으로 최적의 삼각형 메쉬 시퀀스를 생성할 수 있는가?

주요 결과

• 부드러운 함수에 대해 최적의 이방성 삼각형 메쉬를 사용할 경우, L^2 노름에서 근사 오차는 O(N^{-p/2}) 의 속도로 감쇠한다.
• 감쇠율은 가중치가 있는 Sobolev 공간에서의 함수의 정(regularity)에 의해 결정되며, 더 높은 정(regularity)은 더 빠른 수렴을 가능하게 한다.
• 이산 메쉬 최적화 문제는 이방성을 표현하는 리만 메트릭을 사용하여 효과적으로 연속 문제로 일반화할 수 있다.
• 연속적 일반화는 오차 추정의 이론적 근거를 제공하며, 효율적인 정밀화 알고리즘 설계를 안내한다.
• 제안된 계층적 정밀화 알고리즘은 주어진 함수 클래스에 대해 근사 오차 측면에서 증명된 바 근사적으로 최적의 삼각형 메쉬를 생성한다.
• 수치 실험을 통해 이 방법이 최적 수렴 속도를 달성하고, 방향성 특성을 지닌 문제에서 등방성 메쉬보다 정확도를 크게 향상시킴을 확인하였다.