基于一种有限差分方法对复的修正Korteweg-de Vries方程孤子解的数值模拟研究

来源 :中国石油大学(北京) | 被引量 : 0次 | 上传用户:huajinxiu
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
Sasa-Satsuma(SS)方程是高阶非线性薛定谔方程的一种可积情形,可以用来描述飞秒光孤子在光纤中的传播和深水中的内孤波现象,复修正Korteweg-de Vries(MKdV)方程是SS方程的等价变换形式。本文首先基于Crank-Nicolson(CN)格式构造了复MKdV方程的一种守恒数值差分格式,并证明了该格式在时间和空间上都具有二阶精度,同时利用von Neumann方法证明了 CN型有限差分方法满足无条件的线性稳定性。其次,利用CN型数值差分方法对复MKdV方程的单孤子解的传播和二孤子解的碰撞进行了数值模拟。数值实验显示,单峰孤子和双峰孤子之间以及两个多峰孤子之间的相互作用都是标准弹性的,即两者在发生碰撞后均能保持速度和形状不变;双峰孤子和多峰孤子在相互作用之后双峰孤子变成了多峰孤子,但是碰撞前后孤子的质量是保持守恒的,因此这是具有形状改变特性的弹性碰撞。最后,本文从数值角度验证了差分格式的精度和线性稳定性,并得到了与理论分析一致的结论。
其他文献
石墨烯具有丰富的物理、化学性质,在透明导电薄膜、半导体等领域极具应用前景。目前,化学气相沉积(CVD)法由于操作简单、价格低廉并且可以大规模制备出高质量的样品,被认为是最可能实现石墨烯未来应用的制备手段。然而,CVD石墨烯总是存在着褶皱,褶皱作为一种线缺陷,将降解石墨烯的性能,影响薄膜的宏观均匀性。因此,如何彻底消除褶皱,实现无褶皱、超平整石墨烯的大规模制备是该领域的一个重点也是难点。此外,CVD
神经形态计算旨在模拟生物大脑和感知神经系统的信息处理方式,构建一个能够自主学习和具备认知功能的超低功耗智能计算机。突触、神经元和感觉接收器是生物计算系统处理信息的基本单元。因此,神经形态计算硬件实现的第一步就是研制出具有突触、神经元或感觉接收器功能的神经形态器件。新型神经形态器件凭借单个器件就能实现复杂的神经形态功能,其中电解质栅晶体管(Electrolyte Gated Transistors,
偏最小二乘法判别法是被广泛应用的一种高维数据降维分类预测判别方法,但在实际中,该方法通常会产生难以解释的原始变量的线性组合。在本文中把判别纳入到回归的框架之中。提出了一种带有惩罚项的稀疏型偏最小二乘法,结合矩阵的奇异值分解提出了一种秩为一矩阵的近似方法和算法。提出该种方法的目的在于把高维数据的变量筛选和变量的线性组合同时合并到一个步骤之中。尽最大可能地找出原始特征中包含的有用信息,实现数据的降维。
经过近数十年的发展,基于平均场近似的单参考量子化学方法已可对数十乃至上百个原子的一般化学分子进行较精确的结构、能量与性质计算。然而,对于强电子相关体系,如共轭分子和过渡金属配合物,由于存在大量能量近简并的π、d或f轨道,计算时需要对这些轨道组成的活性空间进行精确对角化,其计算量随体系的增大呈指数增长,因此目前的传统多组态和多参考量子化学方法尚难以处理超过18个活性轨道的体系。密度矩阵重正化群(DM
在系统辨识问题中,建立待辨识模型是一个重要的步骤。模块化非线性系统具有模型结构简单的特点,采用不同静态非线性块和动态线性块级联而成,其优点在于将系统的动态静态特性相分离,可以很好描述实际非线性系统。且设计时变模块化非线性系统和切换模块化非线性系统更加全面描述实际过程中的非线性和时变特性。由于模块化非线性系统包含一定的输入非线性或输出非线性,为其辨识带来一定难度,而且复杂噪声环境下,辨识结果存在偏差
研究背景:肾脏纤维化是肾脏衰老以及慢性肾脏病(Chronic kidney disease,CKD)的共有组织病理改变,是肾脏受到持续不断有害刺激过度修复的结果,以肾实质细胞向肌成纤维母细胞转分化(Myofibroblast trans-differentiation,MTD)为特征。肾细胞包括肾小管上皮细胞、血管内皮细胞、成纤维细胞和单核吞噬细胞受损后失去其自身表型并转分化为肌成纤维母细胞,产生
桥梁、道路等许多建筑的振动方程都是偏微分方程。通过有限元方法将这些方程离散为二阶微分方程,再分离变量,就可以得到二次特征值问题(QEP)。本文主要研究的是二次特征值问题中,单特征值以及半单重特征值的一阶导数的误差分析以及高阶导数的数值计算。特征值导数的计算一般来说比较简单,公式是确定的。但是由于特征向量的不唯一性,计算特征向量的导数相对比较复杂。单模态法是计算单特征值所对应的特征向量导数的一种比较
蛋白质作为生物功能的主要承担者,吸引着物理、生物和化学等各领域的广泛关注。而分子动力学模拟因其能同时提供精确的结构和即时的动力学信息,成为了一种极为重要的蛋白质研究方法。然而在现有的计算能力条件下,相对准确的全原子分子动力学模拟很难达到典型蛋白质功能动力学所涉及的时间尺度,这一定程度上限制了它的应用。为了克服这种困难,人们发展出了多种高效采样策略来提高分子动力学模拟的效率,其中包括:1)通过对蛋白
近年来,非局域非线性薛定谔(NNLS)方程在非线性数学物理和可积系统领域受到了密切关注。在本文中,我们首先利用Darboux变换方法和一些极限技巧,获得了自散焦NNLS方程的N阶有理解的行列式结构。然后,通过分析有理解中不同代数项之间的平衡关系,我们提出了改进的渐近分析方法以确定解在|t|→∞时的渐近表达式。基于该方法,我们推得了 1≤N≤4阶有理解的所有渐近孤子的精确表达式,并且发现渐近孤子的中
近十年来,随着团簇化学的蓬勃发展,研究人员开发了多种可控的合成方法来制备组成确定、结构精确的金属团簇。金属团簇的催化性能研究也越来越备受关注,在众多催化反应中表现出高的催化活性和选择性。金属团簇具有明确的组成和结构,可以在原子尺度上建立催化剂结构和催化性能的关系。同时,相比于单金属团簇,双金属团簇不仅具有好的热稳定性和催化活性,而且可以用来研究双金属催化协同作用,指导高性能双金属催化剂的理性设计。