King's College London

Research portal

Lessons for adaptive mesh refinement in numerical relativity

Research output: Contribution to journalReview articlepeer-review

Miren Radia, Ulrich Sperhake, Amelia Drew, Katy Clough, Pau Figueras, Eugene A. Lim, Justin L. Ripley, Josu C. Aurrekoetxea, Tiago França, Thomas Helfer

Original languageEnglish
Article number135006
JournalClassical and Quantum Gravity
Volume39
Issue number13
DOIs
Published7 Jul 2022

Bibliographical note

Funding Information: MR acknowledges support from a Science and Technology Facilities Council (STFC) studentship. KC acknowledges funding from the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation programme (Grant Agreement No. 693024), and an STFC Ernest Rutherford Fellowship. AD is supported by a Junior Research Fellowship (JRF) at Homerton College, University of Cambridge. PF is supported by the European Research Council Grant No. ERC-2014-StG 639022-NewNGR, and by a Royal Society University Research Fellowship Grant Nos. UF140319, RGF∖EA∖180260, URF∖R∖201026 and RF∖ERE∖210291. JCA acknowledges funding from the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation programme (Grant Agreement No. 693024), from the Beecroft Trust and from The Queen’s College via an extraordinary Junior Research Fellowship (eJRF). TF is supported by a Royal Society Enhancement Award (Grant No. RGF∖EA∖180260). TH is supported by NSF Grant Nos. PHY-1912550, AST-2006538, PHY-090003 and PHY-20043, and NASA Grant Nos. 17-ATP17-0225, 19-ATP19-0051 and 20-LPS20-0011. JLR and US are supported by STFC Research Grant No. ST/V005669/1; US is supported by the H2020-ERC-2014-CoG Grant ‘MaGRaTh’ No. 646597. Funding Information: The simulations presented in this paper used DiRAC resources under the projects ACSP218, ACTP186 and ACTP238. The work was performed using the Cambridge Service for Data Driven Discovery (CSD3), part of which is operated by the University of Cambridge Research Computing on behalf of the STFC DiRAC HPC Facility ( www.dirac.ac.uk ). The DiRAC component of CSD3 was funded by BEIS capital funding via STFC capital Grants ST/P002307/1 and ST/cqgac6fa9bib2452/1 and STFC operations Grant ST/cqgac6fa9bib689X/1. In addition we have used the DiRAC at Durham facility managed by the Institute for Computational Cosmology on behalf of the STFC DiRAC HPC Facility ( www.dirac.ac.uk ). The equipment was funded by BEIS capital funding via STFC capital Grants ST/P002293/1 and ST/cqgac6fa9bib2371/1, Durham University and STFC operations Grant ST/cqgac6fa9bib0832/1. DiRAC is part of the National e-Infrastructure. XSEDE resources under Grant No. PHY-090003 were used on the San Diego Supercomputing Center’s clusters Comet and Expanse and the Texas Advanced Supercomputing Center’s (TACC) Stampede2. Furthermore, we acknowledge the use of HPC resources on the TACC Frontera cluster []. PRACE resources under Grant No. 2020225359 were used on the GCS Supercomputer JUWELS at the Jülich Supercomputing Center (JCS) through the John von Neumann Institute for Computing (NIC), funded by the Gauss Center for Supercomputing e.V. ( www.gauss-centre.eu ). The Fawcett supercomputer at the Department of Applied Mathematics and Theoretical Physics (DAMTP), University of Cambridge was also used, funded by STFC Consolidated Grant ST/P000673/1. Publisher Copyright: © 2022 The Author(s). Published by IOP Publishing Ltd

King's Authors

Abstract

We demonstrate the flexibility and utility of the Berger-Rigoutsos adaptive mesh refinement (AMR) algorithm used in the open-source numerical relativity (NR) code GRChombo for generating gravitational waveforms from binary black-hole (BH) inspirals, and for studying other problems involving non-trivial matter configurations. We show that GRChombo can produce high quality binary BH waveforms through a code comparison with the established NR code Lean. We also discuss some of the technical challenges involved in making use of full AMR (as opposed to, e.g. moving box mesh refinement), including the numerical effects caused by using various refinement criteria when regridding. We suggest several ‘rules of thumb’ for when to use different tagging criteria for simulating a variety of physical phenomena. We demonstrate the use of these different criteria through example evolutions of a scalar field theory. Finally, we also review the current status and general capabilities of GRChombo.

View graph of relations

© 2020 King's College London | Strand | London WC2R 2LS | England | United Kingdom | Tel +44 (0)20 7836 5454