渐近保持PIC算法在等离子体粒子模拟中的应用

核聚变与等离子体物理 ›› 2010, Vol. 30 ›› Issue (2): 138-142.

渐近保持PIC算法在等离子体粒子模拟中的应用

孙安邦，毛根旺，Degond P，Deluzet F，Navoret L，Vignal M-H

(1. 西北工业大学航天学院，西安 710072；2. 保尔萨巴蒂埃大学，图卢兹 31062)

收稿日期:2009-05-04 修回日期:2010-03-04 出版日期:2010-06-15 发布日期:2010-06-03
作者简介:孙安邦(1984-)，男，安徽萧县人，中法联合培养博士生，研究领域为电推进技术，航空宇航推进理论与工程。

Application of asymptotic preserving PIC method in the plasma particle simulation

SUN An-bang, MAO Gen-wang, Degond P, Deluzet F, Navoret L, Vignal M-H

(1. College of Astronautics, Northwestern Polytechnical University, Xi’an 710072;2. University Paul Sabatier，Toulouse 31062 )

Received:2009-05-04 Revised:2010-03-04 Online:2010-06-15 Published:2010-06-03

摘要/Abstract

摘要： 针对准中性Vlasov-Poisson系统中多尺度物理参数共存的问题开展研究，采用渐近保持格式推导出与Vlasov方程相耦合的Poisson方程。数值模拟了一维等离子体在真空中的膨胀，比较了渐近保持PIC算法与传统PIC算法模拟结果的异同。结果表明：两种算法在满足多尺度中小尺度量限制的前提下，其结果保持高度一致；且在克服小尺度量束缚时，渐近保持PIC算法依然可以正确地模拟等离子体的物理行为，大大节省了计算时间。

关键词: Vlasov-Poisson, 渐近保持, PIC, 多尺度

Abstract: The existing problem of multi-scale parameters in the quasi-neutral Vlasov-Poisson system was studied, using asymptotic preserving scheme to induce a Poisson equation coupled with Vlasov equation. One dimensional plasma expansion in vacuum was simulated, the asymptotic preserving PIC and classical PIC results were compared. The simulated results show that two methods have a good agreement while satisfying the smaller scale constraint, and the asymptotic preserving scheme still simulates correctly the physics behavior of plasma while free of the smaller scale constraint, the calculated time is saved intensely.

Key words: Vlasov-Poisson, Asymptotic preserving, PIC, Multi-scale

中图分类号:

V439.4

孙安邦，毛根旺，Degond P，Deluzet F，Navoret L，Vignal M-H. 渐近保持PIC算法在等离子体粒子模拟中的应用[J]. 核聚变与等离子体物理, 2010, 30(2): 138-142.

