[논문 리뷰] Simulating Turbulence Using the Astrophysical Discontinuous Galerkin Code TENET

Andreas Bauer, Kevin Schaal외 4명|arXiv (Cornell University)|2016. 02. 29.

Computer Graphics and Visualization Techniques인용 수 4

한 줄 요약

이 논문은 천체물리학적 우주론적 맥락에서 비음속 난류를 시뮬레이션하기 위해 설계된 새로운 비연속 갈레르킨(DG) 코드인 TENET을 소개한다. 고차수 DG 이산화와 밀도 있는 스텐실, 효율적인 병렬 처리를 통해 TENET는 전통적인 유한체적 방법과 유사한 정확도를 달성하면서도 자유도를 1.8배 적게 사용하고, 1.5배 더 빠르게 실행하여 뛰어난 계산 효율성을 입증하며, 에크사스케일 천체물리학 시뮬레이션의 잠재력을 보여준다.

ABSTRACT

In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications.

연구 동기 및 목표

• 비연속 갈레르킨(DG) 방법을 사용하여 천체물리학적 시뮬레이션을 위한 고차수, 효율적이고 확장 가능한 유체역학 해법기 개발.
• 계산 천체물리학의 벤치마크 문제인 비음속 난류 시뮬레이션에서 DG 방법의 성능과 정확도 평가.
• 계산 비용, 정확도, 확장성 측면에서 TENET의 DG 기반 시뮬레이션과 전통적인 2차 유한체적 스킴을 비교.
• 특히 매끄럽고 난류가 있는 유량에 대해 향후 에크사스케일 컴퓨팅에서 DG 방법의 잠재력 평가.
• 정확도를 유지하거나 향상시키면서도 유한체적 방법보다 더 높은 계산 효율성을 달성할 수 있는지 여부 조사.

제안 방법

• 코드는 각 요소 내에서 고차수 다항식 기저 함수를 사용하여 해를 표현하는 비연속 갈레르킨(DG) 이산화를 사용하여 올리어 방정식을 이산화한다.
• DG 방법은 약한 수식을 사용하며, 요소 경계에서 수치적 플럭스를 통해 해를 구해 고차수 정확도와 국소적 보존성을 가능하게 한다.
• 요소 내에서 모드 또는 노드 기반 표현을 사용하여 임의의 정확도 순서를 달성할 수 있으며, 오직 직접 인접 요소 간의 통신만 필요하다.
• 시간 적분은 안정성과 정확도를 유지하기 위해 강안정성 보존 강타워-쿠타 스킴을 사용한다.
• 구조적 직각좌표 메쉬 기반으로 구현되었으며, 적응형 메쉬 정밀도(AMR) 지원 기능을 통해 난류 구조의 효율적 해상도를 가능하게 한다.
• MPI를 통한 도메인 분할을 통한 병렬 처리가 이루어졌으며, SUPERMUC 시스템에서 최대 4096코어까지 뛰어난 강한 스케일링 성능을 보였다.

실험 결과

연구 질문

• RQ1 TENET와 같은 고차수 비연속 갈레르킨 코드가 비음속 난류 시뮬레이션에서 기존 2차 유한체적 스킴과 동일한 정확도를 달성할 수 있는가?
• RQ2 유사한 정확도 수준을 확보할 때, DG 방법의 계산 비용이 난류 유동 시뮬레이션에서 유한체적 방법보다 어떻게 비교되는가?
• RQ3 유한체적 스킴 대비 DG 방법이 주어진 정확도를 확보하기 위해 필요한 자유도의 수를 얼마나 줄일 수 있는가?
• RQ4 병렬 환경에서 해상도와 프로세서 수가 증가함에 따라 DG 방법의 성능은 어떻게 스케일링되는가?
• RQ5 매끄러운 천체물리학적 흐름에서 DG 접근법이 전통적인 유한체적 방법보다 보존 성질(예: 각운동량 또는 div B 제약 조건)에서 우월한가?

주요 결과

• TENET의 DG 구현은 비음속 난류 시뮬레이션에서 2차 유한체적 방법 대비 약 1.8배 적은 자유도로 동일한 정확도를 달성한다.
• 유사한 난류 파wer 스펙트럼을 확보한 경우, 제3차 정확도의 DG 시뮬레이션은 최고 해상도에서 해당되는 2차 정확도의 유한체적 시뮬레이션보다 1.53배 더 빠르게 실행된다.
• DG 방법의 성능 우월성은 해상도가 높아질수록 증가한다: 최고 해상도에서 DG_128_3 실행은 FV_512_2 실행보다 1.53배 더 빠르며, 자유도 감소율은 약 1.86배에 가까워진다.
• DG 방법은 SUPERMUC에서 최대 4096코어까지 뛰어난 강한 병렬 스케일링 성능을 보이며, 대규모 분산 시스템에서의 높은 효율성을 입증한다.
• DG 방법은 매끄러운 흐름 영역에서 각운동량을 명백히 보존하며, 이는 표준 유한체적 스킴이 본질적으로 만족하지 못하는 성질이다.
• DG 접근법은 기저 함수 수준의 제약 조건을 통해 MHD에서 div B = 0 제약 조건을 강력하게 구현할 수 있어, 향상된 MHD 시뮬레이션을 위한 길을 열어준다.