<div dir="ltr"><div><div>Thanks for the reply. My apologies that I&#39;m not completely familiar many of the physical and numerical issues while simulating bbh. I have experience in AMR astrophysical simulations. <br><br></div><div>When I do the bbh simulation of ratio around 1, it evolves perfectly fine. With higher mass ratio, the smaller BH just becomes bigger and blows up. I&#39;m assuming this as a problem of refinement.<br></div><div><br></div>My naive picture was that the grid is created adaptively at every n time step(a parameter) based on lapse values. After locating the horizon, dynamic variables are updated in the outside region only. So, the problem remains with locating the horizon often enough and refining to the appropriate level so that at least 4-5 cells remains across a BH.<br></div><div>Please correct me, if I&#39;m off by any big misunderstanding. <br><br></div><div>Thanks again,<br></div>-rahul<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 6, 2016 at 4:22 PM, Frank Loeffler <span dir="ltr">&lt;<a href="mailto:knarf@cct.lsu.edu" target="_blank">knarf@cct.lsu.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On Mon, Jun 06, 2016 at 12:41:35PM -0400, rahul kashyap wrote:<br>
&gt; I&#39;m working on binary black hole mergers with large mass ratio (20:1 or<br>
&gt; more). By using higher refinement level, I have been able to resolve the<br>
&gt; smaller BH, but it doesn&#39;t remain resolved even for a single time step.<br>
<br>
</span>Hi,<br>
<br>
Thanks for contacting us. I am not that familiar with that particular<br>
type of simulation (bbh with that mass ratio), but I try to give a best<br>
guess. Please everyone: if you disagree, say so. :)<br>
<br>
It would be interesting to know what exactly you mean with &#39;doesn&#39;t<br>
remain resolved&#39;. Assuming that the resolution is high enough, one thing<br>
I would expect a simulation like this to use is some kind of<br>
position-dependent gauge condition. The reason is that one thing the gauge<br>
mostly does is (loosely speaking) to &quot;prevent the grid points from<br>
falling into the black hole, or to come out of it&quot;. In your simulation<br>
you see this as a black hole that gets &#39;larger or smaller&#39; without<br>
changing its mass a lot. The best parameter values that determine how<br>
much the grid points are &quot;pushed out&quot; depend on the black hole mass.<br>
Since you mass ratio is quite large, I would expect the optimal<br>
parameter to be quite different for the two black holes.<br>
<br>
There are groups that did evolve BBHs with even higher mass ratios than<br>
this - not using the exact same code, but the same methods. They should<br>
have much more experience. Anyone - you know whom I mean ...?<br>
<span class="HOEnZb"><font color="#888888"><br>
Frank<br>
<br>
</font></span></blockquote></div><br></div>