Related papers: The Magnetohydrodynamic-Particle-In-Cell Module in…
Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We…
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…
Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview…
We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of…
We describe a new electrostatic Particle-In-Cell (PIC) code in curvilinear geometry called Curvilinear PIC (CPIC). The code models the microscopic (kinetic) evolution of a plasma with the PIC method, coupled with an adaptive computational…
(Abridged) We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics when the light-crossing time is only marginally shorter than dynamical timescales. The transfer…
We present a method for simulating the evolution of HII regions driven by point sources of ionizing radiation in magnetohydrodynamic media, implemented in the three-dimensional Athena MHD code. We compare simulations using our algorithm to…
We have developed a hybrid code GMEC: Gyro-kinetic Magnetohydrodynamics (MHD) Energetic-particle Code that can numerically simulate energetic particle-driven Alfv\'en eigenmodes and energetic particle transport in tokamak plasmas. In order…
When supersonic plasma flows collide, many physical processes contribute to the morphology of the resulting shock. One of these processes is the acceleration of non-thermal ions, which will, eventually, reach relativistic speeds and become…
Collisionless magnetic reconnection typically requires kinetic treatments that are, in general, computationally expensive compared to fluid-based models. In this study, we use the magnetohydrodynamics with adaptively embedded…
We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements…
We implement a local model for a spherical collapsing/expanding gas cloud into the Athena++ magnetohydrodynamic code. This local model consists of a Cartesian periodic box with time-dependent geometry. We present a series of benchmark test…
We describe the implementation of multigrid solvers in the Athena++ adaptive mesh refinement (AMR) framework and their application to the solution of the Poisson equation for self-gravity. The new solvers are built on top of the AMR…
We adapt a modern scheme of smoothed particle hydrodynamics (SPH) to our tree N-body/SPH galactic chemodynamics code GCD+. The applied scheme includes imple- mentations of the artificial viscosity switch and artificial thermal conductivity…
We introduce $\texttt{A}$strophysical $\texttt{H}$ybrid-$\texttt{K}$inetic simulations with the $\texttt{flASH}$ code ($\texttt{AHKASH}$) -- a new Hybrid particle-in-cell (PIC) code developed within the framework of the multi-physics code…
We present a numerical solver for plasma dynamics simulations in Hall magnetohydrodynamic (HMHD) approximation in one, two and three dimensions. We consider both isotropic and anisotropic thermal pressure cases, where a general gyrotropic…
Large scale simulations are a key pillar of modern research and require ever-increasing computational resources. Different novel manycore architectures have emerged in recent years on the way towards the exascale era. Performance…
We describe an initial implementation of an electrostatic Particle-in-Cell (ES-PIC) module with adaptive Cartesian mesh in our Unified Flow Solver framework. Challenges of PIC method with cell-based adaptive mesh refinement (AMR) are…
We perform a suite of two- and three-dimensional magnetohydrodynamic (MHD) simulations with the Athena code of the non-driven Kelvin-Helmholtz instability in the subsonic, weak magnetic field limit. Focusing the analysis on the non-linear…
A brief review of the existing kinetic-magnetohydrodynamic(MHD) hybrid models for the alpha particle physics in burning plasma demonstrates that the pressure-coupling scheme is equivalent to the current-coupling scheme only in a specific…