Related papers: A conservative orbital advection scheme for simula…
In numerical models of thin astrophysical disks that use an Eulerian scheme, gas orbits supersonically through a fixed grid. As a result the time step is sharply limited by the Courant condition. Also, because the mean flow speed with…
We present a description of the adaptive mesh refinement (AMR) implementation of the PLUTO code for solving the equations of classical and special relativistic magnetohydrodynamics (MHD and RMHD). The current release exploits, in addition…
We present an efficient and simple modification of the standard transport algorithm used in explicit eulerian fixed polar grid codes, aimed at getting rid of the average azimuthal velocity when applying the Courant condition. This results…
We describe a newly developed hydrodynamic code for studying accretion disk processes. The numerical method uses a finite volume, nonlinear, Total Variation Diminishing (TVD) scheme to capture shocks and control spurious oscillations. It is…
Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with…
Due to the prevalence of magnetic fields in astrophysical environments, magnetohydrodynamic (MHD) simulation has become a basic tool for studying astrophysical fluid dynamics. To further advance the precision of MHD simulations, we have…
We present a new numerical code, PLUTO, for the solution of hypersonic flows in 1, 2 and 3 spatial dimensions and different systems of coordinates. The code provides a multi-physics, multi-algorithm modular environment particularly oriented…
The transport of energy through radiation is very important in many astrophysical phenomena. In dynamical problems the time-dependent equations of radiation hydrodynamics have to be solved. We present a newly developed…
We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on protoplanetary disks physics and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms…
Our aim is to study the thermal and dynamical evolution of protoplanetary disks in global simulations, including the physics of radiation transfer and magneto-hydrodynamic (MHD) turbulence caused by the magneto-rotational instability. We…
We present a numerical code for radiation hydrodynamics designed as a module for the freely available PLUTO code. We adopt a gray approximation and include radiative transfer following a two-moment approach by imposing the M1 closure to the…
We present a constrained transport (CT) algorithm for solving the 3D ideal magnetohydrodynamic (MHD) equations on a moving mesh, which maintains the divergence-free condition on the magnetic field to machine-precision. Our CT scheme uses an…
We employ the PLUTO code for computational astrophysics to assess and compare the validity of different numerical algorithms on simulations of the magneto-rotational instability in 3D accretion disks. In particular we stress on the…
The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in…
We test a new "hybrid" scheme for simulating dynamical fluid flows in which cylindrical components of the momentum are advected across a rotating Cartesian coordinate mesh. This hybrid scheme allows us to conserve angular momentum to…
Radiation is crucial not only for observing astrophysical objects, but also for transporting energy and momentum. However, accurate on-the-fly radiation transport in astrophysical simulations is challenging and computationally expensive.…
We present a new algorithm for the identification and physical characterization of current sheets and reconnection sites in 2D and 3D large scale relativisticmagnetohydrodynamic numerical simulations. This has been implemented in the PLUTO…
In this paper, we investigate the application of the maximum principle preserving (MPP) parametrized flux limiters to the high order finite volume scheme with Runge-Kutta time discretization for solving convection dominated problems. Such…
We present preliminary performance results of gPLUTO, the new GPU-optimized implementation of the PLUTO code for computational plasma astrophysics. Like its predecessor, gPLUTO employs a finite-volume formulation to numerically solve the…
We describe the implementation of the shearing box approximation for the study of the dynamics of accretion disks in the Athena magnetohydrodynamics (MHD) code. Second-order Crank-Nicholson time differencing is used for the Coriolis and…