Related papers: Exploring the Formation of Resistive Pseudodisks w…
The collapse of singular magnetized toroids (Li & Shu 1996) is a natural representation of an early phase in star formation, bridging the prestellar and protostellar phases of the collapse of molecular cloud cores. We revisit the collapse…
We introduce a new GPU-accelerated general-relativistic magneto-hydrodynamic (GR-MHD) code based on HARM which we call cuHARM. The code is written in CUDA-C and uses OpenMP to parallelize multi-GPU setups. Our code allows us to run high…
We present the development and benchmarking of Poisson solvers for graphics processing units (GPUs). Implemented in the Astaroth platform, the solvers feature high computational efficiency. We present novel combinations of discretizations…
We present a detailed study of the collapse of molecular cloud cores using high resolution 3D adaptive mesh refinement (AMR) numerical simulations. In this first in a series of investigations our initial conditions consists of spherical…
We report results from twelve simulations of the collapse of a molecular cloud core to form one or more protostars, comprising three field strengths (mass-to-flux ratios, {\mu}, of 5, 10, and 20) and four field geometries (with values of…
In the context of star and planet formation, understanding the formation of disks is of fundamental importance. Previous studies found that the magnetic field has a very strong impact on the collapse of a prestellar cloud, particularly in…
We model gas inflow patterns onto circumstellar disks and the evolution of the pseudodisk using three-dimensional resistive MHD simulations. Starting from a prestellar core without turbulence and with a misalignment between the initial…
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…
Numerical methods to improve the treatment of magnetic fields in smoothed field magnetohydrodynamics (SPMHD) are developed and tested. Chapter 2 is a review of SPMHD. In Chapter 3, a mixed hyperbolic/parabolic scheme is developed which…
General-relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a graphics processing unit (GPU) accelerated GRMHD code \hammer{} with multi-faceted optimizations…
Efficient magnetic braking is a formidable obstacle to the formation of rotationally supported disks (RSDs) around protostars in magnetized dense cores. We have previously shown, through 2D (axisymmetric) non-ideal MHD simulations, that…
The numerical study of relativistic magnetohydrodynamics (MHD) plays a crucial role in high-energy astrophysics, but unfortunately is computationally demanding, given the complex physics involved (high Lorentz factor flows, extreme…
We present the first open-source, GPU-based code for complex plasmas. The code, OpenDust, aims to provide researchers both experimenters and theorists user-friendly and high-performance tool for self-consistent calculation forces, acting on…
We study numerically the collapse of rotating, magnetized molecular cloud cores, focusing on rotation and magnetic braking during the main accretion phase of isolated star formation. Motivated by previous numerical work and analytic…
Dynamic Parallelism (DP) is a runtime feature of the GPU programming model that allows GPU threads to execute additional GPU kernels, recursively. Apart from making the programming of parallel hierarchical patterns easier, DP can also…
The pseudospectral method is a highly accurate numerical scheme suitable for turbulence simulations. We have developed an open-source pseudospectral code, \textsc{\textsf{Calliope}}, which adopts the P3DFFT library \citep{Pekurovsky2012} to…
Magnetic fields are critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Dynamo action determines their dynamical structure and strength. Simulations of…
This paper introduces the Sheffield Magnetohydrodynamics Algorithm Using GPUs (SMAUG+), an advanced numerical code for solving magnetohydrodynamic (MHD) problems, using multi-GPU systems. Multi-GPU systems facilitate the development of…
We present a new 2D axisymmetric code, cuDisc, for studying protoplanetary discs, focusing on the self-consistent calculation of dust dynamics, grain size distribution and disc temperature. Self-consistently studying these physical…
Magnetic fields are important contributers to the dynamics of collapsing molecular cloud cores, and can have a major effect on whether collapse results in a single protostar or fragmentation into a binary or multiple protostar system. New…