[논문 리뷰] All-Hex Meshing Strategies For Densely Packed Spheres

Yu-Hsiang Lan, Paul Fischer외 2명|arXiv (Cornell University)|2021. 06. 01.

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

한 줄 요약

이 논문은 구의 중심의 베이저의 분해를 이용하여 조밀하게 배열된 구의 전체 6면체 메시 생성 전략을 제시한다. 이 전략은 베이저 면을 사각형 요소로 테셀레이션하고, 이를 구 표면으로 스윕하여 메시를 생성한다. 이 방법은 고순도 스펙트럼 요소 방법(Spectral Element Method, SEM) 시뮬레이션에 적합한 높은 품질, 낮은 요소 수의 메시를 생성하며, 엑사스케일 시스템에서 최적화된 메시 품질과 솔버 성능을 기반으로 최대 9900만 개의 요소와 505억 개의 격점까지 구현하였다.

ABSTRACT

We develop an all-hex meshing strategy for the interstitial space in beds of densely packed spheres that is tailored to turbulent flow simulations based on the spectral element method (SEM). The SEM achieves resolution through elevated polynomial order N and requires two to three orders of magnitude fewer elements than standard finite element approaches do. These reduced element counts place stringent requirements on mesh quality and conformity. Our meshing algorithm is based on a Voronoi decomposition of the sphere centers. Facets of the Voronoi cells are tessellated into quads that are swept to the sphere surface to generate a high-quality base mesh. Refinements to the algorithm include edge collapse to remove slivers, node insertion to balance resolution, localized refinement in the radial direction about each sphere, and mesh optimization. We demonstrate geometries with 10^2-10^5 spheres using approximately 300 elements per sphere (for three radial layers), along with mesh quality metrics, timings, flow simulations, and solver performance.

연구 동기 및 목표

• 조밀하게 배열된 구로 이루어진 흐름의 난류 시뮬레이션을 위한 스펙트럼 요소 방법(Spectral Element Method, SEM)에 특화된 고품질 전체 6면체 메시 생성 전략을 개발하기 위해.
• 메시 품질을 유지하면서 요소 수를 줄여 좁은 간극 공간에서의 적합한 메시 생성 문제를 해결하기 위해.
• 고다항수 순서(N=7–15)의 대규모 시뮬레이션을 가능하게 하기 위해 요소 수를 최소화하고 솔버 효율을 향상시키기 위해.
• 높은 순서 방법에 대해 높은 요소 수와 낮은 확장성 문제를 야기하는 삼각형-에서-육면체 변환 및 페비드/플라스터링 방법의 한계를 극복하기 위해.
• 엑사스케일 아키텍처에서의 시뮬레이션 준비를 위해 메시 품질, 자코비안 조건 및 압력 푸아송 해법 성능을 최적화하기 위해.

제안 방법

• 메서드는 구의 중심에 대한 베이저 분해로 시작하여, 구 표면에 의해 둘러싸인 볼록 셀로 간극 공간을 분할한다.
• 각 베이저 면은 중간 노드 삽입을 통해 정점 수를 균일하게 유지함으로써, 전적으로 사각형 요소로 이루어진 메시로 테셀레이션된다.
• 사각형 면은 구 표면으로 반경 방향으로 스윕되어 구조적인 층을 형성하는 고품질의 육면체 기반 메시를 생성한다.
• 미세한 요소를 제거하기 위해 모서리 수축을 적용하고, 특히 각 구 주변의 반경 방향에서 해상도를 균형 있게 유지하기 위해 노드 삽입을 사용한다.
• 메시 최적화에는 다중 단계의 라플라시안 스무딩 및 최적화 단계가 포함되어 비율을 개선하고 척도화된 자코비안을 향상시킨다.
• 최종 노드 점들은 스펙트럼 요소 방법(Spectral Element Method, SEM) 호환성을 위해 곡면 구 표면으로 투영되며, 고순도 이산화를 위해 GLL 점이 사용된다.

실험 결과

연구 질문

• RQ1 베이저 기반 전체 6면체 메시 생성 전략이 조밀하게 배열된 구에 대해 요소 수를 크게 줄이고 높은 메시 품질을 달성할 수 있는가?
• RQ2 고순도 스펙트럼 요소 방법(Spectral Element Method, SEM) 시뮬레이션에서 삼각형-에서-육면체 변환 방법과 비교해 본다면, 제안된 방법은 요소 수, 메시 품질 및 솔버 성능 측면에서 어떻게 다른가?
• RQ3 메시 최적화와 반경 방향 정밀 조절을 통해 압력 푸아송 시스템의 조건수는 얼마나 향상될 수 있는가?
• RQ4 10만 개 이상의 구와 십억 개 이상의 격점이 포함된 전체 핵심로드 반응기 시뮬레이션에 대해 이 방법이 확장 가능한가?
• RQ5 점 접촉 특이성을 둥근 다리로 대체할 경우, 솔버 수렴성과 성능에 어떤 영향을 미치는가?

주요 결과

• 이 메서드는 352,000개의 구에 대해 최대 9900만 개의 스펙트럼 요소와 505억 개의 격점까지 생성하여 전체 핵심 반응기 시뮬레이션의 확장성을 입증하였다.
• 최적화 후 메시에는 음수 자코비안이 없었으며, 최소 척도화된 자코비안은 1.93e-2였고, 최대 비율은 30.8이었으며, 이는 높은 품질을 의미한다.
• 4,608개의 노드(27,648개의 V100 GPU)에서 1단계당 0.233초의 월드 타임을 기록했으며, 조정된 다중 격자 솔버를 사용해 단계당 약 6회의 압력 반복을 수행하였다.
• 유사한 해상도에서 삼각형-에서-육면체 방법보다 전체 6면체 접근 방식이 더 낮은 쿠라운트 수를 기록하고 더 빠른 단계를 제공하여 더 큰 시간 간격과 더 낮은 시뮬레이션 비용을 가능하게 하였다.
• 이 메서드는 압력 반복 수를 줄이고 솔버 성능을 향상시켰으며, 점 접촉을 둥근 다리로 대체할 경우 조건수 감소가 기대된다.
• 메시 생성 파이프라인은 MATLAB 루프로 인해 병목 현상이 있었지만, 가장 큰 케이스에서는 10분 이내에 전체 메시 생성을 완료했으며, C 기반의 Mex 함수를 통해 성능 향상이 가능할 것으로 기대된다.