求解哈密顿-雅可比方程和非守恒双曲方程组的高分辨率方法

摘要

本文主要研究高分辨率方法，如高阶有限差分和有限体积加权本质无振荡格式（即：WENO格式）以及间断Gerlerkin限元方法（即：DG方法），用于求解哈密顿-雅可比方程和非守恒双曲方程组，特别是流体力学中的非守恒欧拉方程组；并成功地应用于行人流问题和多相流问题。论文主要分为如下两大部分：
　　 论文的第一大部分，我们推广了求解静态哈密顿-雅可比方程中的快速扫描三阶WENO格式至快速扫描五阶WENO格式，其中包含静态哈密顿-雅可比方程中一类特殊的方程，Eikonal方程。结合快速扫描高阶WENO格式求解Eikonal方程，高阶WENO格式求解双曲守恒律方程，以及Runge-Kutta时间离散，我们对二维宏观动态反应连续行人流模型进行了高阶数值模拟。宏观动态反应连续行人流模型是基于动态反应行人平衡原则，在此原则下行人选择行进路线，使得到达目的地的瞬时行进成本最小。它的控制方程为行人流密度的质量守恒双曲方程以及决定流通量方向的势函数的Eikonal方程的耦合方程组。通过数值模拟比较，高阶格式在粗网格下能够更有效的得出相应低阶格式细网格下的结果。而且我们用高阶高分辨率方法模拟宏观相向行人流模型，得出了明显的相向行人流分层现象，相比于微观模型下的行人流分层现象，宏观模型能更加有效更加快速的得到这样的结果。对于依赖时间的哈密顿—雅可比方程，本质无振荡格式（即：ENO格式）和WENO格式，以及DG方法和局部间断Galerkin有限元方法（即：LDG方法），都是求解此类方程的非常有效的高阶高分辨率方法。本文中，围绕哈密顿-雅可比方程，我们给出了直接求解一维和二维依赖时间的哈密顿—雅可比方程的DG方法和LDG方法的最优L2先验误差估计。
　　 论文的第二大部分，我们用高阶有限体积WENO格式和次区间分辨技巧求解非守恒欧拉方程组，并应用于多相流问题。高分辨率方法如有限差分和有限体积高阶ENO格式和WENO格式，都能够很好的求解单介质的守恒欧拉方程组。但对于多介质的守恒欧拉方程组，高分辨率方法在介质界面处有很强的振荡，这种振荡是目前几乎所有经典格式都本质存在的。因此对于这种多介质问题，如果考虑原始变量的非守恒方程组，它能够更好的模拟多介质界面移动的问题，并能够得到精准单调的解。对于包含非守恒乘积项的方程组，严格的弱解是定义在一组积分路径上的。基于路径积分理论，C．Pares等人发展了一系列的路径守恒格式，但这些格式大部分只应用于求解浅水方程。路径守恒格式主要的问题是如何得到正确的积分路径，并确保数值解能够收敛到正确的解。R．Abgrall和S．Karni在2010年指出这种路径守恒格式在求解非守恒欧拉方程组时存在着问题，路径守恒格式收敛不到正确的激波解。本文中，我们指出大部分本质无振荡格式，如全变差递减格式（即：TVD格式），ENO和WENO格式，它们在间断处都会被抹平。这种抹平，使得格式的解在间断处存在过渡点，并且这种过渡点没有落在真正的间断积分路径上。因此路径守恒格式通过这些过渡点来构造积分曲线，不能够得到正确的积分路径，从而得不到正确的数值解。我们的做法是采用有限体积WENO格式和次区间分辨技巧，尽可能减少间断面的过渡点，来磨尖数值格式在间断面处的解。我们通过相邻区间的近似多项式延伸来得到数值格式在包含间断的区间里间断面处两端的点值，然后求解这组端点值满足的守恒欧拉方程组的黎曼问题。最后我们用该黎曼问题的真解作为包含间断的区间间断面处的积分路径，从而可以得到比较正确的路径守恒格式。一维欧拉方程组的单介质和多介质流问题的数值试验验证了我们新方法的有效性。