Related papers: Hyperbolic Self-Gravity Solver for Large Scale Hyd…
We present a new methodology for simulating self-gravitating general-relativistic fluids. In our approach the fluid is modelled by means of Lagrangian particles in the framework of a general-relativistic (GR) Smooth Particle Hydrodynamics…
We describe a new hybrid N-body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and…
Investigating the evolution of disk galaxies and the dynamics of proto-stellar disks can involve the use of both a hydrodynamical and a Poisson solver. These systems are usually approximated as infinitesimally thin disks using two-…
We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied $f(R)$ gravity models. High resolution simulations for such models are…
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve…
Hydrodynamical simulations are the most accurate way to model structure formation in the universe, but they often involve a large number of astrophysical parameters modeling subgrid physics, in addition to cosmological parameters. This…
We develop a numerical strategy to solve multi-dimensional Poisson equations on dynamically adapted grids for evolutionary problems disclosing propagating fronts. The method is an extension of the multiresolution finite volume scheme used…
The axisymmetric form of the hydrodynamic equations within the smoothed particle hydrodynamics (SPH) formalism is presented and checked using idealized scenarios taken from astrophysics (free fall collapse, implosion and further pulsation…
We describe and demonstrate a method for increasing the resolution locally in a Smoothed Particle Hydrodynamic (SPH) simulation, by splitting particles. We show that in simulations of self-gravitating collapse (of the sort which are…
We propose a novel mathematical method to construct an exact polytropic sphere in self-gravitating hydrostatic equilibrium, improving the non-linear Poisson equation. The central boundary condition for the present equation requires a ratio…
We present three-dimensional simulations of Einstein equations implementing a symmetric hyperbolic system of equations with dynamical lapse. The numerical implementation makes use of techniques that guarantee linear numerical stability for…
We present a simple and effective multigrid-based Poisson solver of second-order accuracy in both gravitational potential and forces in terms of the one, two and infinity norms. The method is especially suitable for numerical simulations…
Hydromagnetic turbulence produced during phase transitions in the early universe can be a powerful source of stochastic gravitational waves (GWs). GWs can be modelled by the linearised spatial part of the Einstein equations sourced by the…
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…
In self-consistent N-body simulations of collisionless systems, gravitational interactions are modified on small scales to remove singularities and simplify the task of numerically integrating the equations of motion. This `gravitational…
Volume-preserving hyperelastic materials are widely used to model near-incompressible materials such as rubber and soft tissues. However, the numerical simulation of volume-preserving hyperelastic materials is notoriously challenging within…
We study the evolution of cosmological perturbations, using a hybrid approximation scheme which upgrades the weak-field limit of Einstein's field equations to account for post-Newtonian scalar and vector metric perturbations and for…
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…
Using a generalization of the Madelung transformation, we derive the hydrodynamic representation of the Klein-Gordon-Einstein equations in the weak field limit. We consider a complex self-interacting scalar field with an arbitrary potential…
Simulating and predicting dynamics of quantum many-body systems is extremely challenging, even for state-of-the-art computational methods, due to the spread of entanglement across the system. However, in the long-wavelength limit, quantum…