Related papers: Binary black hole simulation with an adaptive fini…
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…
Choptuik's critical phenomena in general relativity is revisited in the affine-null metric formulation of Einstein's equations for a massless scalar field in spherical symmetry. Numerical solutions are obtained by evolution of initial data…
There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to…
We present a new computational framework for the Galerkin-collocation method for double domain in the context of ADM 3+1 approach in numerical relativity. This work enables us to perform high resolution calculations for initial sets of two…
For an arbitrary strong, spherically symmetric super-horizon curvature perturbation, we present analytical solutions of the Einstein equations in terms of asymptotic expansion over the ratio of the Hubble radius to the length-scale of the…
Errors due to imperfect boundary conditions in numerical relativity simulations of binary black holes can produce unphysical reflections of gravitational waves which compromise the accuracy of waveform predictions, especially for…
We find a plethora of new analytic black holes and globally regular horizonless spacetimes in three dimensions. The solutions involve a single real scalar field $\phi$ which always admits a magnetic-like expression proportional to the…
We consider the Galerkin boundary element method (BEM) for weakly-singular integral equations of the first-kind in 2D. We analyze some residual-type a posteriori error estimator which provides a lower as well as an upper bound for the…
The well-known $(2+1)$-dimensional Reissner-Nordstrom (BTZ) black hole can be generalized to three dimensional Einstein-nonlinear electromagnetic field, motivated from obtaining a finite value for the self-energy of a pointlike charge.…
We extend a recently developed numerical code to obtain stationary, axisymmetric solutions that describe rotating black hole spacetimes in a wide class of modified theories of gravity. The code utilizes a relaxed Newton-Raphson method to…
A conforming discontinuous Galerkin (DG) finite element method has been introduced in [21] on simplicial meshes, which has the flexibility of using discontinuous approximation and the simplicity in formulation of the classic continuous…
We present a method for computing the evolution of a spacetime containing a massive particle and a black hole. The essential idea is that the gravitational field is evolved using full numerical relativity, with the particle generating a…
The long-term evolution of massive black hole binaries at the centers of galaxies is studied in a variety of physical regimes, with the aim of resolving the ``final parsec problem,'' i.e., how black hole binaries manage to shrink to…
This article is concerned with the numerical solution of convex variational problems. More precisely, we develop an iterative minimisation technique which allows for the successive enrichment of an underlying discrete approximation space in…
Nonlinearly scalarized black holes are investigated in Einstein-scalar-Gauss-Bonnet (EsGB) theory with polynomial coupling functions $\zeta(\phi)$ satisfying $\zeta''(0) = 0$, where $\zeta'(\phi) = 0$ features besides $\phi=0$ solutions…
Allen--Cahn equation with constant and degenerate mobility, and with polynomial and logarithmic energy functionals is discretized using symmetric interior penalty discontinuous Galerkin (SIPG) finite elements in space. We show that the…
In this paper, the numerical approximation of the generalized Burgers'-Huxley equation (GBHE) with weakly singular kernels using non-conforming methods will be presented. Specifically, we discuss two new formulations. The first formulation…
A weak Galerkin (WG) method is introduced and numerically tested for the Helmholtz equation. This method is flexible by using discontinuous piecewise polynomials and retains the mass conservation property. At the same time, the WG finite…
Binary black hole spacetimes with a helical Killing vector, which are discussed as an approximation for the early stage of a binary system, are studied in a projection formalism. In this setting the four dimensional Einstein equations are…
Non-continuous "jumps" of Apparent Horizons occur generically in 3+1 (binary) black hole evolutions. The dynamical trapping horizon framework suggests a spacetime picture in which these "Apparent Horizon jumps" are understood as spatial…