8,342 research outputs found
SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code
Numerical heating in particle-in-cell (PIC) codes currently precludes the
accurate simulation of cold, relativistic plasma over long periods, severely
limiting their applications in astrophysical environments. We present a
spatially higher-order accurate relativistic PIC algorithm in one spatial
dimension, which conserves charge and momentum exactly. We utilize the
smoothness implied by the usage of higher-order interpolation functions to
achieve a spatially higher-order accurate algorithm (up to fifth order). We
validate our algorithm against several test problems -- thermal stability of
stationary plasma, stability of linear plasma waves, and two-stream instability
in the relativistic and non-relativistic regimes. Comparing our simulations to
exact solutions of the dispersion relations, we demonstrate that SHARP can
quantitatively reproduce important kinetic features of the linear regime. Our
simulations have a superior ability to control energy non-conservation and
avoid numerical heating in comparison to common second-order schemes. We
provide a natural definition for convergence of a general PIC algorithm: the
complement of physical modes captured by the simulation, i.e., those that lie
above the Poisson noise, must grow commensurately with the resolution. This
implies that it is necessary to simultaneously increase the number of particles
per cell and decrease the cell size. We demonstrate that traditional ways for
testing for convergence fail, leading to plateauing of the energy error. This
new PIC code enables us to faithfully study the long-term evolution of plasma
problems that require absolute control of the energy and momentum conservation.Comment: 26 pages, 19 figures, discussion about performance is added,
published in Ap
A new hybrid implicit-explicit FDTD method for local subgridding in multiscale 2-D TE scattering problems
The conventional finite-difference time-domain (FDTD) method with staggered Yee scheme does not easily allow including thin material layers, especially so if these layers are highly conductive. This paper proposes a novel subgridding technique for 2-D problems, based on a hybrid implicit-explicit scheme, which efficiently copes with this problem. In the subgrid, the new method collocates field components such that the thin layer boundaries are defined unambiguously. Moreover, aspect ratios of more than a million do not impair the stability of the method and allow for very accurate predictions of the skin effect. The new method retains the Courant limit of the coarse Yee grid and is easily incorporated into existing FDTD codes. A number of illustrative examples, including scattering by a metal grating, demonstrate the accuracy and stability of the new method
Introducing PHAEDRA: a new spectral code for simulations of relativistic magnetospheres
We describe a new scheme for evolving the equations of force-free
electrodynamics, the vanishing-inertia limit of magnetohydrodynamics. This
pseudospectral code uses global orthogonal basis function expansions to take
accurate spatial derivatives, allowing the use of an unstaggered mesh and the
complete force-free current density. The method has low numerical dissipation
and diffusion outside of singular current sheets. We present a range of one-
and two-dimensional tests, and demonstrate convergence to both smooth and
discontinuous analytic solutions. As a first application, we revisit the
aligned rotator problem, obtaining a steady solution with resistivity localised
in the equatorial current sheet outside the light cylinder.Comment: 23 pages, 18 figures, accepted for publication in MNRA
Discontinuous Galerkin method for the spherically reduced BSSN system with second-order operators
We present a high-order accurate discontinuous Galerkin method for evolving
the spherically-reduced Baumgarte-Shapiro-Shibata-Nakamura (BSSN) system
expressed in terms of second-order spatial operators. Our multi-domain method
achieves global spectral accuracy and long-time stability on short
computational domains. We discuss in detail both our scheme for the BSSN system
and its implementation. After a theoretical and computational verification of
the proposed scheme, we conclude with a brief discussion of issues likely to
arise when one considers the full BSSN system.Comment: 35 pages, 6 figures, 1 table, uses revtex4. Revised in response to
referee's repor
Numerical Relativity Using a Generalized Harmonic Decomposition
A new numerical scheme to solve the Einstein field equations based upon the
generalized harmonic decomposition of the Ricci tensor is introduced. The
source functions driving the wave equations that define generalized harmonic
coordinates are treated as independent functions, and encode the coordinate
freedom of solutions. Techniques are discussed to impose particular gauge
conditions through a specification of the source functions. A 3D, free
evolution, finite difference code implementing this system of equations with a
scalar field matter source is described. The second-order-in-space-and-time
partial differential equations are discretized directly without the use first
order auxiliary terms, limiting the number of independent functions to
fifteen--ten metric quantities, four source functions and the scalar field.
This also limits the number of constraint equations, which can only be enforced
to within truncation error in a numerical free evolution, to four. The
coordinate system is compactified to spatial infinity in order to impose
physically motivated, constraint-preserving outer boundary conditions. A
variant of the Cartoon method for efficiently simulating axisymmetric
spacetimes with a Cartesian code is described that does not use interpolation,
and is easier to incorporate into existing adaptive mesh refinement packages.
Preliminary test simulations of vacuum black hole evolution and black hole
formation via scalar field collapse are described, suggesting that this method
may be useful for studying many spacetimes of interest.Comment: 18 pages, 6 figures; updated to coincide with journal version, which
includes some expanded discussions and a new appendix with a stability
analysis of a simplified problem using the same discretization scheme
described in the pape
Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes
We have developed the formalism necessary to employ the
discontinuous-Galerkin approach in general-relativistic hydrodynamics. The
formalism is firstly presented in a general 4-dimensional setting and then
specialized to the case of spherical symmetry within a 3+1 splitting of
spacetime. As a direct application, we have constructed a one-dimensional code,
EDGES, which has been used to asses the viability of these methods via a series
of tests involving highly relativistic flows in strong gravity. Our results
show that discontinuous Galerkin methods are able not only to handle strong
relativistic shock waves but, at the same time, to attain very high orders of
accuracy and exponential convergence rates in smooth regions of the flow. Given
these promising prospects and their affinity with a pseudospectral solution of
the Einstein equations, discontinuous Galerkin methods could represent a new
paradigm for the accurate numerical modelling in relativistic astrophysics.Comment: 24 pages, 19 figures. Small changes; matches version to appear in PR
Formulation of boundary conditions for the multigrid acceleration of the Euler and Navier Stokes equations
An explicit, Multigrid algorithm was written to solve the Euler and Navier-Stokes equations with special consideration given to the coarse mesh boundary conditions. These are formulated in a manner consistent with the interior solution, utilizing forcing terms to prevent coarse-mesh truncation error from affecting the fine-mesh solution. A 4-Stage Hybrid Runge-Kutta Scheme is used to advance the solution in time, and Multigrid convergence is further enhanced by using local time-stepping and implicit residual smoothing. Details of the algorithm are presented along with a description of Jameson's standard Multigrid method and a new approach to formulating the Multigrid equations
The Emergence of Gravitational Wave Science: 100 Years of Development of Mathematical Theory, Detectors, Numerical Algorithms, and Data Analysis Tools
On September 14, 2015, the newly upgraded Laser Interferometer
Gravitational-wave Observatory (LIGO) recorded a loud gravitational-wave (GW)
signal, emitted a billion light-years away by a coalescing binary of two
stellar-mass black holes. The detection was announced in February 2016, in time
for the hundredth anniversary of Einstein's prediction of GWs within the theory
of general relativity (GR). The signal represents the first direct detection of
GWs, the first observation of a black-hole binary, and the first test of GR in
its strong-field, high-velocity, nonlinear regime. In the remainder of its
first observing run, LIGO observed two more signals from black-hole binaries,
one moderately loud, another at the boundary of statistical significance. The
detections mark the end of a decades-long quest, and the beginning of GW
astronomy: finally, we are able to probe the unseen, electromagnetically dark
Universe by listening to it. In this article, we present a short historical
overview of GW science: this young discipline combines GR, arguably the
crowning achievement of classical physics, with record-setting, ultra-low-noise
laser interferometry, and with some of the most powerful developments in the
theory of differential geometry, partial differential equations,
high-performance computation, numerical analysis, signal processing,
statistical inference, and data science. Our emphasis is on the synergy between
these disciplines, and how mathematics, broadly understood, has historically
played, and continues to play, a crucial role in the development of GW science.
We focus on black holes, which are very pure mathematical solutions of
Einstein's gravitational-field equations that are nevertheless realized in
Nature, and that provided the first observed signals.Comment: 41 pages, 5 figures. To appear in Bulletin of the American
Mathematical Societ
Stability of exact force-free electrodynamic solutions and scattering from spacetime curvature
Recently, a family of exact force-free electrodynamic (FFE) solutions was
given by Brennan, Gralla and Jacobson, which generalizes earlier solutions by
Michel, Menon and Dermer, and other authors. These solutions have been proposed
as useful models for describing the outer magnetosphere of conducting stars. As
with any exact analytical solution that aspires to describe actual physical
systems, it is vitally important that the solution possess the necessary
stability. In this paper, we show via fully nonlinear numerical simulations
that the aforementioned FFE solutions, despite being highly special in their
properties, are nonetheless stable under small perturbations. Through this
study, we also introduce a three-dimensional pseudospectral relativistic FFE
code that achieves exponential convergence for smooth test cases, as well as
two additional well-posed FFE evolution systems in the appendix that have
desirable mathematical properties. Furthermore, we provide an explicit analysis
that demonstrates how propagation along degenerate principal null directions of
the spacetime curvature tensor simplifies scattering, thereby providing an
intuitive understanding of why these exact solutions are tractable, i.e. why
they are not backscattered by spacetime curvature.Comment: 33 pages, 21 figures; V2 updated to match published versio
- …