Related papers: Implementation of the Shearing Box Approximation i…
The prevalence and consequences of convection perpendicular to the plane of accretion discs have been discussed for several decades. Recent simulations combining convection and the magnetorotational instability have given fresh impetus to…
In this paper, we present a novel second order in time mixed finite element scheme for the Cahn-Hilliard-Navier-Stokes equations with matched densities. The scheme combines a standard second order Crank-Nicholson method for the…
In this work, a simple fourth-order accurate finite volume semi-discrete scheme is introduced to solve astrophysical magnetohydrodynamics (MHD) problems on Cartesian meshes. Hydrodynamic quantities like density, momentum and energy are…
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must confront the challenge of controlling errors in the discrete divergence of the magnetic field. One approach that has been…
Fully adaptive computations of the resistive magnetohydrodynamic (MHD) equations are presented in two and three space dimensions using a finite volume discretization on locally refined dyadic grids. Divergence cleaning is used to control…
This paper presents applications of weighted meshless scheme for conservation laws to the Euler equations and the equations of ideal magnetohydrodynamics. The divergence constraint of the latter is maintained to the truncation error by a…
The stability of a differentially rotating fluid subject to its own gravity is a problem with applications across wide areas of astrophysics--from protoplanetary discs (PPDs) to entire galaxies. The shearing box formalism offers a…
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…
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.…
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…
In this paper, an implicit nonsymplectic exact energy-preserving integrator is specifically designed for a ten-dimensional phase-space conservative Hamiltonian system with five degrees of freedom. It is based on a suitable…
An algorithm for simulating self-gravitating cosmological astrophysical fluids is presented. The advantages include a large dynamic range, parallelizability, high resolution per grid element and fast execution speed. The code is based on a…
Recently, we developed a pair of meshless finite-volume Lagrangian methods for hydrodynamics: the 'meshless finite mass' (MFM) and 'meshless finite volume' (MFV) methods. These capture advantages of both smoothed-particle hydrodynamics…
We describe implementation and tests of sink particle algorithms in the Eulerian grid-based code Athena. Introduction of sink particles enables long-term evolution of systems in which localized collapse occurs, and it is impractical (or…
In this paper, we study the charged-particle dynamics under strong magnetic field in a toroidal axi-symmetric geometry. Using modulated Fourier expansions of the exact and numerical solutions, the long-term drift motion of the exact…
Port-Hamiltonian systems provide a highly-structured framework for modeling of physical systems. By definition, they encode a balance equation relating energy changes to supplied and dissipated energy. Capturing this energy balance in…
We analyze a suite of global magnetohydrodynamic (MHD) accretion disk simulations in order to determine whether scaling laws for turbulence driven by the magnetorotational instability, discovered via local shearing box studies, are globally…
In this work, we design and analyze semi/fully-discrete virtual element approximations for the time-dependent Navier--Stokes-Cahn--Hilliard equations, modeling the dynamics of two-phase incompressible fluid flows with diffuse interfaces. A…
A new methodology is developed to integrate numerically the equations of motion for classical many-body systems in molecular dynamics simulations. Its distinguishable feature is the possibility to preserve, independently on the size of the…
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…