[论文解读] High-order low-storage explicit Runge-Kutta schemes for equations with quadratic nonlinearities

Pablo D. Mininni, Mininni, Pablo D.+2 more|ArXiv.org|Aug 13, 2008

Computational Fluid Dynamics and Aerodynamics参考文献 21被引用 3

一句话总结

本文通过在JST算法中引入校正项，将高阶、低存储的显式Runge-Kutta格式扩展至具有二次非线性的方程，实现了三阶和四阶精度，同时保持了最小的内存占用。该方法在Burgers方程和不可压Euler方程上得到验证，仅带来适度的计算成本增加，但显著提升了时间精度和稳定性。

ABSTRACT

We show in this paper that third- and fourth-order low storage Runge-Kutta algorithms can be built specifically for quadratic nonlinear operators, at the expense of roughly doubling the time needed for evaluating the temporal derivatives. The resulting algorithms are especially well suited for computational fluid dynamics. Examples are given for the Hénon-Heiles Hamiltonian system and, in one and two space dimensions, for the Burgers equation using both a pseudo-spectral code and a spectral element code, respectively. The scheme is also shown to be practical in three space solving the incompressible Euler equation using a fully parallelized pseudo-spectral code.

研究动机与目标

• 为解决JST算法在非线性二次系统中仅具二阶精度的问题（尽管对线性问题为高阶），提出改进方案。
• 推导适用于JST算法的校正项，使其在右端项包含二次非线性时恢复三阶和四阶精度。
• 保持原始JST格式的低存储特性，确保在内存受限的大规模模拟中具有适用性。
• 在代表性问题（包括一维和二维Burgers方程、多维不可压Euler方程）上验证新格式的性能。
• 证明在高分辨率模拟中，时间误差可能主导空间误差，因此高阶时间积分至关重要。

提出的方法

• 方法从JST算法（一种低存储Runge-Kutta格式）出发，通过引入由精确解的泰勒展开导出的校正项进行改进。
• 校正项专门针对右端项形式为F(u) = L(u) + N(u,u)的系统推导，其中L为线性算子，N为关于其变量对称的二次非线性算子。
• 校正项的施加方式保持了JST算法原有的两级存储结构，从而确保极低的内存占用。
• 该格式在伪谱法和谱元法代码中实现，仅需对现有时间推进框架进行最小修改。
• 该方法应用于Hénon-Heiles系统、一维和二维Burgers方程，以及使用完全并行化伪谱代码求解的三维不可压Euler方程。
• 3+格式（包含三阶校正项）被识别为最具成本效益的实现方式，仅需四次JST迭代。

实验结果

研究问题

• RQ1 能否在不增加存储需求的前提下，对JST算法进行校正，使其在二次非线性系统中实现三阶和四阶精度？
• RQ2 为使JST算法在非线性二次算子下超越二阶精度，所需的特定校正项是什么？
• RQ3 与标准低存储Runge-Kutta方法相比，校正后格式在高分辨率模拟中的精度和稳定性表现如何？
• RQ4 新格式提升的时间精度在大规模模拟（如湍流的直接数值模拟）中是否具有实际效益？
• RQ5 校正后的格式能否在现有高阶谱法和谱元法代码中高效实现，且仅需极少的代码修改？

主要发现

• JST算法虽对线性问题精确，但在处理二次非线性系统时引入三阶误差，因此必须引入校正项以实现更高阶精度。
• 所推导的校正项在保持JST格式原有两级存储结构的同时，成功恢复了二次非线性问题的三阶和四阶精度。
• 3+格式（包含三阶校正项）被证明最具成本效益，仅需四次JST迭代，且全局误差显著低于标准三阶格式。
• Burgers方程的数值结果表明，随着格式阶数提高，截断误差减小，且3+格式优于标准三阶方法。
• 对于三维不可压Euler方程，新格式可在高分辨率模拟中实现精确的长时间积分，避免了低阶时间积分导致的虚假守恒性质。
• 该方法在大规模模拟中具有实际可行性和高效性，包括三维伪谱代码，适用于千万亿次级（petascale）湍流直接数值模拟。