Related papers: Shearingbox-implementation for the central-upwind,…
The nonlinear evolution and the saturation mechanism of the magnetorotational instability (MRI) are investigated using three-dimensional resistive MHD simulations. A local shearing box is used for our numerical analysis and the simulations…
Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of…
We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme. The original scheme has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas,…
Numerical simulations including magnetic fields have become important in many fields of astrophysics. Evolution of magnetic fields by the constrained transport algorithm preserves magnetic divergence to machine precision, and thus…
Most of our knowledge on the nonlinear development of the magneto-rotational instability (MRI) relies on the results of numerical simulations employing the shearing box (SB) approximation. A number of difficulties arising from this approach…
In this paper, we present a new method to perform numerical simulations of astrophysical MHD flows using the Adaptive Mesh Refinement framework and Constrained Transport. The algorithm is based on a previous work in which the MUSCL--Hancock…
This letter investigates the transport properties of MHD turbulence induced by the magnetorotational instability at large Reynolds numbers Re when the magnetic Prandtl number Pm is larger than unity. Three MHD simulations of the…
The local properties of turbulence driven by the magnetorotational instability (MRI) in rotating, shearing flows are studied in the framework of a shearing-box model. Based on numerical simulations, we propose that the MRI-driven turbulence…
We present a new magnetohydrodynamic (MHD) code for the simulation of wave propagation in the solar atmosphere, under the effects of electrical resistivity, but not dominant, and heat transference in a uniform 3D grid. The code is based on…
We describe a novel Godunov-type numerical method for solving the equations of resistive relativistic magnetohydrodynamics. In the proposed approach, the spatial components of both magnetic and electric fields are located at zone interfaces…
A description is given for preserving ${\bmsy\nabla}\cdot{\vec B}=0$ in a magnetohydrodynamic (MHD) code that employs the upwind, Total Variation Diminishing (TVD) scheme and the Strang-type operator splitting for multi-dimensionality. The…
We describe the special relativistic extension of the CRONOS code, which has been used for studies of gamma-ray binaries in recent years. The code was designed to be easily adaptable, allowing the user to easily change existing…
The hydrodynamic response of the inviscid small shearing box model of a midplane section of a rotationally supported astrophysical disk is examined. An energy functional ${\cal E}$ is formulated for the general nonlinear problem. It is…
Self-gravity plays an important role in the evolution of rotationally supported systems such as protoplanetary disks, accretion disks around black holes, or galactic disks, as it can both feed turbulence or lead to gravitational…
Magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability can provide diffusive transport of angular momentum in astrophysical disks, and a widely studied computational model for this process is the ideal, stratified,…
The ideal Magnetohydrodynamics (MHD) equations are challenging because one needs to maintain the divergence free condition, $\nabla \cdot \Bv = 0$. Many numerical methods have been developed to enforce this condition. In this work, we…
We present a new general relativistic magnetohydrodynamics (GRMHD) code integrated into the Athena++ framework. Improving upon the techniques used in most GRMHD codes, ours allows the use of advanced, less diffusive Riemann solvers, in…
In this paper we describe a numerical method designed for modelling different kinds of astrophysical flows in three dimensions. Our method is a standard explicit finite difference method employing the local shearing-box technique. To model…
In this paper, we utilize the maximum-principle-preserving flux limiting technique, originally designed for high order weighted essentially non-oscillatory (WENO) methods for scalar hyperbolic conservation laws, to develop a class of high…
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…