抛物型方程的几种可并行的有限差分方法

论文摘要

抛物型方程是偏微分方程中基本方程之一.在自然科学的众多领域中,许多现象是用抛物型方程或者方程组描述的[1],例如热传导以及其它扩散现象、化学反应、粒子的运输等等.另外在一些问题的数值处理中也经常出现抛物型偏微分方程.在现代数值方法中,最早为人们所注意且理论分析完善的是有限差分法,因此抛物型偏微分方程的有限差分方法一直是人们关心的焦点.随着并行机的问世和发展,传统的有限差分方法在不同方面暴露出各自的弱点.例如,古典显式虽然适合于并行计算,但它是条件稳定的,特别是多维问题中计算步长受到严格的限制;古典隐式和Crank-Nicolson格式是绝对稳定的,但需要求解联立方程组,不便于直接在并行机上应用.因此需要构造具有良好稳定性、并行性和计算精度的新的差分方法.七十年代初,Miranker指出用有限差分逼近偏微分方程时,主要是组织传统差分方法的并行实现,至于设计新算法,推动力是很小的;之后的十几年高阶差分格式方面的研究得到了发展[3-7].八十年代初,上述情况由于Evans和Abdullah的工作[8-12]而发生了变化,他们设计的分组显式方法保证了数值计算的稳定性,同时由于显式求解而使该方法具有很好的并行性质.它是不同类型Saul’yev非对称格式[13]的恰当组合.由于不同的Saul’yev格式的截断误差中某些项绝对值相等,符号相反,在同一时间层和不同时间层上连续交替使用不同的非对称格式,可带来截断误差的部分抵消,从而提高方法的计算精度.这些非对称格式都是隐格式,但由于它们之间的巧妙结合,可以显式求解,这就是Evans-Abdullah分组显式（GE）.这项工作说明了建立满足上述要求的新的差分格式是可能的.但是在将分组显式思想应用于变系数问题时,稳定性的证明遇到了困难.在此基础上,张宝琳等在[14-16]中提出利用Saul’yev非对称格式构造分段隐式的思想,并恰当的使用交替技术建立了多种显-隐式和纯隐式交替并行方法,取得了稳定性和并行兼顾的研究成果.之后又将方法推广到变系数问题,并用能量法证明了方法的绝对稳定性.在数值试验中发现,分段或分块并行计算的结果一般都比原来相应的未加分裂时的结果精确.所以通过分而治之的策略来建立新算法,不但可以用于并行,还可以提高精度.之后涌现出大量的并行差分算法的研究成果,韩臻在[17,18]中详细研究了一类纯显-隐分段和分块交替方法;[19-21]中冯慧等通过不同点的隐式差分格式之间的相互约化来建立新型迭代方法,此方法和Jacobi方法同样具有并行性,却比Jacobi收敛快.[22]中张志跃等给出变系数抛物型问题的分组显式方法,并用能量方法给出稳定性证明.[23-26]中王文洽等针对不同的方程建立了分段的显隐格式,证明方法的稳定性并给出数值算例.上述的方法在并行性和稳定性方面都有其优良的表现,但是他们都存在一个共同的问题,那就是它们都是基于二阶差分格式建立的,这直接影响数值计算中空间的误差精度.近年来,研究人员开始致力于研究高阶紧致差分格式.[27]中Sanjiva K.Lele提出了高阶的紧致差分格式,文中对格式的误差做了Fourier分析,并将它与经典差分格式做了比较.[28-32]中Mark H.Carpenter等针对不同问题提出高阶紧致格式并给出理论分析.将这些高阶差分格式与交替分组思想结合起来,是否可以得到稳定性好,可并行且精度高的数值算法呢?近年来,涌现出大量的高阶交替分组格式的研究工作[44-52].本文作者在王文洽教授的精心指导下,就抛物型问题的几类数学模型利用有限差分方法的技巧,构造了具有良好数值性质和计算效果的迭代方法、分组显式方法和交替分段方法,对方法做了理论分析并给出算例说明方法的适用性.本人拓广了前人的工作,不具有重复性.本文共分为五章.第一章中主要利用[19]中冯慧提出的数值Stencil的概念,将其应用于二维对流扩散方程,建立了比Jacobi迭代收敛快的新型迭代算法.本章首先给出针对对流扩散方程的数值Stencil的概念,经过三次消元过程得到最终的数值Stencil,在此基础上建立了新型迭代算法;通过分析迭代误差证明了方法的收敛性,并与Jacobi迭代比较收敛阶;最后数值试验说明方法的适用性,证实了理论分析的结论.本章内容已被《International Journal of Computer Mathematics》接受.第一章的创新之处在于将[19]中的方法应用到含有时间项的高维抛物型问题中,建立了收敛速度快、具有并行性质的新型迭代格式,通过分析迭代误差证明了方法的收敛性以及与古典迭代法之间收敛阶的比较;最后用实际例子说明了算法的有效性.第二章主要运用[22,24]中的构造思想,将中心差分格式与分组显式思想相结合,针对含有变系数的对流扩散方程建立了分组显式方法,并用能量方法证明了该格式的稳定性.本章首先给出基于Crank-Nicolson差分格式的四种非对称的逼近方程,通过它们的巧妙组合建立交替分组显式方法;由于扩散项为变系数,所以采用能量法证明稳定性;最后数值试验说明方法的适用性.本章内容已投到《International Journal ofComputer Mathematics》.第二章的创新之处在于对于变系数的抛物问题给出和Crank-Nicolson格式相匹配的交替分组显格式,并用能量方法证明了稳定性.接下来的三章内容中,主要借鉴了[14,23-26,46-48]中交替分组（段）格式的思想,将它与高阶差分格式[27-32]相结合,建立了高阶的交替格式;经证明方法都是绝对稳定的,且具有并行性质;在时间步长足够小时,空间的局部截断误差可达到O（h4）第三章中首先给出高阶的显、隐差分格式,在隐格式的基础上构造了四种非对称格式,通过它们之间的巧妙组合建立了交替分段显隐格式;由[33-34]中的Kellogg引理证明了方法的无条件稳定性;得到了方法的局部截断误差;数值算例证实了方法的实用性,并且可以达到O（h4）的误差精度.本章内容已被《计算物理》接受.第三章的创新之处在于将高阶差分格式与交替分段显隐格式的思想相结合,构造出高阶的交替分段显隐格式.方法具有良好的数值稳定性,空间误差阶可以达到O（h4）阶.第四、五章是在第三章的基础上,引入高阶Crank-Nieolson差分格式并适当变形,构造了八个非对称逼近方程;通过交替使用这些差分格式建立了两种不同的数值计算方法.第四章是单独应用八个非对称差分格式构造了交替分组显格式,通过Kellogg引理证明方法的稳定性,通过误差分析得到两层抵消部分误差后的误差可以达到O（Τh）阶;数值试验说明了方法的实用性,并且时间步长充分小的前提下误差对于空间来说可以达到四阶.本章内容发表在《山东大学学报》（理学版）.第五章利用交替分段格式思想,引入了在非对称格式之间插入对称格式的思想,从而建立交替分段Crank-Nicolson格式.本格式同样具有数值稳定性、可并行性质,误差分析时由于插入中心对称格式,所以在这些点处两层之间部分抵消后误差较小,这在数值试验中得到了证实.交替分段Crank-Nicolson格式的结果比相应的交替分组显格式结果要好.本章内容已投稿到《应用数学与力学》.后两章的创新之处在于将高阶差分格式与交替分组、交替分段Crank-Nicolson思想充分的结合,建立与第三章不同的非对称Saul’yev格式,在此基础上构造相应的交替方法.这些方法都具有绝对稳定性、可并行性质;并且在空间上截断误差可以达到O（h4）阶;最后给出数值算例说明方法的适用性.

