<div dir="ltr">Hi Roland,<div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span style="color:rgb(80,0,80)">> Additionally, I find it convenient to smooth all quantities near the<br></span><span style="color:rgb(80,0,80)">> singularity, and to choose to advect both lapse and shift. I don't<br></span><span style="color:rgb(80,0,80)">> know whether the latter is necessary in theory, but I am always using<br></span><span style="color:rgb(80,0,80)">> it.<br></span>Concerning the smoothing it seems that some smoothing was necessary to<br>not crash the code (via the Exact epsilon parameter to transform r -><br>r+epsilon). An alternative may have been to use NoExcision to fill the<br>interior of the AH with smooth data.</blockquote><div><span class="gmail-im" style="color:rgb(80,0,80)"><br></span>I generally suggest shifting the grid so that the closest gridpoint to the origin is (dx/2, dy/2, dz/2), where dx,dy,dz are the resolutions on the finest grid. I've used this trick many times to stabilize evolutions. Also you may want to try the ShiftedKerrSchild thorn (in the ETK), which enables a radial offset that shrinks the BH coordinate size but at the same time kills off enormous constraint violations near r=0. Unlike the existing KerrSchild thorn(s), this one uses the standard Kerr-Schild metric written in spherical coordinates and does the basis xform.<br><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:12.8px">-Zach</div><div style="font-size:12.8px"><br></div><span style="font-size:12.8px">* * *</span><br style="font-size:12.8px"><span style="font-size:12.8px">Zachariah Etienne</span><br style="font-size:12.8px"><span style="font-size:12.8px">Assistant Professor of Mathematics</span><div style="font-size:12.8px">West Virginia University</div><div><font color="#0000ee"><u><a href="https://math.wvu.edu/~zetienne/" target="_blank">https://math.wvu.edu/~zetienne/</a></u></font></div><div><a href="https://blackholesathome.net" style="font-size:12.8px" target="_blank">https://blackholesathome.net</a><br></div></div></div></div></div></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 25, 2019 at 12:30 PM Haas, Roland <<a href="mailto:rhaas@illinois.edu">rhaas@illinois.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello Erik,<br>
<br>
> This is almost possible. For puncture data, one typically evolves<br>
> lapse, shift, and a quantity B (the time derivative of the shift). For<br>
> Kerr-Schild data, one also needs to evolve A, the time derivative of<br>
> the lapse. Otherwise, Kerr-Schild data are not stationary. (One could<br>
> instead add an offset alpha_0 to the evolution equations for K, but<br>
> that is quite non-standard and not implemented in McLachlan.)<br>
Ok so it is possible.<br>
<br>
> Additionally, I find it convenient to smooth all quantities near the<br>
> singularity, and to choose to advect both lapse and shift. I don't<br>
> know whether the latter is necessary in theory, but I am always using<br>
> it.<br>
Concerning the smoothing it seems that some smoothing was necessary to<br>
not crash the code (via the Exact epsilon parameter to transform r -><br>
r+epsilon). An alternative may have been to use NoExcision to fill the<br>
interior of the AH with smooth data.<br>
<br>
> I can send a sample parameter file if that helps.<br>
That would be great if you had something that you could send around<br>
easily.<br>
<br>
Yours,<br>
Roland<br>
<br>
-- <br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://pgp.mit.edu" rel="noreferrer" target="_blank">http://pgp.mit.edu</a> .<br>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote></div>