Related papers: A self-gravity module for the PLUTO code
We obtain the Maxwell-J\"uttner distribution function at first order in the post-Newtonian approximation within the framework of general relativity. Taking into account the aforesaid distribution function, we compute the particle four-flow…
We describe a finite-volume method for solving the Poisson equation on oct-tree adaptive meshes using direct solvers for individual mesh blocks. The method is a modified version of the method presented by Huang and Greengard (2000), which…
We investigate the Gravitational Poissonian Spontaneous Localization (GPSL) model, a hybrid classical-quantum model in which classical Newtonian gravity emerges from stochastic collapses of the mass density operator, and consistently…
We compute the structure of a self-gravitating torus with polytropic equation-of-state (EOS) rotating in an imposed centrifugal potential. The Poisson-solver is based on isotropic multigrid with optimal covering factor (fluid…
A method is presented to include irregular domain boundaries in a geometric multigrid solver. Dirichlet boundary conditions can be imposed on an irregular boundary defined by a level set function. Our implementation employs quadtree/octree…
We generalize the interpolative separable density fitting (ISDF) method, used for compressing the four-index electron repulsion integral (ERI) tensor, to incorporate adaptive real space grids for potentially highly localized single-particle…
In this paper, we use only the equation of motion for an interacting system of gravity, dilaton and antisymmetric tensor to study the soliton solutions by making use of a Poincar\'e invariant ansatz. We show that the system of equations are…
We deal with accelerating the solution of a sequence of large linear systems solved by preconditioned conjugate gradient method (PCG). The sequence originates from time-stepping within a simulation of an unsteady incompressible flow. We…
We describe an implementation to solve Poisson's equation for an isolated system on a unigrid mesh using FFTs. The method solves the equation globally on mesh blocks distributed across multiple processes on a distributed-memory parallel…
We present a numerical code for radiation hydrodynamics designed as a module for the freely available PLUTO code. We adopt a gray approximation and include radiative transfer following a two-moment approach by imposing the M1 closure to the…
We discuss the scalable parallel solution of the Poisson equation within a Particle-In-Cell (PIC) code for the simulation of electron beams in particle accelerators of irregular shape. The problem is discretized by Finite Differences.…
We present a numerical method to solve the equations of general relativistic hydrodynamics in a given external gravitational field. The method is based on a generalization of Roe's approximate Riemann solver for the non relativistic Euler…
We study periodic, two-dimensional, gravity-capillary traveling wave solutions to a viscous shallow water system posed on an inclined plane. While thinking of the Reynolds and Bond numbers as fixed and finite, we vary the speed of the…
In the Klein-Gordon equation, quantum and relativistic parameters are strongly coupled, which poses significant analytical challenges in the derivation and analysis of related classical fluid models. In this paper, starting from the…
A new conservative symmetry-preserving second-order time-accurate PISO-based pressure-velocity coupling for solving the incompressible Navier-Stokes equations on unstructured collocated grids is presented in this paper. This new method for…
We present a polynomial multigrid method for nodal interior penalty and local discontinuous Galerkin formulations of the Poisson equation on Cartesian grids. For smoothing we propose two classes of overlapping Schwarz methods. The first…
We propose in this paper a Proper Generalized Decomposition (PGD) solver for reduced-order modeling of linear elastodynamic problems. It primarily focuses on enhancing the computational efficiency of a previously introduced PGD solver based…
We present a well-balanced nodal discontinuous Galerkin (DG) scheme for compressible Euler equations with gravity. The DG scheme makes use of discontinuous Lagrange basis functions supported at Gauss-Lobatto-Legendre (GLL) nodes together…
We have developed efficient numerical algorithms for solving 3D steady-state Poisson-Nernst-Planck (PNP) equations with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are…
Variable density incompressible flows are governed by parabolic equations. The artificial compressibility method makes these equations hyperbolic-type, which means that they can be solved using techniques developed for compressible flows,…