辛几何算法在计算非磁化等离子体中波传播轨迹时的应用

为了在数值计算中保持哈密顿系统的辛几何结构不变,利用辛几何算法得到了在线性哈密顿系统中射线追踪方程的一般辛差分格式。通过具体算例,利用辛几何算法计算了波在非磁化等离子体中的传播轨迹,并且与传统Runge-Kutta-Fehlberg算法所得结果进行了比较。利用辛几何算法所得传播轨迹与解析解一致,其色散函数值的误差随时间线性增长,能在长时间内保持色散函数值在一个很小的误差范围内。利用传统的Runge-Kutta-Fehlberg算法所得传播轨迹与解析解不一致,其误差随时间做大幅振荡增加。计算结果表明辛几何算法在保持传播轨迹和色散函数值方面具有独特的优势。