Related papers: Comparing Numerical Methods for Isothermal Magneti…
We compare the results of numerical simulations of thin and quasi-spherical (thick) accretion flows with existing analytical solutions. We use a Lagrangian code based on the Smooth Particle Hydrodynamics (SPH) scheme and an Eulerian finite…
Supersonic magnetohydrodynamic (MHD) turbulence is a ubiquitous state for many astrophysical plasmas. However, even the basic statistics for this type of turbulence remains uncertain. We present results from supersonic MHD turbulence…
The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a…
Many systems of current interest in relativistic astrophysics require a knowledge of radiative transfer in a magnetized gas flowing in a strongly-curved, dynamical spacetime. Such systems include coalescing compact binaries containing…
We study numerically the dispersion and dissipation properties of the plane wave virtual element method and the nonconforming Trefftz virtual element method for the Helmholtz problem. Whereas the former method is based on a conforming…
Aims. A realistic model of magnetic linkage between a central object and its accretion disk is a prerequisite for understanding the spin history of stars and stellar remnants. To this end, we aim to provide an analytic model in agreement…
Many problems in physics are inherently of multi-scale nature. The issues of MHD turbulence or magnetic reconnection, namely in the hot and sparse, almost collision-less astrophysical plasmas, can stand as clear examples. The Finite Element…
The Lagrangian vortex method offers an alternative numerical approach for direct numerical simulation of turbulence. The fact that it uses the fast multipole method (FMM)--a hierarchical algorithm for N-body problems with highly scalable…
The fluid dynamics community has found success in explaining both the onset and coherent structure formation in wall-bounded turbulence through examining transient growth and pseudoresonance. Whether similar effects are important in plasmas…
In certain astrophysical systems the commonly employed ideal magnetohydrodynamics (MHD) approximation breaks down. Here, we introduce novel explicit and implicit numerical schemes of ohmic resistivity terms in the moving-mesh code AREPO. We…
We construct a relativistic resistive magneto-hydrodynamic (RRMHD) numerical simulation code for high-energy heavy-ion collisions. We split the system of differential equations into two parts, a non-stiff and a stiff part. For the non-stiff…
(abridged) We apply a second-order Godunov code, Athena, to studies of the magnetorotational instability using unstratified shearing box simulations with a uniform net vertical field and a sinusoidally varying zero net vertical field. The…
The numerical simulation of turbulence in stars has led to a rich set of possibilities regarding stellar pulsations, asteroseismology, thermonuclear yields, and formation of neutron stars and black holes. The breaking of symmetry by…
Most numerical investigations on the role of magnetic fields in turbulent molecular clouds (MCs) are based on ideal magneto-hydrodynamics (MHD). However, MCs are weakly ionized, so that the time scale required for the magnetic field to…
The piece-wise parabolic method (PPM) is applied to simulations of forced isotropic turbulence with Mach numbers $\sim 0.1... 1$. The equation of state is dominated by the Fermi pressure of an electron-degenerate fluid. The dissipation in…
Large eddy simulation (LES) has been increasingly used to tackle vortex-dominated turbulent flows. In LES, the quality of the simulation results hinges upon the quality of the numerical discretizations in both space and time. It is in this…
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 full 2 Pi global 3-D stratified MHD simulations of accretion disks. We interpret our results in the context of proto-planetary disks. We investigate the turbulence driven by the magneto-rotational instability (MRI) using the…
We present a sample of 20 massive galaxy clusters with total virial masses in the range of 6 10^14 M_sol<M(vir)< 2 10^15M_sol, re-simulated with a customized version of the 1.5. ENZO code employing Adaptive Mesh Refinement. This technique…
The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere…