[论文解读] Simulating Turbulence Using the Astrophysical Discontinuous Galerkin Code TENET

Andreas Bauer, Kevin Schaal+4 more|arXiv (Cornell University)|Feb 29, 2016

Computer Graphics and Visualization Techniques被引用 4

一句话总结

本文介绍了TENET，一种专为宇宙学背景下模拟亚音速湍流而设计的新一代天体物理不连续伽辽金（DG）代码。通过采用高阶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的模拟与传统二阶有限体积格式的性能表现。
• 评估DG方法在天体物理学中面向未来百亿亿次计算的潜力，尤其针对平滑的湍流流动。
• 探究DG方法是否能在保持或提升精度的同时，实现优于有限体积方法的计算效率。

提出的方法

• 该代码采用不连续伽辽金（DG）方法对欧拉方程进行离散化，通过每个单元内的高阶多项式基函数表示解。
• DG方法采用方程的弱形式，通过单元界面的数值通量求解，从而实现高阶精度和局部守恒。
• 方法在单元内使用模态或节点基表示，仅需与直接邻近单元通信，即可实现任意阶精度。
• 时间积分采用强稳定性保持的龙格-库塔格式，以维持稳定性和精度。
• 代码在结构化笛卡尔网格上实现，并支持自适应网格细化（AMR），可高效解析湍流结构。
• 通过MPI实现域分解并行化，该方法在SUPERMUC系统上表现出优异的强可扩展性，最高达4096个核心。

实验结果

研究问题

• RQ1 像TENET这样的高阶不连续伽辽金代码能否在模拟亚音速湍流时达到与传统二阶有限体积格式相当的精度？
• RQ2 在实现相同精度水平的湍流模拟中，DG方法的计算成本与有限体积方法相比如何？
• RQ3 与有限体积格式相比，DG方法在达到给定精度时，能将自由度数量减少多少？
• RQ4 在并行环境中，DG方法的性能如何随分辨率和处理器数量的增加而变化？
• RQ5 在平滑的天体物理流动中，DG方法是否在守恒性质（如角动量或div B约束）方面优于传统有限体积方法？

主要发现

• TENET中的DG实现可在亚音速湍流模拟中，以约1.8倍更少的自由度达到与二阶有限体积方法相当的精度。
• 在湍流功率谱相近的模拟中，三阶DG运行比最高分辨率下对应的二阶有限体积运行快1.53倍。
• DG的性能优势随分辨率提升而增强：在最高分辨率下，DG_128_3运行比FV_512_2运行快1.53倍，自由度减少接近1.86倍。
• DG方法在SUPERMUC系统上表现出优异的强并行可扩展性，最高达4096个核心，表明其在大规模分布式系统中具有高效率。
• DG方法在平滑流动区域明显保持角动量守恒，而标准有限体积格式不具备此内在特性。
• DG方法通过基函数层面的约束，能够稳健地实现MHD中div B = 0的约束，为改进MHD模拟提供了可行路径。