Related papers: An Algorithm for Radiation Magnetohydrodynamics Ba…
In this paper we introduce a new method for exact decomposition of propagating, nonlinear magnetohydrodynamic (MHD) disturbances into their component eigenenergies associated with the familiar slow, Alfv\'en, and fast wave eigenmodes, and…
The theory of heat transfer by electromagnetic radiation is based on the radiative transfer equation (RTE) for the radiation intensity, or equivalently on the Boltzmann transport equation (BTE) for the photon distribution. We focus in this…
In this work we develop a class of high-order finite difference weighted essentially non-oscillatory (FD-WENO) schemes for solving the ideal magnetohydrodynamic (MHD) equations in 2D and 3D. The philosophy of this work is to use efficient…
We describe a new hybrid framework to model non-thermal spectral signatures from highly energetic particles embedded in a large-scale classical or relativistic MHD flow. Our method makes use of \textit{Lagrangian} particles moving through…
We present the theory, algorithms and implementation of a parallel finite-volume algorithm for the solution of the incompressible magnetohydrodynamic (MHD) equations using unstructured grids that are applicable for a wide variety of…
Electromagnetic scattering and absorption by material particles is a fundamental physical problem with a broad range of applications, going from laboratory experiments, biology and material sciences, all the way up to environmental studies…
To achieve efficient and accurate long-time integration, we propose a fast, accurate, and stable high-order numerical method for solving fractional-in-space reaction-diffusion equations. The proposed method is explicit in nature and…
The transport of charged particles or photons in a scattering medium can be modelled with a Boltzmann equation. The mathematical treatment for scattering in such scenarios is often simplified if evaluated in a frame where the scattering…
We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial…
We examine the effect of accuracy of high-order spectral element methods, with or without adaptive mesh refinement (AMR), in the context of a classical configuration of magnetic reconnection in two space dimensions, the so-called…
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…
The impressive progress of the kinetic schemes in the solution of gas dynamics problems and the development of effective parallel algorithms for modern high performance parallel computing systems led to the development of advanced methods…
We present a new method for the numerical solution of the radiative-transfer equation (RTE) in multidimensional scenarios commonly encountered in computational astrophysics. The method is based on the direct solution of the Boltzmann…
We present Arepo-MCRT, a novel Monte Carlo radiative transfer (MCRT) radiation-hydrodynamics (RHD) solver for the unstructured moving-mesh code Arepo. Our method is designed for general multiple scattering problems in both optically thin…
Three non-equilibrium photoionisation algorithms for hydrodynamical grid-based simulation codes are compared in terms of accuracy, timestepping criteria, and parallel scaling. Explicit methods with first-order time accuracy for photon…
We present a code for generating synthetic SEDs and intensity maps from Smoothed Particle Hydrodynamics simulation snapshots. The code is based on the Lucy (1999) Monte Carlo Radiative Transfer method, i.e. it follows discrete luminosity…
An algorithm for creating synthetic telescope images of Smoothed Particle Hydrodynamics (SPH) density fields is presented, which utilises the adaptive nature of the SPH formalism in full. The imaging process uses Monte Carlo Radiative…
An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the…
We show that the algorithm based on the weighted essentially nonoscillatory (WENO) scheme with anti-diffusive flux corrections can be used as a solver of the radiative transfer equations. This algorithm is highly stable and robust for…
Modelling the interaction between ionizing photons emitted from massive stars and their environment is essential to further our understanding of galactic ecosystems. We present a hybrid Radiation-Hydrodynamics (RHD) scheme that couples an…