2024/03/19 by Yuxi Chen, Chen, Yuxi, G. Tóth +11
Physics and Astronomy · #FOS: Physical sciences #Ionosphere and magnetosphere dynamics #Magnetic confinement fusion research #Solar and Space Plasma Dynamics #Space Physics (physics.space-ph)
paper · pdf · doi:10.48550/arxiv.2403.12395
openalex publication_date 2024/03/19 · openalex created_date 2024/03/21 · openalex updated_date 2026/08/01
The charge exchange between the interstellar medium (ISM) and the solar wind plasma is crucial for determining the structures of the heliosphere. Since both the neutral-ion and neutral-neutral collision mean free paths are either comparable to or larger than the size of the heliosphere, the neutral phase space distribution can deviate far away from the Maxwellian distribution. A kinetic description for the neutrals is crucial for accurately modeling the heliosphere. It is computationally challenging to run three-dimensional (3D) time-dependent kinetic simulations due to the large number of macro-particles. In this paper, we present the new highly efficient SHIELD-2 model with a kinetic model of neutrals and a magnetohydrodynamic (MHD) model for the ions and electrons. To improve the simulation efficiency, we implement adaptive mesh refinement (AMR) and particle splitting and merging algorithms for the neutral particles to reduce the particle number that is required for an accurate simulation. We present several tests to verify and demonstrate the capabilities of the model.