[논문 리뷰] Support for Non-conformal Meshes in PETSc's DMPlex Interface

Tobin Isaac, Matthew G. Knepley|arXiv (Cornell University)|2015. 08. 11.

Computational Geometry and Mesh Generation참고 문헌 3인용 수 4

한 줄 요약

이 논문은 p4est의 숲 형태의 2진수 트리와 8진수 트리와 같은 비일관 메시를 네이티브로 지원할 수 있도록 PETSc의 DMPlex 인터페이스를 확장한다. 이를 위해 계층적 트리 구조와 참조 트리를 도입하여 부분집합-전체집합 관계와 비일관 메시 인터페이스 유형을 인코딩한다. 이 방법은 떨어진 노드 제약 조건을 자동으로 처리하여, 기존의 일관 메시를 대상으로 한 유한요소 코드가 비일관 메시에서도 그대로 작동하도록 한다.

ABSTRACT

PETSc's DMPlex interface for unstructured meshes has been extended to support non-conformal meshes. The topological construct that DMPlex implements---the CW-complex---is by definition conformal, so representing non- conformal meshes in a way that hides complexity requires careful attention to the interface between DMPlex and numerical methods such as the finite element method. Our approach---which combines a tree structure for subset- superset relationships and a "reference tree" describing the types of non-conformal interfaces---allows finite element code written for conformal meshes to extend automatically: in particular, all "hanging-node" constraint calculations are handled behind the scenes. We give example code demonstrating the use of this extension, and use it to convert forests of quadtrees and forests of octrees from the p4est library to DMPlex meshes.

연구 동기 및 목표

• PETSc의 DMPlex 인터페이스가 2진수 트리와 8진수 트리의 숲과 같은 비일관 계층 메시를 표현할 수 있도록 하기.
• 기존의 일관 메시를 대상으로 한 유한요소 및 유한체적 코드와의 후행 호환성을 유지하기.
• 비일관 인터페이스의 복잡성, 특히 떨어진 노드 제약 조건을 데이터 구조 확장으로써 추상화하기.
• 기본 요소에 대한 피보팅 연산을 통해 비일관 메시 상에서 고급 유한요소 방법(H(curl), H(div))을 지원하기.
• p4est 등의 계층 메시 프레임워크와 PETSc의 광범위한 수치 라이브러리 스택 간의 상호 운용성 촉진하기.

제안 방법

• 기존의 CW-복합체 표현 방식에 계층적 트리 구조를 추가하여 메시 요소 간의 부모-자식 관계를 인코딩하기.
• 비일관 인터페이스 유형(예: 떨어진 노드)의 기하학적 및 위상적 유형을 기술하는 참조 트리를 도입하기.
• 참조 트리를 사용해 떨어진 노드의 제약 행렬을 자동으로 계산하여 사용자 개입 없이 연속성 조건을 강제하기.
• 피보팅 연산(예: 공변 및 반공변 피올라 변환)을 활용해 물리적 요소에서 기준 요소로의 유한요소 함수 매핑을 통해 비일관 인터페이스를 넘어선 연속성 유지하기.
• 메시 계층을 인터페이스 뒤에 감춰서, 유한요소 및 유한체적 솔버의 셀 중심 루프가 일관 메시와 비일관 메시에서 동일한 의미를 가지도록 하기.
• p4est의 2진수 트리 숲과 8진수 트리 숲을 DMPlex 메시로 변환하면서 메시 위상 구조를 유지하고 PETSc의 전체 솔버 및 조절기와의 사용을 가능하게 하기.

실험 결과

연구 질문

• RQ1 비일관 계층 메시는 PETSc의 DMPlex 인터페이스 내에서 어떻게 표현할 수 있는가? 이 인터페이스는 기본적으로 일관 CW-복합체에 기반한다.
• RQ2 기존의 유한요소 및 유한체적 코드와의 호환성을 해치지 않으면서도, 부모-자식 관계와 비일관 인터페이스 유형을 인코딩하기 위해 필요한 데이터 구조는 무엇인가?
• RQ3 참조 트리와 계층 메타데이터를 사용해 적응형 메시 정밀도에서의 떨어진 노드 제약 조건을 자동으로 계산할 수 있는가? 수동 제약 조건 조립이 필요 없어지는가?
• RQ4 기존의 PETSc 유한요소 및 유한체적 솔버가 메시 계층을 추상화함으로써 비일관 메시에서 그대로 작동할 수 있는 정도는 어느 정도인가?
• RQ5 기본 요소에 대한 피보팅 연산을 통해 비일관 메시 상에서 고급 유한요소(H(curl), H(div))를 어떻게 지원할 수 있는가?

주요 결과

• 확장은 DMPlex가 p4est 라이브러리에서 온 2진수 트리 숲과 8진수 트리 숲과 같은 비일관 메시를 성공적으로 표현할 수 있도록 했다.
• 떨어진 노드 제약 조건은 참조 트리와 계층 트리를 사용해 배경에서 자동으로 계산되며, 유한요소 코드에 대한 변경 사항이 전혀 필요 없다.
• 공변 및 반공변 피올라 변환을 요구하는 유한요소 방법도 기준 요소로의 피보팅을 통해 비일관 메시에 적용할 수 있다.
• DMPlex의 DM API와 완전히 호환되며, 하위메시 추출, 메시 분할, PetscFE 및 PetscFV와의 통합 기능을 포함한다.
• 유한체적 방법은 비일관 메시를 지원하도록 확장되고 있으나, 면에 대한 다중 셀 지원과 유량의 중복 계산 문제로 인해 여전히 도전 과제가 남아 있다.
• 이 방법은 다물리스 시뮬레이션에서 유한요소 및 유한체적 필드가 공존할 때 적응형 메시 정밀도를 원활하게 사용할 수 있도록 한다.