<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=big5">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p><br>
</p>
<meta content="text/html; charset=UTF-8">
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size: 12pt; color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif, Helvetica, EmojiFont, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols;">
<p>Dear Roland,</p>
<p>Thank you for the detailed information.</p>
<p>I have tried to simulate unequal mass binary neutron star with new setting of initial position of the two stars.It seems fine at least in first several orbits .I will try the thorn you mentioned to check the difference. <br>
Thank you for your help.</p>
<p>Best regards,</p>
<p>Chia-Hui</p>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>寄件者:</b> Roland Haas <rhaas@illinois.edu><br>
<b>寄件日期:</b> 2018年5月15日 下午 08:16:12<br>
<b>收件者:</b> 林家暉<br>
<b>副本:</b> Einstein Toolkit Users<br>
<b>主旨:</b> Re: [Users] problems of reproducing the simulation of binary neutron star with different mass</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Hello Chia-Hui,<br>
<br>
> I do not quite understand what causes ,in your example,<br>
> <br>
> CarpetRegrid2::position_x_1 = 42.<br>
> CarpetRegrid2::position_x_2 = -17<br>
> <br>
> to become fail.<br>
This has to do with the way that the thorn GRHydro_Analysis that looks<br>
for the location of the star works. It first determines the location on<br>
the grid with the largest density, then it computes a center of mass<br>
integral in a region around that point. This computed center of mass if<br>
used as the location of the star and to move boxes. There is code in<br>
thorn NSTracker that uses this single found location and moves two<br>
refined regions to (x,y,z) and (-x,-y,z) where (x,y,z) is the found<br>
location of the star.<br>
<br>
There are two reasons why this does not work well with unequal masses:<br>
<br>
1. it only looks for a single star and bases that star's location on<br>
the highest density<br>
2. NSTracker makes the hard-coded assumption that the two regions that<br>
cover the star are mirror immages of each other<br>
<br>
In your case since the star at -17 (always assuming you actually have<br>
set up a grid that has no symmetries) is the heavier one (since it is<br>
closer to the center of rotation) and thus has the higher central<br>
density, therefore GRHydro_Analysis will use it as the star to track<br>
and thus will find the star location to be (-17,0,0) when it looks for<br>
the star.<br>
<br>
> In my case, if the initial location of the binary star is (-34.59<br>
> km , 45.41 km) , I should set<br>
> <br>
> CarpetRegrid2::position_x_1 = -34.59/1.48=-23.37<br>
> CarpetRegrid2::position_x_2 = 45.41/1.48=30.68<br>
> <br>
> right? <br>
Yes that is correct.<br>
<br>
> Would this soon be set to (-23.37,23.37)? Because the mass difference<br>
> is not so much, I am wondering whether there exists some maximum value<br>
> of mass ratio which can tolerate small mass difference.<br>
> <br>
> By the way ,the link <a href="https://trac.einsteintoolkit.org/ticket/1974" id="LPlnk178019" previewremoved="true">
https://trac.einsteintoolkit.org/ticket/1974</a><br>
> seems not work. I would like to check it up.<br>
Thank you for pointing this out. It seems that there is an endless<br>
redirection loop in the trac website already for<br>
<a href="https://trac.einsteintoolkit.org" id="LPlnk124250" previewremoved="true">https://trac.einsteintoolkit.org</a><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" id="LPlnk335348" previewremoved="true">
http://pgp.mit.edu</a> .<br>
</div>
</span></font></div>
</body>
</html>