等离子体中Fokker-Planck方程的有限元模拟  

Finite Element Scheme for Fokker-Planck Equation in Plasmas

作　　者：翁苏明[1] 盛政明[1] 

机构地区：[1]中国科学院物理研究所光物理实验室,北京100080

出　　处：《计算物理》2007年第2期134-140,共7页Chinese Journal of Computational Physics

基　　金：国家自然科学基金(10335020;10476033);国家杰出青年基金(10425416);国家863高技术惯性约束聚变主题资助项目

摘　　要：利用有限元方法设计了一套相对简单明了的求解Fokker-Planck方程的方案.这个方案不必严格限制计算格点的步长和时间步长,就可以确保分布函数的非负性和粒子数的守恒.通过一维程序模拟,进一步证实了这个方案的可靠性.对于多维问题的分析和一维问题完全一样,所以非常容易将其推广到多维问题.A numerical scheme for the Fokker-Planck equation is proposed with a finite element method, which is relatively simple for implementation and physically clear. It assures a positive distribution function and particle-number eonservation without specific requirement on grid and time step. The scheme is confirmed in one-dimensional simulations. It can be extended easily to multi-dimensional eases.

关 键 词：FOKKER-PLANCK方程 有限元方法 等离子体 Legendre展开式 

分 类 号：O241.8[理学—计算数学] O53[理学—数学]