Multi-moment advection scheme for Vlasov simulations

被引:11
作者
Minoshima, Takashi [1 ]
Matsumoto, Yosuke [2 ]
Amano, Takanobu [3 ]
机构
[1] Japan Agcy Marine Earth Sci & Technol, Inst Res Earth Evolut, Kanazawa Ku, Yokohama, Kanagawa 2360001, Japan
[2] Nagoya Univ, Solar Terr Environm Lab, Chikusa Ku, Nagoya, Aichi 4648601, Japan
[3] Nagoya Univ, Dept Phys, Chikusa Ku, Nagoya, Aichi 4648602, Japan
关键词
Advection equation; Conservative form; Multi-moment; Information entropy; Vlasov simulations; Magnetized plasmas; CUBIC-POLYNOMIAL INTERPOLATION; SEMI-LAGRANGIAN SCHEME; INCOMPRESSIBLE FLOWS; HYPERBOLIC-EQUATIONS; UNIFIED FORMULATION; CONSERVATIVE SCHEME; UNIVERSAL SOLVER; IMPLEMENTATION; INTEGRATION; SHOCK;
D O I
10.1016/j.jcp.2011.05.010
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We present a new numerical scheme for solving the advection equation and its application to Vlasov simulations. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, for better conservation of the information entropy. We have developed one-and two-dimensional schemes and show that they provide quite accurate solutions within reasonable usage of computational resources compared to other existing schemes. The two-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile for more than hundred rotation periods with little numerical diffusion. This is crucially important for Vlasov simulations of magnetized plasmas. Applications of the one- and two-dimensional schemes to electrostatic and electromagnetic Vlasov simulations are presented with some benchmark tests. (C) 2011 Elsevier Inc. All rights reserved.
引用
收藏
页码:6800 / 6823
页数:24
相关论文
共 40 条
[1]   ELECTRON SHOCK SURFING ACCELERATION IN MULTIDIMENSIONS: TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATION OF COLLISIONLESS PERPENDICULAR SHOCK [J].
Amano, Takanobu ;
Hoshino, Masahiro .
ASTROPHYSICAL JOURNAL, 2009, 690 (01) :244-251
[2]  
[Anonymous], 1991, PLASMA PHYS VIA COMP
[3]   INTEGRATION OF VLASOV EQUATION IN CONFIGURATION SPACE [J].
CHENG, CZ ;
KNORR, G .
JOURNAL OF COMPUTATIONAL PHYSICS, 1976, 22 (03) :330-351
[4]   A forward semi-Lagrangian method for the numerical solution of the Vlasov equation [J].
Crouseilles, Nicolas ;
Respaud, Thomas ;
Sonnendruecker, Eric .
COMPUTER PHYSICS COMMUNICATIONS, 2009, 180 (10) :1730-1745
[5]   Comparison of Eulerian Vlasov solvers [J].
Filbet, F ;
Sonnendrücker, E .
COMPUTER PHYSICS COMMUNICATIONS, 2003, 150 (03) :247-266
[6]   Conservative numerical schemes for the Vlasov equation [J].
Filbet, F ;
Sonnendrücker, E ;
Bertrand, P .
JOURNAL OF COMPUTATIONAL PHYSICS, 2001, 172 (01) :166-187
[7]   A non-periodic 2D semi-Lagrangian Vlasov code for laser-plasma interaction on parallel computer [J].
Ghizzo, A ;
Huot, F ;
Bertrand, P .
JOURNAL OF COMPUTATIONAL PHYSICS, 2003, 186 (01) :47-69
[8]   Total variation diminishing Runge-Kutta schemes [J].
Gottlieb, S ;
Shu, CW .
MATHEMATICS OF COMPUTATION, 1998, 67 (221) :73-85
[9]   THE ELECTROSTATIC EFFECT FOR THE COLLISIONLESS TEARING MODE [J].
HOSHINO, M .
JOURNAL OF GEOPHYSICAL RESEARCH-SPACE PHYSICS, 1987, 92 (A7) :7368-7380
[10]  
HOSHINO M, 1986, THESIS U TOKYO