The muphyII Code: Multiphysics Plasma Simulation on Large HPC Systems
arXiv:2305.01487 · doi:10.1016/j.cpc.2023.109064
Abstract
Collsionless astrophysical and space plasmas cover regions that typically display a separation of scales that exceeds any code's capabilities. To help address this problem, the muphyII code utilizes a hierarchy of models with different inherent scales, unified in an adaptive framework that allows stand-alone use of models as well as a model-based dynamic and adaptive domain decomposition. This requires ensuring excellent conservation properties, careful treatment of inner-domain model boundaries for model coupling, and robust time stepping algorithms, especially with the use of electron subcycling. This multi-physics approach is implemented in the muphyII code, tested on different scenarios of space plasma reconnection and evaluated against space probe data and higher-fidelity simulation results from literature. Adaptive model refinement is highlighted in particular, and a hybrid model with kinetic ions, pressure-tensor fluid electrons, and Maxwell fields is appraised.
References in corpus (24)
- Electron-Scale Dynamics of the Diffusion Region during Symmetric Magnetic Reconnection in Space
- Vlasov methods in space physics and astrophysics
- A Review of the 0.1 Reconnection Rate Problem
- Fully Kinetic Simulation of 3D Kinetic Alfven Turbulence
- How accurately can we measure the reconnection rate for the MMS diffusion region event of 2017-07-11?
- Comparison of multi-fluid moment models with Particle-in-Cell simulations of collisionless magnetic reconnection
- A Multiscale Kinetic-Fluid Solver with Dynamic Localization of Kinetic Effects
- Fully kinetic versus reduced-kinetic modelling of collisionless plasma turbulence
- A Parallel Mesh-Adaptive Framework for Hyperbolic Conservation Laws
- Kinetic Vlasov Simulations of collisionless magnetic Reconnection
- Darwin-Vlasov Simulations of magnetized Plasmas
- A Multi Level Multi Domain Method for Particle In Cell Plasma Simulations
- Comparison of time splitting and backsubstitution methods for integrating Vlasov's equation with magnetic fields
- The Fluid-Kinetic Particle-in-Cell Solver for Plasma Simulations
- Two-way coupling of magnetohydrodynamic simulations with embedded particle-in-cell simulations
- Three-dimensional simulations of solar wind turbulence with the hybrid code CAMELIA
- Coupled Vlasov and two-fluid codes on GPUs
- Temperature gradient driven heat flux closure in fluid simulations of collisionless reconnection
- Magnetohydrodynamic with Adaptively Embedded Particle-in-Cell model: MHD-AEPIC
- Bridging hybrid- and full-kinetic models with Landau-fluid electrons: I. 2D magnetic reconnection
- An Energy Conserving Vlasov Solver That Tolerates Coarse Velocity Space Resolutions: Simulation of MMS Reconnection Events
- A cancellation problem in hybrid particle-in-cell schemes due to finite particle size
- Fluid Simulations of Three-Dimensional Reconnection that Capture the Lower-Hybrid Drift Instability
- Kinetic Modeling of Magnetospheres
Cited by in corpus (4)
- A performance portable implementation of the semi-Lagrangian algorithm in six dimensions
- A Parallel-Kinetic-Perpendicular-Moment Model for Magnetized Plasmas
- Electron neural closure for turbulent magnetosheath simulations: energy channels
- Optimal Landau-type closure parameters for two-fluid simulations of plasma turbulence at kinetic scales