Related papers: Simulation Scenarios in One-Dimensional Self-Gravi…
Quantum computing can be used to speed up the simulation time (more precisely, the number of queries of the algorithm) for physical systems; one such promising approach is the Hamiltonian simulation (HS) algorithm. Recently, it was proven…
We investigate the evolution of the phase-space distribution function around slightly perturbed stationary states and the process of violent relaxation in the context of the dissipationless collapse of an isolated spherical self-gravitating…
A cascadic multigrid method is proposed for eigenvalue problems based on the multilevel correction scheme. With this new scheme, an eigenvalue problem on the finest space can be solved by smoothing steps on a series of multilevel finite…
We describe a method of Monte-Carlo simulations of simplicial quantum gravity coupled to matter fields. We concentrate mainly on the problem of implementing effectively the random, dynamical triangulation and building in a detailed-balance…
The current paradigm for understanding galaxy formation in the universe depends on the existence of self-gravitating collisionless dark matter. Modeling such dark matter systems has been a major focus of astrophysicists, with much of that…
We present a general strategy to solve the stationary Schr\"odinger-Poisson (SP) system of equations for multistates with axial symmetry. The approach allows us to obtain the well known single and multistate solutions with spherical…
We construct an algorithm to determine all stationary axi-symmetric solutions of 3-dimensional Einstein gravity with a minimally coupled self-interacting scalar field. We holographically renormalize the theory and evaluate then the on-shell…
We describe an efficient numerical method for simulating the dynamics and steady states of collective spin systems in the presence of dephasing and decay. The method is based on the Schwinger boson representation of spin operators and uses…
A particle method for reproducing the phase space of collisionless stellar systems is described. The key idea originates in Liouville's theorem which states that the distribution function (DF) at time t can be derived from tracing necessary…
This paper presents a p-adaptive high-order hybridizable discontinuous Galerkin spectral element method (HDG-SEM) for solving the Poisson equation in electrostatic plasma simulations using particle-in-cell (PIC) schemes. This approach…
We construct a pseudospectral method for the solution of time-dependent, non-linear partial differential equations on a three-dimensional spherical shell. The problem we address is the treatment of tensor fields on the sphere. As a test…
High-dimensional simulation optimization is notoriously challenging. We propose a new sampling algorithm that converges to a global optimal solution and suffers minimally from the curse of dimensionality. The algorithm consists of two…
We present a numerical solver for plasma dynamics simulations in Hall magnetohydrodynamic (HMHD) approximation in one, two and three dimensions. We consider both isotropic and anisotropic thermal pressure cases, where a general gyrotropic…
This paper is concerned with the existence of forward self-similar solutions to the three-dimensional Magnetohydrodynamic-Boussinesq (MHD-Boussinesq) system with newtonian gravitational field. By employing a blow-up argument and the…
We present the construction of stationary boson-fermion spherically symmetric configurations governed by Newtonian gravity. Bosons are described in the Gross-Pitaevskii regime and fermions are assumed to obey Euler equations for an inviscid…
The Poisson-Boltzmann equation is widely used to model molecular electrostatics; however, it is usually solved in linearised form because the sinh nonlinearity is challenging, limiting its applicability in highly charged systems such as…
We introduce a novel numerical method to integrate partial differential equations representing the Hamiltonian dynamics of field theories. It is a multi-symplectic integrator that locally conserves the stress-energy tensor with an excellent…
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 develop a universally applicable embedded boundary finite difference method, which results in a symmetric positive definite linear system and does not suffer from small cell stiffness. Our discretization is efficient for the wave, heat…
A well-balanced second order finite volume central scheme for the magnetohydrodynamic (MHD) equations with gravitational source term is developed in this paper. The scheme is an unstaggered central scheme that evolves the numerical solution…