求解弹性波方程的辛RKN格式   总被引：2，自引：2，他引：0       下载免费PDF全文

刘少林  李小凡  汪文帅  朱童 《地球物理学报》2015,58(4):1355-1366

将弹性波方程变换至Hamilton体系,构造适用于弹性波模拟的高效显式二阶辛Runge-Kutta-Nystrm(RKN)格式,运用根数理论得到此格式的阶条件方程组.通过给定系数的限定条件,得到方程的对称解.为了使时间离散误差达到极小,提出数值频率与真实频率比较,通过Taylor展开,得到关于辛系数的限定方程,求解方程组得到最小频散辛RKN格式.对比分析时间演进方程的稳定性,得到使库朗数达到极大值的限定方程,求解方程组得到最稳定辛RKN格式.发现此两种格式为同一格式.新得到的辛RKN格式不依赖于空间离散方法,为了对比的需要,选取有限差分法进行空间离散.在频散、稳定性分析中,与常见辛格式对比,从理论上分析了本文提出的格式在数值频散压制、稳定性提升等方面的优势,数值实验进一步证实了理论分析的正确性.  相似文献   

Behaviour of a new type of Runge–Kutta methods when integrating satellite orbits     

Ana B. González  Pablo Martín  David J. López 《Celestial Mechanics and Dynamical Astronomy》1999,75(1):29-38

Recently, González, Martín and Farto have developed new numerical methods (RKGM methods) of Runge–Kutta type and fixed step size for the numerical integration of perturbed oscillators. Moreover, it seems natural to study the behaviour of these new methods for the accurate integration of orbital problems after the application of linearizing transformation, such us KS or BF due to the fact that in these variables, the structure of the problem is of the form of perturbed oscillators, for which the methods constructed are indicated. In this paper, we check the efficiency of these new methods when integrating the satellite problem. The RKGM methods show a very good behaviour when they compete with other, classical and special, methods. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献