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

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

5次非线性Schrdinger方程的一个线性化4层紧致差分格式

对5次非线性Schrdinger方程提出了一个线性化4层紧致有限差分格式,引入"抬升"技巧,运用标准的能量方法和数学归纳法建立了误差的最优估计,证明数值解在空间和时间2个方向分别具有4阶和2阶精度.数值实验对理论结果进行了验证,并通过对比表明该文格式在保持精度相当的前提下较已有格式具有更高的计算效率.     

一类四阶非线性抛物方程的紧致差分格式

针对一类四阶非线性抛物方程的初边值问题建立紧致差分格式,利用降阶的思想,通过引入中间变量将原四阶问题转化成二阶非线性方程组.对方程中的时间导数项和空间导数项分别采用Crank-Nicolson格式和四阶紧致差分格式进行离散,对非线性项采用外插的方法进行处理,从而得到原问题的三层线性紧致差分格式,其局部截断误差为■.数值算例表明该格式具有良好的计算效果.基于四阶非线性抛物方程在薄膜理论等问题中的重要作用,对此类方程构造高精度的紧致差分格式,可以使该方程在有关工程计算方面得到更好的应用,因此该研究成果具有重要的理论意义和广泛的应用前景.     

KdV方程的一个空间六阶空间精度守恒差分格式

本文对含齐次边界条件的KdV方程的初边值问题进行了数值研究. 通过在时间层进行二阶精度的Crank-Nicolson差分离散、在空间层进行六阶理论精度的外推组合差分离散，本文建立了一个具有六阶空间精度的两层非线性差分格式. 该格式能够合理地模拟原问题的两个守恒量. 然后，本文利用能量方法证明了格式的收敛性和稳定性. 数值算例验证了该方法的有效性.     

求解二维热传导方程的高精度紧致差分方法

基于Richardson外推法提出了一种数值求解二维热传导方程的高阶紧致差分方法.该方法首先利用时间二阶、空间四阶精度的紧致交替方向隐式(ADI)差分格式在不同尺寸的网格上对原方程进行求解,然后利用Richardson外推技术外推一次,最终得到了二维热传导方程时间四阶、空间六阶精度的数值解,数值实验验证了该方法的高阶精度及有效性.     

数值求解一维波动方程的四阶紧致差分方法

针对一维波动方程提出了一种有限差分方法.首先,采用泰勒级数展开公式和原方程代入的方法推导出了第一个时间层未知函数值的四阶紧致差分格式.然后,用四阶紧致差分公式近似空间导数项,采用中心差分格式截断误差余项修正的方法处理时间导数项,推导出了第二个时间层以后未知函数的四阶紧致差分格式.该方法时间和空间具有整体四阶精度.利用Fourier方法分析了所提格式的稳定性.由于本文格式在未知时间层仅涉及3个网格点,因此可采用追赶法求解离散化后所得到的线性方程组.最后,用数值算例验证了本文格式的精确性和稳定性.     

求解Klein-Gordon方程的新型快速紧致时间积分方法

构建了基于Hermite插值的快速紧致时间积分方法求解Klein-Gordon方程.该方法先在空间方向上采用四阶紧致差分格式离散得到了一个半离散格式.然后结合离散正弦变换和常数变易公式给出了半离散格式之解的显示时间积分表示式,并对积分中的非线性源项采用Hermite插值逼近,得到了一个全离散格式.仅需利用前两个时间步的计算结果,就可获得空间和时间方向均为四阶精度的高效算法.数值模拟的结果验证了该方法的有效性.     

一类四阶发展方程的混合有限体积格式

针对一类带对流项的四阶发展方程,运用特征线方法,直接从积分守恒形式出发导出了混合有限体积格式,给出了格式的截断误差,证明了该格式按离散L^2模对时间和空间具有一阶精度．数值算例表明,该格式计算效果良好．     

求解非线性时间分数阶Klein-Gordon方程的谱配置方法

用Jacobi谱配置方法, 数值求解一类非线性时间分数阶导数为Caputo导数的Klein-Gordon方程. 先用Caputo分数阶导数和Riemann-Liouville分数阶积分的关系, 将分数阶Klein-Gordon方程转化为在时间上带奇异核的积分微分方程, 再在时间和空间上采用Jacobi谱配置法, 并用高斯积分公式逼近积分项, 使方程在配置点上 成立, 从而求得其数值解. 数值算例结果表明, 该方法所得数值解很好地逼近了精确解.     

双曲型守恒律的一个4阶精度差分格式

基于TVD限制器函数方法选取数值导数,在空间方向用分段3次多项式进行重构,对时间积分用Simpson求积公式,并用四阶Runge-Kutta NCE方法求中间时间点的值,得到求解一维非线性双曲型守恒律方程的4阶精度差分格式;之后给出2个经典数值算例,以验证格式的高精度高分辨率优点。     

2维Gross-Pitaevskii方程的分裂高阶紧致差分格式

