[논문 리뷰] Function space bases in the dune-functions module

Christian Engwer, Carsten Gräser외 2명|arXiv (Cornell University)|2018. 06. 25.

Advanced Numerical Methods in Computational Mathematics참고 문헌 6인용 수 4

한 줄 요약

이 논문은 Dune-Fem 유한요소 소프트웨어 생태계 내에서 함수 공간 기저를 체계적으로 표현하고 조작하기 위한 프레임워크를 제시한다. 트리 기반 추상화를 사용하여 혼합 및 벡터값 유한요소 공간을 모델링하며, 복잡한 기저 구조를 통해 자유로운 다중 인덱싱을 가능하게 한다. 이는 스토크스 방정식에 대한 테일러–후드 유한요소 방법에 적용되어, Dune 2.6의 dune-functions 모듈을 통해 유형 안정성, 확장 가능성, 효율성 있는 구현을 입증한다.

ABSTRACT

The dune-functions Dune module provides interfaces for functions and function space bases. It forms one abstraction level above grids, shape functions, and linear algebra, and provides infrastructure for full discretization frameworks like dune-pdelab and dune-fem. This document describes the function space bases provided by dune-functions. These are based on an abstract description of bases for product spaces as trees of simpler bases. From this description, many different numberings of degrees of freedom by multi-indices can be derived in a natural way. We describe the abstract concepts, document the programmer interface, and give a complete example program that solves the stationary Stokes equation using Taylor-Hood elements.

연구 동기 및 목표

• Dune 생태계 전반에서 유한요소 함수 공간과 그 기저의 표현을 통합하고 추상화하여 코드 중복과 분열을 줄이기.
• 혼합 및 벡터값 공간에 특히 적합한 복잡한 유한요소 기저를 표현하기 위한 유연하고 확장 가능하며 유형 안정적인 인터페이스를 제공하기.
• 곱공간의 계층적 트리 구조를 반영하는 다중 인덱싱 체계를 체계적으로 구성하여 자유도의 인덱싱을 가능하게 하기.
• 다중 인덱스를 평탄하고 연속적인 국소 번호 체계로 매핑하여 국소 행렬 조립을 효율적으로 지원하기.

제안 방법

• 간단하고 더 이상 분해할 수 없는 기저로 구성된 트리로서 유한요소 공간을 모델링하며, 스칼라 또는 더 이상 분해할 수 없는 공간을 잎으로, 곱 연산을 내부 노드로 사용한다.
• 트리 구조를 반영하는 다중 인덱스 번호 체계를 생성하기 위한 전략 기반 시스템을 정의하여 평탄한 인덱싱과 중첩 인덱싱을 모두 가능하게 한다.
• dune-istl 라이브러리의 중첩 벡터 및 행렬 유형을 사용하여 트리 계층에 대응하는 다차원 데이터 구조를 표현한다.
• 각 격자 요소에 대해 형상 함수 값, 그라디언트 및 자코비안을 통한 변환을 계산하는 국소 기저 인터페이스를 구현한다.
• 행렬 조립을 위해 트리 구조의 다중 인덱스를 평탄한 국소 인덱스로 매핑하기 위해 `localView.tree().child(...).localIndex()` 메서드를 사용한다.
• 스토크스 문제에 대해 테일러–후드 요소를 사용하여 프레임워크를 적용하며, 별도의 속도 및 압력 기저 트리를 사용하고, 적분을 통한 복소 행렬 조립(속도-속도, 속도-압력, 압력-속도 항)을 수행한다.

실험 결과

연구 질문

• RQ1 혼합 및 벡터값 유한요소 공간을 다수의 기저 유형과 인덱싱 체계를 지원하는 방식으로 추상적으로 표현할 수 있는 방법은 무엇인가?
• RQ2 트리 구조 기반 기저에서 다중 인덱스를 평탄한 국소 자유도로 매핑하는 데 있어, 타입 안정성, 확장 가능성, 효율성을 동시에 확보할 수 있는 방법은 무엇인가?
• RQ3 동일한 프레임워크 내에서 노달, 계층적, 비연속 갈레르킨 유한요소와 같은 다양한 유한요소 유형을 지원하기 위해 동일한 추상 기저 인터페이스를 어떻게 활용할 수 있는가?
• RQ4 함수 공간 기저의 트리 기반 추상화를 사용하여 자유도의 다양한 번호 체계를 체계적으로 유도하고 관리할 수 있는가?
• RQ5 결과로 도출된 인터페이스를 dune-pdelab 및 dune-fem과 같은 기존 Dune 모듈에 통합하여 코드 재사용과 마이그레이션을 어떻게 지원할 수 있는가?

주요 결과

• dune-functions 모듈은 트리 기반 표현을 통해 유한요소 함수 공간 기저를 성공적으로 추상화하여 혼합 및 벡터값 공간의 일관되고 확장 가능한 모델링을 가능하게 한다.
• 프레임워크는 노달, 계층적, 단항식 유한요소 기저를 포함한 다수의 기저 유형을 하나의 인터페이스 내에서 지원하며, 동일한 코드베이스 내에서 상호 교환 가능하게 한다.
• 전략 기반 다중 인덱싱 시스템은 트리 구조 기반 기저에서 평탄한 국소 번호 체계로의 매핑을 매우 자유롭게 가능하게 하여 효율적인 행렬 조립에 핵심적이다.
• 테일러–후드 요소를 사용한 스토크스 시스템의 정확하고 효율적인 조립이 성공적으로 구현되었으며, 트리 인터페이스를 통해 속도 및 압력 자유도를 별도로 처리하였다.
• `localView.tree().child(...).localIndex()`의 사용은 구성 요소 간 형상 함수의 유형 안정성과 정확한 인덱싱을 보장하며, 수동으로 인덱스를 기록할 필요가 없어진다.
• 이 방법은 dune-pdelab 및 dune-fem에 대한 원활한 통합을 가능하게 하였으며, 두 프로젝트 모두 dune-functions를 핵심 기저 추상화로 사용하기 위해 마이그레이션을 진행 중이다.