浅水方程阶梯黎曼问题精确解的求解方法研究

浅水方程阶梯黎曼问题精确解的求解方法研究

论文摘要

浅水方程(Shallow Water Equations,SWE)的间断问题被称为黎曼问题(Riemann Problem,RP)。黎曼问题又分为齐次黎曼问题和非齐次黎曼问题。其中含底坡源项的黎曼问题被称为阶梯黎曼问题(Step Riemann Problem),也称非齐次黎曼问题。齐次黎曼问题已经得到了充分的研究,其求解方法有近似黎曼解和精确黎曼解。近似黎曼解有Roe格式,HLL格式以及ASUM格式等,精确黎曼解可以用简单的Newton迭代法来解决,这些都取得了很好的计算结果。而阶梯黎曼问题目前还没有一个好的解决方法。本文尝试利用迭代法解决这一问题。主要工作如下:(1)阶梯黎曼问题是在初始条件的基础上新衍生出两个新的状态,这两个新的状态与两边的初始条件之间存在着数学关系,如果这种关系以激波的形式呈现出来,则新衍生状态的物理量与初始条件的物理量之间满足RH条件;如果这种关系以稀疏波的形式呈现出来,则新衍生状态的物理量与初始条件的物理量之间满足广义黎曼条件。这两种情况的的排列组合构成四种黎曼解的形式。与齐次黎曼问题不同的是,阶梯黎曼问题会衍生出两个新的状态,这两个新衍生状态的物理量之间存在着满足RH条件的关系,会以静态激波的形式呈现出来。这些关系组成方程组,本文的任务就是求解此方程组。(2)研究讨论了使用Newton法求解上述方程组时出现的问题,这些问题体现了 Newton法在求解浅水方程阶梯黎曼问题时的局限性:①Newton迭代法出现不收敛情况;②Newton迭代法收敛至错解处;③Newton法迭代过程中出现雅可比矩阵奇异或接近奇异,迭代难以继续进行下去。针对Newton法的局限性在其他文献中寻找了两种改进算法:中点求积法和连续型修正Newton法来避开Newton法的缺点。其中,中点求积法采用对迭代初值进行优化的方法来达到对Newton法的改进;而连续型修正Newton法则通过对迭代步步长的优化来达到对Newton法的改进。并在几种不同形式的黎曼解下给出了两种算法的算例验证。此外,在算例验证的过程中,采用了简化方程组的方法以加快方程的求解速度。(3)研究讨论了基于Gauss-Newton法的优化型算法——阻尼牛顿法,也即LM算法,并对算法的数值特征进行了研究。介绍了 LM方法阻尼因子的选择及相应配套策略,结合信赖域方法调整自适应因子,从而构成自适应LM算法。详细解释了自适应LM算法的迭代过程;从该数值过程可以看出:通过对阻尼因子的调节,自适应LM方法具有了避开雅克比矩阵近似奇异区的能力,从而能够求解非线性方程组奇异问题。并对算法在较差迭代初值的条件下进行验证。

