[论文解读] Reminiscences about numerical schemes

С. К. Годунов|ArXiv.org|Oct 3, 2008

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

一句话总结

本文呈现了S. K. Godunov对其在1953至1969年间发展双曲守恒律系统Godunov格式的个人反思。内容详述了从带人工粘性的的一阶格式到高阶方法、网格自适应的演变过程，强调了数值粘性在模拟激波间断中的关键作用，最终形成满足熵条件的广义解，并认识到方案依赖的数值耗散效应。

ABSTRACT

This preprint appeared firstly in Russian in 1997. Some truncated versions of this preprint were published in English and French, here a fully translated version is presented. The translation in English was done by O. V. Feodoritova and V. Deledicque to whom I express my gratitude. In the present paper I will describe how the first variant of the Godunov's scheme has been elaborated in 1953-1954 and tell about all modifications realized by myself (until 1969) and the group of scientists from the Institute of Applied Mathematics in Moscow (which has become the M.V.Keldysh Institute of Applied Mathematics).

研究动机与目标

• 重建莫斯科应用数学研究所初创时期Godunov数值格式在双曲守恒律领域的发展历史。
• 解释在激波波前附近稳定二阶格式所面临的挑战，以及由此导致向一阶精度、保守性方法的转变。
• 探讨数值粘性在近似物理耗散中的作用及其对离散化方式的依赖性，进而认识到不同格式对同一问题会产生不同解。
• 建立数值格式、广义解与热力学熵条件之间的联系，尤其通过行波解和小粘性极限进行分析。
• 强调物理学家与数学家之间的跨学科讨论在完善计算流体动力学的理论与计算基础中的影响。

提出的方法

• 该格式通过修改A. I. Zhukov的现有算法而发展而来，后者使用一阶有限差分法，利用截断误差引起的数值粘性来稳定激波轮廓。
• Godunov引入了欧拉坐标系下的保守有限体积格式，采用精确黎曼求解器计算单元界面的通量，后被形式化为Godunov格式。
• 该方法采用一阶微分近似概念分析数值粘性，将截断误差与差分格式中的人工耗散联系起来。
• 采用形式为 $ q_i = q_i((x - \rho t)/\varepsilon) $ 的行波解分析小粘性下的解行为，揭示了解对耗散项形式的依赖性。
• A. V. Zabrodin与G. P. Prokopov提出的移动网格构造与区域分解技术，使得求解具有移动间断的复杂多尺度问题成为可能。
• 对网格生成中出现的椭圆系统进行理论分析，推动了数值线性代数的发展，并提出了矩阵束的谱二分性准则。

实验结果

研究问题

• RQ1 为何二阶格式在激波附近无法产生稳定解？其不稳定性背后的物理或数学原理是什么？
• RQ2 有限差分格式中的数值粘性如何与物理耗散相关联？能否用于近似双曲守恒律的广义（满足熵条件）解？
• RQ3 不同数值格式对同一双曲系统是否会产生不同解？当粘性趋于零时，极限解的唯一性由什么决定？
• RQ4 广义解的结构如何被严格定义并加以证明，特别是在不连续流与熵条件的背景下？
• RQ5 物理直觉以及实验物理学家的批判性反馈在改进数值算法及其理论基础中起到何种作用？

主要发现

• 基于精确黎曼求解器与保守通量的一阶Godunov格式在激波捕捉中表现出鲁棒与稳定，优于产生非物理振荡的二阶格式。
• 有限差分格式中的数值粘性被证明等价于引入小的人工耗散，但所得解依赖于耗散项的形式，而不仅取决于粘性趋于零的极限。
• B. F. Diachenko构造的实例表明，不同数值格式可能对同一双曲系统产生不同解，从而否定了极限解与格式无关的假设。
• 行波解的使用揭示了只有特定形式的粘性才能产生满足熵条件的不连续解，建立了数值格式与热力学原理之间的联系。
• 移动网格与区域分解技术的发展使得模拟具有移动前导波与多重间断的复杂流动成为可能，显著扩展了该方法的应用范围。
• 对网格生成的理论研究促成了谱二分性准则的提出，以及奇异向量计算的高精度算法，对数值线性代数产生了更广泛的影响。