高阶CIP数值方法及其在相关物理问题中的应用

计算物理 ›› 2011, Vol. 28 ›› Issue (2): 259-267.

高阶CIP数值方法及其在相关物理问题中的应用

傅德月¹, 彭晓东²

1. 成都邮电职业技术学院, 四川成都 610061;
2. 核工业西南物理研究院, 四川成都 610041

收稿日期:2010-01-28 修回日期:2010-07-07 出版日期:2011-03-25 发布日期:2011-03-25
通讯作者: 彭晓东,E-mail:pengxd@swip.ac.cn
作者简介:傅德月(1963-),女,重庆市,副教授,硕士,从事微波通信与数值模拟,成都邮电职业技术学院.
基金资助:
国家自然科学基金(10770543);国家磁约束核聚变研究专项(2009GB105002,2010GB106006)资助项目

High-order CIP Numerical Method and Applications in Vlasov-Poisson Equation

FU Deyue¹, PENG Xiaodong²

1. Chengdu Post and Telecommunication College, Chengdu 610061, China;
2. Southwestern Institue of Physics, Chengdu 610041, China

Received:2010-01-28 Revised:2010-07-07 Online:2011-03-25 Published:2011-03-25

摘要/Abstract

摘要： 利用函数的高阶空间导数值构建其高次插值,得到高阶CIP(Constrained Interpolation Profile)数值算法,并在此基础上模拟研究等离子体物理中著名的伏拉索夫-泊松(Vlasov-Poisson)方程相关物理问题.高阶CIP数值方法具有更高数值精度,从而可以在同等精度的情况下减少计算格点数,加速数值计算速度.

关键词: CIP数值方法, 高阶插值, 伏拉索夫-泊松方程

Abstract: A high-order CIP(Constrained Interpolation Profile,HCIP) method is developed by constructing a high-order interpolation function based on high-order spatial derivatives.Numerical errors of both HCIP scheme and a standard CIP(SCIP) scheme proposed by Yabe and et al are studied.With HCIP method,we investigate numerically physical problems in the famous Vlasov-Poisson equation,such as Landau damping and two-stream instability.It shows that HCIP method is a fifth-order accuracy scheme.Its accuracy is higher than that of SCIP method.Dynamical evolution due to Landau damping in an electrostatic plasma simulated by HCIP method agrees well with previous results.By reducing the number of grids calculation speed is increased with same accuracy.

Key words: CIP method, high-order interpolation, Vlasov-Poisson equation

中图分类号:

O241.82

傅德月, 彭晓东. 高阶CIP数值方法及其在相关物理问题中的应用[J]. 计算物理, 2011, 28(2): 259-267.

FU Deyue, PENG Xiaodong. High-order CIP Numerical Method and Applications in Vlasov-Poisson Equation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 28(2): 259-267.