激波捕捉格式与对称性问题.pdfVIP

对称性与激波捕捉问题1 冉 政 (上海大学，上海市应用数学和力学研究所，上海 200072) 摘 要: 差分格式的群不变分析结果与 Lax-Wendroof 格式的数值模拟结果有一种有趣的 对应，即在激波上游 Lax-Wendroff 格式的数值模拟结果出现数值振荡，而在这一区域该格 式并不满足群不变条件；激波下游，Lax-Wendroff 格式数值模拟结果是无振荡的，而此时 的群不变条件是满足的。分析表明：Lax-Wendroff 格式计算无粘 Burgers 方程的问题时， 除了前两个（即：空间平移，时间平移变换）对称条件满足外，无粘 Burgers 方程所具有的 另外两个对称一般无法满足，对称性的丧失可能是导致数值解出现三种异常现象的内在机制。 关键词: Lax-Wendroff 格式， 对称不变性， 激波振荡 PACC: 034G,0220 引 言 流场中激波的数值模拟一直是计算流体力学所研究的重要问题之一。目前使用得最广的 是激波捕捉法。最初人们采用一阶精度格式捕捉激波，可得到过激波的单调解，但因耗散过 大，激波的梯度被抹平，其物理特性失真。二阶精度差分格式所给出的数值模拟解在激波附 近产生非物理振荡，且可能出现非物理弱解。这一直是学术界关心的问题[1,2]。20 世纪 80 年代 TVD 格式的提出使得激波捕捉法有了重大发展，基本上解决了高阶精度差分解在激波附 近的非物理波动问题。在中国[3]，很多学者也深入研究了激波的数值模拟问题，提出了很好 的计算方法，如 NND 格式，MmB 格式，紧致迎风格式，耗散比拟格式，耗散守恒格式等等。 且对 TVD 格式，ENO 格式进行了深入的研究，给出了改进型的 TVD 格式和 ENO 格式。这些方 法都以高分辨率的性能模拟了激波，使激波捕捉法有了新的进展 。目前该问题的研究仍处 于高潮，还存在着有待进一步解决的问题。 差分格式的稳定性分析的本质困难在于方程的非线性，并且解可以具有间断。传统的分 析方法，采用了局部冻结系数的线性化方法，这种线性化方法得到的稳定性和耗散性条件是 必要的但并不充分。大量的数值实验和理论结果表明，由上述方法得到的格式，用于守恒型 问题计算时，仍可能发生如下异常现象：激波振荡；非线性不稳定；非物理解现象。上述的 异常现象是目前激波捕捉法计算中存在的三个主要问题[4,5]。因此，有必要从理论上搞清这 些异常现象产生的原因和条件，找到更为精确的分析方法，才可能构造更具有普遍意义的格 式。在文献[6，7]中，已就一些理论方面问题进行了论述，本文利用 Lie 变换群理论，对有 关著名的 Lax-Wendroff 格式进行非线性行为分析，以期对 Lax-Wendroff 格式所得到的波后 振荡结论提出新的见解。 1 Lax-Wendroff 格式的构造[8] Lax-Wendroff 格式是一种典型的二阶格式。由于它影响，在 60 年代曾经产生过大量的二阶 精度格式，在实际计算中发挥了重要作用。即使是新近的一些研究工作，也会不时地提及 Lax-Wendroff 格式。 对于一维拟线性问题 1 国家自然科学基金资助，项目号 1 ut + f (u) x = 0 （1.1） ?f 引入 a(u) = ,为 f (u) 的 Jacobi 矩阵，由此可以得拟线性问题的 La