求解多维Euler方程的二阶旋转混合型格式-应用数学和力学.PDFVIP

　 ， ３５ ５ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１４年５月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３５，Ｎｏ．５，Ｍａｙ．１５，２０１４ 文章编号：１０００⁃０８８７（２０１４）０５⁃０５４２⁃１２ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 求解多维 Ｅｕｌｅｒ方程的二阶旋转混合型格式∗ 刘友琼，　 封建湖，　 任　 炯，　 龚承启 （长安大学 理学院，西安 ７１００６４） 摘要：　 提出了一个基于旋转Ｒｉｅｍａｎｎ求解器的二阶精度的Ｅｕｌｅｒ（欧拉）通量函数．不同于“网格相 关”的有限体积方法或者维数分裂的有限差分方法，本格式是基于旋转 Ｒｉｅｍａｎｎ求解器将 ＨＬＬＣ 格式与ＨＬＬ格式进行特定结合而得到的一类混合型数值格式．在激波法向采用ＨＬＬ格式从而抑制 红斑现象，在激波方向采用ＨＬＬＣ格式从而避免产生过多的耗散．新的旋转混合型格式具有结构简 单、无红斑、高分辨率等优点．数值算例充分说明了新格式消除Ｅｕｌｅｒ方程激波不稳定现象的有效 性和鲁棒性． 关　 键　 词：　 激波不稳定；　 Ｅｕｌｅｒ方程；　 无红斑；　 旋转混合型格式 中图分类号：　 Ｏ３５４；Ｏ２４１．８２　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１４．０５．００８ 引　 　 言 长期以来，双曲型守恒律方程的计算取得了巨大的进步．以通量差分分裂格式（ＦＤＳ）和 Ｇｏｄｕｎｏｖ格式为代表的流体力学数值方法在双曲守恒律方程的求解中占据了很重要的地位， 这些格式在大多数情况下都表现出很好的鲁棒性．然而，当计算二维问题时，它们还是会出现 ［１］ 一些非物理现象，比如说“红斑（ｃａｒｂｕｎｃｌｅ）”和“奇偶失联（ｏｄｄ⁃ｅｖｅｎ ｄｅｃｏｕｐｌｉｎｇ）”两个现象 ， ［２］ ［３］ 即使是常用的Ｒｏｅ近似Ｒｉｅｍａｎｎ（黎曼）求解器 、ＨＬＬＣ 近似Ｒｉｅｍａｎｎ求解器 也出现了不 同程度的红斑现象．在过去的十几年中，人们进行了一系列研究试图避免非物理现象的产生， 详见文献［４⁃１０］．虽然他们作了详细的解释并提出了不同的格式，但这些格式都没有完全消除 此类现象．众所周知，一条波的Ｒｕｓａｎｏｖ 通量函数及两条波的ＨＬＬ 通量函数都是“无红斑” ［１１］ 的 ．这些通量函数对于涉及激波的非粘性计算是非常稳健的，但对于剪切层特别是 Ｎ⁃Ｓ （Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ）方程边界层，计算时会产生较大耗散． ［１１］ Ｎｉｓｈｉｋａｗａ等提出了将无红斑通量函数与Ｒｏｅ通量函数相结合的方法 ，并利用旋转Ｒｉ⁃ ｅｍａｎｎ求解器［１２⁃１３］将这两个通量函数相结合以达到自动选择迎风方向的目的．在此之前， ［１３］ Ｒｅｎ 已经采取了旋转通量，获得了对激波捕捉的稳健性，但并没有提高精度．受Ｒｅｎ及Ｎｉｓｈ⁃ ｉｋａｗａ等的启示，本文将单元交界面的法向量分解成两个正交方向，一个是垂直于激波的速度 ∗ 收稿日期：　 ２０１３⁃０６⁃０８；修订日期：　 ２０１４⁃０４⁃０３ 基金项目：　 国家自然科学基金；长安大学中央高校基本科研业务费项目 （ＣＨＤ２１０２ＴＤ０１５） 作者简介：　 刘友琼（１９８９—），女，云南人，硕士生（Ｅ⁃ｍａｉｌ：ｙｏｕｑｉｏｎｇｌｉｕ＠１６３．ｃｏｍ）； 封建湖（１９６０—），男，陕西人，教授，博士，博士生导师（通讯作者． Ｅ⁃ｍａｉｌ ：ｊｈｆｅｎｇ＠ｃｈｄ． ｅｄｕ．ｃｎ）．