I have a technical question about numerical relativity, hopefully someone can help.(adsbygoogle = window.adsbygoogle || []).push({});

In the usual 3 + 1 decomposition in NR, the four gauge freedoms are expressed via the lapse, [tex]\alpha[/tex] and three shift components [tex] \beta^i [/tex]. In a finite element numerical scheme, each grid point will have a value for these 4 components.

Now, looking at evolution equations on the other hand, I see terms that want to take the covariant derivative of the lapse function, which I can't understand since it is a scalar, and hence the covariant derivative reduces to the regular derivative.

To be explicit, take for instance the ADM evolution equation for the extrinsic curvature. The textbook I have ("Elements of Numerical Relativity" by Carles Bona and Carlos Palenzuela-Luque), write this down as

[tex] (\partial_t - L_{\beta} ) K_{ij} = -\alpha_{j;i} + \alpha [ ....] [/tex]

I've left the rest of the equation out for simplicity (note that [tex]L_{\beta}[/tex] is the Lie Derivative, I couldn't work out how to make the nice curly L with the tex tags). See that both [tex]\alpha[/tex] and [tex]\alpha^i[/tex] appear which I don't understand. It would make sense if this book was using [tex]\alpha^i[/tex] to denote the shift vector, but as you can see from the LHS (and it made clear in the book) it uses [tex]\beta^i[/tex] for this.

Any ideas? I'm just replacing the Covariant derivative with the regular one for [tex]\alpha[/tex] but maybe the equation is actually telling me something very different?

**Physics Forums - The Fusion of Science and Community**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Numerical Relativity: Components of the Lapse Function?

Loading...

Similar Threads - Numerical Relativity Components | Date |
---|---|

I Non-diagonal metric for testing elliptic PDE solver | Jun 22, 2017 |

B Numerical solution to GRT problems? | May 7, 2017 |

A How to Integrate the geodesic equations numerically? | Mar 7, 2017 |

A Numerical integration of motion | Sep 1, 2016 |

A Numerical simulation the n-body problem under SR model | Apr 23, 2016 |

**Physics Forums - The Fusion of Science and Community**