Code generation for AMReX with applications to numerical relativity

DOI10.1088/1361-6382/AD0B37MaRDI QIDQ6183737zbMATH OpenFDO

Authors Adam J. Peterson, Philipp Mösta

Publication date 4 January 2024

Published in Classical and Quantum Gravity (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2301.08354

zbMATH Keywords

adaptive mesh refinementgravitational wavesnumerical relativityNewman-Penrose formalism

Mathematics Subject Classification ID

Gravitational waves (83C35) Computational methods for problems pertaining to relativity and gravitational theory (83-08)

Abstract: We present a new python/SymPy based code generator for producing executable numerical expressions for partial differential equations in AMReX-based applications. We demonstrate the code generator capabilities for the case of $3+1$ ADM formulations of numerical relativity for the constraint damped, conformal Z4 formulations (Z4c and CCZ4). The generated spacetime solvers are examined for stability and accuracy using a selection of checks from the standard Apples with Apples testbeds for numerical relativity applications. We also explore physically interesting vacuum spacetimes including head-on and inspiraling black hole binary collisions, and investigate the simulated gravitational waveforms from such events with the Newman-Penrose formulation of waveform extraction.

Recommendations

• GRChombo: numerical relativity with adaptive mesh refinement
• ExaGRyPE: numerical general relativity solvers based upon the hyperbolic PDEs solver engine ExaHyPE
• scientific article; zbMATH DE number 5117346
• Numerical relativity as a tool for computational astrophysics
• Kranc: a Mathematica package to generate numerical codes for tensorial evolution equations

Cites work

• A high-order finite-volume method for conservation laws on locally refined grids
• Constraint damping in the Z4 formulation and harmonic gauge
• Implementation of standard testbeds for numerical relativity
• Initial conditions for numerical relativity: introduction to numerical methods for solving elliptic PDEs
• Interaction Energy in Geometrostatics
• Kranc: a Mathematica package to generate numerical codes for tensorial evolution equations

Cited in

(2)

• GYOTO: a new general relativistic ray-tracing code
• ExaGRyPE: numerical general relativity solvers based upon the hyperbolic PDEs solver engine ExaHyPE

This page was built for publication: Code generation for AMReX with applications to numerical relativity

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6183737)