David Fiske
Making the constraint hypersurface an attractor in free evolution
Free evolution, in which the constraints of a system are solved only on the
initial hypersurface and not enforced during a numerical evolution, is a
commonly used technique in numerical relativity. Numerous groups in the
field, however, have seen in their numerical experiments that the exact form
in which the evolution equations are written can greatly influence the
stability of numerical simulations. I will discuss a particular way of
modifying the evolution equations of the Maxwell system such that the
constraint hypersurface (div E = 0) becomes an attractor of the
system of evolution equations.
References:
Back to ECGM6 Program