Hi Maria,
On Tuesday 2016-06-07 10:02, Hamilton, Maria wrote:
Date: Tue, 7 Jun 2016 10:02:36
From: "Hamilton, Maria" <[email protected]>
To: Peter Diener <[email protected]>, Frank Loeffler <[email protected]>
Cc: Einstein Toolkit Users <[email protected]>
Subject: Re: [Users] Resolving two BH
Hello,
Does this paper come with access to a par file that will show how to
set the dynamic shift?
No it doesn't. That paper was based on the Jena groups BAM code and
these features have not yet been implemented in McLachlan.
I need some help as well in setting up a run with two unequal mass
black holes, for an undergraduate visualization project. My
understanding is that the initial conditions are set with the
parameters given in the TwoPunctures thorn. For example, if I want the
first black hole to be equal with the mass of the sun, and the second,
10 times the mass of the sun, I do something like this?
TwoPunctures::par_m_plus = 1.5
TwoPunctures::par_m_minus = 15
The code uses units so G=c=1 so you can get mass you want. But yes, this
will give you a mass ratio of 10. Note however, that if you want this to
be the physical masses, you have to set
TwoPunctures::give_bare_mass = no
as well. Otherwise these will be the bare mass parameters given as input
to the elliptical solve, and the actual physical masses after the solve
might be a little bit different. One comment: I usually preferred to
have the total mass be 1.0 in order to make grid setup as well
as analysis simpler between runs with different mass ratio. So I would
use
TwoPunctures::par_m_plus = 0.09090909090909
TwoPunctures::par_m_minus = 0.9090909090909
But That is a matter of taste of course.
The last thing I need to set is the velocity for each hole, If I want
each of them to move at 1/3 c tangential velocity? I’m assuming I need
to set:
TwoPunctures::par_P_plus [1] = +0.45
TwoPunctures::par_P_minus[1] = -4.5
You most likely want to do the simulation in the rest frame, where the
total momentum is 0. Otherwise your black hole system will drift away
from the system. So the larger black hole will move 10 times slower than
the small black hole. How big a velocity you want depends on exactly
what kind of system you want to evolve. Do you want a hyperbolic flyby,
an elliptical orbit or a quasi-circular orbit?
Is this correct? And one more important thing: if I want the distance
between them to be for example 150, how do I set this? I see a
parameter par_b which gives the x coordinate of the m+ puncture, but
I’m not clear if this is what I need to set, and how.
Yes, par_b gives the x-coordinate of the m+ puncture. By default the m-
puncture will then be at -par_b (on the x-axis). So par_b is 1/2 the
coordinate separation of the punctures. There is in addition a parameter
center_offset that allow you to add an offset to the initial position.
Again the choice of the separation and of the momenta goes hand in hand
to determine what kind of orbit you'll get.
All the other parameters I am assuming are to be kept the way they
are. In Rahul’s parameter file I do not see the PunctureTracker thorn
activated. I do want to know the position of the back holes on the
grid, and the distance between them at each step of the evolution. Do
I need to activate PunctureTracker? What output do I need for this?
As far as I remember you just need to activate PunctureTracker and then
set a few parameters to activate tracking as well as setting the initial
puncture positions consistently with the TwoPuncture parameters.
And as far as I remember PunctureTracker keeps track of the punctures in
Cactus grid scalars and you'd just need to request output for those.
Thank you,
Maria
Cheers,
Peter
On 6/7/16, 10:13 AM, "[email protected] on behalf of Peter Diener"
<[email protected] on behalf of [email protected]> wrote:
Hi,
This issue is discussed for example in
http://arxiv.org/abs/1003.4681
where a non-constant damping in the gamma driver shift is presented to
handle unequal mass binaries. As far as I remember there may also be
later papers with further improvements to the method.
Cheers,
Peter
On Tuesday 2016-06-07 01:00, Frank Loeffler wrote:
Date: Tue, 7 Jun 2016 01:00:59
From: Frank Loeffler <[email protected]>
To: rahul kashyap <[email protected]>
Cc: Einstein Toolkit Users <[email protected]>
Subject: Re: [Users] Resolving two BH
On Mon, Jun 06, 2016 at 06:50:04PM -0400, rahul kashyap wrote:
Thanks for the reply. My apologies that I'm not completely familiar many of
the physical and numerical issues while simulating bbh. I have experience
in AMR astrophysical simulations.
No need to apologize.
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'm assuming this as a problem of refinement.
It is not. It looks like what I described earlier: your gauge
conditions, in particular the shift condition, isn't well suited for the
small black hole. Grid points are "falling in", making it appear growing
on the computational grid. You may try to play with the values such that
the small black hole evolves fine (you can try that with a simulation of
only a small black hole of that size), but then you might find that the
large black hole "vanishes from the grid" because grid points there are
pushed out too much. If this is the case, you might need
position-dependent values.
Frank
_______________________________________________
Users mailing list
[email protected]
http://lists.einsteintoolkit.org/mailman/listinfo/users
_______________________________________________
Users mailing list
[email protected]
http://lists.einsteintoolkit.org/mailman/listinfo/users
