Related papers: Orbital Advection by Interpolation: A Fast and Acc…
Explicit numerical computations of super-fast differentially rotating disks are subject to the time-step constraint imposed by the Courant condition. When the bulk orbital velocity largely exceeds any other wave speed the time step is…
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…
We describe a method for incorporating ambipolar diffusion in the strong coupling approximation into a multidimensional magnetohydrodynamics code based on the total variation diminishing scheme. Contributions from ambipolar diffusion terms…
This work considers flows from an accretion disk corotating with the aligned dipole magnetic field of a rotating star. Ideal magnetohydrodynamics (MHD) is assumed with the pressure and density related as $p \propto \rho^\gamma$ and with…
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…
The study of many astrophysical flows requires computational algorithms that can capture high Mach number flows, while resolving a large dynamic range in spatial and density scales. In this paper we present a novel method, RAM: Rapid…
Accretion disks are three-dimensional, turbulent, often self-gravitating, magnetohydrodynamic flows, which can be modeled in detail with numerical simulations. In this paper, we present a new algorithm that is based on a spectral…
Embedded planets disturb the density structure of the ambient disk and gravitational back-reaction will induce possibly a change in the planet's orbital elements. The accurate determination of the forces acting on the planet requires…
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…
Insoluble surfactants adsorbed at liquid-liquid or gas-liquid interfaces alter interfacial tension, leading to variations in the normal stress jump and generating tangential Marangoni stresses that can dramatically affect the flow dynamics.…
High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $\Delta x$) with only a moderate increase to computational expense. Significant…
Rotationally supported, cold, gaseous disks are ubiquitous in astrophysics and appear in a diverse set of systems, such as protoplanetary disks, accretion disks around black holes, or large spiral galaxies. Capturing the gas dynamics…
The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO…
We present a two-dimensional Cartesian code based on high order discontinuous Galerkin methods, implemented to run in parallel over multiple GPUs. A simple planet-disc setup is used to compare the behaviour of our code against the behaviour…
Implicit time-stepping for advection is applied locally in space and time where Courant numbers are large, but standard explicit time-stepping is used for the remaining solution which is typically the majority. This adaptively implicit…
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…
Magneto-hydrodynamics is one of the foremost models in plasma physics with applications in inertial confinement fusion, astrophysics and elsewhere. Advanced numerical methods are needed to get an insight into the complex physical phenomena.…
Magnetic fields play an important role in many astrophysical systems and a detailed understanding of their impact on the gas dynamics requires robust numerical simulations. Here we present a new method to evolve the ideal…
We study a truncated accretion disk using a well-resolved, semi-global magnetohydrodynamic simulation that is evolved for many dynamical times (6096 inner disk orbits). The spectral properties of hard state black hole binary systems and…