高阶导数方程的数值算法及其应用

发布时间：2017-05-23 13:27

  本文关键词：高阶导数方程的数值算法及其应用，由笔耕文化传播整理发布。

【摘要】：本文主要工作是应用局部间断有限元方法(LDG)求解偏微分方程及利用基于偏微分方程的正则化算法研究生物医学成像中的去噪问题。主要分成五部分。 第一部分主要研究了具有复杂几何结构的量子定向耦合器上的量子传输现象的数值求解问题。这个现象可以由一个静态Schrodinger方程来解释,我们运用最小耗散LDG方法来求解该方程。为了保证数值格式的稳定性,本文在边界数值流通量上添加惩罚项。此外,量子传输本身的物理性质决定了其频率的变化主要在y方向上,因而针对该问题的求解,我们除了应用基于多项式基函数的LDG方法外,还应用了基于指数基函数的LDG方法来减少计算量。数值算例验证了该算法的有效性。 第二部分针对含波动算子的非线性Schrodinger (NLSW)方程给出了守恒的LDG格式。NLSW方程的一个重要性质是能量守恒,而非守恒的数值格式很容易导致数值解爆破,因此我们在空间上使用LDG方法进行离散,而在时间上使用Crank-Nicholson格式进行离散,从而得到一个全离散的守恒数值格式并在理论上给出了能量守恒的严格证明。同时,我们还进一步给出了线性情形下半离散LDG方法最优误差估计的证明。数值结果证明了守恒数值格式在长时间数值模拟中的优越性。 第三部分采用高阶保正LDG格式求解含有爆破解的抛物方程。对于该类方程,如果初始条件及源项都是正的,则根据极大值原理,方程的解也是正的,而不加保正限制器的高阶格式会导致错误的爆破时间和爆破区间。由于该问题是Dirichlet边界条件,为了得到最优收敛阶,我们需要在边界数值流通量上添加惩罚项。数值结果表明,保正LDG格式能够准确捕捉到爆破时间和爆破区间。 第四部分主要研究基于偏微分方程正则化算法的去噪方法,主要应用于乳腺超声波弹性成像问题。传统的临床设备只能得到较为准确的轴向(平行于波)位移测算,而对于侧向(垂直于波)的位移测算质量是较低的。但是在多种临床弹性成像应用中(如模型重建及温度成像)同时获取轴向和侧向的精确超声波测算位移却是很重要的。因此,在该问题的研究中,我们主要使用在临床设备上获取的传统超声回声数据来提升侧向散斑追踪的精确度。该算法已被用于计算机模拟数据、组织模拟仿真及人体数据测试。 第五部分研究如何减少生物医学领域中磁共振成像中的噪音。由于在相位对比磁共振成像中,时域解析的3D(time-resolved3D)相位对比磁共振成像或磁共振造影(PC-MRA/MRI)只用来判定基本的信息,如：临床环境下动脉瘤(aneurysm)的尺寸和流率。而计算误差和潜在的缺陷都会对磁共振的结果产生不利的影响。但是要表征内动脉瘤血液动力学(intra-aneurismal hemodynamics)中的复杂流体,则需要更精确的血管成像。因此这里主要采用基于偏微分方程正则化的去噪算法减少噪音。计算机模拟数据验证了该算法的性能。
【关键词】：局部间断有限元方法 二维静态Schr(o|")dinger方程 量子传输现象 量子定向耦合器 含波动算子的Schr(o|")dinger方程 能量守恒 最优误差估计 保正 爆破解 Diriehlet边界条件 超声 弹性成像 散斑追踪 应变力成像 去噪 相位对比磁共振成像 流体不可压性 正则化
【学位授予单位】：中国科学技术大学
【学位级别】：博士
【学位授予年份】：2015
【分类号】：O241.82
【目录】：