论文目录

  • 摘要
  • Abstract
  • 第一章 引言
  •   1.1 研究背景
  •   1.2 非线性双曲型守恒律方程组近似解法
  •   1.3 非线性双曲型守恒律方程组的精确求解方法
  •     1.3.1 针对不收敛或收敛至错解的改进方法
  •     1.3.2 针对奇异问题的改进方法
  •   1.4 研究内容
  •   1.5 研究思路
  • 第二章 黎曼问题及其解的基本结构
  •   2.1 黎曼精确解
  •     2.1.1 不含源项的浅水方程
  •     2.1.2 含底坡源项的浅水方程
  •     2.1.3 变量之间满足的关系
  •   2.2 主要结果
  •   2.3 缩减变量方法
  •   2.4 本章小结
  • 第三章 两种针对Newton的缺陷的改进算法
  •   3.1 引言
  •   3.2 计算阶梯黎曼问题的几种算法
  •     3.2.1 中点求积法
  •     3.2.2 连续型修正Newton法
  •   3.3 两种改进算法应用
  •     3.3.1 程序正确性验证
  •     3.3.2 求解逆矩阵方法在Newton算法中效率验证
  •     3.3.3 初值条件较差时收敛性比较
  •   3.4 本章小结
  • 第四章 优化型算法
  •   4.1 L-M算法
  •     4.1.1 L-M法基本原理
  •     4.1.2 LM算法的数值特性
  •     4.1.3 结合信赖域算法的自适应LM算法
  •     4.1.4 自适应LM算法迭代数值特点分析
  •   4.2 自适应LM算法应用
  •     4.2.1 程序正确性验证
  •     4.2.2 初值条件较差时收敛性比较
  •     4.2.3 几种算法整体优化效率
  •   4.4 本章小结
  • 第五章 总结与展望
  •   5.1 总结
  •   5.2 展望
  • 参考文献
  • 攻读硕士期间科研成果
  • 致谢
  • 文章来源

    类型: 硕士论文

    作者: 陈伟

    导师: 许仁义,张松

    关键词: 阶梯黎曼问题,精确求解,中点求积法,自适应算法

    来源: 扬州大学

    年度: 2019

    分类: 基础科学,工程科技Ⅱ辑

    专业: 地球物理学,水利水电工程

    单位: 扬州大学

    分类号: TV13

    总页数: 68

    文件大小: 3972K

    下载量: 83

    相关论文文献

    • [1].基于任意拉格朗日—欧拉型移动网格的反应流广义黎曼问题方法(英文)[J]. 计算物理 2020(02)
    • [2].非等熵Chaplygin气体极限黎曼解关于扰动的依赖性[J]. 数学年刊A辑(中文版) 2018(02)
    • [3].二维黎曼问题一些新解[J]. 数学学习与研究(教研版) 2008(11)
    • [4].压差方程的广义黎曼问题格式[J]. 郑州大学学报(理学版) 2012(03)
    • [5].二维线性化压力梯度方程组的黎曼问题[J]. 云南大学学报(自然科学版) 2018(06)
    • [6].气体动力学压力梯度方程组的广义黎曼问题[J]. 昆明学院学报 2008(04)
    • [7].等熵Chaplygin气体动力学系统三片常数的黎曼问题[J]. 安庆师范大学学报(自然科学版) 2018(02)
    • [8].一类解耦非线性双曲守恒律系统的广义黎曼问题[J]. 昆明学院学报 2009(06)
    • [9].非线性波方程含三片常数的黎曼问题[J]. 应用数学学报 2014(04)
    • [10].具有常外力项的零压欧拉方程组初值涉及狄拉克函数的黎曼问题[J]. 数学的实践与认识 2019(07)
    • [11].带有源项的Chaplygin气体非对称Keyfitz-Kranzer方程组含狄拉克初值的广义黎曼问题[J]. 新疆大学学报(自然科学版) 2019(03)
    • [12].初值带δ函数的运输方程的黎曼问题(英文)[J]. 应用数学与计算数学学报 2018(02)
    • [13].2维非凸标量守恒律分3片黎曼问题的数值解[J]. 云南大学学报(自然科学版) 2010(06)
    • [14].一个三角型双曲守恒律系统的狄拉克激波解[J]. 鲁东大学学报(自然科学版) 2017(04)
    • [15].零压相对论欧拉方程组含有狄拉克函数初值的黎曼问题[J]. 数学物理学报 2014(05)
    • [16].一个最简带粘性项的非凸ZND燃烧模型的黎曼问题(英文)[J]. Journal of Shanghai University(English Edition) 2008(05)
    • [17].非线性Schrdinger方程的直接间断Galerkin方法[J]. 计算物理 2012(02)
    • [18].一维等温欧拉方程组的Delta激波解[J]. 上海电机学院学报 2014(02)
    • [19].关于一种微分方程解的研究[J]. 科学中国人 2016(02)
    • [20].带群集耗散项的零压流方程的扰动黎曼问题[J]. 数学物理学报 2020(01)
    • [21].一种带有高精度多维限制器的中心格式[J]. 湘潭大学学报(自然科学版) 2019(01)
    • [22].气体动力学零压流的狄拉克激波和真空相互作用的数值分析[J]. 云南大学学报(自然科学版) 2013(03)

    标签:;  ;  ;  ;  

    浅水方程阶梯黎曼问题精确解的求解方法研究
    下载Doc文档

    猜你喜欢