Related papers: A self-gravity module for the PLUTO code
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 obtain equilibrium solutions for rotating compact stars including the special relativistic effects. The gravity is assumed to be Newtonian, but we used the active mass density, which takes into account all the energies such as motions of…
In ultra-relativistic collisions of heavy ions, the strongly interacting Quark Gluon Plasma (sQGP) is created. The fluid nature of the sQGP was one of the important discoveries of high energy heavy ion physics in the last decades.…
We present a new multilevel method for calculating Poisson's equation, which often arises form electrostatic problems, by using hierarchical loop bases. This method, termed hierarchical Loop basis Poisson Solver (hieLPS), extends previous…
Compressible Euler-Poisson equations are the standard self-gravitating models for stellar dynamics in classical astrophysics. In this article, we construct periodic solutions to the isothermal ($\gamma=1$) Euler-Poisson equations in $R^{2}$…
A new approximation method for inverting the Poisson's equation is presented for a continuously distributed and finite-sized source in an unbound domain. The advantage of this image multipole method arises from its ability to place the…
Using numerical integration, in 1969 Penston [22] and Larson [17] independently discovered a self-similar solution describing the collapse of a self-gravitating asymptotically flat fluid with the isothermal equation of state $p=k\varrho$,…
We present a new method for solving the relativistic Vlasov--Maxwell system of equations, applicable to a wide range of extreme high-energy-density astrophysical and laboratory environments. The method directly discretizes the kinetic…
In recent years several research efforts focused on the development of high-order discontinuous Galerkin (dG) methods for scale resolving simulations of turbulent flows. Nevertheless, in the context of incompressible flow computations, the…
We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a…
Many astrophysical systems can only be accurately modelled when the behaviour of their baryonic gas components is well understood. The residual distribution (RD) family of partial differential equation (PDE) solvers produce approximate…
In this paper, we report the advantages of using a stochastic algorithm in the context of mineral exploration based on gravity measurements. This approach has the advantage over deterministic methods in that it allows one to find the…
We present a matrix-free GPU multigrid preconditioner with algebraically consistent coarsening for solving Poisson equations on adaptive octree grids with irregular domains. Within uniform-resolution regions, the coarsening satisfies the…
An important aspect of computational fluid dynamics is related to the determination of the fluid pressure in isothermal incompressible fluids. In particular this concerns the construction of an exact evolution equation for the fluid…
In this paper we consider Runge-Kutta discontinuous Galerkin (RKDG) schemes for Vlasov-Poisson systems that model collisionless plasmas. One-dimensional systems are emphasized. The RKDG method, originally devised to solve conservation laws,…
We introduce a new code, ECOSMOG, to run N-body simulations for a wide class of modified gravity and dynamical dark energy theories. These theories generally have one or more new dynamical degrees of freedom, the dynamics of which are…
We present an algorithm based on numerical techniques that have become standard for solving nonlinear integral equations: Newton's method, homotopy continuation, the multilevel method and random projection to solve the inversion problem…
We investigate the existence of stable soliton solution in a system described by complex Ginzburg-Landau (CGL) equation with near parity reflection - time reversal ($\mathcal{PT}$) symmetric Rosen-Morse potential. In this study, the…
In this paper, we study the Euler and Euler-Poisson equations in $R^{N}$, with multiple $\gamma$-law for pressure function: \begin{equation} P(\rho)=e^{s}\sum_{j=1}^{m}\rho^{\gamma_{j}}, \end{equation} where all…
Reciprocal space methods for solving Poisson's equation for finite charge distributions are investigated. Improvements to previous proposals are presented, and their performance is compared in the context of a real-space density functional…