Astronomy and Astrophysics – Astrophysics
Scientific paper
2003-12-02
Astrophys.J. 595 (2003) 980-991
Astronomy and Astrophysics
Astrophysics
27 pages, 5 figures
Scientific paper
10.1086/377166
A scheme is presented for accurately propagating the gravitational field constraints in finite difference implementations of numerical relativity. The method is based on similar techniques used in astrophysical magnetohydrodynamics and engineering electromagnetics, and has properties of a finite differential calculus on a four-dimensional manifold. It is motivated by the arguments that 1) an evolutionary scheme that naturally satisfies the Bianchi identities will propagate the constraints, and 2) methods in which temporal and spatial derivatives commute will satisfy the Bianchi identities implicitly. The proposed algorithm exactly propagates the constraints in a local Riemann normal coordinate system; {\it i.e.}, all terms in the Bianchi identities (which all vary as $\partial^3 g$) cancel to machine roundoff accuracy at each time step. In a general coordinate basis, these terms, and those that vary as $\partial g\partial^2 g$, also can be made to cancel, but differences of connection terms, proportional to $(\partial g)^3$, will remain, resulting in a net truncation error. Detailed and complex numerical experiments with four-dimensional staggered grids will be needed to completely examine the stability and convergence properties of this method. If such techniques are successful for finite difference implementations of numerical relativity, other implementations, such as finite element (and eventually pseudo-spectral) techniques, might benefit from schemes that use four-dimensional grids and that have temporal and spatial derivatives that commute.
No associations
LandOfFree
Constrained Transport Algorithms for Numerical Relativity. I. Development of a Finite Difference Scheme does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Constrained Transport Algorithms for Numerical Relativity. I. Development of a Finite Difference Scheme, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Constrained Transport Algorithms for Numerical Relativity. I. Development of a Finite Difference Scheme will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-687507