Three-dimensional general-relativistic hydrodynamic simulations of binary neutron star coalescence and stellar collapse with multipatch grids

C. Reisswig, R. Haas, C. D. Ott, E. Abdikamalov, P. Mösta, D. Pollney, E. Schnetter

Research output: Contribution to journalArticle

29 Citations (Scopus)

Abstract

We present a new three-dimensional, general-relativistic hydrodynamic evolution scheme coupled to dynamical spacetime evolutions which is capable of efficiently simulating stellar collapse, isolated neutron stars, black hole formation, and binary neutron star coalescence. We make use of a set of adapted curvilinear grids (multipatches) coupled with flux-conservative, cell-centered adaptive mesh refinement. This allows us to significantly enlarge our computational domains while still maintaining high resolution in the gravitational wave extraction zone, the exterior layers of a star, or the region of mass ejection in merging neutron stars. The fluid is evolved with a high-resolution, shock-capturing finite volume scheme, while the spacetime geometry is evolved using fourth-order finite differences. We employ a multirate Runge-Kutta time-integration scheme for efficiency, evolving the fluid with second-order integration and the spacetime geometry with fourth-order integration. We validate our code by a number of benchmark problems: a rotating stellar collapse model, an excited neutron star, neutron star collapse to a black hole, and binary neutron star coalescence. The test problems, especially the latter, greatly benefit from higher resolution in the gravitational wave extraction zone, causally disconnected outer boundaries, and application of Cauchy-characteristic gravitational wave extraction. We show that we are able to extract convergent gravitational wave modes up to (ℓ,m)=(6,6). This study paves the way for more realistic and detailed studies of compact objects and stellar collapse in full three dimensions and in large computational domains. The multipatch infrastructure and the improvements to mesh refinement and hydrodynamics codes discussed in this paper will be made available as part of the open-source Einstein Toolkit.

Original languageEnglish
Article number064023
JournalPhysical Review D - Particles, Fields, Gravitation and Cosmology
Volume87
Issue number6
DOIs
Publication statusPublished - Mar 18 2013
Externally publishedYes

Fingerprint

binary stars
coalescing
neutron stars
hydrodynamics
grids
gravitational waves
simulation
high resolution
fluids
geometry
ejection
shock
stars
cells

ASJC Scopus subject areas

  • Nuclear and High Energy Physics

Cite this

Three-dimensional general-relativistic hydrodynamic simulations of binary neutron star coalescence and stellar collapse with multipatch grids. / Reisswig, C.; Haas, R.; Ott, C. D.; Abdikamalov, E.; Mösta, P.; Pollney, D.; Schnetter, E.

In: Physical Review D - Particles, Fields, Gravitation and Cosmology, Vol. 87, No. 6, 064023, 18.03.2013.

Research output: Contribution to journalArticle

@article{19a3567c40dc419d948d949aa06cdd68,
title = "Three-dimensional general-relativistic hydrodynamic simulations of binary neutron star coalescence and stellar collapse with multipatch grids",
abstract = "We present a new three-dimensional, general-relativistic hydrodynamic evolution scheme coupled to dynamical spacetime evolutions which is capable of efficiently simulating stellar collapse, isolated neutron stars, black hole formation, and binary neutron star coalescence. We make use of a set of adapted curvilinear grids (multipatches) coupled with flux-conservative, cell-centered adaptive mesh refinement. This allows us to significantly enlarge our computational domains while still maintaining high resolution in the gravitational wave extraction zone, the exterior layers of a star, or the region of mass ejection in merging neutron stars. The fluid is evolved with a high-resolution, shock-capturing finite volume scheme, while the spacetime geometry is evolved using fourth-order finite differences. We employ a multirate Runge-Kutta time-integration scheme for efficiency, evolving the fluid with second-order integration and the spacetime geometry with fourth-order integration. We validate our code by a number of benchmark problems: a rotating stellar collapse model, an excited neutron star, neutron star collapse to a black hole, and binary neutron star coalescence. The test problems, especially the latter, greatly benefit from higher resolution in the gravitational wave extraction zone, causally disconnected outer boundaries, and application of Cauchy-characteristic gravitational wave extraction. We show that we are able to extract convergent gravitational wave modes up to (ℓ,m)=(6,6). This study paves the way for more realistic and detailed studies of compact objects and stellar collapse in full three dimensions and in large computational domains. The multipatch infrastructure and the improvements to mesh refinement and hydrodynamics codes discussed in this paper will be made available as part of the open-source Einstein Toolkit.",
author = "C. Reisswig and R. Haas and Ott, {C. D.} and E. Abdikamalov and P. M{\"o}sta and D. Pollney and E. Schnetter",
year = "2013",
month = "3",
day = "18",
doi = "10.1103/PhysRevD.87.064023",
language = "English",
volume = "87",
journal = "Physical review D: Particles and fields",
issn = "1550-7998",
publisher = "American Physical Society",
number = "6",

}