该文为带有旋转角动量的Gross-Pitaevskii方程构造了分裂高阶紧致差分格式.首先通过时间分裂将其分为线性方程和非线性方程,非线性方程可以通过质量守恒定律进行精确求解,线性方程通过高阶紧致格式和局部1维方法进行离散,最终得到的格式时间方向2阶收敛和空间方向4阶收敛,并保持质量守恒.最后用数值算例验证了格式的收敛阶以及质量守恒性.     

第三类Dirichlet边界下四阶抛物方程的紧差分格式

应用加权平均和高次Hermite插值等技术,提出逼近四阶导数的几个有用的数值微分公式,并对其截断误差进行分析。在此基础上建立求解第三类Dirichlet边界条件下四阶抛物方程初边值问题的三个高阶紧差分格式,应用Fourier分析方法证明格式的无条件稳定性,并对其进行数值验证。这三个差分格式的差异主要体现在空间导数临近边界处的离散方式不同,所得格式全局精度均达到了时间二阶、空间四阶。     

污染物对流扩散方程的预测校正紧差分格式

结合预报校正线性多步法与高阶紧致差分格式方法的优点,空间导数采用四阶紧致差分格式进行离散之后,对得到的空间半离散格式采用改进的预报校正的线性多步法进行时间推进,得到一种时空方向均为四阶精度的求解非线性对流扩散方程的高精度方法。数值试验表明该格式可以有效求解非线性对流扩散方程,验证了格式的良好性能。     

非线性反应扩散方程的两重网格算法

将两重网格算法和混合有限元方法结合起来,通过对二维非线性反应扩散方程右端的非线性项进行基于粗网格解的泰勒展开,化为细网格上的线性问题,从而为求解该类方程提供了一种有效的数值解法。收敛性分析和数值算例结果表明,该算法与标准有限元方法相比,其优点是在不降低解的精度阶数的条件下,提高了计算速度,同时能够得到精度更高的导数。     

应用载荷增量—最小二乘配点法的圆柱形壳体几何非线性分析

本文采用载荷增量－最小二乘配点分析圆柱壳体和纵、横向加筋圆柱壳体和纵、横向加筋圆柱壳的几何非线性问题。首先使用载荷增量法将两种圆柱形壳体的非线性微分方程线性化。然后使用最小二乘配点法求解这些离散后的线性微分方程。本文对几种典型的边界条件和载荷状况下的壳体编制了通用程度进行计算。结果表明本方法是有效和可靠的。     

一种求解对流扩散反应方程的高阶紧致差分格式

提出一种消除对流扩散反应方程中对流项的处理技巧,结合中心差分格式的新方法与相同节点的迎风差分格式相比具有更好的精度,该方法很容易与Padé格式相结合,构造出具有四阶精度的无条件稳定的高阶差分格式.数值实验表明,新方法具有很好的精度和健壮性,并且可以有效求解对流占优问题.     

立铣加工切削力和振动的计算机仿真与实验

在改进的非线性立铣加工铣削力数学模型和动力学方程基础上,采用变步长数值积分算法(四阶显式Runge-Kutta算法),建立动态铣削加工过程的计算机仿真模型.通过铣削加工动力学实验,对仿真模型输出的铣削力和铣削振动的时域特性预测精度进行验证,进一步对其频域特性进行分析.结果表明,该仿真模型作为高效低耗的研究平台,可较好地用于立铣加工铣削力与铣削振动的预估及其频域特性分析.     

隔水管非线性振动的线性化分析

本文在规则波浪与海流联合作用下对海洋钻井隔水管进行了动力分析,对非线性拖曳力采用改进的线性化方法,推导了有限元线性振动方程。求解这些方程几乎象静力方程一样,只需要少数几次迭代就可求得响应最大值。同时给出一些规则波浪与海流联合作用下的算例,并对由本文方法和由非线性时间域动力分析方法分别求得的两种结果进行了比较。比较表明,两种计算结果吻合。本文的线性化方法及把动力方程变换为静力方程的方法,能节约计算机时。隔水管是海洋工程中较典型的构件,该分析方法便于实际应用。     

一种紧致差分格式及其在新型纺纱研究上的应用

本文导出了一种与涡传输方程的ADI格式相协调的流函数方程的四阶精度五点紧致差分格式,并导出了与此相容的速度差分方程及三阶精度的壁涡公式,以及在非网格边界点上的壁涡公式与相邻内点的涡传输差分方程和流函数差分方程。本文成功地用此方法对某新型纺纱器罩壳内的流动作了数值计算。其结果与流场显示法的实验结果(照片)基本符合,对新型纺纱的理论研究以及罩壳的合理改型提供了可靠的依据。     

抛物型最优控制问题的三次B样条有限元方法

对一类四阶非线性抛物方程最优控制问题提出一种三次B样条有限元方法。状态变量和对偶状态变量用具有更好光滑性的分片三次B样条连续函数进行逼近,控制变量由分片常数函数进行逼近。这样得到的状态变量和对偶状态变量的数值解二阶连续可微。建立最优性系统的全离散格式,并用迭代法进行求解。最后建立数值算例,验证方法的有效性。