双曲守恒律方程降阶Spectral Volume方法研究.docVIP

双曲守恒律方程降阶Spectral Volume方法研究 　　摘 要：文章中提出了降阶的Spectral Volume 方法 （SV方法） 的概念，对方法的公式进行推导，并成功将其应用与一维线性传递方程。对降阶的SV方法的稳定性进行了Fourier分析，并与传统的SV方法进行了对比。实验证明，降阶的SV方法在初值连续情况下，在线性传递方程上可以达到所期望的数值精度，并且在数值色散和数值耗散特性上较传统SV方法有显著提高。 　　关键词：SV方法 一维守恒律 降阶 　　中途分类号：R445 文献标识码：A 文章编号：1674-098X（2014）10（b）-0211-04 　　1 研究背景 　　计算流体力学（CFD）从创立至今得到了飞速发展，在工业界起着不可替代的作用。然而目前工业界所用到的CFD计算通常只能达到二阶精度，远远达不到人们的需要，并且对网格以及内存和CPU的要求超出了当今计算机所能承受的范围。另外，运算区域通常包含复杂单光滑的流体结构，比如涡旋与剪切流层的相互作用区，以及同时包含激波与其他复杂流动的区域等。因此，高精度高分辨率格式的研究一直是CFD领域中非常活跃的一个课题，发展高精度高分辨率的计算格式成为CFD工作者的一大研究方向。为采用高精度格式不但可以降低对网格规模的苛刻要求，而且能够正确分辨其中复杂的流动现象。对网格点数目的要求低，具有减少计算机内存需求和CPU时间消耗的优势。其中比较著名的是Harten提出的高阶精度的ENO（essentially non-oscillatory） 格式，随后Liu，Osher和Chen发展了原有的ENO格式，提出了有显著改进的WENO（weightedessentially non-oscillatory）格式，并获得了迅速的发展，并得到了一系列的研究。ENO/WENO格式善于捕捉物理间断，其应用迅速发展到多相流领域。 　　新型格式发展的目标之一是发展可以应用于非结构网格的可靠的高阶格式。理论和实践都表明，大多数高精度、高分辨率格式的性能都严重地依赖于网格的光滑性，而工业界所需要模拟的物体通常具有复杂的几何外形（图1），CFD应用最多的航空领域所需要的几何外形更是相当复杂，而这给结构网格的生成带来很大困难，一种直接有效的解决办法是采用非结构网格。然而由于稳定性等因素的影响，对于三维有限体积法，非???构算法模板所需单元往往较多，而过多的模板单元会带来边界处理、内存占用、编程复杂等多方面的困难。因此实际应用中大都只使用二阶精度的非结构算法，而这远远满足不了一些工业领域对计算精度的要求。可以应用与非结构网格的高精度、高分辨率、低耗散的实用数值格式的发展成为当今CFD的重要任务之一。以此为目标的一些高阶格式已经发明出来，其中比较著名的是Cockburn发展的间断Galerkin有限元方法（DG），间断 Galerkin有限元方法的插值模板小，各单元之间仅通过界面通量计算相联系，因而比较适合在复杂外形。而另一种成功的格式是Wang在自2002年发展的应用于非结构网格的算法SV方法。SV方法是一种高精度高分辨率低耗散的数值格式。同间断Galerkin有限元方法一样，SV方法可以达到相同的精度。然而相比之下，SV方法比间断Galerkin有限元方法有着更好的稳定性（可以允许更大的CFL数）。与传统的有限体积法相比，对计算机内存和CPU的要求更低。所以SV方法有着非常大的能在将来被应用于工业界的潜力。同时，SV方法的稳定性和可行性还需要进一步验证和改进。 　　2 SV方法介绍 　　SV方法，是比较新提出的另一种可以应用于完全非结构网格的高精度、高分辨率、低耗散的守恒型格式。Wang于2002年提出一维SV方法的概念并随后将其发展到二维和三维。SV方法是经典有限体积法的延伸。它提出了控制体积（Control Volume）的概念，从而避免的复杂模板（stencil）的使用，从而使格式可以轻松应用与非结构网格上。经典的有限体积法中，变量在每个网格中的值被用来进行多项式插值。为了达到相应的精度，由相邻多个单元格组成的模板被用来进行插值。而非结构网格使得模板中相邻网格的选取成为一大难题。不同于经典的有限体积法，SV方法把网格中的每一个单元格进一步分成若干控制体积（CV），对每个控制体积中变量的平均值加以利用来实现多项式插值从而达到所要求的精度。SV方法的优势在于其在保证相同精度的同时免除了对周围网格的依赖，从而降低了对网格的要求。此外，由于不同单元格之间的信息交流少，SV方法非常适宜并行运算的实现。 　　有限体积法（包括ENO和WENO方法）利用相邻多个自由度（模板）内状态变量的值进行多项式构建。而在SV方法中，我们把这个单个的单元格叫做Spectral Volume（SV）