Application of symplectic geometric algorithm in calculating propagation trajectory of wave in unmagnetized plasma

doi:10.16568/j.0254-6086.201801005

Abstract

Abstract:

In order to preserve the symplectic geometric structure of Hamiltonian system, symplectic geometric algorithm is used to solve ray tracing equations. Symplectic geometric algorithm assures that conversion between two steps is symplectic. Symplectic geometric algorithm and Runge-Kutta-Fehlberg method are used to solve ray tracing equations respectively. The accuracy of both methods is fourth order. The results calculated by a simplistic geometric algorithm are in accordance with the analytic solution. In symplectic geometric algorithm, the error of dispersion function value increases with time linearly. The dispersion function value is within high accuracy. The results calculated by Runge-Kutta-Fehlberg method are different with the analytical solution. The calculation results show that the symplectic geometric algorithm has irreplaceable advantages in maintaining propagation trajectory and dispersion function value. Next symplectic geometric algorithm will be applied in magnetized plasma and nonlinear Hamiltonian systems.

Key words: Symplectic geometric algorithm, Ray tracing, Plasma

摘要：

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

关键词: 辛几何算法, 射线追踪, 等离子体

CLC Number:

O534

YAO Kun. Application of symplectic geometric algorithm in calculating propagation trajectory of wave in unmagnetized plasma[J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2018, 2(1): 29-33.

姚琨. 辛几何算法在计算非磁化等离子体中波传播轨迹时的应用[J]. 核聚变与等离子体物理, 2018, 2(1): 29-33.

References

[1] Feng K. Difference schemes for Hamiltonian formalism and symplectic geometry[J]. J. Comput. Math., 1986, 4: 279.

[2] Feng K, Qing M Z. Hamiltonian algorithms for Hamil- tonian systems and a comparative numerical study[J]. Comput. Phys. Comm, 1991, 65: 173.

[3] Feng K. Formal power series and numerical algorithms for dynamical systems[J]. Math. Appl., 1991, 6: 17.

[4] Feng K, Wang D L. Symplectic difference schemes for Hamiltonian systems in general symplectic structure[J]. J. Comput. Math., 1991, 9: 86.

[5] Feng K, Shang Z J. Volume-preserving algorithms for source-free dynamical systems[J]. Numer. Math., 1995, 71: 451.

[6] Feng K. The calculus of generating functions and the formal energy for Hamiltonian algorithms[J]. J. Comput. Math., 1998, 16: 481.

[7] Feng K. The step-transition operators for multi-step methods of ODE's[J]. J. Comput. Math., 1998, 16: 193.

[8] 刘世兴, 刘畅, 常鹏, 等. 辛几何算法在特殊Chaplygin系统中的应用研究[J]. 物理学报, 2011, 60(3): 034501.

[9] 刘学深, 丁培柱. 量子系统保结构计算新进展[J]. 物理学进展, 2004, 24(1): 47-89.

[10] Bernstein I B. Geometric optics in space- and time- varying plasmas[J]. Phys. Fluids, 1975, 18: 320.

[11] Smirnov A P, Harvey R W. Calculations of the current drive in DⅢ-D with the GENRAY ray tracing code[J]. Bull. Am. Phys. Soc, 1995, 40: 35.

[12] Yang C, Bonoli P T, Wright J C, et al. Modelling of the EAST lower-hybrid current drive experiment using GENRAY/CQL3D and TORLH/CQL3D[J]. Plasma Phys. Contr. Fusion, 2014, 56: 125003.

[13] Esterkin A R, Piliya A D. Fast ray tracing code for LHCD simulations[J]. Nucl. Fusion, 1996, 36: 1501.

[14] McVey B D. A ray-tracing analysis of fast-wave heating of tokamaks[J]. Nucl. Fusion, 1979, 19: 461.

[15] 冯康, 秦孟兆. 哈密顿系统的辛几何算法[M]. 杭州: 浙江科学技术出版社, 2003. 189-194.