论文部分内容阅读
浅水方程的数值解法是计算流体学中的一个十分活跃的课题。众所周知,此类方程的一个重要特点是无论初值是否光滑,其解都可能出现间断。解的间断性给数值模拟带来了巨大困难,高于一阶精度的数值方法在间断附近会产生非物理振荡。虽然二阶TVD格式可以很好的控制间断,但是其在临界点退化为一阶格式,使得它对于光滑解也不是高阶的。因此如何克服间断附近的非物理振荡并提高解的分辨率成为首先要解决的问题。本论文基于此,从非物理数值振荡产生的物理原因出发,提出了群速度控制(GVC)格式,我们证明了此方法既能控制振荡,又具有较高的分辨率。针对带源项的非线性双曲守恒律方程,时间上我们利用二阶龙格库塔方法离散,空间上使用群速度控制方法离散,对于源项,采用一个分裂方法进行处理。最后的数值算例求解结果证明了所给方法可以很好的控制非物理数值振荡并具有较高的分辨率,好于一般的二阶TVD方法。此外,这个方法还具有计算量小、形式简单等特点。