Related papers: The 3D MHD code GOEMHD3 for large-Reynolds-number …
Much progress in realistic modeling of core-collapse supernovae has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements are enabling simulations with…
This paper presents the benchmarking and scaling studies of a GPU accelerated three dimensional compressible magnetohydrodynamic code. The code is developed keeping an eye to explain the large and intermediate scale magnetic field…
The objective of this study is to develop a fully compressible magnetohydrodynamic solver for fast simulations of the global dynamo of the Sun using unstructured grids and GPUs. Accurate modeling of the Sun's convective layers is vital to…
We introduce Gradient Particle Magnetohydrodynamics (GPM), a new Lagrangian method for magnetohydrodynamics based on gradients corrected for the locally disordered particle distribution. The development of a numerical code for MHD…
Electromagnetic turbulence is addressed in tokamak and stellarator plasmas with the global gyrokinetic particle-in-cell codes ORB5 [E. Lanti et al, Comp. Phys. Comm, vol. 251, 107072 (2020)] and EUTERPE [V. Kornilov et al, Phys. Plasmas,…
For more than 40 years the quest to understand how large-scale magnetic fields emerge from turbulent flows in rotating astrophysical systems, such as the Sun, has been a major focus of computational astrophysics research. Using a parameter…
Astrophysical turbulent flows display an intrinsically multi-scale nature, making their numerical simulation and the subsequent analyses of simulated data a complex problem. In particular, two fundamental steps in the study of turbulent…
Numerical simulation of plasma turbulence in the Large Plasma Device (LAPD) [Gekelman et al, Rev. Sci. Inst., 62, 2875, 1991] is presented. The model, implemented in the BOUndary Turbulence (BOUT) code [M. Umansky et al, Contrib. Plasma…
Phenomena taking place in capacitively coupled plasmas with large electrodes and driven at very high frequencies are studied numerically utilizing a novel energy- and charge-conserving implicit fully electromagnetic particle-in-cell / Monte…
This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed…
We compute realistic 3D radiative MHD near-surface models of starspots with substantial penumbrae on cool main-sequence stars using the MURaM simulation code. This work is an improvement on the the previous starspot models in a slab…
This paper presents CUBEP3M, a publicly-available high performance cosmological N-body code and describes many utilities and extensions that have been added to the standard package. These include a memory-light runtime SO halo finder, a…
We present 2D MHD numerical simulations of tearing-unstable current sheets coupled to a population of non-thermal test-particles, in order to address the problem of numerical convergence with respect to grid resolution, numerical method and…
We present a computational method for extreme-scale simulations of incompressible turbulent wall flows at high Reynolds numbers. The numerical algorithm extends a popular method for solving second-order finite differences Poisson/Helmholtz…
Realistic numerical simulations, i.e., those that make minimal use of ad hoc modeling, are essential for understanding the complex turbulent dynamics of the interiors and atmospheres of the Sun and other stars and the basic mechanisms of…
iPIC3D is a widely used massively parallel Particle-in-Cell code for the simulation of space plasmas. However, its current implementation does not support execution on multiple GPUs. In this paper, we describe the porting of iPIC3D particle…
In nuclear astrophysics, quantum simulations of large inhomogeneous dense systems as they appear in the crusts of neutron stars present big challenges. The number of particles in a simulation with periodic boundary conditions is strongly…
We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core…
Vortex flows are structures associated with the rotation of the plasma and/or the magnetic field that are present throughout the solar atmosphere. In recent years, their study has become increasingly important, as they are present on a wide…
We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which…