Dear all,
I encountered some problem while carrying out some binary neutron star simulations
and would like to know if anyone can help me out.
I followed the BNS example on the site (nsnstohmns) but adopted the "Hybrid" eos for
the BNS Meudon initial data generated with certain tabular nuclear models.
#####-----#####-----#####-----#####-----#####-----#####-----#####-----#####-----#####
Problem 1:
I understand that I need to assign the following parameters accordingly,
Meudon_Bin_NS::eos_table = "Hybrid"
EOS_Omni::hybrid_gamma_th
EOS_Omni::n_pieces
EOS_Omni::hybrid_k0
EOS_Omni::hybrid_gamma[]
EOS_Omni::hybrid_rho[]
but I am not sure whether it is correct to also set
GRHydro::GRHydro_eos_table = "Hybrid" ?
or should I set
GRHydro::GRHydro_eos_table = "Ideal_Fluid" ?
In the BNS example nsnstohmns, despite using a polytropic BNS model, the parameter file shows
GRHydro::GRHydro_eos_type = "General" instead of "Polytype"
GRHydro::GRHydro_eos_table = "Ideal_Fluid" instead of "2D_Polytrope"
What is the particular reason for this setting ??
#####-----#####-----#####-----#####-----#####-----#####-----#####-----#####-----#####
Problem 2:
Here is one of the error I encountered.
WARNING level 1 from host emc4 process 0
in thorn CarpetLib, file /home/hkhui/Cactus/configs/sim/build/CarpetLib/dh.cc:161:
->
/home/hkhui/Cactus/configs/sim/build/CarpetLib/dh.cc:837:
[ml=0 rl=7 c=0] The following grid structure consistency check failed:
Refinement prolongation: All points must have been received
needrecv.empty()
WARNING level 1 from host emc4 process 0
in thorn CarpetLib, file /home/hkhui/Cactus/configs/sim/build/CarpetLib/dh.cc:161:
->
/home/hkhui/Cactus/configs/sim/build/CarpetLib/dh.cc:983:
[ml=0 rl=7 c=0] The following grid structure consistency check failed:
Synchronisation and boundary prolongation: All points must have been received
needrecv.empty()
It seems this error is associated with the following settings in the parfile
CoordBase::xmin = 0.00
CoordBase::ymin = -396.00
CoordBase::zmin = 0.00
CoordBase::xmax = +396.00
CoordBase::ymax = +396.00
CoordBase::zmax = +396.00
CoordBase::dx = 9.00
CoordBase::dy = 9.00
CoordBase::dz = 9.00
CarpetRegrid2::num_levels_1 = 7
CarpetRegrid2::position_x_1 = 15.1875
CarpetRegrid2::radius_1[1] =240.0
CarpetRegrid2::radius_1[2] =120.0
CarpetRegrid2::radius_1[3] = 60.0
CarpetRegrid2::radius_1[4] = 28.125 # nsnstohmns 26.125
CarpetRegrid2::radius_1[5] = 19.875 # nsnstohmns 17.875
CarpetRegrid2::radius_1[6] = 11.0 # nsnstohmns 13
The initial size of the stars of my model:
Coordinate separation between the two stellar centers : 45 km
Coordinate equatorial radius (phi=0) a1 = 10.16960419148 km
Coordinate equatorial radius (phi=pi/2) a2 = 9.872118702179 km
Coordinate equatorial radius (phi=pi): 10.19931065979 km
Coordinate polar radius a3 = 9.925176860927 km
I used the nsnstohmns parfile while changing the eos part and the resolution only,
for the set of par above, I set [5] and [6] so that they together cover the separation (19.875+11 M > 45km)
and I set [4] so that this refinement region of the 1st star cover some region of the 2nd star.
I would like to know if there is any rule for setting the AMR region (the radius par above),
I varied these radius parameters and end up with either error about the grid structure (see
above err message) or the c2p conversion error, but never been able to carry out simulation to merger.
If the problem is not about the AMR region radius setting, what could be the problem?
#####-----#####-----#####-----#####-----#####-----#####-----#####-----#####-----#####
Problem 3:
I have carried the BNS simulation with stars of slightly larger R initially with my nuclear eos,
and adopted exactly the same parfile setting as nsnstohmns with Time::dtfac = 0.35,
it happens that the NSTracker thorn loses track of the stars.
The rho_max position suddenly moved so fast that I suspect the result is wrong.
And the position of the star is not updated anymore after certain time step.
What could cause such false result??
-----------------------------------------------------------------------------------------------------------------------
Iteration Time | *me_per_hour | GRHYDRO::dens | *ROBASE::rho | *::w_lorentz | *STRAINTS::H | *axrss_mb
| | sum maximum | maximum | maximum | maximum | maximum
-----------------------------------------------------------------------------------------------------------------------
INFO (NSTracker): Found star at (15.1875,0.28125,0)
200 4.922 | 113.9376773 | 0.0001258 0.0027456 | 0.0011040 | 1.3729137 | 0.0031954 | 989
INFO (NSTracker): Found star at (15.1875,0.5625,0)
400 9.844 | 114.5209435 | 0.0001258 0.0026962 | 0.0011407 | 1.5612539 | 0.0050456 | 991
INFO (NSTracker): Found star at (15.1875,0.84375,0)
600 14.766 | 114.7115204 | 0.0001259 0.0014155 | 0.0006785 | 1.6008531 | 0.0038587 | 993
INFO (NSTracker): Found star at (15.1875,1.125,0)
620 15.258 | 114.1468585 | 0.0001259 0.0011319 | 0.0005486 | 1.6009425 | 0.0036545 | 1084
INFO (NSTracker): Found star at (14.9062,1.6875,0)
660 16.242 | 113.7788398 | 0.0001259 0.0007284 | 0.0003032 | 1.6088940 | 0.0033746 | 1084
INFO (NSTracker): Found star at (14.9062,1.96875,0)
668 16.439 | 113.6618707 | 0.0001259 0.0006774 | 0.0002687 | 1.6088572 | 0.0033143 | 1084
INFO (NSTracker): Found star at (14.625,1.96875,0)
672 16.537 | 113.4731182 | 0.0001259 0.0006559 | 0.0002538 | 1.6088969 | 0.0032825 | 1084
INFO (NSTracker): Found star at (14.625,2.25,0)
676 16.636 | 112.2718493 | 0.0001259 0.0006351 | 0.0002415 | 1.6090522 | 0.0032501 | 993
INFO (NSTracker): Found star at (14.625,2.53125,0)
684 16.833 | 112.1253398 | 0.0001259 0.0005958 | 0.0002184 | 1.6093073 | 0.0031844 | 993
INFO (NSTracker): Found star at (14.625,2.8125,0)
688 16.931 | 111.9626309 | 0.0001259 0.0005780 | 0.0002090 | 1.6111194 | 0.0031493 | 993
...
INFO (NSTracker): New star location not set (0.5625,-0.28125,0) because too far from old (13.7812,7.875,0).
884 21.755 | 99.6235324 | 0.0001259 0.0002227 | 0.0000792 | 1.9277784 | 0.0066220 | 1183
INFO (NSTracker): New star location not set (0.84375,0.28125,0) because too far from old (13.7812,7.875,0).
1004 24.708 | 99.9757198 | 0.0001259 0.0001761 | 0.0001082 | 2.0183873 | 0.0088739 | 1183
INFO (NSTracker): New star location not set (1.125,-0.5625,0) because too far from old (13.7812,7.875,0).
1008 24.806 | 100.0596083 | 0.0001259 0.0001782 | 0.0001092 | 2.0194414 | 0.0089929 | 1183
...
INFO (NSTracker): New star location not set (3.9375,-21.375,0) because too far from old (13.7812,7.875,0).
2244 55.223 | 99.3936843 | 0.0001258 0.0000634 | 0.0000347 | 2.5307425 | 0.0051844 | 1189
...
INFO (NSTracker): New star location not set (11.25,-45,33.75) because too far from old (13.7812,7.875,0).
4028 99.127 | 98.5126340 | 0.0001270 0.0000141 | 0.0000075 | 3.0658306 | 0.0034565 | 1190
#####-----#####-----#####-----#####-----#####-----#####-----#####-----#####-----#####
Any help will be greatly appreciated..
Yours,
Hon-ka
--
Open WebMail Project (http://openwebmail.org)
Present: Peter D., Roland H., Federico C., Miguel G., Zach E., Steve B., Beyhan, Alois S., Maria H., Bilal
Chair: Zach E.
Minutes: Miguel G.
* Online ET workshop<https://www.cct.lsu.edu/Einsteintoolkitworkshop>, August 3-7 2020 [SB].
* Schedule (https://docs.google.com/spreadsheets/d/11waMmS73k6mxBxQVIC6eZH4hFjLsxVpH_tg…)
Steve: Trying to get the agenda together, does not have anything on Friday but a discussion panel could fit.
Zach: suggested future directions in numerical relativity (code development like gpu support) as topic.
Steve: needs titles and abstracts (Peter D. will send it today)
Steve is setting up the new workshop tutorial server, similar to the ETK tutorial. Steve needs volunteers to try the new server (people interested should provide him their email address).
Steve: Speakers should send him a list of their technical needs (python packages, thorns, etc.) and all the extra material needed for their tutorials (jupyter notebooks, etc.).
Roland: asked if there’s a Workshop tutorial repo… Steve said No but he can do it, he’s still making changes and running tests, eventually will require some cleaning.
Steve: Students thinks students should run the workshop tutorial in advance.
Steve: Recording the session would be great.
Steve: 30 people signed up already. He will send a second announcement to the mailing list soon.
* Gallery testing: tasks<https://bitbucket.org/einsteintoolkit/tickets/issues?kind=task&milestone=ET…> [RH]
* BNS: In progress, needs to be uploaded, some differences post-merger.
Roland: No movie for BNS, hopefully the person in charge will have it soon.
* BBH: Peter Schaffarczyk (Kiel), Atul Kedia (Notre Dame) In progress.
Peter: editing webpage where they have collected all the material, should be finished soon.
* gcc-10 / gfortran-10 [RH]
About backporting process:
Roland: asked Zach if there any pending backporting for WVUThorns , Boundary fix should go to Master or not? Zach will make a decision and try to make time to patch in a reasonable time frame.
* Open tickets
No new open tickets.
Zach’s student while reading documentation found that the web browser (safari and Chrome) were not rendering the latex equations correctly. The issue could cause missing math symbols . The problem is the Latex to HTML conversion process, done automatically by Jenkins server… Two possible solutions 1) go back to old version, 2) update Jenkins server. Federico C. could help.
* Extra questions:
Bilal Asked: “Can I solve Einstein's constraints by conformal method using the ETK?”
Zach: Elliptic solvers in the toolkit like CT_MultiLevel thorn could be useful.
Beyhan——Rotation profile in BNS remnants?
Roland: PostCactus wiki page (https://docs.einsteintoolkit.org/et-docs/Analysis_and_post-processing) could be helpful.
* Next minute taker and call chair:
07/30/2020:
Chair: Roland H.
Minute taker: Steven B.
* No call on the Workshop week.
Call ended.
———————————————
Miguel Gracia,
Postdoctoral Researcher,
Center for Relativistic Astrophysics,
Boggs Bldg. 1-35, Georgia Institute of Technology, 770 State Street, Atlanta, GA, 30313, USA.
Hello,
Please consider joining the weekly Einstein Toolkit phone call at
9:00 am US central time on Thursdays. For details on how to connect
and what agenda items are to be discussed, use the link below.
https://docs.einsteintoolkit.org/et-docs/Main_Page#Weekly_Users_Call
--The Maintainers
Respected sir/maam
im facing issue while running simple "Hello world command. im running it on
online sever. Please help me facing this issue. I am attaching the
screenshot for your reference.
Dear colleagues,
I am looking for a postdoctoral researcher working in the area of
computational astrophysics.
The main focus of the planned research is on gamma ray bursts central
engines, and their outflows, including chemical enrichment of the
interstellar plasma by these explosions and heavy element radioactivities.
The post-doc will also collaborate with other team members on the
development of the GR MHD numerical tools and large scale jet modeling.
The position will be for 3 years. The starting date is negotiable but could
be as early as October 2020.
The formal job advertisement has been put on AAS Job register and on the
EAS website.
Please forward this email to anybody that you think might be interested!
Contact address for further questions: agnes @ cft.edu.pl
regards
Agnieszka Janiuk
--
--
dr hab. Agnieszka Janiuk (prof. CFT PAN) | tel. +48 22 116 20 06
Center for Theoretical Physics | fax. +48 22 843 13 69
Polish Academy of Sciences | www:
http://www.cft.edu.pl/astrofizyka
Al. Lotnikow 32/46 |
02-668 Warsaw, Poland |
We are planning to have an online EinsteinTookit Workshop at LSU this
August. You can register for this event.
Registration is free, but it helps us to know how many students to
expect so that we can properly provision our tutorial servers.
Planned Speakers include: Peter Diener, Zachary Etienne, Bruno
Giacomazzo, Roland Haas, Eric Schnetter, Helvi Witek, and Yosef Zlochower.
Students will be able to make short presentations if they are
interested. So far, these include Terrence Pierre Jacques and Patrick
Nelson. More are wellcome.
https://www.cct.lsu.edu/Einsteintoolkitworkshop
Present: Peter D, Bill G, Roland H, Federico C, Yosef Z, Miguel G, Atul K, Zach E, Steve B
Chair: Peter D
Minutes: Bill G
* Online ET Workshop, https://www.cct.lsu.edu/Einsteintoolkitworkshop , August 3-7 2020 [SB]
** Speakers for August meeting https://docs.google.com/spreadsheets/d/11waMmS73k6mxBxQVIC6eZH4hFjLsxVpH_tg…
Use above link to signup for presenting talks. Looking for more workshop-like, or if interested in hearing some topics...lecture
and tutorials. Maybe also lightning talks for small projects, or for students. Have 7 speakers listed (as of now). No registration
for attendance, so that count is unclear. Format is four days and only mornings, 9am to noon US Central Time, so about half of
the usual time...more school than workshop, mabe. Might be YouTube live stream, or other. Want more social interaction, the
"coffee break"-like time. Expect to transfer earlier, in-person, registration to the online ones.
Zach, helped with a Zoom workshop recently, instead of asking for contributions, they asked people directly "please give a talk"---Steve
is doing this. Yosef Z could give a NL wave equation tutorial...hands-on would be one hour. Frank Loeffler talking about Neutron Stars.
See above link for others.
Attendees? Steve will add a registration to capture attendees. Will advertise on ETK outlets. We all should encourage students and colleagues to attend.
Zach, tutorials hosted on Jupyter Hub at NCSA? Steve, no it will not work. Have one method at LSU that Steve knows will work and
another that may. Zach, wonders about a docker image with compiled toolkit on it. Steve, should have that because tutorial server
uses image. For other tutorials we have trimmed the thorn list to match the par file. Roland, it was Perl and also on Python. But it
still makes sense to have it pre-compiled.
Steve, switch from OpenMPI to MPIch because running multiple images on the same host is problematic with OpenMPI.
Yosef, used VirtualBox for his class he taught, and that can also be difficult.
* test.cactuscode.orghttps://github.com/EinsteinToolkit/cactuscode.org/issues?q=is:open+is:issue… [SB]
Ongoing effort away from our supporting another server, want to move to GitHub. Steve needs help for this...and it is important to do.
Student with summer project for web development would be nice. Jekyll is a static web page generation tool, powerful, with built-in loops,
not fully dynamic, but might update/change from time to time. Ian H was using these features to do news feeds and table of contents.
GH pages can do some of this functionality. Liquid macro language.
* gallery testing: https://bitbucket.org/einsteintoolkit/tickets/issues?kind=task&milestone=ET… [RH]
** BNS: Shawn Rosofsky (UIUC). In progress, needs to be uploaded, some differences post-merger.
Roland has to poke Shawn again.
** BBH: Peter Schaffarczyk (Kiel), Atul Kedia (Notre Dame) In progress.
Atul, Contacted Wardell for the movies, and Bentivigna for the plots.
** Poisson gallery example, Bill was able to run Visit via a container from GitHub, on data using the psi.session file and plots (3D with slice of Psi
function) are the same as before.
* ET release manager search [RH]
Manager should also have assistants. Is a bit of work, especially reminders to others to finish their work, reviews, etc.
* reviewers for proposed inclusions [RH]
Three inclusions below that are left for inclusion. Want into the Toolkit so that they are downloaded sooner rather than later, so that users
could try this. Look to see if the Thorn breaks anything, and if not, let us include soon. This is different than a positive review...it not a
positive review we can always remove it from the list.
** ReadInterpolate https://github.com/rhaas80/ReadInterpolate.git thorn ticket 2416 https://bitbucket.org/einsteintoolkit/tickets/issues/2416 no reviewer [RH]
** Power https://git.ncsa.illinois.edu/elihu/Gravitational_Waveform_Extractor ticket 2413 https://bitbucket.org/einsteintoolkit/tickets/issues/2413 Zach is reviewer waveform extractor (or its post-summer version) [RH]
Zach, Power Package, not a thorn per se, more of a tool. He would like it as an option and not a separate download. Is really helps the integration of the
toolkit, takes output and makes is useful for GW science. Need to tell user how to access Power. Can this be added to Jenkins? Doc test? Written in Python.
Roland, Anything that is testable by a shell script can be tested by Jenkins. Roland, will try to make a release repo without the sample data (100 GB of so), and
find a way for it to be downloaded. It will be part of the toolkit, like GetComponents. Will be in the Thornlist and it will have statements to download the Power
script. Would probably put into bin directory as link.
** NRPYPn https://nbviewer.jupyter.org/github/zachetienne/nrpytutorial/blob/master/NR… ticket 2417 https://bitbucket.org/einsteintoolkit/tickets/issues/2417 [RH]
* gcc-10/ gfortran-10 [RH]
Roland, Needs review! Peter D has reviewed some parts. Will push into Master some of this. Backporting is a little more touchy, and needs separate review
by someone other than Roland. This is affecting people, hearing from users that this hits them.
** have patches to make all tests run with gcc-10/gfortran-10 with no more failures directly tracked down to gcc-10/gfortran-10 documented changes
** https://bitbucket.org/einsteintoolkit/tickets/issues/2431/review-changes-re…
** should be backported along with other fixes
** urgently needs a reviewer
unanswered question on mailing list
https://www.einsteintoolkit.org/tools/unanswered.php
open tickets sorted by update time
https://bitbucket.org/einsteintoolkit/tickets/issues?status=open&status=new…
Ticket 2436, Using Comet for production run is having some problems, but with System Topology helps with this. Peter D's example that all
16 threads ran on the same core. Roland, a caveat hwloc is important, and hwloc 1 and hwloc 2 are not API compatible. Not without risks
but worth the effort.
tickets ready for review
https://bitbucket.org/einsteintoolkit/tickets/issues?status=new&status=open…
Peter, Other items?
None
--
=====================================
William Gabella
Research Assistant Professor
Department of Physics and Astronomy
Vanderbilt University
Nashville, TN USA
b.gabella(a)vanderbilt.edu<mailto:b.gabella@vanderbilt.edu>
(o) 615-343-2713
Hello,
Please consider joining the weekly Einstein Toolkit phone call at
9:00 am US central time on Thursdays. For details on how to connect
and what agenda items are to be discussed, use the link below.
https://docs.einsteintoolkit.org/et-docs/Main_Page#Weekly_Users_Call
--The Maintainers