Dear Roberto,
With respect to the question, I wanted to state that the issue with the GRHydro thorn which gave "retval" was overcomed by increasing resolution to dx=0.375 as you have suggested. But simulation stopped with NANs again at 1300 simulation time. Just as in the previous simulation, AHFinderDirect has searched for horizon after finding it for the first time at 825 but could not find again. I do not know how this problem could be solved.
Output and error file are attached.
Best regards,
Beyhan.
On Fri, Apr 17, 2020 at 6:45 PM Beyhan Karakaş beyhannkarakas@gmail.com wrote:
Dear Roberto,
After increasing resolution to dx=0.375 there has not been any error in the merger and afterwards so far (around 1077 simulation time) But AHFinderDirect has been searching for horizon after finding it at 825. Output file is attached. Doesn't it have to find the same horizon and how it is changing during the evolution?
I am looking forward to hearing from you.
Best regards,
Beyhan.
On Tue, Apr 14, 2020 at 9:10 PM Beyhan Karakaş beyhannkarakas@gmail.com wrote:
Dear Roberto,
Thank you very much for your consideration and reply. When I change the resolution to dx=0.375 it requires 16 gb of ram and simulation becomes 9 times slower and it would take 2.5 days to reach the 1000 simulation time. I will try to do this and , if possible, inform you about the result.
Best regards,
Beyhan.
On Tue, Apr 14, 2020 at 8:32 PM Roberto De Pietri < roberto.depietri@unipr.it> wrote:
Dear Beyjann
From the output, it seems that the two stars just merged and then a prompt collapse to Black Hole. An Apparent Horizon was found.
INFO (AHFinderDirect): proc 0: searching for horizon 1/1 INFO (AHFinderDirect): AH 1/1: r=2.18345 at (-0.000000,0.000000,0.000000) INFO (AHFinderDirect): AH 1/1: area=376.8755139 m_irreducible=2.738192857 5952 558.000 | 129.2759814 | 0.0003459 | 0.0000439 | 8.9127188 | 8.4903683 | 151.5039356 | 301.5059049
Dear Beyhan;
Quite likely just resolution.
But at the next iteration, AH finder was not able to find it again, and indeed the interior of the black hole was not excised and conservative to primitive failed close to the singularity.
Resolution of dx=0.75 it is not enough for simulation, and leas dx=0.375 is needed for having running simulation when black hole form and at least 0.25 to have some useful data.
Roberto
On 14 Apr 2020, at 17:22, Beyhan Karakaş beyhannkarakas@gmail.com wrote:
Dear to whom it may concern,
I am using Mayer release on Fedora 31. I have made some changes in the parameter file that was created by De Pietri, Feo, Maione and Loeffler in 2016 (which could be downloaded from http://www.fis.unipr.it/numrel/Research/BNS2016.html https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.fis.unipr.it%2Fnumrel%2FResearch%2FBNS2016.html&data=02%7C01%7Croberto.depietri%40unipr.it%7C9856f8be4a614c36092808d7e087c54c%7Cbb064bc5b7a841ecbabed7beb3faeb1c%7C0%7C0%7C637224746075422862&sdata=dwMfDo1WyBsGfCHoyLC%2BxsgD8e7LywtPY%2BCajZwPheA%3D&reserved=0). Some of the Parma Groups' thorns were used in the parameter file added to the end of the mail.
Initial data for equal mass (M_1,2=1.7 M_sun) binary neutron stars was generated on LORENE using SLy EOS. Table of it was attached. Initial data could be downloaded from https://bitbucket.org/beyhank/bns https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fbitbucket.org%2Fbeyhank%2Fbns&data=02%7C01%7Croberto.depietri%40unipr.it%7C9856f8be4a614c36092808d7e087c54c%7Cbb064bc5b7a841ecbabed7beb3faeb1c%7C0%7C0%7C637224746075432863&sdata=%2FihD6LLxM2FTPBZpMQUevJecdN1yJf6Y0%2B%2FEN7xnJv4%3D&reserved=0 .
GRHydro thorn starts to give retval after 504 simulation time and stops the simulation at 756 with NANs (BH formed at 558). I suppose the error persists either "EinsteinEvolve/GRHydro" or "PRThorns/GRHydro" is used.
I do not know whether this error is resulted from low resolution or something else. Same setting does not seem to give any error for the configuration M_1,2=1.2 M_sun.
Output, error and parameter files are attached.
Best regards,
Beyhan.
# The Parma Thorns !TARGET = $ARR !TYPE = git !URL = https://bitbucket.org/GravityPR/prthorns.git https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fbitbucket.org%2FGravityPR%2Fprthorns.git&data=02%7C01%7Croberto.depietri%40unipr.it%7C9856f8be4a614c36092808d7e087c54c%7Cbb064bc5b7a841ecbabed7beb3faeb1c%7C0%7C0%7C637224746075432863&sdata=PeD%2BZDceaHgO2mnDg8kKI2dmfc34Y7YmB%2BqOG383Gco%3D&reserved=0 !REPO_PATH= $2 !CHECKOUT = PRThorns/TempPool PRThorns/Meudon_Bin_NS PRThorns/Hydro_RNSID PRThorns/GRHydro PRThorns/GRHydro_Analysis PRThorns/GRHydro_InitData <tnisan12.par><tnisan12.err><tnisan12.out><eos_sly4_pp.d> _______________________________________________ Users mailing list Users@einsteintoolkit.org
https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.einst...