SUN An-bang, MAO Gen-wang, Degond P, Deluzet F, Navoret L, Vignal M-H. Application of asymptotic preserving PIC method in the plasma particle simulation[J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2010, 30(2): 138-142.

参考文献

[1] Birdsall C K, Langdon A B. Plasma physics via computer simulation [M]. Bristol: Adam Hilger, 1991.

[2] Hockney R W, Eastwood J W. Computer simulation using particles [M]. Institute of Physics Publishing, 1988.

[3] Cohen B I, Langdon A B, A Friedman. Implicit time integration for plasma simulation [J]. J. Comput. Phys., 1982, 46: 15.

[4] Mason R J. Implicit moment PIC-hybrid simulation of collisional plasma [J]. J. Comput. Phys., 1983, 51: 484.

[5] Grismayer T, Mora P, Adam J C, et al. Electron kinetic effects in plasma expansion and ion acceleration [J]. Phys. Review E, 2008, 77: 066407.

[6] Crispel P, Degond P, Vignal M-H. An asymptotic preserving scheme for the two-fluid Euler-Poisson model in the quasi-neutral limit [J]. J. Comput. Phys., 2007, 223: 204.

[7] Crispel P, Degond P, Vignal M-H. Quasi-neutral fluid models for current-carrying plasmas [J]. J. Comput. Phys., 2005, 205: 408.

[8] Degond P, Deluzet F, Navoret L. An asymptotic stable particle-in-cell (PIC) scheme for collisionless plasma simulations near quasi-neutrality [J]. C. R. Acad. Sci. Paris., Ser. I, 2006, 343: 613.

[9] Friedman A, Parker S E, Ray S L et al. Multi-scale particle-in-cell plasma simulation [J]. J. Comput. Phys., 1991, 96: 54.

[10] Larsen E W, Morel J E, Jr W F M. Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes [J]. J. Comput. Phys., 1987, 69: 283-324.

[11] Larsen E W, Morel J E. Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes Ⅱ [J]. J. Comput. Phys., 1989, 83: 212-236.

[12] Golse F, Jin S, Levermore C D. The convergence of numerical transfer schemes in diffusive regimes I: the discrete-ordinate method [J]. SIAM J. Num. Anal., 1999, 36: 1333–1369.

[13] Jin S. Efficient asymptotic-preserving (AP) schemes for some multi-scale kinetic equations [J]. SIAM J. Sci. Comp., 1999, 21: 441–454.

[1]	王泽, 胡纯栋, 谢亚红, 刘伟, 崔庆龙, 赵远哲, 许永建, 班挺, 段练. 射频负离子源中性束关键参数监控报警系统研制[J]. 核工业西南物理研究院, 2024, 44(2): 183-190.
[2]	张正浩, 赵杰, 李平川, 唐德礼. 双阳极霍尔推力器效率及推力的模拟研究[J]. 核工业西南物理研究院, 2023, 43(3): 347-351.
[3]	袁啸林1, 2，陈跃1, 李建刚1，胡建生1，周跃1，侯吉磊1, 2，黄明1. 基于实验物理与工业控制系统的EAST 真空抽气及加料控制系统设计[J]. 核聚变与等离子体物理, 2019, 39(1): 34-40.
[4]	许婕，董方正，田培红，颉延风，唐芳群，蔡萧，赵丽，魏凌峰，杨青巍. 实验物理与工业控制系统平台在HL-2A主机测控系统中的应用[J]. 核聚变与等离子体物理, 2018, 2(1): 87-92.
[5]	李欣怡1，夏凡1，杨超2，潘莉1，张刚1，陈燎原1. HL-2A 装置EPICS 过程变量管理自动化系统[J]. 核聚变与等离子体物理, 2016, 36(2): 167-173.
[6]	张刚，夏凡，李欣怡，陈燎原. 基于EPICS 的HL-2A 装置实时监测和报警系统[J]. 核聚变与等离子体物理, 2016, 36(1): 83-87.
[7]	李宏斌，唐德礼，聂军伟，耿少飞，张帆. 圆柱形霍尔推进器的三维刻蚀模拟与实验研究[J]. 核聚变与等离子体物理, 2015, 35(2): 181-186.
[8]	王川，夏凡，潘莉，潘卫，赵丽，宋啸，李欣怡，陈燎原. EPICS 在HL-2A 控制系统PLC 上的部署[J]. 核聚变与等离子体物理, 2013, 33(2): 142-146.
[9]	杨莹1，付东儒1，张连珠2. N₂射频放电带电粒子(e, N₂⁺, N⁺)的PIC-MC模拟[J]. 核聚变与等离子体物理, 2012, 32(3): 206-212.
[10]	孙倩，杨莹，赵海涛，郝莹莹，张连珠，赵国明. N₂-Ar射频放电表面氮化的等离子体特性模拟[J]. 核聚变与等离子体物理, 2012, 32(2): 128-132.
[11]	姚福宝, 郝莹莹, 张连珠, 赵海涛. 氢气甚高频与高频容性耦合放电等离子体特征的比较[J]. 核聚变与等离子体物理, 2011, 31(4): 320-326.
[12]	邱清泉, 励庆孚, 苏静静, JIAO Yu , Finley Jim. 工作参数对平面直流磁控溅射放电特性的影响[J]. 核聚变与等离子体物理, 2009, 29(2): 182-188.
[13]	袁忠才，时家明，黄勇，马柳. 低温等离子体数值模拟方法的分析比较[J]. 核聚变与等离子体物理, 2008, 28(3): 278-284.
[14]	袁忠才, 时家明. PICPMCC 数值模拟四阳极直流辉光放电[J]. 核聚变与等离子体物理, 2005, 25(4): 305-310.