【摘 要】
:
提出不同类型的增量未知元用于构造有限差分数值格式。本文主要考虑增量未知元以下方面。首先,通过增量未知元方法建立适合三维偏微分方程的清晰矩阵框架,在三维情形提出增量未知元的多层格式,通过数值试验,证实包括求解泊松方程,增量未知元的分层预处理形式可以用于更一般的方程。其次,将增量未知元方法和一些现代迭代方法相结合,如MR,GCR,Orthomin(k),Bi-CGSTAB,HSS,BTSS等,并求解由
论文部分内容阅读
提出不同类型的增量未知元用于构造有限差分数值格式。本文主要考虑增量未知元以下方面。首先,通过增量未知元方法建立适合三维偏微分方程的清晰矩阵框架,在三维情形提出增量未知元的多层格式,通过数值试验,证实包括求解泊松方程,增量未知元的分层预处理形式可以用于更一般的方程。其次,将增量未知元方法和一些现代迭代方法相结合,如MR,GCR,Orthomin(k),Bi-CGSTAB,HSS,BTSS等,并求解由三维对流扩散方程多层离散生成的非对称正定线性系统。通过理论分析,估计出与静态对流扩散方程相关的增量未知元矩阵条件数,以及MR,GCR,Orthomin(k)迭代法求解数值解所需的迭代步数。给出数值试验,通过两层逼近证明增量未知元是一个好的预处理子,尤其在结合一些迭代方法时。第三,在二维情形,使用两层增量未知元的θ格式求解依赖时间的对流扩散方程。证实关于二阶增量未知元的惯性流形多重网格算法是收敛的。增量未知元θ格式在0≤θ<1/2时是条件稳定的,在1/2≤θ≤1时是无条件稳定的。在时间层每一步迭代使用GMRES方法求解依赖时间的对流扩散方程,数值试验说明增量未知元方法能控制计算中带来的扰动。第四,提出一个修正的,基于一边差分的Crank-Nicolson格式,用于求解二维依赖时间的对流占优扩散方程。二阶一边差分逼近用于对流项的离散,二阶中心差分逼近用于扩散项的离散。该数值格式是相容的,无条件稳定的。对全离散格式,推导了先验误差估计。使用增量未知元预处理子之后,数值结果证实了修正的Crank-Nicolson格式的稳定性和有效性。接着,本文提出一个隐式差分格式求解上述类似问题。对常系数的对流占优扩散方程而言,基于一边有限差分和中心差分的隐格式是相容的并且无条件稳定,对应L2误差估计被推导。而且Burgers类型的方程使用类似半隐格式,其中只是对非线性项线性化,其时间上采用显式,这种半隐格式是条件稳定的。通过增量未知元方法的预处理技术,这种隐格式或半隐格式是一种有效格式,它能避免数值扰动,在稳定性和精度上表现和理论结果一致。同时,提出一种新的类小波增量未知元,它具有一些好的性质。一方面,近似解空间能被类小波增量未知元分解成两个L2正交子空间,这能使耦合系统自动消掉一些项,简化计算。另一方面,从类小波增量未知元到节点未知元的变化矩阵是正交矩阵,更方便计算。当应用于惯性流形多网格算法和非线性伽略金方法时,可证明类小波增量未知元方法的收敛性。最后,提出二维和三维类小波增量未知元方法,用于求解二维和三维具有多项式非线性项的反应扩散方程,推导类小波增量未知元的先验估计,证实如同所希望的那样,它是比较小的量,并且类小波增量未知元分解后的空间具有正交性。本文推导了基于类小波增量未知元方法的显格式和半隐格式的稳定性充分条件,这要优于相应标准算法的稳定性,即类小波增量未知元改进了算法稳定性条件。在二维情形,提供数值算例指出类小波增量未知元方法的有效性。并且,三维情形类小波增量未知元方法的稳定性可以类似二维情形进行推导。
其他文献
测度链上动力方程理论不但可以统一微分方程和差分方程,更好地洞察二者之间的本质差异,而且还可以更精确地描述那些有时随时间连续出现而有时又离散发生的现象.在探讨测度链上的动力方程的动力学行为时人们所熟悉的基本工具诸如Fermat定理,Rolle定理以及介值定理等不再成立,同时很难找到适应不同测度链的模拟程序,这些在给测度链理论研究带来诸多困难的同时,也更引起了广大学者的兴趣.本文首先考虑了测度链上一阶
双光子激光共聚焦荧光显微镜(Two-photon laser scanning fluorescencemicroscopy,TPM)因其具有独特的优势而被广泛应用于生命科学及医学领域。然而由于缺乏真正意义的双光子荧光探针,使双光子荧光显微镜的应用受到了很大程度上的限制。利用双光子荧光显微镜,长时间高空间分辨地活体观测核酸的生命活动是非常有意义的工作,然而目前尚无同时满足优良双光子特性和活体染色要
这篇博士论文着重研究了在吸收集不具有紧性的条件下,Banach空间上指数吸引子的存在性问题.并讨论了含任意p次多项式增长的非线性反应扩散方程指数吸引子的存在性.设{Sn}n=1∞为定义在Banach空间X上的离散半群,(?)为其吸引子.假设S在集合B∈0((?))上是C1的,并且S在集合B∈0((?))上任意一点的线性化算子可以分解为紧算子K与压缩算子C(||C||<λ<1)之和,即L=K+C.则
结构活性/性质关系方法(Structure Activity-Property Relationship,SAR/SPR)是目前国际上一个相当活跃的研究领域,近些年人们对该领域研究的投入呈现逐年递增的趋势。SAR/SPR方法的研究对象主要包含物质各种各样的物理化学性质参数,生物活性,毒性,以及药物的生物利用度等等,研究领域涉及化学、生物学、药学以及环境化学等诸多学科。该方法主要是从化合物的分子结构
在本博士论文中,首先我们在无界区域上考察了下面非自治反应扩散方程解的渐近行为这里(?)是一个N×N的实矩阵,并且具有正的对称项(?)(a + a*)≥βI,β> 0, a*表示a的转置,u = u(x,t) = (u1,...,uN),g=g(x,t)=(g1,...,gN),f=f(u,t)=(f1,...fN)..我们假定外力项g = g(x, t)∈Lb2(R; H),非线性项f = f(u
在这篇论文中,通过运用无穷维动力系统关于吸引子理论的最新研究成果并且结合一些能量估计技巧,我们研究了两类方程:具有衰退记忆的非经典扩散方程和具有衰退记忆的半线性热方程.我们对其弱解和强解的长时间动力学行为进行了深入的讨论,并且证明了以上方程对应动力系统的全局吸引子或一致吸引子的存在性.首先,我们研究了在自治外力项作用下具有衰退记忆的非经典扩散方程ut—△ut—△u—∫0∞k(s)△u(t—s)ds
众所周知,对于多个非奇异矩阵乘积的逆有如下的反序律成立:然而,当矩阵乘积A1A2…Am奇异时(此时,矩阵Ai可为奇异矩阵或长方形矩阵),这种所谓的反序律对于广义逆就不一定成立了.如何给出广义逆反序律成立的充要条件是矩阵广义逆理论中一个重要而又有趣的问题.假设Ai∈Cli×li+1,i=1,…,m为任意的m个复矩阵.本文利用广义Schur补的极大极小秩这一途径研究矩阵广义逆如下反序律成立的充要条件,
高等植物中,胚胎的发育关系到种子的形成,而种子的萌发又关系到物种的延续,这两个发育过程都起到至关重要的作用。拟南芥作为一种模式植物,是研究胚胎发育和种子萌发的理想材料。我们的研究采用分析拟南芥突变体表型的方法,将拟南芥肌动蛋白相关蛋白ARP6定位到这两个发育过程中。拟南芥的胚胎发育是一个复杂的基因调控过程,到目前为止,已经发现了不下于200个基因在胚胎发育过程中起作用。胚胎极性的建立,细胞分裂,母
本文主要研究了图的自同态幺半群的性质和结构,全文共分七章.第三章研究了路的补图的自同态幺半群,证明了路的补图的自同态幺半群是一个纯整半群,同时,解决了自同态幺半群的一些有关计数问题.特别地,确定了路的补图的自同态谱和自同态型.第四章研究了分裂图的联图的自同态幺半群,给出了此类图是自同态正则和自同态纯整图的充分必要条件,并且证明了分裂图的联图的自同态幺半群不可能是逆半群;同时,刻画了分裂图的联图的半
在数学生态学中,一个生态系统共存态的存在性及各种群的长时间行为是种群动态模型研究的主要内容.近年来,随着对种群动态模型研究的深入,同时考虑了扩散、自扩散和交错扩散作用的强耦合反应扩散方程组受到越来越多的关注.在这篇博士学位论文中,我们主要考虑种群动力学中的一类具有代表性的强耦合非线性反应扩散系统,通过对此类系统解的整体性态以及非常数稳态解存在性的研究,得到了一系列新的结果.整篇论文由五章组成.第一