Helmholtz传输输特征值问题的高效计算方法

来源 :贵州师范大学 | 被引量 : 1次 | 上传用户:Water_E
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
传输特征值问题出现在针对一个非均匀的介质的逆散射理论中。它有着广泛的物理背景,比如:使用它们可以估计出散射对象材料的性质,它在逆散射理论的唯一性和重构方面也有重要理论意义。传输特征值问题是一个非线性二次特征值问题,其研究难度无论从高效计算方法还是方法的理论分析来看均明显大于线性特征值问题。由于该问题的重要性,该问题在被提出后就吸引众多相关领域的专家学者的关注。他们的研究工作主要集中在探讨非均匀介质中传输特征值的存在性和性质,折射率的上下界等理论问题上。近年来,传输特征值问题的数值计算方法是工程和计算数学界的热点课题。2010年,Colton等学者首次提出了传输特征值的三种数值方法(Argynis法,连续法,散度协调法),后来他们的工作引起的计算数学界的高度重视,许多著名学者提出了许多高效计算方法,包括迭代法,多网格方法,谱方法,混合有限元法等。其中迭代法需要一个迭代初值,这就需要确定特征值的一个近似值,此方法适合求解实传输特征值,其理论结果在单特征值前提下得到验证。现有的多网格方法是以迭代法为基础的,理论上具有同样类似的问题。而谱方法对区域的要求比较严格,仅适用于圆、矩形、圆柱等可分离的区域。现有的混合有限元方法虽然高效但缺乏理论分析作其支撑。本文致力于研究特征值问题的高效数值求解方法,其主要内容包括:第一章,介绍传输特征值有限元方法的预备知识,我们回顾了Helmholtz传输特征值问题的经典弱形式,并将其线性化为一个等价的线性非对称特征值问题。线性弱公式是我们研究该问题高效数值算法的基础。以线性弱形式作基础,我们进一步介绍了求解该问题的一个H2-协调元离散,该离散与2010年Colton等学者提出的Argyris元离散在代数上在等价。第二章,在协调元离散下,我们使用多水平校正技术建立了一个新的多网格方案,该方案将传输特征值在细网格上的解简化为一系列非常粗的网格上的解和相应的边值问题在一系列多水平网格上的解。最后我们给出方法的理论分析和数值结果。第三章,参考Verf¨urth建立的四阶问题后验误差估计理论,我们研究了传输特征值问题的协调元后验误差估计和自适应算法,这是一个新的工作,而且这也是高效计算方法中主流方向;我们首先给出源特征函数和共轭特征函数误差的后验估计子,并证明它的可靠性和有效性。根据特征值和数值特征值的基本关系式,我们将给出特征值误差的后验估计子。根据后验误差估计子我们将设计一个高效的自适应算法。这个算法可以用于求解复和重的特征值。在数值实验中我们将采用著名的Argyris元实现我们的自适应算法,并验证后验估计子的可靠性和有效性。从目前我们已经完成的数值例子来看,采用这个自适应算法,在很多情况下我们将获得具有最优收敛阶的数值特征值,即使对一个非凸的区域上的问题仍然有如此的最优收敛阶。第四章,我们提出了一个多维区域上的Hm-协调谱元法(m≥1),并将H2-协调谱元法用于高效求解传输特征值。我们将按下面的思路进行研究:(1)首先考虑一维标准单元[-1,1]上基函数的构造。由广义雅可比多项式的性质,我们可以利用勒让德多项式或切比雪夫多项式来构造[-1,1]上的泡泡函数,这些泡泡函数在端点-1和1上函数值连同一些低阶导数值均为0。为了保证谱元离散的Hm-协调性,我构造了节点函数值和一些低阶导数值基函数,此外,我们构造了一维标准区间上的谱元插值算子,并证明了相应的插值误差估计。(2)对于多维区域上的标准单元,我们采用一维基函数的张量积作为基函数,和一维插值算子的张量积作为高维单元上的插值算子。利用一维区间上的插值误差估计,我们证明了高维标准单元上的插值误差估计。在高维单元上,我们容易构造标准单元到一般单元的仿射变换。以此做铺垫,我们采用比例论证的方式证明了一般矩形单元上的插值误差估计。(3)我们将构造的H2-协调谱元离散用到传输特征值问题上,并证明该离散下特征值和特征函数的误差估计。最后我们将用数值试验验证该谱元法的高效性。第五章,作为现有文献中的H1型非协调元方法的补充,我们提出了L2型的非协调元方法,使方法覆盖不连续的非协调元空间。我们在分片光滑的Sobolev空间中使用Babuska-Osborn谱逼近理论,证明了该方法的最优收敛阶。在数值试验部分,我们呈现了使用三次四面体元求解该问题的几个数值例子。数值结果表明使用本文的离散和三次四面体元-线性元可以捕捉到比使用现有文献中的Morley-Zienkiewicz元方法更高精度的特征值。第六章,介绍在现有有限元离散下的求解传输特征值的三种其它数值算法和数值试验结果,包括二网格离散算法、H1型非协调元离散算法和Ciarlet-Raviart谱混合离散算法,并用数值试验证明了用前两种算法求解特征值可达到渐进最优收敛阶,使用Ciarlet-Raviart谱混合离散算法数值特征值可以达到谱精度。
其他文献
近年来,互联网金融飞速发展,客户需求的多元化和差异化日益凸显,在此背景下,商业银行网点改革亟待开展,商业银行逐步朝数字化和智能化方向转型。该文主要结合互联网的定义及特征,基于互联网视角,就互联网给商业银行主动业务、被动业务以及中间业务等所造成的影响展开探讨,进一步分析互联网背景下商业银行在基础设施、营销服务、业务流程、产品开发、经营等方面创新变革的必要性,并有针对性地提出商业银行进行网点转型的策略
学位
随着我国互联网金融的不断发展,中小商业银行传统的金融业务如贷款、支付结算和金融管理等都受到严重冲击,迫切需要转变其经营模式。本文通过对商业银行发展情况及转型状况进行综合分析,以徽商银行为例,探讨互联网金融背景下安徽省中小商业银行网点建设与互联网金融对接的优质途径。
推理是数学的基本思维方式,也是人们学习和生活中经常使用的思维方式,推理能力也是问题解决能力的核心,具有良好的数学推理能力对于学生今后进一步就业和工作有着重要的作用,学生只有“具有良好的推理能力,才能够形成有条理、有逻辑、有论据的良好思维习惯,从而提高探究事物本源的能力”,但“工欲善其事,必先利其器”。故研究在遵循一致性原则、完备性原则、本土化原则的基础上,拟建构高中生数学推理能力测评模型,力求为提
矩阵优化问题是数值代数和数值优化研究的重要课题之一.它在信号处理、图像处理、机器学习、数据分析、金融工程、量子计算、系统与控制理论和高维统计等科学与工程领域有着广泛的应用.本论文系统研究了如下几类具有广泛应用背景的矩阵优化问题的数值求解方法.为了提高聚类的精度,首次将Q-加权范数应用于非负矩阵分解,从而考虑了Q-加权范数下的非负矩阵分解问题(?)基于Q-加权范数的可加性表示方法,将该问题转化为一类
天基组网雷达是目标检测和空间探测的重要发展趋势之一,科学研究、经济社会、军事价值十分巨大,得到世界主要航天大国、强国的高度重视。但其实现技术难度大,仍需要解决航天器轨道高精度控制和相对位置高精度感知等重大难题。因此,论文开展天基簇飞行组网雷达覆盖、定位的随机几何优化研究,以期解决天基簇飞行组网雷达一些关键数学问题,为构建天基簇飞行组网雷达打下理论基础,提供技术支撑。论文主要进行了如下工作:(1)提
作为分布式空间系统的一种创新性结构,不像传统的星座,编队,集群卫星队形间需要保持紧控制,分离航天器由多个自由飞行,可物理分割的称作模块的小卫星组成,模块可以是同构的,也可以是异构的.通过将质量和故障风险进行分散,分离模块的可维护性和灵活性显著提高;模块间通过无线组网实现空间中信息交互和资源共享,大大提高了系统的可靠性.这种"功能分解,结构分离,无线连接"的技术特征,使得分离航天器网络具有快速响应,
阅读的好处:1.书籍是全世界的营养品。2.读书有助于孩子形成良好的品格和健全的人格。3.读书可以给我们打拼的勇气和战胜困难的力量。4.读课外书有助于积累词汇,提高写作能力。5.提高孩子的学习兴趣和积极性。总之,可以用培根的一句话来高度概括读书的好处:
期刊
C0内部惩罚Galerkin(C0IPG)方法是最近十年才发展起来的Galerkin方法的一个新类.该方法融合了连续Glerkin方法,不连续Glerkin方法和稳定性技巧的思想.对于四阶微分算子,它是一种不连的续Galerkin方法.与传统的有限元相比较,C0IPG方法采用Lagrange基函数,因而基函数容易构造,并且高阶Lagrange基函数能够高效地捕获光滑解,而且不破坏原问题的正定性.而