Related papers: A robust numerical scheme for highly compressible …
Numerical solutions of the incompressible magnetohydrodynamic (MHD) equations are reported for the interior of a rotating, perfectly-conducting, rigid spherical shell that is insulator-coated on the inside. A previously-reported spectral…
Context: Calculating stellar pulsations requires a sufficient accuracy to match the quality of the observations. Many current pulsation codes apply a second order finite-difference scheme, combined with Richardson extrapolation to reach…
We present a detailed study of localised magnetohydrodynamical (MHD) instabilities occuring in two--dimensional magnetized accretion disks. We model axisymmetric MHD disk tori, and solve the equations governing a two--dimensional magnetized…
Multidimensional shock-capturing numerical schemes for special relativistic hydrodynamics (RHD) are computationally more expensive than their correspondent Euler versions, due to the nonlinear relations between conservative and primitive…
In this paper, we present exact divergence-free spectral method for solving the incompressible and resistive magneto-hydrodynamic (MHD) equations in two and three dimensions, as well as the efficient solution algorithm and unconditionally…
This paper is denoted to the study of dynamical behavior near explicit finite time blowup solutions for three dimensional incompressible Magnetohydrodynamics (MHD) equations. More precisely, we find a family of explicit finite time blowup…
We present the theory, algorithms and implementation of a parallel finite-volume algorithm for the solution of the incompressible magnetohydrodynamic (MHD) equations using unstructured grids that are applicable for a wide variety of…
The HLLEM scheme is a popular contact and shear preserving approximate Riemann solver for cheap and accurate computation of high speed gasdynamical flows. Unfortunately this scheme is known to be plagued by various forms of numerical shock…
We present a thorough numerical study on the MRI using the smoothed particle magnetohydrodynamics method (SPMHD) with the geometric density average force expression (GDSPH). We perform shearing box simulations with different initial setups…
Compressible magnetohydrodynamic (MHD) turbulence is ubiquitous in astrophysical phenomena ranging from the intergalactic to the stellar scales. In studying them, numerical simulations are nearly inescapable, due to the large degree of…
We present a finite-volume, genuinely 4th-order accurate numerical method for solving the equations of resistive relativistic magnetohydrodynamics (Res-RMHD) in Cartesian coordinates. In our formulation, the magnetic field is evolved in…
We propose a new Harten-Lax-van Leer discontinuities (HLLD) approximate Riemann solver to improve the stability of shocks and the accuracy of low-speed flows in multidimensional magnetohydrodynamic (MHD) simulations. Stringent benchmark…
In this paper, we develop a new method for magnetohydrodynamics (MHD) using smoothed particle hydrodynamics (SPH). To describe MHD shocks accurately, the Godunov method is applied to SPH instead of artificial dissipation terms. In the…
Numerically solving magnetohydrodynamic (MHD) equations faces many challenges: avoiding divergence error, maintaining positivity, and satisfying entropy conditions. Among discontinuous Galerkin (DG) schemes, there has been a modal version…
Performing accurate large eddy simulations in compressible, turbulent magnetohydrodynamics is more challenging than in non-magnetized fluids due to the complex interplay between kinetic, magnetic and internal energy at different scales.…
Magnetic reconnection and non-thermal particle distributions associated with current-driven instabilities are investigated by means of resistive magnetohydrodynamics (MHD) simulations combined with relativistic test particle methods. We…
We propose and analyse a hybrid high-order (HHO) scheme for stationary incompressible magnetohydrodynamics equations. The scheme has an arbitrary order of accuracy and is applicable on generic polyhedral meshes. For sources that are small…
We present numerical simulations and explore scalings and anisotropy of compressible magnetohydrodynamic (MHD) turbulence. Our study covers both gas pressure dominated (high beta) and magnetically dominated (low beta) plasmas at different…
High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its…
Due to increase in computing power, high-order Eulerian schemes will likely become instrumental for the simulations of turbulence and magnetic field amplification in astrophysical fluids in the next years. We present the implementation of a…