Related papers: Binary black hole simulation with an adaptive fini…
In this work we address the analysis of the stationary generalized Burgers-Huxley equation (a nonlinear elliptic problem with anomalous advection) and propose conforming, nonconforming and discontinuous Galerkin finite element methods for…
In order to perform model-dependent tests of general relativity with gravitational wave observations, we must have access to numerical relativity binary black hole waveforms in theories beyond general relativity (GR). In this study, we…
We present a new scheme for constructing initial data for the Einstein field equations using the conformal thin-sandwich formulation that does not assume conformal flatness or approximate Killing vectors. This includes a method for…
The Einstein field equations of gravitation are characterized by cross-scale, high-order nonlinear terms, representing a challenge for numerical modeling. In an exact spectral decomposition, high-order nonlinearities correspond to a…
In this paper, we aim to develop a hybridizable discontinuous Galerkin (HDG) method for the indefinite time-harmonic Maxwell equations with the perfectly conducting boundary in the three-dimensional space. First, we derive the wavenumber…
We consider the modified Einstein equations obtained in the framework of effective spherically symmetric polymer models inspired by Loop Quantum Gravity. When one takes into account the anomaly free point-wise holonomy quantum corrections,…
We investigate scalarized black holes in the Einstein-minimally coupled scalar theory with a negative potential $V(\phi)=-\alpha^2\phi^6$. The tachyonic instability is absent from analyzing the linearized scalar equation, which could not…
In this paper, we define new unfitted finite element methods for numerically approximating the solution of surface partial differential equations using bulk finite elements. The key idea is that the $n$-dimensional hypersurface, $\Gamma…
We study the dynamics of black holes in Einstein-scalar-Gauss-Bonnet theories that exhibit spontaneous black hole scalarization using recently introduced methods for solving the full, non-perturbative equations of motion. For one sign of…
We present a new numerical code developed for the evolution of binary black-hole spacetimes using different initial data and evolution techniques. The code is demonstrated to produce state-of-the-art simulations of orbiting and inspiralling…
We present an analytic, perturbative solution to the Einstein equations with a scalar field that describes dynamical black holes in a slow-roll inflationary cosmology. We show that the metric evolves quasi-statically through a sequence of…
Recent numerical simulations have found that the Cauchy horizon inside spherical charged black holes, when perturbed nonlinearly by a self-gravitating, minimally-coupled, massless, spherically-symmetric scalar field, turns into a null weak…
Spherically symmetric (1D) black-hole spacetimes are considered as a test for numerical relativity. A finite difference code, based in the hyperbolic structure of Einstein's equations with the harmonic slicing condition is presented.…
We calculate the gravitational perturbations produced by a small mass in eccentric orbit about a much more massive Schwarzschild black hole and use the numerically computed perturbations to solve for the metric. The calculations are…
The study of the interaction of a massive black hole binary with its gaseous environment is crucial in order to be able to predict merger rates and possible electromagnetic counterparts of gravitational wave signals. The evolution of the…
We develop a new $C^{\,0}$-continuous Petrov-Galerkin spectral element method for one-dimensional fractional elliptic problems of the form ${}_{0}{\mathcal{D}}_{x}^{\alpha} u(x) - \lambda u(x) = f(x)$, $\alpha \in (1,2]$, subject to…
We study here the evolution of a massless scalar field in a spacetime, developing from a regular initial spacelike surface. The Einstein equations and regularity and boundary conditions governing the same are specified. Both homogeneous and…
A new weak Galerkin (WG) method is introduced and analyzed for the second order elliptic equation formulated as a system of two first order linear equations. This method, called WG-MFEM, is designed by using discontinuous piecewise…
We describe the first axisymmetric numerical code based on the generalized harmonic formulation of the Einstein equations which is regular at the axis. We test the code by investigating gravitational collapse of distributions of complex…
This work extends the high-resolution isogeometric analysis approach established for scalar transport equations to the equations of gas dynamics. The group finite element formulation is adopted to obtain an efficient assembly procedure for…