Related papers: An implicit numerical algorithm for solving the ge…
We review the issue of steady spherically symmetric accretion onto a renormalization group improved Schwarzschild space-time which is solution to an asymptotically safe theory (AS) containing high-derivative terms. We use a Hamiltonian…
In this paper, we propose a robust solver for the finite element discrete problem of the stationary incompressible magnetohydrodynamic (MHD) equations in three dimensions. By the mixed finite element method, both the velocity and the…
The high-order numerical solution of the non-linear shallow water equations (and of hyperbolic systems in general) is susceptible to unphysical Gibbs oscillations that form in the proximity of strong gradients. The solution to this problem…
In this work, we design an entropy stable, finite volume approximation for the shallow water magnetohydrodynamics (SWMHD) equations. The method is novel as we design an affordable analytical expression of the numerical interface flux…
We re-derive the equation for vertical hydrodynamical equilibrium of stationary thin and slim accretion disks in the Kerr spacetime. All the previous derivations have been unsatisfactory, yielding unphysical singularities. Our equation is…
We present a novel high-order nodal artificial viscosity approach designed for solving Magnetohydrodynamics (MHD) equations. Unlike conventional methods, our approach eliminates the need for ad hoc parameters. The viscosity is…
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…
General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black hole can be divided into inflow and outflow parts, according to the result of the competition between the black hole gravity and…
A technique is proposed for constructing three-dimensional solutions comlplying with the self-consistent magnetohydrostatic (MHS) equations and with observations along the line of sight of the magnetic field at the photosphere. The…
A numerical scheme is developed for systems of conservation laws on manifolds which arise in high speed aerodynamics and magneto-aerodynamics. The systems are presented in an arbitrary coordinate system on the manifold and involve source…
The characteristic decomposition for GRMHD in the comoving frame of the fluid has been known for a long time. However, it has not been known in the coordinate frame of the simulation and in terms of the conserved variables evolved in…
In two previous papers (Price & Monaghan 2004a,b) (papers I,II) we have described an algorithm for solving the equations of Magnetohydrodynamics (MHD) using the Smoothed Particle Hydrodynamics (SPH) method. The algorithm uses dissipative…
Recently, compressible magnetohydrodynamics (MHD) has been elegantly formulated in terms of Lie derivatives. This paper exploits the geometrical properties of the Lie bracket to give new insights into the properties of compressible MHD…
Many astrophysical systems of interest to numerical relativity-such as rapidly rotating stars, black hole accretion disks, and core-collapse supernovae-exhibit near-symmetries. These systems generally consist of a strongly gravitating…
A fully implicit finite difference scheme has been developed to solve the hydrodynamic equations coupled with radiation transport. Solution of the time dependent radiation transport equation is obtained using the discrete ordinates method…
The radiation hydrodynamics equations for smoothed particle hydrodynamics are derived by operator splitting the radiation and hydrodynamics terms, including necessary terms for material motion, and discretizing each of the sets of equations…
This paper presents a general positivity-preserving algorithm for implicit high-order finite volume schemes solving Euler and Navier-Stokes equations. Previous positivity-preserving algorithms are mainly based on mathematical analyses,…
We present a detailed three dimensional magnetohydrodynamic (MHD) simulation describing the inner region of a disk accreting onto a black hole. To avoid the technical complications of general relativity, the dynamics are treated in…
The relativistic method of moments is one of the most successful approaches to extract second order viscous hydrodynamics from a kinetic underlying background. The equations can be systematically improved to higher order, and they have…
In this paper, we present the stationary axisymmetric configuration of a resistive magnetised thick accretion disc in the vicinity of external gravity and intrinsic dipolar magnetic field of a slowly rotating black hole. The plasma is…