Related papers: ExaGRyPE: Numerical General Relativity Solvers Bas…
By applying the method of moving frames modelling one and two dimensional local anisotropies we construct new solutions of Einstein equations on pseudo-Riemannian spacetimes. The first class of solutions describes non-trivial deformations…
Iterative solvers are widely used to accurately simulate physical systems. These solvers require initial guesses to generate a sequence of improving approximate solutions. In this contribution, we introduce a novel method to accelerate…
We present the development and benchmarking of Poisson solvers for graphics processing units (GPUs). Implemented in the Astaroth platform, the solvers feature high computational efficiency. We present novel combinations of discretizations…
We numerically investigate the dynamics near black hole formation of solutions to the Einstein--Vlasov system in axisymmetry. Our results are obtained using a particle-in-cell and finite difference code based on the $(2+1)+1$ formulation of…
Rapid and reliable solvers for parametric partial differential equations (PDEs) are needed in many scientific and engineering disciplines. For example, there is a growing demand for composites and architected materials with heterogeneous…
We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems…
We present a new class of near-horizon geometries which solve Einstein's vacuum equations, including a negative cosmological constant, in all even dimensions greater than four. Spatial sections of the horizon are inhomogeneous S^2-bundles…
In this paper we propose an extension of the generalized Lagrangian multiplier method (GLM) of Munz et al. (JCP 2000, JCP 2002), which was originally conceived for the numerical solution of the Maxwell and MHD equations with divergence-type…
GRTresna is a multigrid solver designed to solve the constraint equations for the initial data required in numerical relativity simulations. In particular, it is focussed on scenarios with fundamental fields around black holes and…
Coarse Grid Projection (CGP) methodology is used to accelerate the computations of sets of decoupled nonlinear evolutionary and linear static equations. In CGP, the linear equations are solved on a coarsened mesh compared to the nonlinear…
Historically, the interpolation of large geophysical datasets has been tackled using methods like Optimal Interpolation (OI) or model-based data assimilation schemes. However, the recent connection between Stochastic Partial Differential…
In this study, we present a novel computational framework that integrates the finite volume method with graph neural networks to address the challenges in Physics-Informed Neural Networks(PINNs). Our approach leverages the flexibility of…
We present and discuss the characteristics and performances, both in term of computational speed and precision, of a numerical code which numerically integrates the equation of motions of N 'particles' interacting via Newtonian gravitation…
Physics informed neural networks (PINNs) have emerged as a powerful tool to provide robust and accurate approximations of solutions to partial differential equations (PDEs). However, PINNs face serious difficulties and challenges when…
The Event Horizon Telescope (EHT) has recently released high-resolution images of accretion flows onto two supermassive black holes. Our physical understanding of these images depends on accuracy and precision of numerical models of plasma…
Previous math word problem solvers following the encoder-decoder paradigm fail to explicitly incorporate essential math symbolic constraints, leading to unexplainable and unreasonable predictions. Herein, we propose Neural-Symbolic Solver…
We have developed a gravity solver based on combining the well developed Particle-Mesh (PM) method and TREE methods. It is designed for and has been implemented on parallel computer architectures. The new code can deal with tens of millions…
Elliptic partial differential equations (PDEs) frequently arise in continuum descriptions of physical processes relevant to science and engineering. Multilevel preconditioners represent a family of scalable techniques for solving discrete…
The 'anholonomic frame' method (see gr-qc/0005025, gr-qc/0001060 and hep-th/0110250) is applied for constructing new classes of exact solutions of vacuum Einstein equations with off-diagonal metrics in 4D and 5D gravity. We examine several…
Generalized Disjunctive Programming (GDP) provides a natural framework for optimization models that combine logical decisions with nonlinear constraints. The Hull Reformulation (HR) is attractive because it yields tight continuous…