[논문 리뷰] A wavelet-adaptive method for multiscale simulation of turbulent flows in flying insects

Thomas Engels, Kai Schneider외 2명|arXiv (Cornell University)|2019. 12. 11.

Fluid Dynamics and Turbulent Flows인용 수 26

한 줄 요약

이 논문은 인공 압축성과 체적 벌점 기법을 사용하여 복잡하고 시간에 따라 변화하는 기하구조에서 난류성 불압축성 유동을 시뮬레이션하기 위한 웨이브렛 적응형 블록 기반 수치 방법을 제시한다. 특히 나비의 비행을 대상으로 한다. 이 방법은 이중직교 웨이브렛 임계처리를 통해 동적으로 격자를 적응시킴으로써 기존의 푸리에 준스펙트럴 방법 대비 최대 160배의 속도 향상을 달성하여 최적의 오차 제어와 계산 효율성을 확보한 고성능 대량 병렬 시뮬레이션을 가능하게 한다.

ABSTRACT

We present a wavelet-based adaptive method for computing 3D multiscale flows in complex, time-dependent geometries, implemented on massively parallel computers. While our focus is on simulations of flapping insects, it can be used for other flow problems, including turbulence, as well. The incompressible fluid is modeled with an artificial compressibility approach in order to avoid solving elliptical problems. No-slip and in/outflow boundary conditions are imposed using volume penalization. The governing equations are discretized on a locally uniform Cartesian grid with centered finite differences, and integrated in time with a Runge--Kutta scheme, both of 4th order. The domain is partitioned into cubic blocks with equidistant grids with different resolution and, for each block, biorthogonal interpolating wavelets are used as refinement indicators and prediction operators. Thresholding the wavelet coefficients allows to generate dynamically evolving grids, and an adaption strategy tracks the solution in both space and scale. Blocks are distributed among MPI processes and the global topology of the grid is encoded using a tree-like data structure. Analyzing the different physical and numerical parameters allows balancing their individual error contributions and thus ensures optimal convergence while minimizing computational effort. Different validation tests score accuracy and performance of our new open source code, WABBIT (Wavelet Adaptive Block-Based solver for Interactions with Turbulence), on massively parallel computers using fully adaptive grids. Flow simulations of flapping insects demonstrate its applicability to complex, bio-inspired problems.

연구 동기 및 목표

• 복잡하고 시간에 따라 변화하는 기하구조(예: 나비의 비행)에서 난류 흐름을 시뮬레이션하기 위한 고성능 적응형 수치 방법을 개발하는 것.
• 웨이브렛 기반 정밀화를 통해 공간과 스케일에 따라 동적으로 계산 격자를 적응시킴으로써 정확하고 효율적인 대규모 시뮬레이션을 가능하게 하는 것.
• 압축성, 벌점, 이산화, 임계처리 등의 다수의 수치 오차를 균형 잡고 계산 비용을 최소화하는 것.
• 과학적 활용을 위해 개방소스 코드(WABBIT)에 구현하여 생물학적 영감을 받은 유체역학 분야에서 널리 사용할 수 있도록 하는 것.
• 나무 기반의 블록 구조적 데이터 레이아웃을 사용하여 수천 개의 CPU 코어에서의 확장성과 효율성을 입증하는 것.

제안 방법

• 타원형 압력 푸아송 방정식을 피하기 위해 인공 압축성을 사용하여 명시적 시간 적분이 가능하도록 한다.
• 복잡하고 움직이는 기하구조에서의 슬립 금지 및 유입/유출 경계 조건을 몸체에 맞는 격자 없이도 체적 벌점 기법을 적용하여 구현한다.
• 유체역학적 방정식은 국소적으로 균일한 직각좌표 격자 위에서 4차 중심 유한차분법으로 이산화하며, 블록 구조적 정밀화를 적용한다.
• 이중직교 보간 웨이브렛은 정밀화 지표와 예측 연산자로 사용되어 계수의 임계처리를 통해 동적 격자 적응을 가능하게 한다.
• 나무 구조의 데이터 구조가 계층적 격자 구조를 인코딩하며, 블록은 MPI 프로세스 간에 분배되어 고스트 노드 교환을 통한 병렬 계산을 수행한다.
• 블록 크기는 압축 효율성과 현대 아키텍처에서의 CPU 성능을 균형 잡기 위한 수치적 파rameter로 조정된다.

실험 결과

연구 질문

• RQ1 웨이브렛 기반의 적응형 정밀화가 복잡하고 시간에 따라 변화하는 기하구조에서 난류 흐름에 대해 최적의 수렴성과 높은 계산 효율성을 달성할 수 있는가?
• RQ2 압축성, 벌점, 이산화, 임계처리 등의 다수의 수치 오차를 어떻게 균형 잡아 정확성을 유지하면서 계산 비용을 최소화할 수 있는가?
• RQ3 블록 기반의 웨이브렛 적응형 접근 방식이 기존의 균일 격자 또는 준스펙트럴 방법에 비해 속도와 확장성 측면에서 어느 정도 뛰어난가?
• RQ4 이 방법은 고해상도와 동적 격자 적응을 동시에 구현하면서도 병렬 HPC 환경에서 실제 나비의 비행을 시뮬레이션할 수 있는가?
• RQ5 블록 크기의 선택이 웨이브렛 적응 프레임워크 내에서 성능과 압축 효율성에 어떤 영향을 미치는가?

주요 결과

• 웨이브렛 적응형 방법은 동일한 정확도에서 고해상도 꿀벌 시뮬레이션에 대해 균일 격자 기반의 푸리에 준스펙트럴 코드(FLUSI) 대비 최대 160배의 빠른 계산 성능을 달성한다.
• 이 방법은 수천 개의 CPU 코어를 통해 높은 병렬 효율성을 유지하며, 격자 포인트 및 시간 단계당 성능이 최적화된 준스펙트럴 코드 FLUSI와 유사하다.
• 웨이브렛 계수의 임계처리를 통해 솔루션의 특징이 공간과 스케일 모두에서 추적되는 동적 격자 적응이 가능하여 정확도를 유지하면서도 계산 비용을 감소시킨다.
• 블록 기반의 구현은 현대 아키텍처에서 캐시 활용도를 향상시키고 성능을 향상시켰으며, 블록 내 일부 비의미 있는 점이 포함되더라도 여전히 유리하다.
• 이 방법은 나비의 비행과 함께 분수형 나무를 연계한 시뮬레이션을 성공적으로 수행하여, 표준 준스펙트럴 방법을 초월한 복잡하고 다스러진 환경 상호작용의 가능성을 입증한다.
• 개방소스 WABBIT 코드는 이동하는 복잡한 기하구조에서 난류 흐름의 재현 가능하고 확장 가능한 시뮬레이션을 가능하게 하며, 특히 생물학적 영감을 받은 유체역학 분야에 적합하다.