1999/12/01 by Walter Landry, Saul A. Teukolsky, Landry, Walter +1
Engineering · Physics and Astronomy · #FOS: Physical sciences #Gamma-ray bursts and supernovae #General Relativity and Quantum Cosmology (gr-qc) #Pulsars and Gravitational Waves Research #Reservoir Engineering and Simulation Methods #gr-qc
paper · pdf · doi:10.48550/arxiv.gr-qc/9912004
10 pages, 24 figures, fixed some typos, added a paragraph about the diffusion method and time stepping
openalex publication_date 1999/12/01 · arxiv created 2000/05/31 · arxiv updated 2009/11/30 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
The merger of two neutron stars has been proposed as a source of gamma-ray bursts, r-process elements, and detectable gravitational waves. Extracting information from observations of these phenomena requires fully relativistic simulations. Unfortunately, the only demonstrated method for stably evolving neutron stars requires solving elliptic equations at each time step, adding substantially to the computational resources required. In this paper we present a simpler, more efficient method. The key insight is in how we apply numerical diffusion. We perform a number of tests to validate the method and our implementation. We also carry out a very rough simulation of coalescence and extraction of the gravitational waves to show that the method is viable if realistic initial data are provided.