基于ADI和IIM的界面问题算法研究

来源 :中南大学 | 被引量 : 3次 | 上传用户:yingying0615
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
界面问题是自然界中一种常见的现象,对界面问题的数值方法研究在工业、生物、军事等方面有着重要的理论意义和实际应用价值,近些年一直受到学者们的广泛关注,也成为计算数学研究的前沿问题之一.  首先,本文对界面问题数值方法的研究现状进行了综述,基于IIM方法,通过对传统的LOD差分格式在非正则点处的差分方程右端加上由跳跃条件确定的修正项,使得在非正则点处的局部截断误差为O(h),而正则点处的局部截断误差仍为O(h2),构造了二维和三维热传导界面问题的LOD-IIM差分格式,所得差分格式在时间和空间方向以无穷范数均二阶收敛,并用三个数值实验验证了格式的收敛阶和稳定性,  其次,利用IIM方法的思想,通过对高阶紧致ADI差分格式在非正则点处的差分方程右端加上由跳跃条件确定的修正项,使得在非正则点处的局部截断误差为O(h3),而正则点处的局部截断误差仍为O(h4),构造了二维热传导界面问题的高阶紧致ADI差分格式,所得的差分格式以无穷范数在时间和空间方向分别为二阶和四阶收敛且无条件稳定.基于固定界面的热传导问题的解关于时间的连续性,采用Richardson外推方法将时间方向的精度提高到四阶,从而得到二维热传导界面问题在时间和空间方向均四阶收敛的差分格式.基于所构造的差分格式,对两个分片光滑的数值算例进行求解,验证了格式的收敛阶和无条件稳定性.  然后,采用Adams-Bashforth方法处理非线性对流项,通过对传统的ADI差分格式在非正则点处的差分方程右端加上由跳跃条件确定的修正项,使得在非正则点的局部截断误差为O(h),而正则点处的局部截断误差仍为O(h2),构造了二维非线性对流扩散界面问题的ADI-IIM格式,所得格式在时间和空间方向以无穷范数均二阶收敛且无条件稳定.利用修正项关于时间的连续性,将当前时间层的修正项用前两个时间层上的修正项通过Adams-Bashforth方法逼近得到.对三个分片光滑的数值算例进行求解,验证了格式的收敛阶和稳定性以及通过Adams-Bashfoth方法逼近得到的修正项的可靠性.  最后,首次将Bernstein多项式用于椭圆界面问题数值方法的研究,在每个子区域上将方程的解表示为Bernstein多项式的线性组合,采用Galerkin方法或者配置方法计算Bernstein多项式的展开系数,得到一维椭圆界面问题和二维直线界面椭圆界面问题的高精度数值方法,数值实验表明该方法具有谱精度,对于二维曲线界面的椭圆界面问题,先通过一个非线性变换将其转化为直线界面问题并求解,再通过逆变换得到曲线界面问题的解,数值实验表明,该方法的逼近精度依赖于界面曲线的复杂程度.
其他文献
在传统的可靠性理论中,是将系统的寿命理解为随机变量.但在无法获得统计所需的大量数据的情况下,传统可靠性理论不能很好的解释现实生活中的不确定现象,使得研究结果与实际情况
考虑带有Hardy-Sobolev临界指数项和奇异项的Kirchhoff方程(此处公式省略)  其中,Ω是R3中的一个有界光滑区域且0∈Ω,a>0,b>0,0< s<1.  本文中,我们将在f(x,u)满足不同的条件下,
由泉州市美术家协会、泉州画院主办,桑莲居艺术馆和中国海西艺术网承办的《天地间的感悟——郭宁油画、水彩画作品集》首发式于2014年5月24日在桑莲居艺术馆举行。同时展出郭
数据收集是指一个传感器网络中,传感器采集环境中的相关数据,进行简单处理后,传输回网络基站的过程。数据收集是传感网络的基础功能。实际上,网络的大多数应用依赖于其数据收集的
设α为d次代数整数,它的极小多项式为P(χ)=χd+b1χd-1+…+bd-1χ+bd,其中bi∈Z,α1=α,α2,…,αd为α的所有共轭根.如果α的所有共轭根都是正实数,则称α是全实正的.α的所
学位
上世纪20年代,芬兰数学家Nevanlinna创立了值分布理论,通常为了纪念他而被称为Nevanlinna理论.这个理论包括了两个基本定理,我们把它们称之为第一基本定理和第二基本定理,该理论
学位
当系统所在区间不规则时,将变量简化为一维的方法便不再适用.因此,对于此情况需要建立一个新的模型以解决其边界控制问题.  本文主要研究了曲边矩形和曲顶柱体上反应扩散方