Two-way coupling of magnetohydrodynamic simulations with embedded particle-in-cell simulations

Abstract

We describe a method for coupling an embedded domain in a magnetohydrodynamic (MHD) simulation with a particle-in-cell (PIC) method. In this two-way coupling we follow the work of Daldorff et al. in which the PIC domain receives its initial and boundary conditions from MHD variables (MHD to PIC coupling) while the MHD simulation is updated based on the PIC variables (PIC to MHD coupling). This method can be useful for simulating large plasma systems, where kinetic effects captured by particle-in-cell simulations are localized but affect global dynamics. We describe the numerical implementation of this coupling, its time-stepping algorithm, and its parallelization strategy, emphasizing the novel aspects of it. We test the stability and energy/momentum conservation of this method by simulating a steady-state plasma. We test the dynamics of this coupling by propagating plasma waves through the embedded PIC domain. Coupling with MHD shows satisfactory results for the fast magnetosonic wave, but significant distortion for the circularly polarized Alfvén wave. Coupling with Hall-MHD shows excellent coupling for the whistler wave. We also apply this methodology to simulate a Geospace Environmental Modeling (GEM) challenge type of reconnection with the diffusion region simulated by PIC coupled to larger scales with MHD and Hall-MHD. In both these cases we see the expected signatures of kinetic reconnection in the PIC domain, implying that this method can be used for reconnection studies.

arXiv commentTo be published in Computer Physics Communication

References in corpus (10)

• Turbulent Heating in Galaxy Clusters Brightest in X-rays
• MPI-AMRVAC for Solar and Astrophysics
• Solar wind turbulent spectrum at plasma kinetic scales
• Adaptive mesh refinement for conservative systems: multi-dimensional efficiency evaluation
• Energy dynamics and current sheet structure in fluid and kinetic simulations of decaying magnetohydrodynamic turbulence
• Reconnection and particle acceleration in interacting flux ropes I. Magnetohydrodynamics and test particles in 2.5D
• Evolution of a Current Sheet in a Solar Flare
• Coupled Vlasov and two-fluid codes on GPUs
• On the Electron Agyrotropy during Rapid Asymmetric Magnetic Island Coalescence in Presence of a Guide Field
• Missing bits of the solar jigsaw puzzle: small-scale, kinetic effects in coronal studies

Cited by in corpus (13)

• MPI-AMRVAC 3.0: updates to an open-source simulation framework
• PIC methods in astrophysics: Simulations of relativistic jets and kinetic physics in astrophysical systems
• A fully self-consistent model for solar flares
• Combining MHD and kinetic modelling of solar flares
• Magnetohydrodynamic with Adaptively Embedded Particle-in-Cell model: MHD-AEPIC
• Ideal MHD instabilities for coronal mass ejections
• The muphyII Code: Multiphysics Plasma Simulation on Large HPC Systems
• Toward Realistic Solar Flare Models: An explicit Particle-In-Cell solver in the DISPATCH framework
• Black hole outflows initiated by a large-scale magnetic field
• The PLUTO Code on GPUs: Offloading Lagrangian Particle Methods
• Multi-scale simulations of particle acceleration in astrophysical systems
• Multiphysics simulations of collisionless plasmas
• Optimal Landau-type closure parameters for two-fluid simulations of plasma turbulence at kinetic scales