分数阶Klein-Gordon-Schrödinger方程的保能量方法

该文先将分数阶Klein-Gordon-Schrödinger方程转化成辛结构的哈密尔顿系统,利用傅里叶拟谱方法对Riesz空间分数阶导数进行近似离散,得到分数阶Klein-Gordon-Schrödinger方程有限维哈密尔顿系统; 再利用2阶平均向量场方法对有限维哈密尔顿系统离散,得到分数阶Klein-Gordon-Schrödinger方程新的保能量格式; 最后利用新的保能量格式数值模拟方程孤立波的演化行为,并分析新格式的保能量守恒特性.     

分数阶Klein-Gordon-Zakharov方程新的保能量格式

将分数阶Klein-Gordon-Zakharov方程转化成多辛结构的偏微分方程,利用傅里叶拟谱方法对方程Riesz空间分数阶导数进行近似离散,得到有限维的常微分方程组,再利用二阶平均向量场方法对常微分方程组离散,得到方程新的保能量格式,最后利用新格式数值模拟分数阶Klein-Gordon-Zakharov方程孤立波的...     

空间分数阶扩散方程的有限元方法

考虑空间分数阶扩散方程的数值解,利用有限元的思想构造了一个半离散数值格式,并严格证明格式的收敛性分析,数值例子支持理论分析的结果.     

变时间分数阶扩散方程的数值模拟

考虑变时间分数阶扩散方程。首先利用分段线性插值法结合对一阶时间导数的一个二阶近似离散Coimbra变时间分数阶导数,然后利用Richardson外推法改进精度,最后用数值例子来验证提出的数值方法,从而说明数值方法的有效性。     

空间分数阶Edwards-Wilkinson方程的数值研究

为了探讨空间分数阶随机生长模型的动力学标度行为,利用Grümwald-Letnikov分数阶导数定义方法求解空间分数阶Edwards-Wilkinson（SFEW）方程在1＋1维情况下的数值解,得到了在不同分数阶导数值时的生长指数、粗糙度指数、动力学指数和局域粗糙度指数,这些结果与标度分析得到的结果是一致的。研究结果表明SFEW模型没有出现奇异动力学行为,仍然遵守Family-Vicsek正常标度规律。同时结果也显示,非局域相互作用对SFEW方程的动力学标度行为有着显著的影响。     

分数阶电报方程的近似解析解与数值解

研究了分数阶电报方程的近似解析解与数值解.首先用Adomian拆分法讨论了它的近似解析解;其次用差分法求解它的数值解,构造出隐式差分格式;最后给出数值例子,把近似解析解、数值解与精确解进行了比较,显示方法是有效的.     

非线性分数阶演化方程的新解

通过使用改进的分数阶sub-equation 方法寻求一些非线性分数阶演化方程的精确解, 如分数阶Burgers 方程、耦合分数阶Burgers 方程与非线性分数阶Klein-Gordon 方程等, 并得到了这些非线性分数阶演化方程的新解.     

时间分数阶薛定谔方程的数值方法

结合非标准有限差分格式给出了求解分数阶薛定谔方程的一种数值解法,对时间导数离散后的分母构造了一个关于时间步长的函数来近似,证明了该差分格式是无条件收敛和稳定的.数值算例表明该方法不仅有非常好的收敛性和稳定性,还有较高的精度,因此该方法是有效的.     

分数阶非线性Schrodinger方程的守恒算法

利用傅里叶谱方法对空间分数阶非线性Schrodinger方程进行数值求解,并证明该格式保持了能量和质量的守恒性且无条件稳定。该方法在空间方向具有谱精度,在时间方向具有二阶精度。还对该格式进行误差分析及收敛性分析。最后通过数值实验验证了该算法的守恒性、准确性和有效性。     

多重网格法求解原始变量形式的Navier-Stokes方程

采用多重网格法 ,求解原始变量形式的Navier Stokes方程 .控制方程在交错网格上离散 ,并采用SIMPLE算法计算 .为加速收敛 ,结合使用适合非线性方程的全近似格式的多重网格法FAS .文章计算了雷诺数Re =10 0 ,10 0 0时不可压粘性流体绕圆柱的流动 .计算结果与实验结果和国外的计算结果基本吻合 ,说明了该方法的准确性和有效性     

复修正KdV方程的高阶保能量方法

利用4阶平均向量场方法和拟谱方法构造了复修正KdV方程的高阶保能量平均向量场格式,并利用构造的高阶保能量格式数值模拟了方程孤立波的演化行为.数值结果表明:构造的4阶格式具有好的稳定性,可以很好地模拟孤立波的演化行为,并且精确保持方程的能量守恒特性.     

乙烯裂解炉管中时变应力场和蠕变损伤分数数值模拟

通过对乙烯裂解炉管典型工况分析,建立了炉管在温度、内压、渗碳和蠕变等交互作用下的有限元计算模型,并对HK40奥氏体耐热钢制成的乙烯裂解炉管在1173K下的应力场和蠕变损伤分数进行了模拟;结果表明渗碳是引起炉管失效的主要原因。通过计算得到的炉管失效方式、位置和炉管寿命均与实际运行结果基本吻合。     

Banach空间中分数阶微分方程初值问题解的性质

讨论了Banach空间中的分数阶微分方程解的性质,利用Schauder不动点定理及Gronwall不等式证明了初值问题解的存在唯一性.当右端函数f(t,u)关于u线性增长时,得到了解的整体存在性.进一步讨论了分数阶方程的解对初值和阶数的连续相依性.     

用无网格法求解不同Re下圆柱绕流问题

采用无网格伽辽金法对不同Re下二维不可压粘性流场圆柱绕流问题进行了数值模拟。该计算不仅预测了不同条件下的流动特性,而且预测了流体作用在圆柱上的曳力和升力。这些力通过积分圆柱表面的压力和粘性力获得。结果显示,当Re大于80时,圆柱所受曳力和升力开始振荡。计算所得的脱落频率(Sr)与前人的结果吻合良好。     

基于波能守恒方程模拟近岸波流场

基于波能守恒方程对平直斜坡上近岸波浪场以及近岸波浪斜向入射破碎后所产生的沿岸流场进行数值模拟.鉴于平直斜坡地形较简单,便于简单计算,利用基于波能守恒方程建立的波浪场数学模型计算得到波高的变化,从而利用Longuet-Higgins计算辐射应力.在此基础上,利用近岸波流场数学模型对实验波况下的波流场进行模拟,并将数值模拟...     

一阶线性变系数对流方程的稳定性分析

采用能量不等式方法并结合冻结系数法,对一阶线性变系数对流方程迎风格式的稳定性进行研究,得到迎风格式稳定性的条件,并通过数值实验验证稳定性分析理论的正确性。     

长直椭柱导体表面静电场的数值计算

依据最低能量原理,将椭柱导体表面的电荷离散成有限电荷微元,通过电荷微元的坐标轮动,寻求体系的最低能量位形,得到了体系达到静电平衡时的电荷密度分布,进而依据高斯定理描述了导体表面的电场分布。