偏微分方程数值计算及增量未知元方法研究

来源 :兰州大学 | 被引量 : 1次 | 上传用户:liupengo0308
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文基于微分方程的有限差分技术以及一致网格增量未知元方法,分别对一维和二维具有时间依赖系数的热方程以及一类一般的三维对流扩散方程进行了不同的研究。由于一致网格增量未知元方法可以很好地降低矩阵条件数,所以该方法的优越性在我们的理论分析和数值实验中都很好地体现了出来。非一致网格作为一种更为灵活的形式,对于许多问题,特别是边界层问题的求解,有着一致网格所无法比拟的优势。相应地非一致网格上的增量未知元方法便自然地引起了我们的注意。所以在本文的后半部分,我们从理论方面,以Dirichlet问题为例,对一维和二维非一致网格增量未知元方法下的系数矩阵条件数进行了详细得分析,并用数值算例对我们的理论分析进行了验证。对于具有时间依赖系数的一维热方程,我们提出了一类增量未知元型半隐θ格式,仔细分析了这一新格式的稳定性,并对其误差进行了估计。结果显示,当θ趋近于1/2时,格式的稳定性条件会明显改善.而当θ趋近于零时,我们得到了一个能帮助我们恢复初始格式误差的条件。对于二维情形,我们依然用有限差分进行离散,构造了一类交替方向增量未知元型半隐(ADIUSI)格式。并在傅立叶方法的帮助下,对格式的稳定性进行了详细的分析。数值试验验证了理论分析的正确性,其结果表明这一新的格式,在某些问题上,会比经典的交替方向格式更为有效。对于一类一般形式的三维对流扩散方程,在有限差分和增量未知元方法下,可以得到一个增量未知元型正定但非对称的线性方程组。其系数矩阵条件数要远远优于不用增量未知元方法的情形。考虑到该方法的这一优点,我们在文中将其与几种经典的迭代方法相结合,来求解上述线性系统.并从理论上对该系统的增量未知元型系数矩阵条件数进行了估计,然后通过数值试验验证了这几种增量未知元型迭代方法的有效性。我们注意到上述的差分离散和增量未知元方法都是在一致网格上进行的,然而对于许多问题,例如边界层问题、流体力学问题等,一致网格上的差分离散已不能满足它们求解精度的需要。很自然地,要考虑非一致网格上的差分离散以及相应地非一致网格增量未知元方法。随之而来的问题是,这种非一致网格增量未知元方法能否也像一致网格增量未知元方法那样,能有效降低系数矩阵条件数?我们在本文后半部分对于该方法的一维和二维情形都进行了详细分析。理论结果表明该方法依然可以很有效地降低矩阵条件数,数值试验结果与我们的理论分析完全吻合。
其他文献
本研究首次以生长于乌鲁木齐河源区高海拔、强辐射、温度波动大的极端环境下的地衣为研究对象,通过对该地区地衣种类统计、生活型分析、区系比较、不同海拔的地理分布以及分子进化的研究,得到以下几点结论:1.该地区地衣分布有35种,隶属于14属,10科,其中,中国新记录种1种,即Rhizoplaca haydenii,揭示了极端环境下地衣分布数量较之其他地区较少,作为“先锋植物”的地衣,在高海拔、强辐射、温度
学位
双非线性抛物型方程在金属材料合成中有着重要的应用价值,但人们对其解的形态以及动力学行为了解甚少.在这篇博士学位论文中考虑如下形式的双非线性抛物型方程:其中α和f是任意次多项式增长的非线性函数,Ω是RN中的有界光滑区域.由于非线性项α的存在,给该问题的研究带来了很大的困难.首先,我们采用Faedo-Galerkin逼近方法,对更一般形式的双非线性P-Laplacian方程,证明了解的存在性.其次,利
学位
低温环境是一个独特的生态系统,它被认为是天然的“菌种保藏中心”。近年来冷环境下微生物的研究取得了长足的发展,尤其是微生物生态学的研究。天山地区具有寒冷、中纬度和高海拔等特点,是研究低温微生物的理想生境。本文选取了天山2米深的雪坑、1.6米深的活动层冻土和3.0米深的永冻土作为材料,应用PCR-DGGE技术研究了其中微生物的群落结构和分布,分析和讨论了它们之间的关系及其与温室气体(CH4和N2O)排
本工作采用能量为190keV-450keV的不同价态的Arq+(q=15~18)离子轰击原子序数Z=4~79的金属靶,对Ar离子的X射线发射机制进行研究。文中依次对Ar17+离子,Ar16+离子,Ar18+离子和Ar15+离子进行了分析。主要的实验工作及结论如下:对Ar17+离子的K壳层X射线分支比进行分析,研究了Kβ/Kα分支比随电荷态与能量的变化规律以及Kα和Kβ射线的平均能量。在级联跃迁模型
学位
本工作利用兰州大学2×1.7MV串列加速器产生的能量为17.9-120 keV/amu的C1+、30-323keV/amu的C2+、71.7-438 keV/amu的C3+、287-480 keWamu的C4+、15.9-80.4 keV/amu的O1+、30-254 keV/amu的O2+及63.5-321.6 keV/amu的O3+与惰性气体原子He、Ne、Ar进行碰撞,对碰撞中的转移电离过程
学位
本文以单调动力系统(monotone systems)为研究对象,在放松K-型序下的拟单调条件限制的同时考虑了斜积半流的动力学行为.分为自治和非自治动力系统两部分.在自治动力系统(autonomous systems)方面,已有的拟单调条件在许多应用中不能满足.一个典型的例子就是在神经元之间具有抑制和刺激的联络时的时滞Hopfield-型神经网络模型.受此启发,我们引入了弱拟单调条件(WQM),建
在这篇博士学位论文中,我们主要考虑非自治无穷维动力系统的拉回吸引子的存在性问题,针对拉回吸引子存在的关键性条件——(?)-拉回渐近紧的验证,提出了两种有效的验证方法,并将这两种方法应用到具体的非自治的无穷维动力系统中,得到了一系列新的深刻的结果.全文共分五章:第一章,介绍无穷维动力系统的背景,拉回吸引子的发展及研究进展情况,详细介绍了本文所讨论的主要问题和研究思想.第二章,给出了本文用到的一些基础
学位
由于兼具毛细管电泳(CE)超高分离效率和激光诱导荧光检测(LIF)超高灵敏度的优势,CE-LIF联用技术已经成为测定多组分化合物的一种有力工具,在药物及生命分析等领域得到了广泛的应用。但是,很多化合物本身并没有荧光,即使有些化合物自身具有荧光,其激发波长和所用激光光源波长的不一致也会限制CE-LIF方法的应用。为了克服这一缺点,分析工作者通常采用衍生反应修饰分析物,使之转变成具有适合光学性质的衍生
学位
作为匹配和拟阵交的共同推广,Cunningham和Geelen在1996年引入了图的路匹配的概念.他们指出许多领域的问题都可以转化为路匹配问题,也就是说,利用路匹配可以解决例如匹配、拟阵、多面体以及代数等很多方面的问题.作为路匹配的应用,他们给出了可匹配集合多面体的强多项式分离算法,并证明了最大路匹配的值就等于给定图所确定的匹配拟阵中顶点集合的秩,同时也等于Tutte矩阵的秩等等.本文共分为六章,
在这篇博士学位论文中,我们主要研究如下的两类反应扩散方程和解的长时间行为,主要是全局吸引子存在性和局部几何结构问题.对于第一类方程,我们从方程弱解的存在性出发,应用强弱连续半群的概念以及相关的判断吸引子存在性的方法,在f(u)是任意次多项式增长且λ>0是任意常数的情况下,得到方程在空间Lq(Ω)和H01(Ω)全局吸引子的存在性.而后,对全局吸引子的维数下界做出估计.从理论上说,应用Z2指标理论,我
学位