[论文解读] Strong Stability Preserving Multistep Runge-Kutta Methods

Christopher Bresten, Sigal Gottlieb+4 more|arXiv (Cornell University)|Jul 30, 2013

Computational Fluid Dynamics and Aerodynamics参考文献 25被引用 12

一句话总结

本文针对具有单调性和稳定性特性的双曲PDE，开发了高阶显式强稳定保持（SSP）多步Runge-Kutta方法。推导了阶次与单调性条件，证明了二阶及以上方法SSP系数的上界，并通过数值优化构造了最高达10阶、8级、5步的方法，其有效SSP系数显著高于标准方法。

ABSTRACT

High-order spatial discretizations with strong stability properties (such as monotonicity) are desirable for the solution of hyperbolic PDEs. Methods may be compared in terms of the strong stability preserving (SSP) time-step. We prove an upper bound on the SSP coefficient of explicit multistep Runge--Kutta methods of order two and above. Order conditions and monotonicity conditions for such methods are worked out in terms of the method coefficients. Numerical optimization is used to find optimized explicit methods of up to five steps, eight stages, and tenth order. These methods are tested on the advection and Buckley-Leverett equations, and the results for the observed total variation diminishing and positivity preserving time-step are presented.

研究动机与目标

• 开发高阶显式多步Runge-Kutta方法，以保持强稳定性特性，如总变差递减（TVD）和保正性。
• 基于方法系数，推导多步Runge-Kutta方法所需的阶次与单调性条件。
• 建立显式多步Runge-Kutta方法在二阶及以上的SSP系数的上界。
• 利用数值算法对最高达10阶、8级、5步的方法进行优化，以最大化有效SSP系数。
• 通过迎风方程和Buckley-Leverett方程的数值实验，验证理论上的SSP特性。

提出的方法

• 利用方法系数，推导显式多步Runge-Kutta方法的阶次条件与单调性条件。
• 将强稳定保持（SSP）特性表述为前向欧拉步骤的凸组合，确保在减小的时间步长下保持稳定。
• 将SSP系数C定义为相对于前向欧拉步骤的最大允许时间步长，有效SSP系数定义为C/s，其中s为步骤数。
• 通过数值优化，寻找在2至10阶、最多5步和8级下具有最大有效SSP系数的方法。
• 使用保守的空间离散化方法，结合Koren限制器处理Buckley-Leverett方程，采用半拉格朗日格式处理迎风方程。
• 将观测到的TVD与保正性保持的时间步长与理论预测值进行比较，以验证SSP特性。

实验结果

研究问题

• RQ1 显式多步Runge-Kutta方法在二阶及以上的理论SSP系数上界是什么？
• RQ2 如何系统地推导多步Runge-Kutta方法的阶次与单调性条件？
• RQ3 在有限级数与步数限制下，高阶显式多步Runge-Kutta方法可达到的最大有效SSP系数是多少？
• RQ4 在数值实验中，观测到的TVD与保正性保持的时间步长与理论预测值相比如何？
• RQ5 优化后的多步Runge-Kutta方法在双曲PDE的稳定性与效率方面是否优于标准Runge-Kutta方法？

主要发现

• 证明了显式多步Runge-Kutta方法在二阶及以上的SSP系数存在上界，限制了在凸组合框架下可达到的最大稳定性。
• 对于5步8级的10阶优化方法，有效SSP系数达到0.603，显著高于非SSP的RK4,4方法的0.250。
• 在迎风方程中，优化后的(7,3,3)方法观测到的TVD时间步长为理论值的1.113倍，表明其稳定保持性能超出理论界限。
• 在Buckley-Leverett问题中，观测到的保正性保持时间步长始终大于TVD时间步长，表明所选稳定性特性会影响实际观测到的稳定极限。
• 对于(20,3,10)方法，观测到的TVD时间步长为0.356Δx，对应有效SSP系数为0.110，表明增加步数可提升高阶方法的稳定性。
• 观测到的TVD与保正性保持时间步长随级数增加而增大，证实了理论预期：更多级数可增强稳定性特性。