【摘 要】
：
本文主要对Schro（?）dinger方程保持质量守恒的DDG（directdiscontinuousGalerkin）方法进行了研究.在第三章,对线性一维及二维Schro（?）dinger方程的守恒数值流进行了选取,对于其半离散格式和全离散格式进行了守恒性分析.其中,在空间离散上采用DDG方法,在时间离散上采用Crank-Nicolson格式.数值实验表明,对于k次元逼近能得到k+1阶精度,而且
论文部分内容阅读
本文主要对Schro（?）dinger方程保持质量守恒的DDG（directdiscontinuousGalerkin）方法进行了研究.在第三章,对线性一维及二维Schro（?）dinger方程的守恒数值流进行了选取,对于其半离散格式和全离散格式进行了守恒性分析.其中,在空间离散上采用DDG方法,在时间离散上采用Crank-Nicolson格式.数值实验表明,对于k次元逼近能得到k+1阶精度,而且在长时间的计算后,同样能得到最优收敛阶.实验结果显示,对于取不同的时间,不同的剖分进行计算,都能得到相同的质量,这表明质量是守恒的.在第四章,讨论了势函数取两种情况时的非线性Schro（?）dinger方程,对其守恒数值流进行了选取,并对于其半离散格式和全离散格式进行了守恒性分析.其中,和线性方程不同的是,在空间离散上还是采用DDG方法,而在时间离散上采用Strangsplitting方法.数值实验表明,对于k次元逼近同样能得到k+1阶精度,而且在长时间的计算后,同样能得到最优收敛阶.值得注意的是,在二维情况下,在用偶数次多项式逼近的时候,数值流中的β0[u]/h起到了稳定性的作用.而用奇数k次多项式逼近的时候,即使β0=0,也能达到最优收敛阶k+1阶.最后实验结果显示,质量是守恒的.在第五章,对于线性Schro（?）dinger方程半离散格式进行了误差分析.首先提出了一种全局投影的定义,引进全局投影主要的想法,就是因为界面条件是受到数值流选择的限制的,要引进的全局投影,它不仅能满足相应的界面条件,而且能使得在节点处的“麻烦”项,能够被消除或者被控制住.接着对于全局投影的存在性进行了分析,表明满足一定条件的全局投影是唯一确定的.最后给出了线性Schro（?）dinger方程半离散格式的误差估计式.而对于二维情况的误差分析,如果继续采用一维的方法,则有些“麻烦”项不能被消除或者被控制住,这需要考虑另外的解决方法,这也是正在继续的工作.和传统的Schro（?）dinger方程的保守恒格式相比,我们的方法不仅具有质量守恒性,高精度性和长时间的稳定性,而且不用将方程进行变形求解,大大缩短了计算时间.
其他文献
热障涂层（Thermal Barrier Coatings,TBCs）因为耐高温、高隔热、抗腐蚀等优异性能,已成为航空发动机涡轮叶片等热端部件不可缺少的热防护材料。然而,在含有杂质颗粒的高温、高速燃气环境下服役时,会因为界面氧化、颗粒冲蚀、CMAS腐蚀（钙镁铝硅等金属氧化物的混合物,简称CMAS）等多种因素导致涂层剥落。其中,CMAS腐蚀指高温下钙镁铝硅混合氧化物熔融、渗透到涂层中,引起涂层结构、
本文利用谱配置方法求解一类带弱奇异核的和非线性的Volterra型积分微分方程,并且构造高精度算法,着重分析该方法的误差估计和收敛性,并进行数值实验,验证所给方法的有效性。第二章,用Legendre谱配置方法对非线性Volterra-Fredholm-Hammerstein积分方程进行求解。并且对于所给的数值格式和误差进行了分析,即当核函数充分光滑时,计算所得数值格式的L2范数和L∞范数误差呈指数
半拉格朗日（半拉氏）方法没有CFL稳定性条件限制且比欧拉方法更加稳定。在文献当中,水平集重新初始化方程的半拉氏方法是复杂的,这是没有必要的。因为重新初始化过程是辅助的,我们提出了耦合有改善紧邻交界面网格点处精度的投影技术的一阶半拉氏格式。标准的二阶半拉氏格式被用于推进水平集对流方程。实现是简单的,包括在块结构自适应网格上。半拉氏方法的有效性通过大量的数值算例被证实,其中包括在给定速度场中做被动对流
复杂网络已经吸引了来自科学技术不同领域研究者的大量关注。继复杂网络中的小世界特性与无标度性质之后,自相似性已经成为复杂网络的第三大基本特征并且在近年来得到了广泛的研究。分形分析可以有效地揭示一些分形对象的自相似性。然而,重分形分析是一个用来系统性地刻画理论和实验复杂分形对象空间异质性的更加强大而有效的工具。尽管已有一些研究复杂网络自相似性的分形和重分形分析算法,但是这些算法的效率并不高而且已有的算
本文主要研究了阶数为α∈ （0,1）的具有Caputo导数的时间分数阶Navier-Stokes方程,这类方程可以用来模拟分形介质中的反常扩散现象.我们主要讨论它的适定性.本文的工作主要分成五个部分:在第二章中我们首先通过把Helmholtz投影子作用到目标方程上消去压强项,将其转化为抽象形式的发展方程;然后根据分数阶发展方程理论给出其相应的积分方程并定义适度解;利用半群理论和不动理论得到方程的全
本论文主要研究了M-模糊化凸结构理论,（L,M）-模糊凸结构理论,及（L,M）-模糊凸结构（子）范畴理论.在M-模糊化凸结构理论中,首先,我们引入并刻画了M-模糊化（弱）JHC凸空间、M-模糊化几何区间空间、M-模糊化（弱）Peano、Pasch、沙漏区间空间,并讨论了它们之间的关系。其次,我们引入了M-模糊化基点序的概念,并用其刻画了 M-模糊化几何性质、M-模糊化Peano性质、M-模糊化Pa
对于非线性扩散占优Volterra型偏泛函微分方程的时间离散化,以往最常用的方法是Volterra泛函微分方程（VFDE）的隐式Euler方法及梯形方法,有时候也使用隐式中点法及Bellen和Zennaro所建议的二级Lobatto IIIC Runge-Kutta方法。为了避免出现Bellen和Zennaro所指出的“order failure”和“stability failure”现象,本文
分数微分方程（Fractional Differential Equations）在许多学科领域有广泛的应用,这些学科领域的许多数学模型都是用分数微分方程来描述的.近十年来,分数微分方程边值问题得到了迅速的发展.随着这一方向研究的进一步深入,研究内容和研究方法不断得到丰富,非线性分数微分方程边值问题获得了许多研究成果.本文主要研究非线性分数微分方程边值问题解的存在性.本文由四章组成,主要内容如下：
本文研究二维各向异性Delaunay三角化网格的生成与优化算法及其应用.当物体在不同方向上的变化速率不同时,例如树木向上生长与横向生长的速率是不同的,我们称其具有各向异性性质.各向异性问题在各向异性晶体材料、多孔介质、石油工程、粘性流模拟、图像处理等科学工程计算中有着广泛的应用.给出一个各向异性问题,进行有限元数值求解时,合适的网格剖分至关重要,选择符合各向异性特征的三角化网格是自然而然的想法.各
有限体积法由于能保持某些物理量（如质量、能量等）的局部守恒性,它已成为一种广泛应用于科学与工程计算领域的重要数值方法.有限体元（FVE）法是一种重要的有限体积法,虽然对它已有大量的研究工作,但仍存在许多需要研究的问题.ICF问题描述的是一种聚变等离子体流体力学问题,其中三温辐射热传导方程的数值求解在该问题的数值模拟中占有主要工作量,因此,对该类问题的有限体元格式及其快速算法进行研究,是一项具有重要