[Users] "Spikes" in density for FLRWSolver single-mode simulation

[Cupp, Samuel D.]

Hello Hayley,
   Thanks for letting us know about this issue. Since this was not a problem in earlier versions of the Einstein Toolkit, could you check out older ET versions and find when this change appears? It would be helpful if we knew that this was introduced in the changes from ET_2020_05 to RT_2020_11, for example. Then, we can focus on how the thorns changed between those two releases.

Also, are you using the FRLWSolver in the Toolkit, or a newer in-house version?

   Samuel Cupp
   Postdoctoral Researcher
   Department of Physics
   University of Idaho
________________________________
From: Users <users-bounces at einsteintoolkit.org> on behalf of Hayley Macpherson <hayleyjmacpherson at gmail.com>
Sent: Thursday, November 9, 2023 3:28 PM
To: users at einsteintoolkit.org <users at einsteintoolkit.org>
Subject: [Users] "Spikes" in density for FLRWSolver single-mode simulation

Hi all,

I’ve noticed an issue with the single-mode setting of FLRWSolver (cosmological initial data thorn) when evolving with GRHydro.
This set up is three intersecting sine-waves in 3D, see the attached 1D slice for the relative density perturbation (delta = rho / rho.mean - 1) on the initial slice (first plot attached).

To do cosmological simulations which are dark-matter dominated we want something as close to P=0 as possible. GRHydro does not allow this so instead we choose a pressure such that P<<rho. Specifically, we choose a Polytropic EOS with poly_gamma=2 and set poly_k to be very small to achieve this. This has worked in the past, but recently I’ve noticed something strange in the perturbations themselves. I have previously used this setup back in 2017 and didn’t notice this issue, so I’m unsure if it’s something in a recent version of the ET or something different in the parameter file I’m setting.
Also just to note - I am aware that Eloisa’s code allows for P=0, however, having a small amount of pressure in the simulations is actually beneficial at late times when the structures are nonlinear (it stops the code from crashing too early and allows us to get further into the nonlinear regime). So for this reason I prefer to continue using GRHydro for now.

I’ve attached a parameter file which will reproduce the problem I’m seeing (which initialises the figure shown below).
The second plot I’ve attached shows the same simulation 5 iterations later (the 1D slice is 1/4 of the way through the domain, at index 8 in this example), where at the base of the sine wave you can see the “spike” starting to grow. This continues to grow throughout the simulation to be very large (the last plot is the same simulation after 400 iterations, for example). Also to note, this perturbation is small enough such that for the entire simulation it *should* remain very closely approximated by a sine wave.

I’ve pinpointed this issue to be sensitive to the value of poly_k that I choose for the simulations. Changing the value of poly_k changes the time in the simulation at which this occurs (it happens later on for larger values of poly_k), but then we reach a point at which the P<<rho condition we require is not satisfied anymore for large values of poly_k. I’ve chosen an extreme value of poly_k in the attached parameter file so that the problem happens early on and it might be easier to diagnose.
It seems like this is occurring at the minima and maxima of the sine wave, so I thought it might be something related to the fact the derivatives are very small at these points with a small poly_k such as we’re choosing.

I thought that instead of trying to find the “sweet spot” where that these spikes don’t occur *and* we can satisfy P<<rho, there might actually be a fix from someone who knows GRHydro much better than I do.

If anyone is able to help me I would very much appreciate it!

Best wishes,
Hayley


[cid:59880ece-9242-414c-8d61-d6789b11e8ce at prod.exchangelabs.com]
[cid:0a0007b2-4e60-4bfc-a0bf-87121af0aec0 at prod.exchangelabs.com]

[cid:e00bd1cc-55a9-41df-8c70-d68ac41191c5 at prod.exchangelabs.com]




----

Hayley Macpherson | NASA Einstein Fellow

Email: hjmacpherson at uchicago.edu<mailto:hjmacpherson at uchicago.edu>
Pronouns: she/her/hers

Office: ERC 479
Kavli Institute for Cosmological Physics
Eckhardt Research Center
5640 South Ellis Avenue
Chicago, IL, 60637

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20231116/92a420d2/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot 2023-11-09 at 13.11.29.png
Type: image/png
Size: 28991 bytes
Desc: Screenshot 2023-11-09 at 13.11.29.png
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20231116/92a420d2/attachment-0003.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot 2023-11-09 at 13.11.52.png
Type: image/png
Size: 29324 bytes
Desc: Screenshot 2023-11-09 at 13.11.52.png
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20231116/92a420d2/attachment-0004.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot 2023-11-09 at 13.17.32.png
Type: image/png
Size: 36297 bytes
Desc: Screenshot 2023-11-09 at 13.17.32.png
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20231116/92a420d2/attachment-0005.png>