强健一维守恒型间断跟踪法.pdfVIP

摘 要 本文进一步分析、推广和提高文[29]中提出的守恒型间断跟踪法。该跟踪法是 条件。在[29]中，当将该方法运用到方程组的情况时只有一阶精度。这是因为在方程 组的情况，在被跟踪的间断处会有其它特征域的波，并且它们可能会穿过被跟踪的 间断，这会影响算法的高精度实现。文[29]中所设计的技巧只能保证一阶精度。在本 文中，我们设计了该算法在一维方程组情况的一种高精度实现，即通过合适地修正 间断两侧的数值解来将间断一侧的其它特征域的波传到被跟踪间断的另一侧，以此 来实现被跟踪间断线上的移动边界条件。我们通过在间断处解Riemannf可题得到那 些需要穿过被跟踪间断的其它特征域的波，并且采用了高阶的重构、高阶Lagrangc 插值和高阶的数值积分以及有限制的波分解等，来修正间断两侧的数值解。我们还 对算法作了严格的局部截断误差分析，证明了该算法是高阶精度的。 我们还讨论了间断相互作用可通过求解Riemann问题来实现，并讨论了1Rie— rilailrl问题分解出来的需要被跟踪的间断所对应的网格上的数值解如何来确定。我 们还设计了一些规则来处理多个间断的情况，把间断的移动和相互作用的无穷多种 情况进行了有限化，使得我们的算法可以做成一个强键的算法，它可以处理任意多 个间断的移动和相互作用。 最后，我们将此算法应用到一维的Euler方程组上，其中处理了激波在固壁上 的反射，自发生激波，推广到多介质流体等，并对算法用Fortran90语言按面向对象 的编程方式进行了程序实现。从而完成了一个几乎可以处理任何情况的和捕捉法几 乎一样强健的一维跟踪法。 关键词： 守恒律，间断跟踪法，高阶精度，间断网格 Abstract Inthis further and the paper．we conservativefront— analyze．improvedevelop method in themethod al*e tracking developed tracked 129]．In in【29]，discontinuities the conservationof PDEs by the ratherthan theRankine- enforcing properties by using condition．Themethodin firstorder when Hugoniot only accurate to 129]is applied ofconservationlaws．Thisis inthe because neartraked systems ease disconti． system nuitiesthererewa-∞tsinothercharacteristicfields and 8∞：rossthe