[논문 리뷰] Universal Meshes: A new paradigm for computing with nonconforming triangulations

Ramsharan Rangarajan, Adrián J. Lew|arXiv (Cornell University)|2012. 01. 24.

Lattice Boltzmann Simulation Studies참고 문헌 25인용 수 20

한 줄 요약

이 논문은 고차원 유한요소 계산을 위한 새로운 프레임워크인 유니버설 메쉬를 제안한다. 이는 변화하는 곡면 영역에서 고정 배경 메쉬를 사용하여, 가장 가까운 점 투영과 정점 리밸런싱을 통해 배경 메쉬 삼각형을 곡선 요소로 매핑함으로써 메쉬 재구성 없이 정확하게 일치하는 고차원 이산화를 가능하게 한다. 이로 인해 최적의 수렴 속도를 달성하고, 단일 유니버설 메쉬를 사용하여 이동 경계 문제를 효율적으로 시뮬레이션할 수 있다.

ABSTRACT

We describe a method for discretizing planar C2-regular domains immersed in non-conforming triangulations. The method consists in constructing mappings from triangles in a background mesh to curvilinear ones that conform exactly to the immersed domain. Constructing such a map relies on a novel way of parameterizing the immersed boundary over a collection of nearby edges with its closest point projection. By interpolating the mappings to curvilinear triangles at select points, we recover isoparametric mappings for the immersed domain defined over the background mesh. Indeed, interpolating the constructed mappings just at the vertices of the background mesh yields a fast meshing algorithm that involves only perturbing a few vertices near the boundary. For the discretization of a curved domain to be robust, we have to impose restrictions on the background mesh. Conversely, these restrictions define a family of domains that can be discretized with a given background mesh. We then say that the background mesh is a universal mesh for such a family of domains. The notion of universal meshes is particularly useful in free/moving boundary problems because the same background mesh can serve as the universal mesh for the evolving domain for time intervals that are independent of the time step. Hence it facilitates a framework for finite element calculations over evolving domains while using a fixed background mesh. Furthermore, since the evolving geometry can be approximated with any desired order, numerical solutions can be computed with high-order accuracy. We demonstrate these ideas with various numerical examples.

연구 동기 및 목표

• 비일치하는 메쉬에 잠기거나 둘러싸인 평면 C²-정규 영역에 대해 강건하고 고차원 이산화 방법을 개발하는 것.
• 시간에 따라 변화하는 영역 문제에서 메쉬 재구성의 필요성을 제거하기 위해 단일 배경 메쉬를 유니버설 메쉬로 사용하는 것.
• 침입 영역 위에 정확하게 일치하는 곡선 요소를 구성함으로써 유한요소 해에서 최적의 수렴 속도를 달성하는 것.
• 판 굽힘과 같은 경계 표현에 민감한 문제에서 고차원 정확도를 확보하기 위해 곡선 경계를 정확히 근사하는 것.
• 배경 메쉬의 변형만을 사용하여 등가형 매핑과 고차원 곡선 유한요소를 체계적으로 구성하는 프레임워크를 제공하는 것.

제안 방법

• 경계 선분의 가장 가까운 점 투영을 통해 배경 메쉬 삼각형을 침입 영역에 정확히 일치하는 곡선 요소로 매핑하는 방법을 구성한다.
• 영역 내에 적어도 한 정점이 포함된 삼각형을 식별하고, 양 끝점이 영역 외부에 있는 선분을 가장 가까운 점 투영을 통해 경계로 매핑한다.
• 영역 내부 및 경계 근처의 정점을 경계에서 멀리 이동시켜 요소 품질을 향상시키고 기형 삼각형을 방지한다.
• 선택된 점들(예: 정점)에서 구성된 곡선 매핑을 보간하여 배경 메쉬 전체에 등가형 매핑을 복원한다.
• 경계 근처의 좁은 영역 내 정점을 편향시키기 위해 局부 변형 지도 𝔪ₕ를 사용하며, 이때 변형 반경 r은 메쉬 크기의 몇 배로 선택한다.
• 메쉬의 기하 조건(예: 예리한 조건 각도, 충분한 해상도)을 만족시켜 호메오모르피즘과 잘 정의된 투영의 존재를 보장한다.

실험 결과

연구 질문

• RQ1 고정된 배경 메쉬를 사용하여 메쉬 재구성 없이 변화하는 곡면 영역에 대해 고차원 정확하게 일치하는 이산화를 생성할 수 있는가?
• RQ2 비일치하는 메쉬에서 곡선 요소를 어떻게 구성하여 최적의 수렴 속도를 달성할 수 있는가?
• RQ3 배경 메쉬가 침입 영역의 일군의 클래스를 위한 유니버설 메쉬로 기능하기 위해 만족해야 할 기하학적 및 메쉬 조건은 무엇인가?
• RQ4 침입 유한요소 방법에서 고차원 해의 정확도는 경계 표현의 정밀도에 어떻게 의존하는가?
• RQ5 경계 근처에서 정점 리밸런싱을 통해 요소 품질 악화를 방지하면서도 기하학적 일치성을 유지할 수 있는가?

주요 결과

• 이 방법은 비일치하는 배경 메쉬 위에서 정점과 선분의 변형만을 사용하여 침입 영역 위에 정확하게 일치하는 곡선 요소를 구성함으로써 최적의 수렴 속도를 달성한다.
• 메쉬 정점에서 곡선 매핑을 보간함으로써 동일한 배경 메쉬 위에서 고차원 유한요소 방법을 지원하는 등가형 매핑을 복원한다.
• 동일한 배경 메쉬는 시간 간격에 관계없이 변화하는 영역에 대해 유니버설 메쉬로 사용될 수 있으며, 메쉬 재구성의 필요성을 제거한다.
• 경계에서 멀리 떨어진 정점 리밸런싱은 특히 정점들이 경계에 가까이 있을 경우 열악한 형태의 또는 기형 요소가 생성되는 것을 방지한다.
• 고정된 메쉬 구조를 유지함으로써 희소 행렬의 구조와 연결성 패턴을 유지하여 계산 효율성을 크게 향상시킨다.
• 경계 정확도가 해 품질에 결정적으로 영향을 주는 문제, 예를 들어 단순 지지된 원형 판의 굽힘 문제에서 이 방법은 특히 효과적이다.