1,984 research outputs found
Advantages of modified ADM formulation: constraint propagation analysis of Baumgarte-Shapiro-Shibata-Nakamura system
Several numerical relativity groups are using a modified ADM formulation for
their simulations, which was developed by Nakamura et al (and widely cited as
Baumgarte-Shapiro-Shibata-Nakamura system). This so-called BSSN formulation is
shown to be more stable than the standard ADM formulation in many cases, and
there have been many attempts to explain why this re-formulation has such an
advantage. We try to explain the background mechanism of the BSSN equations by
using eigenvalue analysis of constraint propagation equations. This analysis
has been applied and has succeeded in explaining other systems in our series of
works. We derive the full set of the constraint propagation equations, and
study it in the flat background space-time. We carefully examine how the
replacements and adjustments in the equations change the propagation structure
of the constraints, i.e. whether violation of constraints (if it exists) will
decay or propagate away. We conclude that the better stability of the BSSN
system is obtained by their adjustments in the equations, and that the
combination of the adjustments is in a good balance, i.e. a lack of their
adjustments might fail to obtain the present stability. We further propose
other adjustments to the equations, which may offer more stable features than
the current BSSN equations.Comment: 10 pages, RevTeX4, added related discussion to gr-qc/0209106, the
version to appear in Phys. Rev.
Numerical experiments of adjusted BSSN systems for controlling constraint violations
We present our numerical comparisons between the BSSN formulation widely used
in numerical relativity today and its adjusted versions using constraints. We
performed three testbeds: gauge-wave, linear wave, and Gowdy-wave tests,
proposed by the Mexico workshop on the formulation problem of the Einstein
equations. We tried three kinds of adjustments, which were previously proposed
from the analysis of the constraint propagation equations, and investigated how
they improve the accuracy and stability of evolutions. We observed that the
signature of the proposed Lagrange multipliers are always right and the
adjustments improve the convergence and stability of the simulations. When the
original BSSN system already shows satisfactory good evolutions (e.g., linear
wave test), the adjusted versions also coincide with those evolutions; while in
some cases (e.g., gauge-wave or Gowdy-wave tests) the simulations using the
adjusted systems last 10 times as long as those using the original BSSN
equations. Our demonstrations imply a potential to construct a robust evolution
system against constraint violations even in highly dynamical situations.Comment: to be published in PR
Constraint propagation in the family of ADM systems
The current important issue in numerical relativity is to determine which
formulation of the Einstein equations provides us with stable and accurate
simulations. Based on our previous work on "asymptotically constrained"
systems, we here present constraint propagation equations and their eigenvalues
for the Arnowitt-Deser-Misner (ADM) evolution equations with additional
constraint terms (adjusted terms) on the right hand side. We conjecture that
the system is robust against violation of constraints if the amplification
factors (eigenvalues of Fourier-component of the constraint propagation
equations) are negative or pure-imaginary. We show such a system can be
obtained by choosing multipliers of adjusted terms. Our discussion covers
Detweiler's proposal (1987) and Frittelli's analysis (1997), and we also
mention the so-called conformal-traceless ADM systems.Comment: 11 pages, RevTeX, 2 eps figure
Adjusted ADM systems and their expected stability properties: constraint propagation analysis in Schwarzschild spacetime
In order to find a way to have a better formulation for numerical evolution
of the Einstein equations, we study the propagation equations of the
constraints based on the Arnowitt-Deser-Misner formulation. By adjusting
constraint terms in the evolution equations, we try to construct an
"asymptotically constrained system" which is expected to be robust against
violation of the constraints, and to enable a long-term stable and accurate
numerical simulation. We first provide useful expressions for analyzing
constraint propagation in a general spacetime, then apply it to Schwarzschild
spacetime. We search when and where the negative real or non-zero imaginary
eigenvalues of the homogenized constraint propagation matrix appear, and how
they depend on the choice of coordinate system and adjustments. Our analysis
includes the proposal of Detweiler (1987), which is still the best one
according to our conjecture but has a growing mode of error near the horizon.
Some examples are snapshots of a maximally sliced Schwarzschild black hole. The
predictions here may help the community to make further improvements.Comment: 23 pages, RevTeX4, many figures. Revised version. Added subtitle,
reduced figures, rephrased introduction, and a native checked. :-
Development of laser pulse sampled DC voltammetry and its application to the determination of glucose
ArticleBUNSEKI KAGAKU. 57(1): 61-65 (2008)journal articl
Constructing hyperbolic systems in the Ashtekar formulation of general relativity
Hyperbolic formulations of the equations of motion are essential technique
for proving the well-posedness of the Cauchy problem of a system, and are also
helpful for implementing stable long time evolution in numerical applications.
We, here, present three kinds of hyperbolic systems in the Ashtekar formulation
of general relativity for Lorentzian vacuum spacetime. We exhibit several (I)
weakly hyperbolic, (II) diagonalizable hyperbolic, and (III) symmetric
hyperbolic systems, with each their eigenvalues. We demonstrate that Ashtekar's
original equations form a weakly hyperbolic system. We discuss how gauge
conditions and reality conditions are constrained during each step toward
constructing a symmetric hyperbolic system.Comment: 15 pages, RevTeX, minor changes in Introduction. published as Int. J.
Mod. Phys. D 9 (2000) 1
Hyperbolic formulations and numerical relativity II: Asymptotically constrained systems of the Einstein equations
We study asymptotically constrained systems for numerical integration of the
Einstein equations, which are intended to be robust against perturbative errors
for the free evolution of the initial data. First, we examine the previously
proposed "-system", which introduces artificial flows to constraint
surfaces based on the symmetric hyperbolic formulation. We show that this
system works as expected for the wave propagation problem in the Maxwell system
and in general relativity using Ashtekar's connection formulation. Second, we
propose a new mechanism to control the stability, which we call the ``adjusted
system". This is simply obtained by adding constraint terms in the dynamical
equations and adjusting its multipliers. We explain why a particular choice of
multiplier reduces the numerical errors from non-positive or pure-imaginary
eigenvalues of the adjusted constraint propagation equations. This ``adjusted
system" is also tested in the Maxwell system and in the Ashtekar's system. This
mechanism affects more than the system's symmetric hyperbolicity.Comment: 16 pages, RevTeX, 9 eps figures, added Appendix B and minor changes,
to appear in Class. Quant. Gra
- …