TY - JOUR

T1 - Three-dimensional general-relativistic hydrodynamic simulations of binary neutron star coalescence and stellar collapse with multipatch grids

AU - Reisswig, C.

AU - Haas, R.

AU - Ott, C. D.

AU - Abdikamalov, E.

AU - Mösta, P.

AU - Pollney, D.

AU - Schnetter, E.

PY - 2013/3/18

Y1 - 2013/3/18

N2 - We present a new three-dimensional, general-relativistic hydrodynamic evolution scheme coupled to dynamical spacetime evolutions which is capable of efficiently simulating stellar collapse, isolated neutron stars, black hole formation, and binary neutron star coalescence. We make use of a set of adapted curvilinear grids (multipatches) coupled with flux-conservative, cell-centered adaptive mesh refinement. This allows us to significantly enlarge our computational domains while still maintaining high resolution in the gravitational wave extraction zone, the exterior layers of a star, or the region of mass ejection in merging neutron stars. The fluid is evolved with a high-resolution, shock-capturing finite volume scheme, while the spacetime geometry is evolved using fourth-order finite differences. We employ a multirate Runge-Kutta time-integration scheme for efficiency, evolving the fluid with second-order integration and the spacetime geometry with fourth-order integration. We validate our code by a number of benchmark problems: a rotating stellar collapse model, an excited neutron star, neutron star collapse to a black hole, and binary neutron star coalescence. The test problems, especially the latter, greatly benefit from higher resolution in the gravitational wave extraction zone, causally disconnected outer boundaries, and application of Cauchy-characteristic gravitational wave extraction. We show that we are able to extract convergent gravitational wave modes up to (ℓ,m)=(6,6). This study paves the way for more realistic and detailed studies of compact objects and stellar collapse in full three dimensions and in large computational domains. The multipatch infrastructure and the improvements to mesh refinement and hydrodynamics codes discussed in this paper will be made available as part of the open-source Einstein Toolkit.

AB - We present a new three-dimensional, general-relativistic hydrodynamic evolution scheme coupled to dynamical spacetime evolutions which is capable of efficiently simulating stellar collapse, isolated neutron stars, black hole formation, and binary neutron star coalescence. We make use of a set of adapted curvilinear grids (multipatches) coupled with flux-conservative, cell-centered adaptive mesh refinement. This allows us to significantly enlarge our computational domains while still maintaining high resolution in the gravitational wave extraction zone, the exterior layers of a star, or the region of mass ejection in merging neutron stars. The fluid is evolved with a high-resolution, shock-capturing finite volume scheme, while the spacetime geometry is evolved using fourth-order finite differences. We employ a multirate Runge-Kutta time-integration scheme for efficiency, evolving the fluid with second-order integration and the spacetime geometry with fourth-order integration. We validate our code by a number of benchmark problems: a rotating stellar collapse model, an excited neutron star, neutron star collapse to a black hole, and binary neutron star coalescence. The test problems, especially the latter, greatly benefit from higher resolution in the gravitational wave extraction zone, causally disconnected outer boundaries, and application of Cauchy-characteristic gravitational wave extraction. We show that we are able to extract convergent gravitational wave modes up to (ℓ,m)=(6,6). This study paves the way for more realistic and detailed studies of compact objects and stellar collapse in full three dimensions and in large computational domains. The multipatch infrastructure and the improvements to mesh refinement and hydrodynamics codes discussed in this paper will be made available as part of the open-source Einstein Toolkit.

UR - http://www.scopus.com/inward/record.url?scp=84875478943&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84875478943&partnerID=8YFLogxK

U2 - 10.1103/PhysRevD.87.064023

DO - 10.1103/PhysRevD.87.064023

M3 - Article

VL - 87

JO - Physical review D: Particles and fields

JF - Physical review D: Particles and fields

SN - 1550-7998

IS - 6

M1 - 064023

ER -