Numerical relativity solves Einstein's equations as an initial-value problem, usually after decomposing spacetime into spatial hypersurfaces and evolution in time.
A 3+1 decomposition writes in terms of lapse, shift, and spatial metric.
The lapse function measures proper-time separation between neighboring spatial hypersurfaces along their unit normal.
The shift vector describes how spatial coordinates move tangentially between neighboring hypersurfaces.
With future unit normal , the convention measures how the hypersurface is embedded in spacetime.
For with spatial , current conservation gives under the stated extrinsic-curvature convention.
A harmonic coordinate satisfies , equivalently or .
Harmonic slicing imposes the harmonic-coordinate condition only on time. In 3+1 form it obeys .
The Bona--Masso family is . Harmonic slicing is the choice .
In harmonic gauge, the vacuum Einstein equations have principal part , forming a quasilinear wave system supplemented by constraints.

Articles by others on the same topic (1)

Numerical relativity is a subfield of computational physics that focuses on solving the equations of general relativity using numerical methods. General relativity, formulated by Albert Einstein, describes the gravitational interaction as a curvature of spacetime caused by mass and energy. The equations governing this curvature, known as the Einstein field equations, are highly complex and often impossible to solve analytically in realistic scenarios, especially in dynamic situations like the collision of black holes or neutron stars.