Mixed Hyperbolic - Second-Order Parabolic Formulations of General Relativity

Two new formulations of general relativity are introduced. The first one is a parabolization of the Arnowitt, Deser, Misner (ADM) formulation and is derived by addition of combinations of the constraints and their derivatives to the right-hand-side of the ADM evolution equations. The desirable property of this modification is that it turns the surface of constraints into a local attractor because the constraint propagation equations become second-order parabolic independently of the gauge conditions employed. This system may be classified as mixed hyperbolic - second-order parabolic. The second formulation is a parabolization of the Kidder, Scheel, Teukolsky formulation and is a manifestly mixed strongly hyperbolic - second-order parabolic set of equations, bearing thus resemblance to the compressible Navier-Stokes equations. As a first test, a stability analysis of flat space is carried out and it is shown that the first modification exponentially damps and smoothes all constraint violating modes. These systems provide a new basis for constructing schemes for long-term and stable numerical integration of the Einstein field equations.Comment: 19 pages, two column, references added, two proofs of well-posedness added, content changed to agree with submitted version to PR

Persistent holes in a fluid

We observe stable holes in a vertically oscillated 0.5 cm deep aqueous suspension of cornstarch for accelerations a above 10g. Holes appear only if a finite perturbation is applied to the layer. Holes are circular and approximately 0.5 cm wide, and can persist for more than 10^5 cycles. Above a = 17g the rim of the hole becomes unstable producing finger-like protrusions or hole division. At higher acceleration, the hole delocalizes, growing to cover the entire surface with erratic undulations. We find similar behavior in an aqueous suspension of glass microspheres.Comment: 4 pages, 6 figure

Numerical performance of the parabolized ADM (PADM) formulation of General Relativity

In a recent paper the first coauthor presented a new parabolic extension (PADM) of the standard 3+1 Arnowitt, Deser, Misner formulation of the equations of general relativity. By parabolizing first-order ADM in a certain way, the PADM formulation turns it into a mixed hyperbolic - second-order parabolic, well-posed system. The surface of constraints of PADM becomes a local attractor for all solutions and all possible well-posed gauge conditions. This paper describes a numerical implementation of PADM and studies its accuracy and stability in a series of standard numerical tests. Numerical properties of PADM are compared with those of standard ADM and its hyperbolic Kidder, Scheel, Teukolsky (KST) extension. The PADM scheme is numerically stable, convergent and second-order accurate. The new formulation has better control of the constraint-violating modes than ADM and KST.Comment: 20 two column pages, 20 figures, submitted to PRD, two typos correcte

Difference methods for initial-value problems

Includes bibliography.Mode of access: Internet

Los Alamos National Laboratory Report LA-1914 (Del.)

The initial growth of irregularities on an interface between two compressible fluids is studied for impulsive (i.e., shock) acceleration. It was found that the ultimate rate of growth is roughly the same as that given by the incompressible theory, if the initial compression of the irregularities and of the fluids is taken into account