35,554 research outputs found
Adaptive Mesh Refinement for Characteristic Grids
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 & Lehner (J. Comp. Phys. 198 (2004), 10), recurses on
individual "diamond" characteristic grid cells. This leads to the use of
fine-grained memory management, with individual grid cells kept in
2-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 \emph{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 2nd and 4th 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.Comment: 37 pages, 15 figures (40 eps figure files, 8 of them color; all are
viewable ok in black-and-white), 1 mpeg movie, uses Springer-Verlag svjour3
document class, includes C++ source code. Changes from v1: revised in
response to referee comments: many references added, new figure added to
better explain the algorithm, other small changes, C++ code updated to latest
versio
The remapped particle-mesh advection scheme
We describe the remapped particle-mesh method, a new mass-conserving method
for solving the density equation which is suitable for combining with
semi-Lagrangian methods for compressible flow applied to numerical weather
prediction. In addition to the conservation property, the remapped
particle-mesh method is computationally efficient and at least as accurate as
current semi-Lagrangian methods based on cubic interpolation. We provide
results of tests of the method in the plane, results from incorporating the
advection method into a semi-Lagrangian method for the rotating shallow-water
equations in planar geometry, and results from extending the method to the
surface of a sphere
Generalized, energy-conserving numerical simulations of particles in general relativity. II. Test particles in electromagnetic fields and GRMHD
Direct observations of compact objects, in the form of radiation spectra,
gravitational waves from VIRGO/LIGO, and forthcoming direct imaging, are
currently one of the primary source of information on the physics of plasmas in
extreme astrophysical environments. The modeling of such physical phenomena
requires numerical methods that allow for the simulation of microscopic plasma
dynamics in presence of both strong gravity and electromagnetic fields. In
Bacchini et al. (2018) we presented a detailed study on numerical techniques
for the integration of free geodesic motion. Here we extend the study by
introducing electromagnetic forces in the simulation of charged particles in
curved spacetimes. We extend the Hamiltonian energy-conserving method presented
in Bacchini et al. (2018) to include the Lorentz force and we test its
performance compared to that of standard explicit Runge-Kutta and implicit
midpoint rule schemes against analytic solutions. Then, we show the application
of the numerical schemes to the integration of test particle trajectories in
general relativistic magnetohydrodynamic (GRMHD) simulations, by modifying the
algorithms to handle grid-based electromagnetic fields. We test this approach
by simulating ensembles of charged particles in a static GRMHD configuration
obtained with the Black Hole Accretion Code (BHAC)
- …