Related papers: Finite differencing second order systems describin…
We introduce a second-order time discretization method for stiff kinetic equations. The method is asymptotic-preserving (AP) -- can capture the Euler limit without numerically resolving the small Knudsen number; and positivity-preserving --…
We perform a statistical analysis of the binary black hole problem in the post-Newtonian approximation by systematically sampling and evolving the parameter space of initial configurations for quasi-circular inspirals. Through a principal…
We present a new derivation of the equations governing the oscillations of slowly rotating relativistic stars. Previous investigations have been mostly carried out in the Regge-Wheeler gauge. However, in this gauge the process of…
In this paper, we revisit the stability of power-law models, focusing on an alternative approach that differs significantly from the standard approaches used in studying power-law models. In the standard approach, stability is studied by…
Accurate calculation of the motion of a compact object in a background spacetime induced by a supermassive black hole is required for the future detection of such binary systems by the gravitational-wave detector LISA. Reaching the desired…
We study the linear perturbations about nonrotating black holes in the context of degenerate higher-order scalar-tensor (DHOST) theories, using a systematic approach that extracts the asymptotic behaviour of perturbations (at spatial…
This is the first of a series of papers describing a numerical implementation of the conformally rescaled Einstein equation, an implementation designed to calculate asymptotically flat spacetimes, especially spacetimes containing black…
In order to derive the precise gravitational waveforms for extreme mass ratio inspirals (EMRI), we develop a formulation for the second order metric perturbations produced by a point particle moving in the Schwarzschild spacetime. The…
We study second-order cosmological perturbations in scalar-tensor models of dark energy that satisfy local gravity constraints, including f(R) gravity. We derive equations for matter fluctuations under a sub-horizon approximation and…
Various higher-dimensional black holes have been shown to be unstable by studying linearized gravitational perturbations. A simpler method for demonstrating instability is to find initial data that describes a small perturbation of the…
A wide range of implicit time integration methods, including multi-step, implicit Runge-Kutta, and Galerkin finite-time element schemes, is evaluated in the context of chaotic dynamical systems. The schemes are applied to solve the Lorenz…
We give a general scheme for finite-differencing partial differential equations in flux-conservative form to second order, with a stencil that can be arbitrarily tilted with respect to the numerical grid, parameterized by a "tilt" vector…
In the paper, only Static Spherically Symmetric space-times in four dimensions are considered within modified gravity models. The non-singular static metrics, including black holes not admitting a de Sitter core in the center and…
We propose an algorithm to numerically determined whether a second-order linear PDE problem satisfying a Garding inequality is well-posed. This algorithm further provides a lower bound to the inf-sup constant of the weak formulation, which…
This article concerns second-order time discretization of subdiffusion equations with time-dependent diffusion coefficients. High-order differentiability and regularity estimates are established for subdiffusion equations with…
We develop a general framework for designing conservative numerical methods based on summation by parts operators and split forms in space, combined with relaxation Runge-Kutta methods in time. We apply this framework to create new classes…
We study the even-parity $\ell=2$ perturbations of a Schwarzschild black hole to second order. The Einstein equations can be reduced to a single linear wave equation with a potential and a source term. The source term is quadratic in terms…
In this article, we propose high-order finite-difference entropy stable schemes for the two-fluid relativistic plasma flow equations. This is achieved by exploiting the structure of the equations, which consists of three independent flux…
Using a world-sheet covariant formalism, we derive the equations of motion for second order perturbations of a generic macroscopic string, thus generalizing previous results for first order perturbations. We give the explicit results for…
Strong Stability Preserving (SSP) time integration schemes maintain stability of the forward Euler method for any initial value problem. However, only a small subset of Runge-Kutta (RK) methods are SSP, and many efficient high-order time…