JpGU-AGU Joint Meeting 2017

講演情報

[JJ] ポスター発表

セッション記号 P (宇宙惑星科学) » P-EM 太陽地球系科学・宇宙電磁気学・宇宙環境

[P-EM21] [JJ] 宇宙プラズマ理論・シミュレーション

2017年5月25日(木) 13:45 〜 15:15 ポスター会場 (国際展示場 7ホール)

コンビーナ:梅田 隆行(名古屋大学 宇宙地球環境研究所)、成行 泰裕(富山大学人間発達科学部)、三宅 洋平(神戸大学計算科学教育センター)、中村 匡(福井県立大学)

[PEM21-P03] Evaluation of Numerical Properties of Constrained-Transport-Type Schemes for Hybrid Simulations

*天野 孝伸1 (1.東京大学 地球惑星科学専攻)

キーワード:マクスウェル方程式、数値シミュレーション、プラズマ

Numerical simulations of space and astrophysical plasmas need an accurate method for solving Maxwell's equations. The divergence-free property of the magnetic field is the fundamental constraint in the system that must be satisfied in a numerical solution because otherwise, the simulation will become unstable. The Constrained-Transport (CT) scheme, which exactly preserves the discrete divergence-free property, has been quite successful in numerical Magnetohydrodynamics (MHD). In recent years, the CT scheme for Maxwell's equations has been ingeniously combined with an HLL-type (Harten-Lax-Van Leer) Riemann solver for the hydrodynamics part in a consistent fashion. The scheme known as the HLL-UCT shows excellent performance in numerical MHD as well as two-fluid plasma simulations.

It is straightforward to apply the same technique to kinetic Particle-in-Cell (PIC) type simulation method. However, the numerical properties of the scheme as applied to kinetic simulations are not known very well. In fact, artifacts arising from numerical noise inherent in the PIC method (which is absent in a grid-based fluid code) should carefully be analyzed.

In this study, we apply the HLL-UCT scheme to a quasi-neutral plasma hybrid code in which ions are treated as kinetic macroparticles whereas electrons are assumed to be a fluid. We found that naive application of HLL-UCT to a hybrid code may lead to artificial heating and/or cooling of ions, presumably because of excessive dissipation in the HLL-UCT scheme. We thus quantify the numerical artifact by extensive numerical experiments with varying mesh size, the number of particles per cell, plasma beta, etc. We found that the numerical heating/cooling may be explained by absorption (or dissipation) of spontaneous emission of waves arising from a discrete particle effect. Practical workarounds to minimize the numerical artifact for long time simulations will be discussed.