共查询到20条相似文献,搜索用时 15 毫秒
1.
提出由第三类生成函数法构造高阶Schroedinger方程δu/δt=i(-1)^nδ^2mu/δx^2m的高精度辛格式.首先,给出它的典则Hamilton方程组;然后,成功地克服了本质上是困难的高阶变分导数的计算,并利用第三类生成函数法得到在时间方向具有任意阶精度的半离散方程,进而得到原始方程相关的修正方程的离散形式,最后得到各种精度的辛格式.数值结果表明该格式是有效的,具有高精度及良好的长时间数值行为等特性. 相似文献
2.
提出一种基于WENO重构的高阶(至少三阶)移动网格动理学格式.利用流体力学方程的积分形式得到移动网格上离散格式,再利用自适应移动网格方法移动网格,进而得到网格速度,利用WENO重构得到高阶插值多项式,最后使用时间方向上精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高精度、高分辨率的特点. 相似文献
3.
4.
给出对流扩散方程的一种高精度交替分组8点格式,可以用于并行计算,且无条件稳定.数值实验证实此格式具有高阶精度. 相似文献
5.
提出了一种新的算法——高阶辛时域有限差分法(SFDTD(3, 4): symplectic finite-difference time-domain)求解含时薛定谔方程.在时间上采用三阶辛积分格式离散, 空间上采用四阶精度的同位差分格式离散, 建立了求解含时薛定谔方程的高阶离散辛框架;探讨了高阶辛算法的稳定性及数值色散性.通过理论上的分析及数值算例表明:当空间采用高阶同位差分格式时, 辛积分可提高算法的稳定度;SFDTD(3, 4)法和FDTD(2, 4)法较传统的FDTD(2, 2)法数值色散性明显改善.对二维量子阱和谐振子的仿真结果表明: SFDTD(3, 4)法较传统的FDTD(2, 2)法及高阶FDTD(2, 4)法有着更好的计算精度和收敛性, 且SFDTD(3, 4)法能够保持量子系统的能量守恒, 适用于长时间仿真. 相似文献
6.
双曲型守恒律的一种高精度TVD差分格式 总被引:3,自引:0,他引:3
构造了一维双曲型守恒律方程的一个高精度高分辨率的守恒型TVD差分格式.其主要思想是:首先将计算区域划分为互不重叠的小单元,且每个小单元再根据希望的精度阶数分为细小单元;其次,根据流动方向将通量分裂为正、负通量,并通过小单元上的高阶插值逼近得到了细小单元边界上的正、负数值通量,为避免由高阶插值产生的数值振荡,进一步根据流向对其进行TVD校正;再利用高阶Runge KuttaTVD离散方法对时间进行离散,得到了高阶全离散方法.进一步推广到一维方程组情形.最后对一维欧拉方程组计算了几个算例. 相似文献
7.
提出了一种新的算法—高阶辛时域有限差分法(SFDTD(3,4):symplectic finite-difference time-domain)求解含时薛定谔方程.在时间上采用三阶辛积分格式离散,空间上采用四阶精度的同位差分格式离散,建立了求解含时薛定谔方程的高阶离散辛框架;探讨了高阶辛算法的稳定性及数值色散性.通过理论上的分析及数值算例表明:当空间采用高阶同位差分格式时,辛积分可提高算法的稳定度;SFDTD(3,4)法和FDTD(2,4)法较传统的FDTD(2,2)法数值色散性明显改善.对二维量子阱和谐振子的仿真结果表明:SFDTD(3,4)法较传统的FDTD(2,2)法及高阶FDTD(2,4)法有着更好的计算精度和收敛性,且SFDTD(3,4)法能够保持量子系统的能量守恒,适用于长时间仿真. 相似文献
8.
给出了柱坐标下Euler方程数值边界条件的一种处理方法.径向第一个网格点设在距离中心半点位置上.根据相应物理量的特性,在中心附近进行边界延拓,使得内点的高精度差分格式可以同样应用在网格中心附近,从而无需单侧差分格式,保持了一致的高阶精度.对于周向边界,也建立了一种周期延拓方法,使得在周向所有节点处都能够采用同样的高精度格式离散,并进行了数值试验. 相似文献
9.
从一阶麦克斯韦旋度方程出发,研究一种区域分解时域有限元目的——高阶间断伽辽金时域有限元目的.其中对时间的离散采用Crank-Nicolson差分格式,电场和磁场采用相同阶数的高阶矢量基函数展开.分析三维谐振腔问题,数值结果表明,目的 中时间步长的选取可以摆脱CFL稳定性条件的限制;此外,与基于常用Whitney矢量基函数的目的 相比,采用高阶矢量基函数可以明显地提高计算精度及计算效率. 相似文献
10.
提出了一种新的算法一高阶辛时域有限差分法(SFDTD(3,4):symplectic finite—difference time-domain)求解含时薛定谔方程.在时间上采用三阶辛积分格式离散,空间上采用四阶精度的同位差分格式离散,建立了求解含时薛定谔方程的高阶离散辛框架;探讨了高阶辛算法的稳定性及数值色散性.通过理论上的分析及数值算例表明:当空间采用高阶同位差分格式时,辛积分可提高算法的稳定度;SFDTD(3,4)法和FDTD(2,4)法较传统的FDTD(2,2)法数值色散性明显改善.对二维量子阱和谐振子的仿真结果表明:SFDTD(3,4)法较传统的FDTD(2,2)法及高阶FDTD(2,4)法有着更好的计算精度和收敛性,且SFDTD(3,4)法能够保持量子系统的能量守恒,适用于长时间仿真. 相似文献
11.
12.
13.
提出一种求解二维功能梯度材料(FGMs)稳态热传导问题的重心Lagrange插值配点法.基于Chebyshev节点构造二维重心Lagrange插值函数及其偏导数,然后基于配点法将其直接代入FGMs热传导问题的控制方程和边界条件,得到系统离散方程.重心Lagrange插值配点法是一种真正的无网格方法,很好地融合了重心Lagrange插值和配点格式的优势,具有高效、稳定、高精度和易于数值实现的优点.采用重心Lagrange插值配点法分别对指数型、二次型和三角型FGMs热传导问题进行数值模拟.结果表明:该方法具有较高的计算效率和计算精度,对材料梯度参数的变化不敏感.可以进一步拓展到FGMs瞬态问题和FGMs的热力耦合分析. 相似文献
14.
15.
16.
高阶CIP数值方法及其在相关物理问题中的应用 总被引:1,自引:0,他引:1
利用函数的高阶空间导数值构建其高次插值,得到高阶CIP(Constrained Interpolation Profile)数值算法,并在此基础上模拟研究等离子体物理中著名的伏拉索夫-泊松(Vlasov-Poisson)方程相关物理问题.高阶CIP数值方法具有更高数值精度,从而可以在同等精度的情况下减少计算格点数,加速数值计算速度. 相似文献
17.
假设爆轰波阵面的法向速度是曲率的线性函数,在非结构四边形网格上采用水平集方法模拟爆轰波阵面的运动过程.水平集方程的曲率无关项采用正格式离散,曲率项采用伽辽金等参有限元方法空间离散,时间离散采用半隐格式.在笛卡儿网格和随机网格上,含曲率的水平集方程的离散格式为强一阶精度,重新初始化方程的离散格式精度为近似一阶精度.曲率收缩的不光滑界面和多个爆轰波阵面相互作用的算例说明格式可有效地模拟爆轰波与曲率相关的运动. 相似文献
18.
重心Lagrange插值配点法求解二维双曲电报方程 总被引:1,自引:0,他引:1
提出一种求解二维双曲电报方程的高精度重心Lagrange插值配点法.采用重心Lagrange插值构造包含时间和空间变量的近似函数.在给定Chebyshev-Gauss-Lobatto节点上,将多变量重心Lagrange插值近似函数代入双曲电报方程及其定解条件,得到离散代数方程组.包含狄里克雷和诺依曼边界条件的数值算例表明,本文方法程序实现方便并具有高精度,可应用于求解高维问题. 相似文献
19.
基于欧拉框架下ADER格式,构造一维守恒只有一个时间步的、高精度中心型拉格朗日ADER(LADER)格式.构造r阶LADER格式包括:从欧拉方程出发推导拉格朗日框架下积分形式的方程、采用WENO方法高精度重构节点处守恒量和从1阶到r-1阶的空间导数、求拉氏框架下这些变量的Godunov值,并计算1阶到r-1阶的时间全导数,最后高精度离散积分形式的流通量函数.对光滑流场的模拟表明,LADER格式达到设计的精度;对含强间断的流场模拟表明,数值解在间断附近基本无振荡. 相似文献
20.
采用间断有限元法(discontinuous finite element method,DFEM)求解非规则形状介质内的辐射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热问题的高精度数值结果.和传统连续型有限元方法不同,DFEM将计算区域划分成相互独立的离散单元,形函数的构造、未知量的加权近似以及控制方程的求解均在每一个离散单元上进行.通过在单元之间施加迎风格式的数值通量,DFEM保证了整个计算区域的连续性,因此这种方法兼具良好的几何灵活性和局部守恒性.推导了辐射传输方程和能量扩散方程的射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热的高精度数值结果. 相似文献