Adaptive mesh refinement for characteristic grids

Abstract

I consider techniques for Berger-Oliger adaptive mesh refinement (AMR) when numerically solving partial differential equations with wave-like solutions, using characteristic (double-null) grids. Such AMR algorithms are naturally recursive, and the best-known past Berger-Oliger characteristic AMR algorithm, that of Pretorius and Lehner (J Comp Phys 198:10, 2004), recurses on individual "diamond" characteristic grid cells. This leads to the use of fine-grained memory management, with individual grid cells kept in two-dimensional linked lists at each refinement level. This complicates the implementation and adds overhead in both space and time. Here I describe a Berger-Oliger characteristic AMR algorithm which instead recurses on null slices. This algorithm is very similar to the usual Cauchy Berger-Oliger algorithm, and uses relatively coarse-grained memory management, allowing entire null slices to be stored in contiguous arrays in memory. The algorithm is very efficient in both space and time. I describe discretizations yielding both second and fourth order global accuracy. My code implementing the algorithm described here is included in the electronic supplementary materials accompanying this paper, and is freely available to other researchers under the terms of the GNU general public license. This paper is dedicated to the memory of Thomas Radke, my late friend and colleague in many computational adventures. Electronic supplementary materialThe online version of this article (

Hidden momentum and black hole kicks

“…where the indices N, S, W, E refer to north, south, west, east on the grid as in reference [20,21], h is the grid spacing, and V c is approximated by V c ≈ (1/2)(V E + V W ) and V = (1 − 2M/r) ( + 1)/(4r 2 ). We evolve equation (23) with M = 1 on a grid of u ∈ [u 0 , u max ] and v ∈ [0, v max ] with initial data 1/r for = 1, 2, 3.…”

Black hole factory: A review of double-null formalism

“…A local error indicator related to the variables and changing its value significantly in the adequate region is the function ∆r/r with the difference in r, ∆r, calculated along the u-coordinate. 9,10 In general, one can also apply the adaptive mesh for both of u and v directions. 11 In addition, it is also equivalently useful to adaptively choose the gauge degrees of freedom for several circumstances.…”

Time-parallel gravitational collapse simulation

“…In particular, adaptive mesh refinement [5] is used to resolve the black hole formation physics. In [22] the same investigation is carried out in double null or characteristic coordinates (τ, ρ) without mesh refinement (see, however, [41,45]). Finally, in [31] the effect of quantum gravity modifications on the collapse is studied in adjusted characteristic coordinates.…”