Source term method for binary neutron stars initial data
From MaRDI portal
Publication:5163391
Iterative numerical methods for linear systems (65F10) Galactic and stellar dynamics (85A05) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Euler-Poisson-Darboux equations (35Q05) Boundary value problems on manifolds (58J32) Two-body problems (70F05) Propagation of singularities; initial value problems on manifolds (58J47)
Abstract: The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: a) the boundary is not known a-priori, but constitutes part of the solution of the problem; b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver COCAL, we demonstrate agreement to approximately . Our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.
Recommendations
- A new numerical method to construct binary neutron star initial data
- A multi-domain spectral method for initial data of arbitrary binaries in general relativity
- Initial data for high-compactness black hole-neutron star binaries
- Helically reduced wave equations and binary neutron stars
- Models of helically symmetric binary systems
Cites work
- A new numerical method to construct binary neutron star initial data
- A source term method for Poisson problems on irregular domains
- Finite difference methods for approximating Heaviside functions
- Level set methods and dynamic implicit surfaces
- Two methods for discretizing a delta function supported on a level set
Cited in
(5)- A source term method for Poisson problems with a discontinuous diffusion coefficient
- Methods for relativistic self-gravitating fluids: from binary neutron stars to black hole-disks and magnetized rotating neutron stars
- A new numerical method to construct binary neutron star initial data
- Accurate evolutions of inspiralling neutron-star binaries: assessment of the truncation error
- Realistic binary neutron star initial data with \texttt{Elliptica}
This page was built for publication: Source term method for binary neutron stars initial data
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5163391)
