[论文解读] A general approach to enhance slope limiters on non-uniform grids

Xianyi Zeng|arXiv (Cornell University)|Jan 6, 2013

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

一句话总结

本文提出了一种通用框架，通过推导保证二阶精度、TVD稳定性及对称性保持的充分条件，以增强非均匀网格上高分辨率格式（特别是MUSCL方法）的斜率限制函数。该方法将经典重构与Harten的稳定性理论扩展至非均匀一维网格，并通过一维和二维数值算例验证，恢复了传统限制器在非规则网格上损失的二阶收敛率。

ABSTRACT

Abstract. A general approach to study and enhance the slope limiter functions on non-uniform grids is presented. Slope limiters are preferred in high-resolutions schemes in general and MUSCL in particular to solve hyperbolic conservation laws. However, most 1D limiters are developed assuming uniform meshes in space, which are shown to be inadequate on non-uniform grids. Especially, second-order convergence is shown to be lost when the conventional limiters are applied on irregular grids in the case of smooth solutions. A methodology based on the classical reconstruct-evolve-project approach and Harten’s stability theory is presented to study the slope limiters on 1D non-uniform computational grids. Sufficient conditions for the limiters to lead to formal second-order spatial accuracy, total-variational-diminishing stability and symmetry-preserving property are derived. The analysis and results extend naturally to cell-centered finite volume methods in multiple dimensions. Several most widely used conventional limiters are enhanced to satisfy these conditions, and their performances are illustrated by various 1D and 2D numerical examples.

研究动机与目标

• 解决在高分辨率格式中，将传统斜率限制器应用于非均匀网格时导致的二阶空间精度损失问题。
• 基于重构-演化-投影框架与Harten的稳定性理论，为非均匀一维网格上的斜率限制器分析，开发系统化的方法论。
• 推导确保非均匀离散化中形式二阶精度、总变差递减（TVD）稳定性及对称性保持性质的充分条件。
• 将分析扩展至多维的单元中心有限体积方法。
• 增强广泛使用的传统限制器，使其满足推导出的条件，并通过数值方法验证其性能。

提出的方法

• 采用重构-演化-投影方法作为分析非均匀网格上斜率限制器的基础。
• 应用Harten的稳定性理论，推导非均匀网格背景下TVD稳定性与对称性保持的充分条件。
• 制定限制函数的数学约束，以确保在非均匀网格上实现形式二阶空间精度。
• 通过修改其限制函数，将传统限制器（如minmod、van Leer、Superbee）进行变换，使其满足推导出的条件。
• 在一维和二维有限体积格式中实现增强后的限制器，以在不规则网格上测试其精度与稳定性。
• 使用光滑解量化收敛率，并验证二阶精度的恢复。

实验结果

研究问题

• RQ1 斜率限制器函数在非均匀网格上的一维双曲守恒律中，必须满足何种条件才能保持二阶空间精度？
• RQ2 Harten的稳定性理论如何适配用于分析非均匀计算网格上的斜率限制器？
• RQ3 对传统限制器需进行何种修改，才能确保在非均匀网格上实现TVD稳定性与对称性保持？
• RQ4 增强后的限制器在不规则网格上对光滑解的二阶收敛性恢复程度如何？
• RQ5 所提出的框架能否扩展至多维单元中心有限体积方法？

主要发现

• 即使在光滑解情况下，传统斜率限制器在应用于非均匀网格时仍会损失二阶空间精度。
• 推导出的充分条件可确保非均匀一维网格上斜率限制器的形式二阶精度、TVD稳定性及对称性保持。
• 广泛使用的限制器（如minmod、van Leer、Superbee）的增强版本在非均匀网格上成功恢复了二阶收敛性。
• 该方法自然推广至多维单元中心有限体积方法。
• 一维与二维数值算例证实，增强后的限制器在保持高分辨率特性的同时，可在不规则网格上实现最优精度。
• 该方法为实际计算流体动力学应用中非均匀网格上设计鲁棒且精确的限制器提供了一个系统化框架。