论文目录

• 中文部分
• 符号
• 中文摘要
• 英文摘要
• 第一章 二维对流扩散方程有限差分方法的Stencil以及新迭代格式
• §1.1 引言
• §1.2 预备知识
• §1.3 数值Stencil和Stencil消元
• §1.4 新迭代法和收敛性
• §1.5 数值算例
• 第二章 变系数对流占优扩散问题的交替分组显格式
• §2.1 引言
• §2.2 预备知识
• §2.2.1 差分格式
• §2.2.2 基本引理
• §2.3 稳定性
• §2.4数值试验
• 第三章 对流扩散问题的高精度交替分组显隐格式
• §3.1 引言
• §3.2 求解公式
• §3.3 新交替分段显隐格式
• §3.3.1 新格式
• §3.3.2 数值分析
• §3.4 数值试验
• 第四章 对流扩散方程的高精度交替分组显格式
• §4.1 引言
• §4.2 求解格式
• §4.3 交替分组八点格式
• §4.3.1 格式构造
• §4.3.2 数值分析
• §4.4 数值算例
• 第五章 对流扩散方程的新型高阶交替分段Crank-Nicolson格式
• §5.1 引言
• §5.2 差分格式
• §5.3 新型交替分段Crank-Nicolson格式
• §5.3.1 方法构造
• §5.3.2 数值分析
• §5.4 数值算例
• 参考文献
• 致谢
• 攻读博士学位期间完成论文情况
• 作者简介
• 学位论文评阅及答辩情况表
• 英文部分
• Notation
• ABSTRACT
• CHINESE ABSTRACT
• Chapter 1 A Stencil of the Finite Difference Method for the 2D Convection Diffusion Equation and Its New Iterative Scheme
• §1.1 Introduction
• §1.2 Preliminaries
• §1.3 Numerical stencil and stencil elimination
• §1.4 The new iteration and its convergence
• §1.5 Numerical experiments
• Chapter 2 The Alternating 4opoint Group scheme for Solving Convection-Dominated Diffusion Problem
• §2.1 Introduction
• §2.2 Preliminaries
• §2.2.1 Tarabicsubsection Difference scheme
• §2.2.2 Trarabicsubsection Basic Lemmas
• §2.3 Stability
• §2.4 Numerical experiments
• Chapter 3 High-Order Accurate New Alternating Segment Explicit-Implicit Schemes for The Convection Diffusion Problem
• §3.1 Introduction
• §3.2 Solving formula
• §3.3 The NASE-I scheme
• §3.3.1 Tarabicsubsection The scheme
• §3.3.2 Tarabicsubsection Numerical Analysis
• §3.4 Numerical Simulations
• Chapter 4 High-order Accurate Alternating Group 8-Point Schemes for The Convection Diffusion Problem
• §4.1 Introduction
• §4.2 Solving formula
• §4.3 The AG 8-point scheme
• §4.3.1 Tarabicsubsection The scheme
• §4.3.2 Tarabicsubsection Numerical Analysis
• §4.4 Numerical Simulations
• Chapter 5 A New High-Order Accurate Alternating Segment Crank-Nicolson Schemes for The Convection Diffusion Problem
• §5.1 Introduction
• §5.2 Solving formula
• §5.3 The nASC-N scheme
• §5.3.1 Tarabicsubsection The scheme
• §5.3.2 Tarabicsubsection Numerical Analysis
• §5.4 Numerical Simulations
• Bibliography
• Acknowledgement
• List of Publications during Study for the Doctorate
• Curriculum Vitae
• 学位论文评阅及答辩情况表