<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)">&gt; Additionally, I find it convenient to smooth all quantities near the<br></span><span style="color:rgb(80,0,80)">&gt; singularity, and to choose to advect both lapse and shift. I don&#39;t<br></span><span style="color:rgb(80,0,80)">&gt; know whether the latter is necessary in theory, but I am always using<br></span><span style="color:rgb(80,0,80)">&gt; 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 -&gt;<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&#39;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 &lt;<a href="mailto:rhaas@illinois.edu">rhaas@illinois.edu</a>&gt; 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>
&gt; This is almost possible. For puncture data, one typically evolves<br>
&gt; lapse, shift, and a quantity B (the time derivative of the shift). For<br>
&gt; Kerr-Schild data, one also needs to evolve A, the time derivative of<br>
&gt; the lapse. Otherwise, Kerr-Schild data are not stationary. (One could<br>
&gt; instead add an offset alpha_0 to the evolution equations for K, but<br>
&gt; that is quite non-standard and not implemented in McLachlan.)<br>
Ok so it is possible.<br>
<br>
&gt; Additionally, I find it convenient to smooth all quantities near the<br>
&gt; singularity, and to choose to advect both lapse and shift. I don&#39;t<br>
&gt; know whether the latter is necessary in theory, but I am always using<br>
&gt; 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 -&gt;<br>
r+epsilon). An alternative may have been to use NoExcision to fill the<br>
interior of the AH with smooth data.<br>
<br>
&gt; 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>