[ET Trac] [Einstein Toolkit] #1974: Add WVUThorns_Diagnostics to the ET
Einstein Toolkit
trac-noreply at einsteintoolkit.org
Mon Oct 24 10:28:08 CDT 2016
#1974: Add WVUThorns_Diagnostics to the ET
--------------------------------+-------------------------------------------
Reporter: zachetie@… | Owner:
Type: enhancement | Status: new
Priority: unset | Milestone:
Component: Other | Version: development version
Keywords: |
--------------------------------+-------------------------------------------
I have open-sourced a very useful set of diagnostic routines for GRMHD,
GRHD, and even vacuum simulations.
The thorns are available here:
https://bitbucket.org/zach_etienne/wvuthorns_diagnostics
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 x^i = v^i 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 b^i^, b^2^, and three spatial components of
Poynting flux. It also computes (-1-u,,0,,), which is useful for tracking
unbound matter.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/1974>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
More information about the Trac
mailing list