<html><head></head><body><div>
<meta charset="utf-8">
<div id="compose-body-wrapper" dir="auto"><div dir="auto"><div dir="auto">Hi Erik,</div><div dir="auto"><br></div><div dir="auto">Sorry for the delay, I have written below the answers.</div><div dir="auto"><br></div><div dir="auto">1. What is the size of the simulation domain?</div><div dir="auto"><br></div><div dir="auto">Simulation domain along x and y direction is [-1000,1000] and along z direction is [0,1000].</div><div dir="auto"><br></div><div dir="auto">2. At which radius do you evaluate the angular momentum?</div><div dir="auto"><br></div><div dir="auto">There are two extraction radii, 200 and 300. The plot that I sent was for 300.</div><div dir="auto"><br></div><div dir="auto">3. What is the spatial resolution there?</div><div dir="auto"><br></div><div dir="auto">Spatial resolution is 5.</div><div dir="auto"><br></div><div dir="auto">4. What is the time resolution there?</div><div dir="auto"><br></div><div dir="auto">The time resolution is 0.75.</div><div dir="auto"><br></div><div dir="auto">5. Is there any matter nearby?</div><div dir="auto"><br></div><div dir="auto">There is no matter during the inspiral. </div><div dir="auto"><br></div><div dir="auto">6. What happens if you change the evaluation radius? What if you make it</div><div dir="auto">smaller or larger?</div><div dir="auto"><br></div><div dir="auto">For smaller radius the amplitude of oscillations is comparatively smaller.</div><div dir="auto"><br></div><div dir="auto">7. How close is this to the outer boundary?</div><div dir="auto"><br></div><div dir="auto">The extraction radius is far from the outer boundary.</div><div dir="auto"><br></div><div dir="auto">8. Do you have a visualisation of the geometry (e.g. the imaginary part</div><div dir="auto">of the Weyl scalar Psi_2)?</div><div dir="auto"><br></div><div dir="auto">I am afraid that I do not have Psi_2 in the output. If necessary then I can output it.</div><div dir="auto"><br></div><div dir="auto">9. In which direction does the angular momentum point?</div><div dir="auto">Does this change over time?</div><div dir="auto"><br></div><div dir="auto">The plot I sent was for the angular momentum along z direction.</div><div dir="auto"><br></div><div dir="auto">10. What is the magnitude of the angular momentum? How does this compare</div><div dir="auto">to the mass?</div><div dir="auto"><br></div><div dir="auto">The initial magnitude of the angular momentum is 7.3 and total gravitational mass of the system is 2.75 M_Sun.</div><div dir="auto"><br></div><div dir="auto">I am leaving below some information from the parameter file that maybe relevant.</div><div dir="auto"><br></div><div dir="auto">CoordBase::xmin = -1000</div><div dir="auto">CoordBase::xmax = 1000</div><div dir="auto">CoordBase::ymin = -1000</div><div dir="auto">CoordBase::ymax = 1000</div><div dir="auto">CoordBase::zmin = 0</div><div dir="auto">CoordBase::zmax = 1000 </div><div dir="auto">CoordBase::spacing = "numcells"</div><div dir="auto">CoordBase::ncells_x = 200</div><div dir="auto">CoordBase::ncells_y = 200</div><div dir="auto">CoordBase::ncells_z = 100</div><div dir="auto"><br></div><div dir="auto">Carpet::max_refinement_levels = 8</div><div dir="auto">Carpet::time_refinement_factors = "[1,1,2,4,8,16,32,64,128,256,512]"</div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">CarpetRegrid2::num_centres = 3</div><div dir="auto">CarpetRegrid2::regrid_every = 128</div><div dir="auto">CarpetRegrid2::snap_to_coarse = "yes"</div><div dir="auto">CarpetRegrid2::freeze_unaligned_levels = "yes"</div><div dir="auto">CarpetRegrid2::freeze_unaligned_parent_levels = "yes"</div><div dir="auto"><br></div><div dir="auto"># ----------- Region 1 --------------------</div><div dir="auto">CarpetRegrid2::active_1 = "yes"</div><div dir="auto">CarpetRegrid2::num_levels_1 = 7</div><div dir="auto">CarpetRegrid2::position_x_1 = 15.2365094091</div><div dir="auto">CarpetRegrid2::position_y_1 = 0.0</div><div dir="auto">CarpetRegrid2::radius_1[1] = 300.0</div><div dir="auto">CarpetRegrid2::radius_1[2] = 160.0</div><div dir="auto">CarpetRegrid2::radius_1[3] = 80.0</div><div dir="auto">CarpetRegrid2::radius_1[4] = 40.0</div><div dir="auto">CarpetRegrid2::radius_1[5] = 20.0</div><div dir="auto">CarpetRegrid2::radius_1[6] = 10.0</div><div dir="auto"><br></div><div dir="auto"># ----------- Region 2 --------------------</div><div dir="auto">CarpetRegrid2::active_2 = "yes"</div><div dir="auto">CarpetRegrid2::num_levels_2 = 7</div><div dir="auto">CarpetRegrid2::position_x_2 = -15.2365094091</div><div dir="auto">CarpetRegrid2::position_y_2 = -0.0</div><div dir="auto">CarpetRegrid2::radius_2[1] = 300.0</div><div dir="auto">CarpetRegrid2::radius_2[2] = 160.0</div><div dir="auto">CarpetRegrid2::radius_2[3] = 80.0</div><div dir="auto">CarpetRegrid2::radius_2[4] = 40.0</div><div dir="auto">CarpetRegrid2::radius_2[5] = 20.0</div><div dir="auto">CarpetRegrid2::radius_2[6] = 10.0 </div><div dir="auto"><br></div><div dir="auto"># ----------- Region 3 --------------------</div><div dir="auto">CarpetRegrid2::active_3 = "yes"</div><div dir="auto">CarpetRegrid2::num_levels_3 = 3</div><div dir="auto">CarpetRegrid2::position_x_3 = 0 </div><div dir="auto">CarpetRegrid2::radius_3[1] = 320.0</div><div dir="auto">CarpetRegrid2::radius_3[2] = 160.0</div><div dir="auto">CarpetRegrid2::radius_3[3] = 80.0</div><div dir="auto">CarpetRegrid2::radius_3[4] = 40.0</div><div dir="auto">CarpetRegrid2::radius_x_3[5] = 30.0</div><div dir="auto">CarpetRegrid2::radius_y_3[5] = 30.0</div><div dir="auto">CarpetRegrid2::radius_z_3[5] = 25.0</div><div dir="auto">CarpetRegrid2::radius_3[6] = 20.0</div><div dir="auto">CarpetRegrid2::radius_3[7] = 5.0</div><div dir="auto"># -----------------------------------------</div><div dir="auto"><br></div><div dir="auto">Thanks for your time. I would look forward to hearing from you.</div></div><div dir="auto"><br></div><div dir="auto" id="tmjah_g_1299"><div><div>Thanks,</div><div>Sudipta</div><div><br></div><div>Sudipta Hensh</div><div>JSPS Postdoctoral Researcher</div><div>Pronouns: He/Him</div><div>Graduate School of Science</div><div>Osaka University</div><div>1-1 Machikaneyama-cho, Toyonaka</div><div>Osaka 560-0043, Japan</div></div></div><br></div><div class="replyHeader" dir="auto">Erik Schnetter wrote:</div><br><br><div><blockquote cite="mid:CADKQjjebmiND1uJZ7u8HN4a=9qyw3xV++CJGfwUkMEhLZuzrYQ@mail.gmail.com" type="cite" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Sudipta</div><div><br></div><div>What is the size of the simulation domain?</div><div>At which radius do you evaluate the angular momentum?</div><div>What is the spatial resolution there?</div><div>What is the time resolution there?</div><div>Is there any matter nearby?</div><div>What happens if you change the evaluation radius? What if you make it</div><div>smaller or larger?</div><div>How close is this to the outer boundary?</div><div>Do you have a visualisation of the geometry (e.g. the imaginary part</div><div>of the Weyl scalar Psi_2)?</div><div>In which direction does the angular momentum point?</div><div>Does this change over time?</div><div>What is the magnitude of the angular momentum? How does this compare</div><div>to the mass?</div><div><br></div><div>-erik</div><div><br></div><div><br></div><div>On Tue, Mar 19, 2024 at 9:35 AM Sudipta Hensh <hensh@astro-osaka.jp> wrote:</div><br><blockquote class="gmail_quote" type="cite" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br><div> Dear all,</div><br><div> I hope this email finds you well.</div><br><div> I am performing simulations of BNS mergers and utilizing the QuasiLocalMeasures thorn to plot the total angular momentum of the spacetime. However, I have noticed high-frequency oscillations in the plotted data, as shown in the attached plot. I was wondering if we should expect such oscillations. I would also like to know the reason behind such oscillations. I would appreciate if you kindly help me to understand this issue.</div><br><div> I am looking forward to hearing from you. Many thanks in advance!</div><br><div> Regards,</div><div> Sudipta Hensh</div><div> JSPS Postdoctoral Researcher</div><div> Pronouns: He/Him</div><div> Graduate School of Science</div><div> Osaka University</div><div> 1-1 Machikaneyama-cho, Toyonaka</div><div> Osaka 560-0043, Japan</div><div> _______________________________________________</div><div> Users mailing list</div><div> Users@einsteintoolkit.org</div><br></blockquote><div> http://lists.einsteintoolkit.org/mailman/listinfo/users</div><div><br></div><div><br></div><div><br></div><div><div>-- </div><div>Erik Schnetter <schnetter@gmail.com></div><div>http://www.perimeterinstitute.ca/personal/eschnetter/</div><br><br></div></blockquote></div>
</div></body></html>