几类偏微分方程的保结构算法研究

来源 :国防科学技术大学 | 被引量 : 1次 | 上传用户:zjs999
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
许多重要的数学物理方程都可以表示成Hamilton系统的形式,Hamilton系统内在具有能量守恒特性和辛几何结构。现代计算方法的基本原则是尽可能保持原问题的本质特征。因此,研究Hamilton系统框架下的偏微分方程,以及能够保持其能量守恒性及辛几何结构特征的数值方法是非常有意义的。本文对一些重要的非线性偏微分方程进行数值研究,并构造了一系列的保结构算法,同时给出了这些算法的离散守恒性质、收敛特性以及数值稳定性等必要的理论分析,保证了算法在长时间数值模拟中的可靠性。本研究主要内容包括:  ⑴对二维Zakharov-Kuznetsov方程构造了一种显式的多辛离散算法。该算法在空间方向采用多辛Fourier拟谱格式进行离散,时间方向采用辛Euler格式进行离散。证明了该多辛离散框架具有相应的离散多辛守恒律。数值实验表明该方法能够清晰地模拟孤立波的传播与碰撞,具有较高数值精度、不变量保持特性以及长时间数值模拟能力,相较于隐式方法,计算时间大大缩短。  ⑵对一系列的一维和二维Schr?dinger方程构造了若干高效的保结构算法。首先,对一种强耦合非线性Schr?dinger方程提出了一种多辛小波分裂算法,并对弹性碰撞和非弹性碰撞问题进行了数值模拟,结果表明了该算法针对强耦合问题在长时间数值计算中的优越性;其次,针对一维和二维变系数Schr?dinger方程,包括一维立方非线性Schr?dinger方程,Gross-Pitaevskii方程,二维时变Schr?dinger方程,结合分裂技术、小波配点离散和Fourier拟谱离散的思想,提出了两种新颖的守恒格式,并分析了算法的误差精度,严格证明了方法能够保持电荷守恒,并且在一定条件下能够保持动量和能量守恒;再次,针对二维非线性Schr?dinger方程,在不考虑边界条件的情况下,提出了一种局部能量守恒算法和一种局部动量守恒算法,证明了算法理论上可以在任何局部时空区域精确保持离散的局部能量(动量),在周期边界条件下,算法仍然保持电荷守恒和全局守恒律,并分析了算法的守恒性、稳定性及误差估计;最后,针对多变量耦合非线性Schr?dinger系统,提出了一种构造高精度保能量算法的简单框架,它们分别在空间方向基于高阶紧致差分方法、Fourier拟谱方法以及小波配点方法进行离散,严格证明了该框架能够保持离散电荷和能量守恒律。  ⑶对量子物理中非线性Dirac模型,构造了一种保能量的小波算法和一种多辛小波算法,证明了两种算法分别可以保持相应的离散守恒特性和几何结构。为了提高计算效率,我们利用分裂技术和显式加速方法对算法进行改造。数值结果表明所有算法都在长时间计算过程中都比较稳定,经过加速改造后的算法效率更高。  ⑷基于Degasperis-Procesi方程的双曲—椭圆耦合的特征,构造出一种双曲—椭圆分裂步算法。把Degasperis-Procesi方程分裂成一个Burgers方程和一个Benjamin-Bona-Mahony方程,对于Burgers方程,利用经典有限体积型本征非震荡格式进行计算,对于Benjamin-Bona-Mahony方程,利用多辛Fourier拟谱方法进行离散。数值结果表明该算法能够较好地处理间断,且具有长时间数值模拟的能力。
其他文献
设G是由有限维复单李代数£与其伴随表示确定的伴随型Chevalley-Demazure群概形([9],[10]),G(Z)是整数环Z上的Chevalley群.设E(Z)是G(Z)的初等子群,U是E(Z)的幺幂子群.该文的目
该论文的第一项工作是,在大量测试研究的基础上,发现Huang提出的EMD(Empirical Mode Decomposition)方法存在尺度分辨率的限制.对于一个信号来讲,这种限制直接导致在判断产生
通信的目的是要把对方不知道的信息及时可靠的传递给对方.所以,一个通信系统传输消息必须要可靠及快速.但是可靠与快速却又往往相互矛盾,纠错码理论就是在解决这对矛盾的过程中
James Stirling于1730年首先提出了Stirling数对的概念,Thiele和Nielsen(1904)正式运用了这一名称.第一类Stirling数s(n,k)和第二类Stirling数S(n,k)是Bell多项式的特殊形式.
该文仅考虑Hilebert空间L(R)及其闭子空间中的小波变换和小波框架等问题.该文主要考虑的问题是:L(R)上的连续小波变换,正尺度小波变换和s-进小波变换,以及L(R)中的小波框架,
在这篇论文中,我们讨论了Bell多项式中的递归关系和同余问题,给出了广义Stirling数对(GSN对)的一类递归关系式,提出了扩展的GSN对G(n,k|w,u)和(Gn,k|u,w)的概念,统一了二项式