差分格式的群不变特性和激波捕捉问题的研究.pdfVIP

差分格式的群不变特性与 激波捕捉问题的研究1 冉政，胡文沛 (国家计算流体力学实验室，北京100083) 摘要 结合非线性分析中的无穷小Lie群分析方法．发展典型激波捕捉格式的群不变分析方法，以 期对目前激波捕捉法所面I临的一些挑战问题，进行深入系统的研究，揭示一些新的差分格式内在的非线 性本质特征：在此基础上建立一些新的差分格式非线性稳定性原则，以及新的基本无波动的激波捕捉 差分计算方法。 关键词差分格式；变换群；波动解 引言 气体动力学中的激波捕捉问题研究一直是计算流体力学的核心问题[1-9]，尽管经过 四十多年的努力，但目前仍然有许多有待解决的问题[8-14]，这构成了对计算流体力学学 科的挑战。发展更为精确的分析方法，建构更具有普遍性的算法，己成为学术界的关注 点。计算含有激波的流动问题是计算方法一直面临的挑战，主要因为激波的流动问题是强 非线性的，计算方法要求尽可能地体现激波的流动问题的物理本质。当20世纪80年代高 分辨率差分格式研究趋于成熟时．便有学者指出[8]，计算流体力学已有的方法已经能够 解决所有问题。但是，一系列的研究表明，已有的计算方法计算含有激波的流动问题是还 面临一系列无法解决的困难，可以列举下面一些问题： 慢激波穿越网格界面的守恒问题：动量峰值问题：激波后的数值振荡问题：膨胀波问题：负内 能问题：红斑问题：壁面过热问题：马赫杆扭曲问题：奇偶网格点失耦问题：启动误差问题． 以上大部分问题在Yee[11]和Quirk[10]著名的文献中已有讨论，之后有些问题有了一些 新的进展(例见[123，[13]，[14])，但离彻底的理解还有一段距离。 与此同时，有限元中的激波捕捉法，较有限差分方法更不完善(例见[15]，[16])， 也存在一些问题和困难，例如： ●由于它离散复杂，不易采用有限差分方法中关于处理间断的高分辨的思想，如迎风 通量分裂和总变差衰减等； ●在有限差分方法中行之有效的多维分裂技术也不易移植用于有限元法； ●对存在激波等强间断的流场的分辨率较差。 因此，研究激波捕捉法中差分计算更为本质的属性，发展更为精确的分析方法，成为学术 焦点。在这种背景下，我们认为：目前计算流体力学的大多数理论结论都缘于线性化的分 析结果，并多为一维问题，这种传统的分析方法所得到的稳定性和耗散性条件一般是必要 国家自然科学基金，国防预研基金资助项目 4 的但并不充分。从根本上讲，上述的大多数挑战问题可能都与问题的非线性本质有关．因 此，从理论上探讨激波捕捉法中差分格式的非线性行为。将有助于建构更具普遍意义的算 法。就目前的学科发展来看．非线性数学也已为我们提供了可以使用的工具(例见[171， f181，[19】)，因而使得本项研究成为了可能。 同时．我们也受益于计算数学发展的启发(特别是冯康先生发展辛算法思想的启 发)。冯康先生曾指,m[20l：当代计算方法研究的一条不成文的基本法则是，问题原型的 基本特征在离散后应尽可能地得到满足，而为了达到这一效果则离散化应尽可能在问题原 型的同一形式框架中进行。根据这一思想，激波捕捉法中差分格式除满足传统的稳定性和 耗散性原则外，一定的对称性要求是必要的，而对称性的分析方法是无穷小Lie交换群。 因而具有一定的可操作性。 非线性分析方法中，Lie变换群的方法具有独特意义．有关Lie变换群的基本知识可参 见文献[17．191，Lie群很有可能成为数值解问题强有力的分析工具则部分地可以从下歹0一 系列工作中反映出： ●1974年。Yanenko和Shokin首先提出了差分格式的群分类问题【21】； ●1983年，Shokin将对称群方法用于了数值分析【22】； ●1987年，Maeda Shigeru研究了Lie变换群方法对差分方程的推r[23b ●1991年，Levi和Wintemitz讨论了离散方程的连续对称性问题[24】； ●1992年，Gear讨论了群不变性与常微分方程的数值积分问题[25】： 12000年，我们首先得到了Lax-Wendrofr格式的群不变分析结果与激波振荡问题的理论 五个对称中的两个，给出了一种有关Lax-Wenoff格式出现激波振荡问题的理论说明 【26】； 12001年，我们首先得到了TVD格式的群不变分析结果。对一