Opened 2 years ago

Last modified 11 months ago

#1974 assigned enhancement

Add WVUThorns_Diagnostics to the ET

Reported by: zachetie@… Owned by: zachetie
Priority: unset Milestone: ET_2018_08
Component: Other Version: development version
Keywords: Cc:


I have open-sourced a very useful set of diagnostic routines for GRMHD, GRHD, and even vacuum simulations.

The thorns are available here:

Here's a brief description of each thorn:
Brief description of included thorns

Seed_Magnetic_Fields-modified Extended Seed_Magnetic_Fields thorn for binary neutron stars. Supercedes Seed_Magnetic_Fields thorn.

Meudon_Bin_NS-modified Modifications to Meudon BNS initial data thorn to disable the overwriting of initial lapse/shift, which acts to significantly reduce coordinate eccentricity. Supercedes Meudon_Bin_NS thorn.

VolumeIntegrals_GRMHD Nice GRMHD volume integration thorn, currently depends on IllinoisGRMHD and Carpet. Performs volume integrals on arbitrary "Swiss-cheese"-like topologies, and even interoperates with Carpet to track NS centers of mass.

VolumeIntegrals_vacuum Nice GRMHD volume integration thorn, currently depends on ML_BSSN. There is a bit of code duplication and duplicated functionality between VI_GRMHD and VI_vacuum, to ensure that VI_vacuum can be used without enabling a GRMHD code. There is probably a better way of doing this, but I haven't had the time to think deeply about this.

particle_tracerET Solves the ODE \partial_t xi = vi for typically thousands of tracer particles, using an RK4 integration atop the current timestepping. E.g., one RK4 substep in the particle integration might occur every 16 RK4 substeps in the GRMHD evolution. These tracer particle positions are quite useful for visualizing magnetic field lines in a consistent way from frame-to-frame in a movie (recall that in the GRMHD approximation, the magnetic field lines stay attached the fluid elements they thread Flux Freezing?). Note that the velocity must be consistent with the velocity appearing in the GRMHD induction equation. This thorn reads in the HydroBase vel[] vector gridfunction, which assumes the Valencia formalism, and converts it into the induction equation velocity.

smallbPoynET Computes bi, b2, and three spatial components of Poynting flux. It also computes (-1-u0), which is useful for tracking unbound matter.

Attachments (0)

Change History (5)

comment:1 Changed 20 months ago by Frank Löffler


comment:2 Changed 20 months ago by Frank Löffler


Can you see in how far the changes described in maybe do the same as what you intended here?

comment:3 Changed 12 months ago by Ian Hinder

Milestone: ET_2018_02

comment:4 Changed 12 months ago by Ian Hinder

Milestone: ET_2018_02ET_2018_08

comment:5 Changed 11 months ago by Steven R. Brandt

Owner: set to zachetie
Status: newassigned

Modify Ticket

Change Properties
Set your email in Preferences
as assigned The owner will remain zachetie.
Next status will be 'review'.
as The resolution will be set.
to The owner will be changed from zachetie to the specified user.
The owner will be changed from zachetie to anonymous.

Add Comment

E-mail address and name can be saved in the Preferences.

Note: See TracTickets for help on using tickets.