论文部分内容阅读
本文我们主要考虑两类微分方程:时间依赖的不可压缩Navier-Stokes方程和二维空间Riemann-Liouville分数阶方程.对于Navier-Stokes方程,首先,我们引入一个辅助变量来分离扩散算子使原始的高阶方程转化为一阶系统,从而降低解决高阶方程的难度.其次,通过变分、精心地选择数值流、添加罚项,我们设计了一个稳定对称的局部间断Galerkin(LDG)格式.这样我们完成了方程的空间离散,但是我们并没有离散非线性对流项.我们知道求解Navier-Stokes方程最大的困难有两个:一是如何处理非线性对流项,另一个是如何处理方程中的压强函数.由于特征线方法在求解对流占优问题上有很大的优势,所以我们考虑特征线方法同时离散非线性对流项和时间导数项.从而,我们成功地解决了非线性问题,并且得到了一个稳定的对称的特征线局部间断Galerkin格式(CLDG).在稳定性推导的过程中,我们发现特征线方法使得理论推导简易明了.特别是应用格式的对称性,使得我们消去了一些比较难处理的算子.由于在NavierStokes方程中速度函数和压强函数没有紧密的联系,所以当我们得到速度的误差估计时,不是很容易得到压强的误差估计.为了解决这个困难,我们应用经典的速度和压强的连续上下界条件,建立速度误差和压强误差的联系.实际上,在间断有限元空间我们通常利用速度和压强的离散上下界条件,但是在本文,我们将连续上下界条件和离散空间结合起来得到压强的误差估计,这也是第三章的一个亮点.最后,我们给出四个不同的数值例子来验证理论结果,并且看到数值结果不仅达到了预期的效果而且比预期的更好.对于二维空间Riemann-Liouville分数阶方程,首先,我们引进两个辅助变量来分离Riemann-Liouville分数阶导数.由于Riemann-Liouville分数阶导数本身具有奇异性,Riemann-Liouville分数阶积分没有奇异性,所以一个辅助变量用来代替函数的梯度项,另一个辅助变量用来代替Riemann-Liouville分数阶积分.这样我们成功地分离了Riemann-Liouville分数阶导数并且把高阶导数方程降低为一阶系统.然后,通过变分、精确地选择数值流、添加罚项,我们设计了一个混合的间断Galerkin(HDG)格式,从而完成了方程的半离散.最后,我们给出三种时间离散方法,针对Riemann-Liouville分数阶扩散问题,我们应用一般的差分方法进行时间离散,由于方法的简易性,我们没有给出相应的离散过程.对于Riemann-Liouville分数阶对流扩散问题,我们应用一阶和二阶特征线方法同时离散时间导数项和对流项,并且给出相应的离散格式、稳定性分析、误差估计.在研究二维空间Riemann-Liouville分数阶方程时,我们发现如何找到一种有效的方法进行理论分析是比较困难的,特别是误差估计.为了解决这一困难,我们没有考虑有效的分析工具而是从格式设计出发,设计有效的数值格式.基于这样的考虑我们得到了不同的数值格式,不仅将一维问题推广到二维而且得到较好的数值结果.在第四章最后一节,我们分别应用三个数值例子验证Riemann-Liouville分数阶方程和一阶二阶HDG格式.总之,本篇论文我们成功地将间断Galerkin方法和特征线方法结合起来求解不可压缩的Navier-Stokes方程和二维空间分数阶微分方程.我们可以看到数值实验结果与理论结果是一致的并且是有效的.因此,这些方法可进一步改进且应用到其他问题中.