Multi-moment advection scheme for Vlasov simulations

JOURNAL OF COMPUTATIONAL PHYSICS

卷:230

Article

Minoshima, Takashi¹  Matsumoto, Yosuke²  Amano, Takanobu³ 

[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;   

DOI  :  10.1016/j.jcp.2011.05.010

来源: Elsevier

【 摘 要 】

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.

【 授权许可】

Free