分数阶微积分的预估-校正算法及其应用

选用分数阶微分方程的预估-校正数值算法,对Chen混沌系统进行仿真研究．首先,讨论分数阶Chen混沌系统在一定的初始条件下,系统为混沌的并且仍然呈现出丰富和复杂的分数阶混沌动力学行为;然后,利用预估-校正数值计算方法,对分数阶Chen混沌系统方程进行离散化处理,得到系统方程组的离散化式;最后通过MATLAB软件进行计算,得到分数阶Chen混沌系统的仿真相图．根据初始状态变量的不同,得到相应混沌系统的仿真图．证明了分数阶预估-校正法可以很好地对分数阶系统方程进行数值稳定分析．     

非结构网格下两步压力校正算法的潮流模拟

针对非结构网格下的潮流数值计算,采用有限体积法离散基本方程,对流-扩散项离散采用幂率格式,交叉扩散项引入超松弛校正方法,水位校正方程应用Rhie-Chow动量插值思想和SIMPLEC类算法导出.采用两步压力校正法将交叉扩散项单独考虑,以提高模型对高度畸变网格的适应能力.计算结果表明,该处理方法可提高压力欠松弛系数,并能有效改善复杂区域潮流模拟的健壮性.通过对黄浦江感潮河段潮流的模拟,表明该模型的预测结果与实测资料吻合良好,模型对潮流的模拟结果是可靠的.     

一种改进型刚性常微分方程数值解法——T—R预估—校正法

Treanor显式方法是一种求解刚性常微分方程的数值积分法,用它来解刚性方程时仍然存在一些问题。本文对Treanor方法的稳定性进行了分析和讨论,以Treanor方法为基础,采用Richardson外推技巧,构造出一种改进型方法—T-R预估-校正法。分析和计算都表明,与Treanor方法相比,改进的方法具有更大的稳定域,且提高了数值解的精度。它更适合于刚性方程的求解,是数字仿真中一种较好的数值解法。     

求非线性方程分歧点的一种新的预估-校正法

提出了求非线性方程简单分歧点的一种新的预估-校正法.该法采用Seydel方法的第一步做预估,其次在与其正交方向上用Newton迭代法进行校正,数值算例表明了该方法的有效性.     

飞行器操纵面嗡鸣的非结构网格并行计算方法

为了数值模拟飞行器操纵面的嗡鸣现象,在集群计算机MPI并行计算环境下建立基于三维非定常欧拉方程耦合结构运动方程的嗡鸣计算方法.气动流场求解采用基于非结构网格的中心有限体积法进行空间离散,时间推进采用双时间方法,结构运动方程采用Adams预估校正方法求解.针对翼面与操纵面缝隙间存在的网格运动问题,在非结构网格系统上采用Delaunay图映射方法实现网格的运动变形.最后,使用飞行器操纵面标准嗡鸣计算模型对计算方法进行验证,结果表明:所建立的并行计算方法正确,程序具有很好的计算效率,能够对飞行器操纵面嗡鸣进行高效的数值分析.     

求解紊流流动方程组的预估-校正法

SIMPLE类方法是一种有效的求解Navier-Stokes方程组的原始变量法．然而，对非恒定流，尽管SIMPLE类方法采用隐式差分离散，允许较大时间步长，但对每一时间步都需要反复的迭代计算，求解相当费时，从分步法的概念出发，本文讨论了用预估-校正法求解隐式离散的Navier-Stokes方程组和k-ε紊流方程的过程，并与SIMPLE方法进行了比较．     

K-∈湍流模型下门桥水阻力的数值计算

采用RANS方程和RNGκ-ε二方程湍流模型及VOF算法对绕门桥的自由表面流动进行了数值计算。数值计算时，使用有限体积法离散控制方程，对流项采用二阶迎风差分格式，扩散项采用中心差分格式。对门桥模型进行了一系列速度下的数值模拟；计算结果与试验测量数据相比较，该计算方法具有较好的计算精度，可用于门桥性能分析和优化。     

基于有限差分方法的水中爆炸三维数值模拟

针对水中爆炸的三维数值模拟问题,采用基于有限差分的拉格朗日与欧拉相结合的方法,在矩形网格上离散差分基本方程组;在欧拉步中引入模糊方法处理界面,计算各输运量;编写了数值模拟程序,进行了炸药在水中爆炸的三维数值模拟.爆炸场中考虑了有刚性障碍物及有反射有限边界计算域.结果表明,模拟结果与实际情况基本吻合,特别是冲击波在水中传播的模拟效果较好,证明所采用的三维数值模拟方法是有效和可行的.     

基于重心插值的多体系统动力学离散变分方法

针对多体系统动力学仿真欧拉—拉格朗日方程的数值求解,基于哈密顿原理和离散变分原理的方法,以平面双连杆为数值算例,使用插值方法和数值积分得到其离散的欧拉—拉格朗日方程,对该离散方程进行求解。离散过程采用重心拉格朗日插值提高算法稳定性,插值节点分别选取等距节点和非均匀节点,其中非均匀节点包括第一类、第二类Chebyshev节点,数值积分采用精度较高的高斯勒让德积分。数值结果表明,该方法在步长较大时相比传统采用的龙格库塔法得到较好的结果,并且具有更高的效率,适用于长时间仿真。     

气垫船兴波破冰问题的数值计算

将求解船舶兴波问题的Rankine源数值计算方法推广应用到水面边界条件为碎冰面、纯冰面和半冰面-半水面等气垫船破冰中不同类型的兴波问题.对不同类型兴波问题采用统一形式的基本方程,通过面元法与有限差分法相结合的方法对方程进行数值离散.重点对气垫船在全冰面、半冰面-半水面上运动时的兴波进行了数值计算,兴波阻力和波形等计算结果与相关文献符合较好.