Related papers: Enzo: An Adaptive Mesh Refinement Code for Astroph…

Hydrodynamic cosmological simulations at present usually employ either the Lagrangian SPH technique, or Eulerian hydrodynamics on a Cartesian mesh with adaptive mesh refinement. Both of these methods have disadvantages that negatively…

We present a novel shape-approximating anisotropic re-meshing algorithm as a geometric generalization of the adaptive moving mesh method. Conventional moving mesh methods reduce the interpolation error of a mesh that discretizes a given…

We present the extension of the differentiable hydrodynamics code, diffhydro, enabling scalable PDE-constrained inference and integrated hybrid physics-ML models for a wide range of astrophysical applications. New physics additions include…

We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step…

We present a new numerical code, ECHO, based on an Eulerian Conservative High Order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a…

In certain astrophysical systems the commonly employed ideal magnetohydrodynamics (MHD) approximation breaks down. Here, we introduce novel explicit and implicit numerical schemes of ohmic resistivity terms in the moving-mesh code AREPO. We…

Highly accurate simulation of plasma transport is needed for the successful design and operation of magnetically confined fusion reactors. Unfortunately, the extreme anisotropy present in magnetized plasmas results in thin boundary layers…

Many systems of current interest in relativistic astrophysics require a knowledge of radiative transfer in a magnetized gas flowing in a strongly-curved, dynamical spacetime. Such systems include coalescing compact binaries containing…

Collsionless astrophysical and space plasmas cover regions that typically display a separation of scales that exceeds any code's capabilities. To help address this problem, the muphyII code utilizes a hierarchy of models with different…

There has been a lot of research interest in modified gravity theories which utilise the Vainshtein mechanism to recover standard general relativity in regions with high matter density, such as the Dvali-Gabadadze-Porrati and Galileon…

Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in…

Epsilon-near-zero (ENZ) metamaterials represent a powerful toolkit for selectively transmitting and localizing light through cavity resonances, enabling the study of mesoscopic phenomena and facilitating the design of photonic devices. In…

We study dendritic microstructure evolution using an adaptive grid, finite element method applied to a phase-field model. The computational complexity of our algorithm, per unit time, scales linearly with system size, rather than the…

Computational fluid dynamics is a crucial tool to theoretically explore the cosmos. In the last decade, we have seen a substantial methodological diversification with a number of cross-fertilizations between originally different methods.…

We present a new three-dimensional radiative transfer (RT) code, RADAMESH, based on a ray-tracing, photon-conserving and adaptive (in space and time) scheme. RADAMESH uses a novel Monte Carlo approach to sample the radiation field within…

The tremendous challenge of comparing our theoretical models with the gravitational-wave observations in the new era of multimessenger astronomy requires accurate and fast numerical simulations of complicated physical systems described by…

We present Rico, a code designed to compute the ionization fraction of the Universe during the epoch of hydrogen and helium recombination with an unprecedented combination of speed and accuracy. This is accomplished by training the machine…

We present a new numerical code which solves the general relativistic magneto-hydrodynamics (GRMHD) equations coupled to the Einstein equations for the evolution of a dynamical spacetime within the conformally-flat approximation. This code…

Numerical methods have become a powerful tool for research in astrophysics, but their utility depends critically on the availability of suitable simulation codes. This calls for continuous efforts in code development, which is necessitated…

We present the first implementation of the Active Flux method on adaptively refined Cartesian grids. The Active Flux method is a third order accurate finite volume method for hyperbolic conservation laws, which is based on the use of point…