Premium
Extended magnetohydrodynamics with embedded particle‐in‐cell simulation of Ganymede's magnetosphere
Author(s) -
Tóth Gábor,
Jia Xianzhe,
Markidis Stefano,
Peng Ivy Bo,
Chen Yuxi,
Daldorff Lars K. S.,
Tenishev Valeriy M.,
Borovikov Dmitry,
Haiducek John D.,
Gombosi Tamas I.,
Glocer Alex,
Dorelli John C.
Publication year - 2016
Publication title -
journal of geophysical research: space physics
Language(s) - English
Resource type - Journals
eISSN - 2169-9402
pISSN - 2169-9380
DOI - 10.1002/2015ja021997
Subject(s) - magnetohydrodynamics , magnetosphere , physics , magnetohydrodynamic drive , computational physics , magnetic reconnection , mechanics , astrophysics , magnetic field , quantum mechanics
We have recently developed a new modeling capability to embed the implicit particle‐in‐cell (PIC) model iPIC3D into the Block‐Adaptive‐Tree‐Solarwind‐Roe‐Upwind‐Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHD‐EPIC) algorithm is a two‐way coupled kinetic‐fluid model. As one of the very first applications of the MHD‐EPIC algorithm, we simulate the interaction between Jupiter's magnetospheric plasma and Ganymede's magnetosphere. We compare the MHD‐EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the importance of kinetic effects in controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHD‐EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular, the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD‐EPIC model. The MHD‐EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3‐D structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHD‐EPIC simulation was only about 4 times more than that of the Hall MHD simulation.