论文部分内容阅读
本文考虑的一维双曲守恒律方程(组)。我们给出了一种设计高分辨率格式的方法。此格式是Godunov型的,用的是分片线性重构。与传统的Godunov型格式不同的是此格式在计算过程中不仅计算数值解还计算了数值熵。在每个网格中线性重构函数的斜率是根据熵耗散得到的,它要求此网格上重构函数熵的网格平均应与此网格上的数值熵相等。数值解是根据有限体积法(finite-volume)求得的,同时数值熵计算的时候格式中有一熵耗散项,它在计算过程中耗散熵.通过这种方式为数值计算引入了适当的粘性,稳定了计算。
在[25]中我们对于特殊的熵函数分析了格式的精度。本文中对于一般的熵函数,我们分析了格式的精度,格式在远离极值点附近为二阶精度,在极值点附近为一阶精度。因为格式与传统的守恒型格式不同,所以我们给出了它的相容性定义和Lax-Wendroff定理,定理说如果用我们的数值格式求得的数值解收敛,则它一定收敛到方程的满足熵条件的弱解。
设计这样一种格式的一个重要动机是期望用此类格式来克服守恒型方程组的数值模拟中,线性特征场上的数值耗散问题。为此对于线性传输方程,我们用不带熵耗散函数的格式进行了一些数值实验,以研究数值熵对消除线性耗散问题的作用。在此研究的基础上,我们设计了两种格式,一种格式中我们仍给线性方程或方程组的线性特征场以一定的熵耗散以稳定计算。另一种格式中我们保持线性方程或方程组的线性特征场上的熵守恒,同时为数值解的重构设计了一种所谓的“极大值减少极小值增”(Minimums-Increase-and-Maximums-Decrease或MIMD)斜率控制因子以消除间断附近的非物理振荡。
最后我们分别给出了用带线性熵耗散的格式和不带线性熵耗散但在数值解的重构中运用了MIMD斜率控制因子的格式进行数值实验,包括单个方程的和方程组的。从中我们可看出熵耗散因子是如何抑制非物理的振荡的,以及格式对计算的有效性。