论文部分内容阅读
在计算流体力学中，对同时包含间断和复杂流动现象的流场进行数值模拟是一个迫切的并且也是非常困难的任务。近年来，针对此类问题的捕捉激波的高精度高分辨率计算方法得到了迅速的发展。而在流体力学数值算法中，有限差分方法历史最悠久也最为成熟，并且适宜于构造高精度的格式。本论文主要研究流体力学高精度高分辨率差分格式，对格式的精度和稳定性进行了分析和讨论。论文内容主要包含以下几个方面：1. WENO格式的精度以及光滑因子的研究给出了一个方便使用的五点WENO格式达到最高精度的充分条件，以及一般情况下的WENO格式的精度公式。另外，对Jiang和Shu的WENO格式的光滑因子进行了分析，给出了两个低耗散的光滑因子，并分析了Liu等人的WENO格式的光滑因子，通过Borges等人的WENO权因子计算方法保证收敛精度。数值算例表明，使用这些光滑因子的WENO格式对连续流场的捕捉能力有所提高。2.紧致格式及紧致WENO混合格式的研究构造了一个使用对称网格模板的五阶守恒型迎风紧致格式，与经典的五阶迎风紧致格式相比，该格式耗散较低，且拥有更大的稳定性范围，并且给出了一种显式的保证整体精度为五阶的边界格式。另外，对Ren等人的紧致WENO混合格式做出了改进，对欧拉方程采用了新的处理，从而避免了块三对角方程组的求解，大大降低了计算量。3.简化了半离散差分格式稳定性条件傅立叶分析方法得到的稳定性条件一般为余弦函数的不等式，不易于使用。本文将余弦函数转化为多项式形式，同时证明了当差分格式每增加两阶收敛精度，该多项式就可以提取一个因子，这使得该多项式不等式的求解大为简化。该方法对判断半离散差分格式的稳定性十分有效。而在构造高精度优化格式时，该方法得到的稳定性条件可以表示为差分格式系数的有限个约束的形式，这使得在优化过程中加入稳定性限制变得可行。4.优化型格式的研究通过将一个给定模板上的最高阶精度格式与一种优化格式进行加权平均，权系数由本文提出的局部波数指示子确定，得到了一种保精度优化（Maximum order preserving optimized,MOPO）格式。然后将这种格式与六点对称WENO格式融合，得到了一种MOPOWENO格式。数值结果表明，MOPO格式与MOPOWENO格式不仅保持了最高阶精度，并且在所测试问题中比相应的两种子格式的相位误差更小。