[Users] problems of reproducing the simulation of binary neutron star with different mass

Mon May 14 13:34:26 CDT 2018

Dear Roland,

Yes, I did use the parameter file from gallery page and use unequal mass initial data from Lorene.

I just check the initial location of the binary star is (-34.59 km , 45.41 km).However you said it makes the hard-coded assumption that if one of the stars is at (x,y,z) the other one is at (-x,-y,z) .How should I set CarpetRegrid2::position_x_1 and CarpetRegrid2::position_x_2 ?

Thank you.

Best regards,

Chia-Hui

________________________________
寄件者: Roland Haas <rhaas at illinois.edu>
寄件日期: 2018年5月15日 上午 02:13:06
收件者: 林家暉
副本: Einstein Toolkit Users
主旨: Re: [Users] problems of reproducing the simulation of binary neutron star with different mass

Hello Chia-Hui,

thank you for the attached file.

Yes, the error is related to the grid structure:

WARNING level 0 from host nid05920 process 0
  while executing schedule bin (none), routine (no thorn)::(no routine)
  in thorn CarpetLib, file /global/u1/c/chlin/Cactus/arrangements/Carpet/CarpetLib/src/gh.cc:178:
  -> The refinement hierarchy is not properly nested.

This is for the run where you used the parameter file from the gallery
page but used unequal mass initial data from LORENE?

I am asking because, given that the parameter file specifies
Rotating180 symmetry, the evolution will be as if equal mass NS were
given in the initial data file (it will duplicate the right hand NS and
ignore the left hand one). So the error that you saw should not be due
to the unequal mass per-se. It could be if the location of the stars
changed initially (eg for unequal mass the more massive one will be
closer to the origin than the other one). This still should not give
you this particular error (it would be wrong in other ways though)
since the parameter file uses CarpetRegrid2::ensure_proper_nesting   =
"yes" which should avoid this particular error.

Looking at the plot from your initial email it may be that there was an
error in the initial setups as it looks to me as if the neutron stars
have actually drifted apart. At least the gallery examples has the NS
initially at ~+/-15 units from the origin but in you plot it looks like
much more.

Wen you added the new ID you must make sure to adjust the initial
location in the parfile in the lines reading:

CarpetRegrid2::position_x_1 = 15.1875
CarpetRegrid2::position_x_2 = -15.1875

The binary black hole example does not have the restriction of
assuming 180 degree symmetry, no.

Yours,
Roland

> Dear Roland,
> Thank you for your reply.
> The attached file is .err file corresponding to the simulation of binary star with different central enthalpy.  Can you have a look at it ?Thank you.
> For not able to track multiple stars, does it occur in binary blackhole case(http://einsteintoolkit.org/gallery/bbh/index.html) ?
> Best regards,
> Chia-Hui
> ________________________________
> 寄件者: Roland Haas <rhaas at illinois.edu>
> 寄件日期: 2018年5月13日 下午 09:28:51
> 收件者: 林家暉
> 副本: Einstein Toolkit Users
> 主旨: Re: [Users] problems of reproducing the simulation of binary neutron star with different mass
>
> Hello Chia-Hui Lin,
>
> changing the central enthalpy seems like the right thing to do to
> change masses in LORENE since it prescribes the central enthalpy of
> each star.
>
> The backtrace that you showed and the error message is unfortunately
> not detailed enough to make a guess as to what could have gone wrong.
> This output is produced by a call to CCTK_Error() however the
> CCTK_Error function takes an error string as an argument which is not
> present in the lines that you have copied. It is likely in the *.err
> file that the simulation produced and you should try and include that
> file as an attachment in your email.
>
> If you have non-equal mass stars then you have to remove the 180 degree
> rotation symmetry from the parameter files. Ie you have to change the
> lines (these are scattered throughout throughout the par file)
>
> ActiveThorns = "RotatingSymmetry180 ReflectionSymmetry"
> CoordBase::xmin =    0.00
> CoordBase::boundary_shiftout_x_lower = 1
> CarpetRegrid2::symmetry_rotating180    = "yes"
> Hydro_Analysis::Hydro_Analysis_rho_max_loc_only_positive_x = "true"
> Hydro_Analysis::Hydro_Analysis_average_multiple_maxima_locations = "yes"
>
> to
>
> ActiveThorns = "ReflectionSymmetry"
> CoordBase::xmin = 396.00
> CoordBase::boundary_shiftout_x_lower = 0
> CarpetRegrid2::symmetry_rotating180    = "no"
> Hydro_Analysis::Hydro_Analysis_rho_max_loc_only_positive_x = "false"
> Hydro_Analysis::Hydro_Analysis_average_multiple_maxima_locations = "no"
>
> Please note though that currently the thorn NSTracker and HydroAnalysis *cannot* properly track multiple stars. It makes the hard-coded assumption that if one of the stars is at (x,y,z) the other one is at (-x,-y,z) which is only true for equal mass neutron stars. This will mean that one of the stars will eventually leave the box that is supposed to track it.
>
> Yours,
> Roland
>
> > Dear whom it may concern,
> >
> > I am trying to reproduce the simulation of binary neutron star (http://einsteintoolkit.org/gallery/bns/index.html) with different mass. I use Lorene code to produce the initial data . The parameter files I used is locate in /Lorene/Codes/Bin_star/Parameters/Polytrope_irrot/Test_GR . And I changed  central enthalpy of each star in par_init.d to get different mass (Is there other way I can get different mass?). And then I simulate with the ET. However I met some problems:
> >
> > 1.I requested 48 cores with 10 hours to run the code , however the code terminated after running ~4 hours. And the end of .err file reads:
> >
> >
> >
> > cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.
> >
> > cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.
> >
> > Rank 0 with PID 47930 received signal 6
> >
> > Writing backtrace to nsnstohmns/backtrace.0.txt
> >
> > Rank 1 with PID 47931 received signal 6
> >
> > Writing backtrace to nsnstohmns/backtrace.1.txt
> >
> > cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.
> >
> > cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.
> >
> > Rank 2 with PID 13239 received signal 6
> >
> > Rank 3 with PID 13240 received signal 6
> >
> > Writing backtrace to nsnstohmns/backtrace.2.txt
> >
> > Writing backtrace to nsnstohmns/backtrace.3.txt
> >
> > srun: error: nid05930: tasks 2-3: Aborted
> >
> > srun: Terminating job step 9192104.0
> >
> > srun: error: nid05920: tasks 0-1: Aborted
> >
> >
> > The code is aborted without completing the whole process until the merger. How can I perform the complete process?
> >
> > 2. When I use visit to do the visualization ,I found the file rho.xy.h5 only has half of the whole simulation as :
> >
> > [cid:e66a53d5-d406-4599-997f-00e90a864dc0]
> >
> >
> >
> >
> >
> > The tutorial(scripts.tar.gz<http://einsteintoolkit.org/gallery/bns/scripts.tar.gz>) tells that I should use symmetric reflection to perform the whole visualization. However I cannot use this method since I simulated binary star with different mass . So there is no symmetry. How can I get the whole picture instead of only half ?
> >
> >
> > Thank you.
> >
> > Best regards,
> >
> > Chia-Hui Lin
>
>
>
> --
> My email is as private as my paper mail. I therefore support encrypting
> and signing email messages. Get my PGP key from http://pgp.mit.edu .

--
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://pgp.mit.edu .
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20180514/37c54a14/attachment-0001.html 