Related papers: Regularization of spherically symmetric evolution …
Several interesting astrophysical phenomena are symmetric with respect to the rotation axis, like the head-on collision of compact bodies, the collapse and/or accretion of fields with a large variety of geometries, or some forms of…
In this pedagogically structured article, we describe a generalized harmonic formulation of the Einstein equations in spherical symmetry which is regular at the origin. The generalized harmonic approach has attracted significant attention…
Instabilities in finite difference codes due to the singularity of spherical coordinates at the center are studied. In typical Numerical Relativity applications, standard regularization techniques by themselves do not ensure long term…
One method for the numerical treatment of future null-infinity is to decouple coordinates from the tensor basis and choose each in a careful manner. This dual-frame approach is hampered by logarithmically divergent terms that appear in a…
Computational techniques which establish the stability of an evolution-boundary algorithm for a model wave equation with shift are incorporated into a well-posed version of the initial-boundary value problem for gravitational theory in…
Brown has recently introduced a covariant formulation of the BSSN equations which is well suited for curvilinear coordinate systems. This is particularly desirable as many astrophysical phenomena are symmetric with respect to the rotation…
We investigate how the accuracy and stability of numerical relativity simulations of 1D colliding plane waves depends on choices of equation formulations, gauge conditions, boundary conditions, and numerical methods, all in the context of a…
We review recent efforts to re-formulate the Einstein equations for fully relativistic numerical simulations. The so-called numerical relativity (computational simulations in general relativity) is a promising research field matching with…
We consider the hyperboloidal initial value problem in numerical relativity, motivated by the goal to evolve radiating compact objects such as black hole binaries with a numerical grid that includes null infinity. Unconstrained evolution…
We perform a numerical free evolution of a selfgravitating, spherically symmetric scalar field satisfying the wave equation. The evolution equations can be written in a very simple form and are symmetric hyperbolic in Eddington-Finkelstein…
We investigate the issue of regularization/renormalization in the presence of a nontrivial background in the case of 1+1-(supersymmetric) solitons. In particular we study and compare the commonly employed regularization methods (mode-…
This work presents a new evolutionary optimization algorithm in theoretical mathematics with important applications in scientific computing. The use of the evolutionary algorithm is justified by the difficulty of the study of the…
This paper presents a new Bayesian framework for quantifying discretization errors in numerical solutions of ordinary differential equations. By modelling the errors as random variables, we impose a monotonicity constraint on the variances,…
We present a new formulation of Einstein's equations for an axisymmetric spacetime with vanishing twist in vacuum. We propose a fully constrained scheme and use spherical polar coordinates. A general problem for this choice is the…
Taking the continuum limit is essential for extracting physical observables from quantum simulations of lattice gauge theories. Achieving the correct continuum limit requires careful control of all systematic uncertainties, including those…
Despite the strong focus of regularization on ill-posed problems, the general construction of such methods has not been fully explored. Moreover, many previous studies cannot be clearly adapted to handle more complex scenarios, albeit the…
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…
In the absence of symmetry assumptions most numerical relativity simulations adopt Cartesian coordinates. While Cartesian coordinates have some desirable properties, spherical polar coordinates appear better suited for certain applications,…
We present the first numerically stable nonlinear evolution for the leading-order gravitational effective field theory (Quadratic Gravity) in the spherically-symmetric sector. The formulation relies on (i) harmonic gauge to cast the…
We propose and analyze a perturbative regularization method to approximate quadratic optimization problems with finite-dimensional degeneracy. The original problem is first approximated by a regularized problem depending on a small positive…