Related papers: Binary black hole simulation with an adaptive fini…
This paper is concerned with a numerical solution to the scattering of a time-harmonic electromagnetic wave by a bounded and impenetrable obstacle in three dimensions. The electromagnetic wave propagation is modeled by a boundary value…
An approximate analytical and non-linear solution of the Einstein field equations is derived for a system of multiple non-rotating black holes. The associated space-time has the same asymptotic structure as the Brill-Lindquist initial data…
A finite element method using B-splines is presented and compared with a conventional finite element method of Lagrangian type. The efficiency of both methods has been investigated at the example of a coupled non-linear system of Dirac…
Keeping Einstein's equations in second order form can be appealing for computational efficiency, because of the reduced number of variables and constraints. Stability issues emerge, however, which are not present in first order…
Approximate solutions to the Einstein field equations are a valuable tool to investigate gravitational phenomena. An important aspect of any approximation is to investigate and quantify its regime of validity. We present a study that…
The description of extreme-mass-ratio binary systems in the inspiral phase is a challenging problem in gravitational wave physics with significant relevance for the space interferometer LISA. The main difficulty lies in the evaluation of…
This paper introduces a novel a posteriori error estimation framework for the enriched Galerkin (EG) finite element method applied to linear parabolic equations. While the EG method has been recognized for its local conservation property…
We present and analyse a new conforming space-time Galerkin discretisation of a semi-linear wave equation, based on a variational formulation derived from De Giorgi's elliptic regularisation viewpoint of the wave equation in second-order…
We developed a numerical code which evolves the semiclassical Einstein's equation (with the quantum stress-energy contribution added as a source term) for the spherically symmetric metric inside an evaporating semiclassical charged black…
The weak Galerkin (WG) finite element method is an effective and flexible general numerical technique for solving partial differential equations. It is a natural extension of the classic conforming finite element method for discontinuous…
This paper presents a primal-dual weak Galerkin (PD-WG) finite element method for a class of second order elliptic equations of Fokker-Planck type. The method is based on a variational form where all the derivatives are applied to the test…
In this work we develop and analyze an adaptive finite element method for efficiently solving electrical impedance tomography -- a severely ill-posed nonlinear inverse problem for recovering the conductivity from boundary voltage…
We consider a linear elliptic partial differential equation (PDE) with a generic uniformly bounded parametric coefficient. The solution to this PDE problem is approximated in the framework of stochastic Galerkin finite element methods. We…
Einstein's equations are known to lead to the formation of black holes and spacetime singularities. This appears to be a manifestation of the mathematical phenomenon of finite-time blowup: a formation of singularities from regular initial…
We derive the BPS type of first order differential equations for the rotating black hole solutions in the three-dimensional Einstein gravity coupled minimally with a self-interacting scalar field, using fake supersymmetry formalism. It…
Fourier solvers have become efficient tools to establish structure-property relations in heterogeneous materials. Introduced as an alternative to the Finite Element (FE) method, they are based on fixed-point solutions of the…
Novel applications of Numerical Relativity demand for more flexible algorithms and tools. In this paper, I develop and test a multigrid solver, based on the infrastructure provided by the Einstein Toolkit, for elliptic partial differential…
A generalized Black-Scholes equation is considered on the semi-axis. It is transformed on the interval (0,1) in order to make the computational domain finite. The new parabolic operator degenerates at the both ends of the interval and we…
Consider the scattering of an elastic plane wave by a rigid obstacle, which is immersed in a homogeneous and isotropic elastic medium in two dimensions. Based on a Dirichlet-to-Neumann (DtN) operator, an exact transparent boundary condition…
In this article, interior penalty discontinuous Galerkin methods using immersed finite element functions are employed to solve parabolic interface problems. Typical semi-discrete and fully discrete schemes are presented and analyzed.…