[论文解读] Efficient Implementation of ADER Schemes for Euler and Magnetohydrodynamical Flows on Structured Meshes -- Comparison with Runge-Kutta Methods

Dinshaw S. Balsara, Chad Meyer+3 more|arXiv (Cornell University)|Jun 10, 2010

Computational Fluid Dynamics and Aerodynamics参考文献 47被引用 5

一句话总结

本文提出了一种在结构化网格上求解欧拉方程和磁流体动力学（MHD）方程的ADERS（任意导数）格式的实用且高效实现，提供了一种单步时间积分方法，其性能优于传统的强稳定保模（SSP）龙格-库塔格式。作者在节点空间和模态空间中提供了详细的公式推导，针对散度自由MHD优化了数值通量和电场计算，并证明ADERS格式在所有测试的精度阶次下，计算速度几乎比龙格-库塔方法快一倍。

ABSTRACT

ADER (Arbitrary DERivative in space and time) methods for the time-evolution of hyperbolic conservation laws have recently generated a fair bit of interest. The ADER time update can be carried out in a single step, which is desirable in many applications. However, prior papers have focused on the theory while downplaying implementation details. The purpose of the present paper is to make ADER schemes accessible by providing two useful formulations of the method as well as their implementation details on three-dimensional structured meshes. We therefore provide a detailed formulation of ADER schemes for conservation laws with non-stiff source terms in nodal as well as modal space along with useful implementation-related detail. We also provide details for the efficient use of ADER schemes in obtaining the numerical flux for conservation laws as well as electric fields for divergence-free magnetohydrodynamics. An efficient WENO-based strategy for obtaining zone-averaged magnetic fields from face-centered magnetic fields in MHD is also presented. The schemes catalogued here have been implemented in the first author's RIEMANN code. The speed of ADER schemes is shown to be almost twice as fast as that of strong stability preserving Runge-Kutta time stepping schemes for all the orders of accuracy that we tested.

研究动机与目标

• 弥合理论ADERS格式与计算流体动力学中实际实现之间的差距。
• 为结构化网格上的ADERS方法在节点空间和模态空间中提供详细且可直接使用的公式推导。
• 针对散度自由MHD模拟，优化数值通量和电场计算。
• 提出一种基于WENO的策略，用于从面中心数据计算单元平均磁感应强度。
• 证明ADERS格式在计算速度和精度方面相对于强稳定保模龙格-库塔格式具有显著优势。

提出的方法

• ADERS格式采用单步时间更新实现，避免了龙格-库塔方法的逐阶段推进。
• 该方法在三维结构化网格上，采用节点和模态公式求解具有非刚性源项的双曲守恒律。
• 利用为守恒律量身定制的高阶精度黎曼求解器，高效计算数值通量。
• 在MHD中，电场通过通量积分直接计算，以保持磁感应强度的散度自由。
• 采用基于WENO的重构策略，从面中心的磁感应强度数据计算单元平均磁感应强度。
• 所有格式均在RIEMANN代码中实现并验证，支持高阶精度和性能评估。

实验结果

研究问题

• RQ1 如何在计算流体动力学的结构化网格中高效实现ADERS格式？
• RQ2 ADERS格式实现高性能的关键细节是什么？
• RQ3 ADERS格式的计算成本与强稳定保模龙格-库塔格式在不同精度阶次下的比较如何？
• RQ4 在基于ADERS的MHD模拟中，如何确保磁感应强度的散度自由？
• RQ5 ADERS格式能否在保持高阶精度的同时，实现计算速度接近龙格-库塔方法两倍的提升？

主要发现

• ADERS格式在所有测试的精度阶次下，计算速度几乎比强稳定保模龙格-库塔格式快一倍。
• 在节点空间和模态空间中实现ADERS格式，使得在结构化网格上实现高效且高阶精度的求解成为可能。
• 采用基于WENO的重构策略，可确保从面中心数据准确且稳定地获得单元平均磁感应强度。
• 高效计算数值通量和电场对于保持MHD模拟的精度和散度自由约束至关重要。
• ADERS方法的单步时间积分机制相较于多阶段龙格-库塔方法具有显著的性能优势。
• RIEMANN代码的实现验证了所提出的ADERS格式在复杂流体动力学问题中的鲁棒性和可扩展性。