• 摘要5-7
• ABSTRACT7-9
• 目录9-16
• 第1章 绪论16-24
• 1.1 间断有限元方法介绍18-20
• 1.2 时间离散方法20
• 1.3 变分方法20-22
• 1.4 本文的工作22-24
• 第2章 二维量子传输现象的局部间断有限元方法24-38
• 2.1 引言24-25
• 2.2 模型问题以及数值方法25-27
• 2.2.1 模型问题25
• 2.2.2 模型问题的数值方法25-27
• 2.3 典型的定向量子耦合器上的量子传输现象27-30
• 2.3.1 量子传输现象的二维静态Schrodinger方程28
• 2.3.2 量子传输现象的LDG方法28-30
• 2.4 数值结果30-35
• 2.5 本章小结35-38
• 第3章 含波动算子的非线性Schrodinger方程38-66
• 3.1 引言38-39
• 3.2 含波动算子的非线性Schrodinger方程的局部间断有限元方法39-52
• 3.2.1 复数域上的符号表示、内积和范数39-40
• 3.2.2 局部间断有限元方法40-41
• 3.2.3 能量守恒41-43
• 3.2.4 线性方程的误差估计43-52
• 3.3 时间离散52-55
• 3.3.1 差分算子52-53
• 3.3.2 时间离散53-55
• 3.3.3 全离散格式的实现55
• 3.4 数值结果55-64
• 3.4.1 一维情形55-59
• 3.4.2 二维情形59-62
• 3.4.3 三维情形62-64
• 3.5 本章小结64-66
• 第4章 具有爆破解的抛物方程66-88
• 4.1 引言66-68
• 4.2 保正高阶的局部间断有限元方法68-78
• 4.2.1 一维情形68-74
• 4.2.2 二维情形74-78
• 4.3 数值结果78-85
• 4.3.1 一维情形78-82
• 4.3.2 二维情形82-85
• 4.4 本章小结85-88
• 第5章 乳腺超声波弹性成像去噪88-102
• 5.1 引言88-90
• 5.2 材料和方法90-93
• 5.2.1 理论背景90
• 5.2.2 去噪算法的实现90-92
• 5.2.3 应变测算92
• 5.2.4 实验验证92-93
• 5.2.5 数据分析93
• 5.3 结果93-95
• 5.3.1 数值模型结果93-94
• 5.3.2 组织模拟模型数据结果94
• 5.3.3 体内乳腺组织数据结果94-95
• 5.4 讨论95-100
• 5.5 本章小结100-102
• 第6章 磁共振成像去噪102-110
• 6.1 引言102-103
• 6.2 材料和方法103-108
• 6.2.1 PC-MRA数据获取104
• 6.2.2 线下数据进程104-105
• 6.2.3 基于偏微分方程去噪算法的理论背景105-106
• 6.2.4 计算机合成数据的产生106-107
• 6.2.5 壁剪应力及相对血压估计107
• 6.2.6 数据分析107-108
• 6.3 结果108
• 6.4 本章小结108-110
• 第7章 总结与展望110-112
• 7.1 总结110
• 7.2 展望110-112
• 附录112-118
• A.1 透明边界112-113
• A.2 二维情形下最优问题的解113-115
• A.3 三维情形下最优问题的解115-117
• A.4 壁剪应力估计和相对血压问题的计算117-118
• A.4.1 壁剪应力估计117
• A.4.2 相对血压问题的计算117-118
• 参考文献118-132
• 致谢132-134
• 攻读博士学位期间发表的学术论文134

【参考文献】

中国期刊全文数据库 前3条

1 ;Global Existence and Uniqueness of Weak Solution to Nonlinear Viscoelastic Full Marguerre-von Kármán Shallow Shell Equations[J];Acta Mathematica Sinica(English Series);2009年12期

2 ;Parabolic Equations with VMO Coefficients in Generalized Morrey Spaces[J];Acta Mathematica Sinica(English Series);2010年01期

3 ;MULTISYMPLECTIC FOURIER PSEUDOSPECTRAL METHOD FOR THE NONLINEAR SCHR■DINGER EQUATIONS WITH WAVE OPERATOR[J];Journal of Computational Mathematics;2007年01期

  本文关键词：高阶导数方程的数值算法及其应用，由笔耕文化传播整理发布。

，

本文编号：388110