MFCAV近似Riemann解法器在相容拉氏方法中的熵条件分析.pdfVIP

2015年 8月 计 算 数 学 第 37卷第 3期 Aug．，2015 M ATHEMATICA NUMERICA SINICA Vo1．37，No．3 MFCAV近似 Riemann解法器在相容拉氏 方法中的熵条件分析木1) 刘 妍 田保林 申卫东 (北京应用物理与计算数学研究所，计算物理重点实验室，北京100094) 茅德康 (上海大学数学系，上海200436) 摘 要 对基于MFCAV(MultiFluidChannelonAveragedVolume)近似 Riemann解法器的相 容拉氏方法的熵条件进行了分析．结果表明与满足声学形式 Riemann解法器的熵不同，前者只 能在每个网格边界左、右两侧网格的熵随时间变化的和保证大于零，即能保证整体熵增，但不保 证传统意义上的在每个网格中的熵增；而后者不仅保证整体熵增，而且还满足传统意义上的熵增． 因此 MFCAV 的熵增相对声学形式解法器而言要弱一些，由此表明其熵增可能要小些，使得格式 的耗散可能要小些．数值算例也验证了分析的正确性． 关键词：熵条件；MFCAV近似 Riemann解法器；相容拉氏方法 MR (2000)主题分类：35L65 1．引 言 拉氏方法 (Lagrangian)是可压缩多介质流动的一种重要数值方法．该类方法的计算网格 随流体运动，数值耗散小，界面清晰，且易于处理 内爆压缩问题中的自由面，在天体物理、惯 性约束聚变 (ICF)等领域的数值模拟中是主要使用的计算方法．但 由于网格跟随流体运动， 因此会受到因网格虚假变形而导致计算中止的困扰 [1-4J，难以处理大变形流动问题．多年来， 如何抑制网格虚假变形，并从根本上提高拉 氏方法对流场数值模拟的分辨率和置信度，是一 直受到关注的研究内容．如 2【]中分析了多维中心型拉氏方法产生网格虚假变形的原因，并提 出了一种求解涡量输运方程的办法来抑制非物理变形 ．2jj但该方法过于复杂，据作者所知 目 前尚未见到推广应用．P．H．Maire等人提出了一种相容有限体积中心型拉氏方法 }6，J．与传 统中心型拉氏方法不同的是算法满足几何守恒律 (GCL)，即保证角点速度计算与网格边上数 值流相容，大大提高了网格的抗变形能力．该方法跨过各网格边界的数值流并不平衡，而是各 网格角点周围的数值流平衡，且算法满足熵条件．该方法的一个缺陷是只适应于声学形式的 Riemann解法器，任意形式的Riemann解法器并不能直接应用于其中．在 [9]中我们提出了一 种方式，成功地将不具有声学形式的体平均多流管MFCAV(MultiFluidChannelonAveraged Volume)Riemann解法器 IS】应用到相容方法 。【J中．我们之所以考虑选择MFCAV解法器，是 因为该解法器与 f6]中的声学形式解法器相 比考虑 了力学平衡弛豫时间的影响L8J．对纯流体 2015年 1月 16日收到． )基金项目： 国家自然科学基金 11472059、中国工程物理研究院联合基金 、 中国工程物理研究院科学技术发展基金 (2012A0202010)和上海市重点科学建设项 目(J5OlO1)． 3期 刘妍 等：MFCAV近似 Riemann解法器在相容拉氏方法中的熵条件分析 287 力学问题，由于可以忽略达到力学平衡所需的弛豫时间，则 MFCAV 中的考虑力学平衡弛豫 时间的处理方式对计算结果的影响不大，但是所计算的问题中还包含有其它物理过程，如化学 反应过程，由于它的弛豫时间比力学平衡弛豫时间小的多，则用 [6]中的声学形式解法器计算 时，就有可能得到定性上不准确的计算结果．因此 MFCAV在很多实际问题，特别是许多内爆 问题的数值模拟中发挥了很大作用 _5j． 熵耗散对应格式的数值粘