无界区域问题和外部问题的区域分解谱方法及其应用

来源 :上海师范大学 | 被引量 : 1次 | 上传用户:bkguo2008
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
近三十年来,谱方法蓬勃地发展起来,为数值求解偏微分方程提供了又一个强有力的工具。谱方法的主要优点是高精度,从而被广泛应用于计算流体力学、数值天气预报、海洋工程、化学反应数值模拟和量子力学计算等。已有的谱方法一般适用于周期问题和有界区域问题,然而,科学和工程中的许多问题可归结为无界区域问题或外部问题。求解这类问题的方法之一是设定一个人工边界,加上人工边界条件,然后在有限子区域中用通常的方法计算,例如差分法、有限元方法或者有界区域上的谱方法等。但是,这种截断区域的办法会导致相应的误差。因此需要研究直接计算无界区域和外部区域问题的高精度算法。 本文研究无界区域和外部区域问题的区域分解谱方法,其理论基础是以无界区域中的正交多项式或正交函数系为基底的正交和插值逼近理论,并伴随区域分解法。 本文由以下四部分组成: 第一章简单回顾已有的无界区域问题的数值方法,特别是无界区域和外部问题谱方法研究的进展情况。概述了本文研究工作的动机、困难和所得到的主要结果,以及本文的结构。 在第二章中介绍一些一维Laguerre正交逼近和一维Legendre正交逼近的基本结果。将以带参数β的Laguerre正交函数和Legendre正交多项式为基函数,来逼近无界区域和外部区域问题的解,并通过参数β的适当选择,使得数值解更好地吻合精确解的渐近行为。 在第三章中,我们研究无穷带状区域上Fokker-Planck方程的高精度数值方法,它是非标准类型的偏微分方程。在已有的文献中,对一维问题采用Hermite谱方法,对二维问题采用Hermite谱-差分混合方法,但限制了解的精度,见[12,60]。一个赋有挑战性的问题是如何应用谱方法计算Fokker-Planck方程,从而得到高精度数值解。这个问题主要有以下几方面的困难:第一,它在一个自变量方向类似抛物型方程,而在另一个自变量方向类似双曲型方程。因此,我们不能用通常的方法构造合理的谱格式和分析数值误差;第二,其解在不同的子区域上满足不同类型的边界条件,因此不能应用统一的基函数;第三,方程中的某些系数在原点退化,而在无穷远处趋向无穷大,这给实际计算和数值分析带来许多困难。为克服上述这些困难,我们首先建立一维区域上以Laguerre函数为基底的拟正交逼近和Legendre拟正交逼近理论,然后在各个子区域上建立一些特殊的Laguerre-Legendre混合谱逼近理论。它们结合区域分解方法有机地组成了无穷带状区域上Laguerre-Legendre混合的谱逼近理论。在这些结果的基础上,我们构造了Fokker-Planck方程的Laguerre-Legendre混合谱格式,证明了它的收敛性,并在空间方向上有谱精度。数值例子显示了该算法的优越性,并和理论分析相吻合。为了简化计算,我们还考虑相应的拟谱方法。我们发展了一维区域上的以Laguerre函数为基底的Laguerre-Gauss-Radau插值逼近和Legendre-Gauss-Radau插值逼近理论,并在此基础上建立各个子区域及整个无穷带状区域上Laguerre-Legendre混合插值逼近理论。应用这些结果,我们构造了Fokker-Planck方程的Laguerre-Legendre混合拟谱格式,证明了它的收敛性及其在空间方向上的谱精度。数值例子同样表明了该算法的有效性和高精度。与相应的谱方法相比较,拟谱方法在实际计算中更简单、更节省计算时间。 在第四章中,我们研究二维外部区域问题的区域分解谱和拟谱方法。文献[25,56,57]给出了有关一维问题区域分解谱方法的一些结果。在文献[27,29,36]中,研究了圆外或球外问题谱方法,然而更困难也更有实际意义的问题是带有多角形障碍物的外部问题,至今尚未见到任何结果。计算这个问题的困难之处在于:如何在各子区域上合理逼近所讨论的问题?如何在两个相邻子区域的公共边界上匹配数值解?如何在空间方向上保持解的整体谱精度?如何估计整体数值误差?为解决这些困难,我们先研究一维区域上的有关Laguerre拟正交逼近和Legendre拟正交逼近,然后在不同的子区域上分别建立一些特殊的Laguerre-Legendre混合逼近和二维Laguerre逼近方法,并在此基础上,建立整个外部区域上的Laguerre-Legendre混合逼近理论。应用这些结果,我们构造两个模型问题的区域分解Laguerre-Legendre混合谱格式,证明了该方法的收敛性及其在空间方向上的谱精度。在实际计算中,我们构造了两类特殊的基函数,简化了计算。数值例子显示了该算法的有效性并与理论分析相一致。我们还发展了一维区域上以Laguerre函数为基底的Laguerre-Gauss-Radau插值逼近和Legendre-Gauss-Lobatto插值逼近理论,并给出各个不同子区域及整个外部区域上的混合插值逼近方法。在这些结果的基础上,我们构造了两个模型问题的Laguerre-Legendre混合拟谱格式,证明了它的收敛性及其在空间方向上的谱精度。数值例子同样表明了该算法的有效性和高精度,并比谱方法更简单、更节省计算时间。 本文有关逼近结果丰富与发展了有关谱方法的理论及其应用,特别是区域分解谱方法。文中采用的技巧同样适用于无界区域上的其它偏微分方程和外部问题的数值求解。
其他文献
不适定问题广泛的存在于地球物理、生物医学、材料科学、金融工程和工程控制等各个领域中,许多应用科学和技术都有这个问题,譬如:电容层析成像图像重建、地应力测量、晴空卫星红
在2011年出版的义务教育《美术课程标准》中提到:积极探索有效地学习方法,明确以“学生为本”的教学设计的指导思想;引导学生采用自主的合作、探究学习,帮助他们会学、爱学、
期刊
近些年来,许多学者对于Logistic方程解的存在性等问题进行了深入了研究,其中对Logistic方程的两个重要参数r,k是周期函数、渐近周期函数、概周期函数等方面(参见文献[1~7]),近期
对未知蛋白的功能注释是蛋白质组学的主要目标,其中一个关键的注释是对蛋白质亚细胞定位的识别。蛋白质必须处于合适的亚细胞位置才能发挥其功能,因此蛋白质的亚细胞定位信息
本文研究的主要内容包括两个方面:孤立子方程的可积系统与Darboux变换。其中主要从以下两个方面研究了孤立子方程的可积系统:即孤立子方程族的生成及其可积性和可积方程族的可
期刊
跟踪性概念是D.V.Anosov于1976年在讨论具有负曲率的微分流形上的测地流的时候首先提出来的.它在微分动力系统稳定性理论中有着重要的研究价值,也在计算数学中有着广泛的应用.199
记得儿时看过一部电影——《从奴隶到将军》,喜欢这部电影有诸多原因,但最吸引我的无疑是罗霄身上散发出来的独特个人魅力和他对于理想信念执着的追求.回顾自己十几年的从教
本文旨在研制基于平板探测器的动态DSA(数字减影血管造影术)软件系统,及X线人体自动融合拼接组件。DSA是显示血管的标准技术,在血管病的诊断及介入治疗中具有重要意义;X线人体