[논문 리뷰] Parallel two-scale finite element implementation of a system with varying microstructures

Omar Richardson, Omar Lakkis외 2명|arXiv (Cornell University)|2021. 03. 31.

Advanced Mathematical Modeling in Engineering참고 문헌 44인용 수 4

한 줄 요약

이 논문은 영역 미분형과 풀백 연산자를 사용하여 비균질 미세구조를 기준 영역으로 매핑함으로써 다스케일 PDE 시스템의 스케일링 가능하고 정확한 시뮬레이션을 가능하게 하는 병렬 이중 스케일 유한요소법을 제시한다. 이는 식물 조직 내 영양소 이동에 대한 사례에서 검증된 수렴성과 뛰어난 병렬 속도 향상을 보여준다.

ABSTRACT

We propose a two-scale finite element method designed for heterogeneous microstructures. Our approach exploits domain diffeomorphisms between the microscopic structures to gain computational efficiency. By using a conveniently constructed pullback operator, we are able to model the different microscopic domains as macroscopically dependent deformations of a reference domain. This allows for a relatively simple finite element framework to approximate the underlying PDE system with a parallel computational structure. We apply this technique to a model problem where we focus on transport in plant tissues. We illustrate the accuracy of the implementation with convergence benchmarks and show satisfactory parallelization speed-ups. We further highlight the effect of the heterogeneous microscopic structure on the output of the two-scale systems. Our implementation (publicly available on GitHub) builds on the deal.II FEM library. Application of this technique allows for an increased capacity of microscopic detail in multiscale modeling, while keeping running costs manageable.

연구 동기 및 목표

• 다양하고 비균질적인 미세구조를 가진 다스케일 PDE 시스템을 위한 계산 효율적인 유한요소 프레임워크를 개발하기 위해.
• 각 미세영역에 대해 반복적인 메쉬 생성을 피하기 위해 미세영역을 기준 영역으로 매핑하는 미분형을 사용하여 복잡한 미세구조 기하구조를 해소하는 데서 계산 비용을 감소시키기 위해.
• 도메인 분할과 다중 스레딩을 통해 병렬 계산을 가능하게 하여 HPC 자원에서의 확장성을 확보하기 위해.
• 결합된 타원형 PDE를 포함한 영양소 이동 모델에 대해 제안된 방법의 정확성과 효율성을 입증하기 위해.
• 생물학 및 재료 과학 분야의 다스케일 모델링을 위한 공개 가능하고 확장 가능한 프레임워크를 제공하기 위해, deal.II 기반으로 구축한다.

제안 방법

• 다양한 미세영역을 단일 기준 영역으로 매핑하기 위해 풀백 연산자를 사용하여 반복적인 메쉬 생성을 방지한다.
• 강한 결합이 아닌 인터페이스 항을 통해 약한 결합된 방식으로 매크로스케일 및 마이크로스케일 문제를 동시에 해결하는 이중 스케일 유한요소 설정을 사용한다.
• 매크로스케일 좌표에 따라 의존하는 미분형을 사용하여 변형된 기준 영역에서 마이크로스케일 PDE를 해결한다.
• 매크로스케일 도메인 분할과 각 하위도메인 내 다중 스레딩을 통해 병렬화를 지원하여 효율적인 HPC 배포를 가능하게 한다.
• 이 방법은 식물 조직 내 영양소 이동, 소비 및 생성을 모델링하는 결합된 타원형 PDE 시스템에 적용된다.
• 수렴 기준 테스트를 통해 유한요소 근사의 기대 이론적 수렴 순서가 확인된다.

실험 결과

연구 질문

• RQ1 각각의 미세영역에 대해 개별 메쉬를 생성하지 않고도 다스케일 PDE 시스템 내 비균질 미세구조를 효율적으로 모델링할 수 있는 방법은 무엇인가?
• RQ2 변동하는 미세구조 기하구조가 수송 문제에서 매크로스케일 및 마이크로스케일 해의 프로파일에 미치는 영향은 무엇인가?
• RQ3 이중 스케일 유한요소법이 최적의 수렴 속도를 달성하면서도 계산 확장성을 유지할 수 있는가?
• RQ4 다양한 프로세서에서의 속도 향상과 로드 밸런싱 측면에서 제안된 병렬 구현의 효과성은 어떠한가?
• RQ5 이 프레임워크는 식물 영양소 이동 외의 다른 다스케일 문제로 일반화될 수 있는가?

주요 결과

• 유한요소 근사는 기대 이론적 수렴 순서를 보이며, 이는 방법의 정확성을 확인한다.
• 병렬 구현은 만족스러운 속도 향상을 달성하여 테스트한 노드 수에서 양호한 확장성을 보여준다.
• 비균질적인 미세구조 기하구조는 수치 결과에서 매크로스케일 및 마이크로스케일 해의 프로파일에 상당한 영향을 미친다.
• 미분형과 풀백 연산자의 사용은 각 도메인에 대해 명시적인 메쉬 생성 없이 다양한 미세구조를 효율적으로 모델링할 수 있도록 한다.
• 프레임워크는 GitHub에 공개되어 있으며, 재현 가능성과 확장성을 보장하기 위해 deal.II 유한요소 라이브러리 기반으로 구축되어 있다.
• 복잡하고 공간적으로 변화하는 미세구조가 존재하더라도 계산 효율성이 유지되어 높은 해상도의 모델링을 관리 가능한 비용으로 수행할 수 있다.