GRay2

A General Purpose Geodesic Integrator for Kerr Spacetimes

Chi Kwan Chan, Lia Medeiros, Feryal Ozel, Dimitrios Psaltis

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

Fast and accurate integration of geodesics in Kerr spacetimes is an important tool in modeling the orbits of stars and the transport of radiation in the vicinities of black holes. Most existing integration algorithms employ Boyer-Lindquist (BL) coordinates, which have coordinate singularities at the event horizon and along the poles. Handling the singularities requires special numerical treatment in these regions, often slows down the calculations, and may lead to inaccurate geodesics. We present here a new general-purpose geodesic integrator, GRay2, that overcomes these issues by employing the Cartesian form of Kerr-Schild (KS) coordinates. By performing particular mathematical manipulations of the geodesic equations and several optimizations, we develop an implementation of the Cartesian KS coordinates that outperforms calculations that use the seemingly simpler equations in BL coordinates. We also employ the OpenCL framework, which allows GRay2 to run on multicore CPUs as well as on a wide range of graphics processing units hardware accelerators, making the algorithm more versatile. We report numerous convergence tests and benchmark results for GRay2 for both time-like (particle) and null (photon) geodesics.

Original languageEnglish (US)
Article number59
JournalAstrophysical Journal
Volume867
Issue number1
DOIs
StatePublished - Nov 1 2018

Fingerprint

integrators
event horizon
hardware
manipulators
accelerators
poles
orbits
stars
optimization
photons
radiation
modeling
calculation

Keywords

  • black hole physics
  • gravitation
  • methods: numerical

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

GRay2 : A General Purpose Geodesic Integrator for Kerr Spacetimes. / Chan, Chi Kwan; Medeiros, Lia; Ozel, Feryal; Psaltis, Dimitrios.

In: Astrophysical Journal, Vol. 867, No. 1, 59, 01.11.2018.

Research output: Contribution to journalArticle

@article{ea01b905193f453c8a6171532af4d87e,
title = "GRay2: A General Purpose Geodesic Integrator for Kerr Spacetimes",
abstract = "Fast and accurate integration of geodesics in Kerr spacetimes is an important tool in modeling the orbits of stars and the transport of radiation in the vicinities of black holes. Most existing integration algorithms employ Boyer-Lindquist (BL) coordinates, which have coordinate singularities at the event horizon and along the poles. Handling the singularities requires special numerical treatment in these regions, often slows down the calculations, and may lead to inaccurate geodesics. We present here a new general-purpose geodesic integrator, GRay2, that overcomes these issues by employing the Cartesian form of Kerr-Schild (KS) coordinates. By performing particular mathematical manipulations of the geodesic equations and several optimizations, we develop an implementation of the Cartesian KS coordinates that outperforms calculations that use the seemingly simpler equations in BL coordinates. We also employ the OpenCL framework, which allows GRay2 to run on multicore CPUs as well as on a wide range of graphics processing units hardware accelerators, making the algorithm more versatile. We report numerous convergence tests and benchmark results for GRay2 for both time-like (particle) and null (photon) geodesics.",
keywords = "black hole physics, gravitation, methods: numerical",
author = "Chan, {Chi Kwan} and Lia Medeiros and Feryal Ozel and Dimitrios Psaltis",
year = "2018",
month = "11",
day = "1",
doi = "10.3847/1538-4357/aadfe5",
language = "English (US)",
volume = "867",
journal = "Astrophysical Journal",
issn = "0004-637X",
publisher = "IOP Publishing Ltd.",
number = "1",

}

TY - JOUR

T1 - GRay2

T2 - A General Purpose Geodesic Integrator for Kerr Spacetimes

AU - Chan, Chi Kwan

AU - Medeiros, Lia

AU - Ozel, Feryal

AU - Psaltis, Dimitrios

PY - 2018/11/1

Y1 - 2018/11/1

N2 - Fast and accurate integration of geodesics in Kerr spacetimes is an important tool in modeling the orbits of stars and the transport of radiation in the vicinities of black holes. Most existing integration algorithms employ Boyer-Lindquist (BL) coordinates, which have coordinate singularities at the event horizon and along the poles. Handling the singularities requires special numerical treatment in these regions, often slows down the calculations, and may lead to inaccurate geodesics. We present here a new general-purpose geodesic integrator, GRay2, that overcomes these issues by employing the Cartesian form of Kerr-Schild (KS) coordinates. By performing particular mathematical manipulations of the geodesic equations and several optimizations, we develop an implementation of the Cartesian KS coordinates that outperforms calculations that use the seemingly simpler equations in BL coordinates. We also employ the OpenCL framework, which allows GRay2 to run on multicore CPUs as well as on a wide range of graphics processing units hardware accelerators, making the algorithm more versatile. We report numerous convergence tests and benchmark results for GRay2 for both time-like (particle) and null (photon) geodesics.

AB - Fast and accurate integration of geodesics in Kerr spacetimes is an important tool in modeling the orbits of stars and the transport of radiation in the vicinities of black holes. Most existing integration algorithms employ Boyer-Lindquist (BL) coordinates, which have coordinate singularities at the event horizon and along the poles. Handling the singularities requires special numerical treatment in these regions, often slows down the calculations, and may lead to inaccurate geodesics. We present here a new general-purpose geodesic integrator, GRay2, that overcomes these issues by employing the Cartesian form of Kerr-Schild (KS) coordinates. By performing particular mathematical manipulations of the geodesic equations and several optimizations, we develop an implementation of the Cartesian KS coordinates that outperforms calculations that use the seemingly simpler equations in BL coordinates. We also employ the OpenCL framework, which allows GRay2 to run on multicore CPUs as well as on a wide range of graphics processing units hardware accelerators, making the algorithm more versatile. We report numerous convergence tests and benchmark results for GRay2 for both time-like (particle) and null (photon) geodesics.

KW - black hole physics

KW - gravitation

KW - methods: numerical

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

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

U2 - 10.3847/1538-4357/aadfe5

DO - 10.3847/1538-4357/aadfe5

M3 - Article

VL - 867

JO - Astrophysical Journal

JF - Astrophysical Journal

SN - 0004-637X

IS - 1

M1 - 59

ER -