带有限制器的cpr方法和交通流模型研究计算数学专业论文.docxVIP

研究生优秀毕业论文 中国科学技术大学学位论文原创性声明本人声明所呈交的学位论文，是本人在导师指导下进行研究工作所取得的成 中国科学技术大学学位论文原创性声明 本人声明所呈交的学位论文，是本人在导师指导下进行研究工作所取得的成 果。除己特别加以标注和致谢的地方外，论文中不包含任何他人已经发表或撰 写过的研究成果。与我一同工作的同志对本研究所做的贡献均己在论文中作了 明确的说明。 1 作者签名： 牵名 签字日期： 沙|s．S。站 中国科学技术大学学位论文授权使用声明 作为申请学位的条件之一，学位论文著作权拥有者授权中国科学技术大学 拥有学位论文的部分使用权，即：学校有权按有关规定向国家有关部门或机构 送交论文的复印件和电子版，允许论文被查阅和借阅，可以将学位论文编入 《中国学位论文全文数据库》等有关数据库进行检索，可以采用影印、缩印或扫 描等复制手段保存、汇编学位论文。本人提交的电子文档的内容和纸质论文的 内容相一致。 保密的学位论文在解密后也遵守此规定。 翻公开 口保密——年 ， 签字日期：盘!盘．皇!羔J 签字日期 塑年．一 作者签名： 导师签名 万方数据 摘要摘要 摘要 摘要 本文关注守恒律和哈密顿．雅可比这两类偏微分方程的数值解法及其在交通 流模型中的应用。我们可以大致地将论文分为以下两个部分： 在第一部分里，我们研究用于求解双曲守恒律的重构修正方法(correction procedure via reconstruction，简称CPR)。由于CPR方法为一个高阶线性方法， 其在处理解中含有较强间断的问题时可能会产生数值震荡。因此，我们 分别针对结构网格和非结构网格(直边或曲边)，将一个原本为间断Galerkin (discontinuous Galerkin，简称DG)格式设计的简单的加权本质无震荡(weighted essentially non．oscillatory，简称WENO)限制器推广到了CPR格式上。其目的是 能够控制CPR格式的解在间断处的数值震荡，同时又能保持解在光滑区域内的 原有高精度。我们使用的WENO限制器十分简单，其不会破坏CPR格式本身的 守恒性。此外，WENO限制器作用在目标单元上时只需用到该单元及其直接相 邻的几个单元上的信息，因此它可以维持CPR格式本身的模板紧凑性。标量守 恒律的熵解的一个重要性质是其满足最大值原理。特别地，当初值为正时解在 后续时刻都将保持为正。在实际中一些物理量应该为正数，例如可压缩气体动 力学里的欧拉方程中的密度及压力。由于带有wEN0限制器的CPR格式在上述 情况里不能自动地保持数值解为正数，我们还将一些原本为DG格式设计的保正 限制器推广到了CPR格式中。论文中给出了一维及二维上的数值算例来证明这 些限制器的有效性。 在第二部分里，我们分别针对各向同性和各向异性这两种不同的情形，研 究了动态交通流问题的建模及数值解法。对于各向同性问题，Jiang等人[541提 出了一个预测型连续动态用户平衡模型。由于模型的路径选择策略中的一些问 题，其模型在数学上是不适定的。因此，我们重新考虑了该问题，提出了一个 新的路径选择策略并构建了一个改进的模型，用于处理含有单个商业中心的 任意形状的密集城市内的交通问题。对于各向异性问题，Hoogendoorn和Bovy [45]构建了一个用于解决行人用户最优动态分配问题的方法。虽然这个模型的 适用性非常广泛，但其文章中只给出了各向同性情形下的数值算例。我们指出， 其模型中的哈密顿．雅可比一贝尔曼方程在各向异性条件下较难进行数值求解。 为了克服这一困难，我们针对含有多个商业中心的密集城市重新构造了他们的 模型。在我们的模型中，哈密顿一雅可比．贝尔曼方程被简化为了一个更容易计算 的哈密顿一雅可比方程。本文构造的各向同性和各向异性模型均由一个守恒律和 一个依赖时间的哈密顿．雅可比方程组成。守恒律用来控制交通流的密度，该方 I 万方数据 摘要程中的交通流方向由所构造的路径选择策略来决定。哈密顿．雅可比方程则用来 摘要 程中的交通流方向由所构造的路径选择策略来决定。哈密顿．雅可比方程则用来 求解实际总交通花费。我们在论文中运用稳健的数值格式来求解这些守恒律和 哈密顿．雅可比方程。对于一般的由两个方程组合而成的模型而言，给定的都是 初始时刻的值，因而我们可以在时间上同时求解这两个方程。然而本论文的模 型中两个方程在时间上的求解方向不同，因此没法同时进行计算。事实上，同 时满足这两个方程的问题可以看做一个不动点问题。我们在论文中构造了一个 自适应的连续平均法来求解此不动点问题。该方法可以利用最小二乘法来自动 地寻找连续平均法的最佳步长。在各向异性的模型中，我们需要求解一个最小 值问题。在本文中我们构造了一个简单的方法来数值求解此最小值问题。