计算物理 ›› 2015, Vol. 32 ›› Issue (3): 253-263.

• 论文 •    下一篇

多尺度等离子体动力学数值模拟中最初的无方程投影积分技术

李百文1, 石黑静儿2, M. M. Skoric2   

  1. 1. 北京应用物理与计算数学研究所, 北京 100094;
    2. 日本核融合科学研究所, 土歧 日本
  • 收稿日期:2014-04-11 修回日期:2014-07-22 出版日期:2015-05-25 发布日期:2015-05-25
  • 作者简介:李百文(1968-),男,辽宁昌图人,博士,研究员,从事等离子体物理、带电粒子输运等研究,E-mail:li_baiwen@iapcm.ac.cn
  • 基金资助:
    国家自然科学基金(11375032,11175030);中物院重点基金(2011A0102008);973课题(6131940301)资助项目

A Primal Equation Free Projective Integration Scheme in Multi-Scale Plasma Dynamics

LI Baiwen1, SEIJI Ishiguro2, M. M. Skoric2   

  1. 1. Institute of Applied Physics and Computational Mathematics;
    2. National Institute for Fusion Science, Toki, Japan
  • Received:2014-04-11 Revised:2014-07-22 Online:2015-05-25 Published:2015-05-25

摘要: 发展一种针对等离子体多尺度问题的"最初的无方程投影积分(p-EFPI)"方法并成功应用于离子声波的模拟.该方法包括时间推进和时间投影两部分.电子速度通过概率分布函数或累积分布函数投影外推,电子位置坐标则通过投影外推的电子电荷密度重新调出;与电子相比,假设粗糙的离子运动可以很好地描述宏观尺度等离子体的演化行为,通过追踪和投影外推每个离子轨道保持全部离子运动学效应.尝试和检验p-EFPI方法在模拟等离子体物理现象中的应用.以离子声波为范例,数值计算表明:在离子相空间出现离子捕获效应前,p-EFPI模拟能够很好地再现PIC的结果,在离子捕获发生后,p-EFPI模拟仍然可以抓住除了离子捕获特征之外的其它主要物理特征.

关键词: 多尺度模拟, 最初的无方程投影积分, 非线性等离子体, 粒子网格法

Abstract: A coarse projective method for simulation of multi-dimensional macroscopic systems,Primal Equation Free Projective Integration(p-EFPI)method,was developed and applied to simulation of plasma ion acoustic wave.The method consists of two parts: Time stepper and time projection.For election velocities, probability distribution function or cumulative distrihution function are recorded in time to project,while,corresponding electron positions are reloaded(restricted)via electron charge density.Ion motion is assumed inherently coarse-grained as compared to electron dynamics,to describe possibly macro scale plasma evolution sufficiehtly well. To preserve full ion kinetic effects,individual ion orbits are tracked in time and simply extrapolate to pmject.This is one attempt to test method in p-EFPI framework on plasma physics phenomenon.Our calculations for ion acoustic wave paradigm indicate that P-EFPI simulation follows PIC results nice before ion trapping in phase space appears.However,as ion trapping takes place,main characteristic still keeps,except for ion trapping characteristic.

Key words: multi-scale simulations, p-FFPI, nonlinear plasmas, particle-in-cell(PIC)

中图分类号: