论文部分内容阅读
本文为数值求解依赖时间的偏微分方程提出两类基于特征思想的高分辨率格式。从而我们主要考虑两部分内容。首先基于CIP方法和高阶紧致方法，我们提出一种新型特征插值高阶紧致差分方法，并将其运用于数值模拟非线性波动方程。然后我们对守恒型双曲方程提出一种基于特征思想的有限体积方法，并应用于数值求解Eulelr方程组和浅水方程组。
论文的第一大部分，我们对非线性波动方程提出一种非守恒Semi—Lagrangian方法。这种格式的主要构造思想如下：根据两相邻节点的节点值和导数值构造三次多项式，然后运用Semi—Lagrangian的方法求解节点沿特征线移动的位置，由此三次多项式和节点位置得到时间上的推进。不同于传统的CIP格式，我们利用最初由Lele提出的高阶紧致格式直接利用节点值求解一阶导数值。最后我们将这种格式运用于数值模拟Burgers方程和KdV方程验证格式的高分辨率的性质。通过一系列的数值试验验证格式的有效性。最后将这种格式推广到一维的Euler方程组。
论文的第二大部分，我们对双曲守恒方程组构造有限体积格式，由于Euler方程组在空气动力学中的的重要性，我们主要基于Euler方程组构造格式。这种格式的空间导数通过有限体积离散，时间上采用Simpson公式离散。其中有限体积格式中的通量值通过单元边界点的值得到，而该单元边界点的值是通过中心加权重构得到。求解此节点值的主要思想如下：首先运用三阶或四阶Runge—Kutta方法求解特征方程，从而得到节点沿特征线的位置。然后采用CWENO重构得到多项式，我们分别构造三阶和五阶CWENO格式构造多项式，最后利用节点位置和此多项式得到节点值。这种高分辨率的有限体积格式结合特征思想和中心加权格式的优点。然后我们将这种基于特征思想的有限体积格式应用于一维的标量和Euler守恒方程，通过经典的算例验证格式的高分率和收敛的性质。
最后将这种基于特征思想的高阶有限体积格式应用于一些浅水问题，从而验证格式，如一维溃坝问题，临界左稀疏波问题，两稀疏波中间几乎为干底问题，干溃坝问题，生成干溃坝问题等。另外，我们同时将这种格式按照分裂方法推广到二维守恒律情形。并通过对二维的溃坝问题的应用验证格式的有效性，通过与WAF格式以及其它高阶数值格式相比较可以看到该格式得到非常好的数值结果。