Shock Finding

New in version 2.1.

Source: Grid_FindShocks.C

Shock finding is accomplished using one of four methods. The primary method uses a coordinate-unsplit temperature jump (method 1), as described in Skillman et. al. 2008 with the exception that instead of searching across multiple grids for the pre- and post-shock cells, we terminate the search at the edge of the ghost zones within each grid.

Shock finding is controlled by the ShockMethod parameter, which can take the following values:

0 - Off

1 - Unsplit Temperature Jumps

2 - Dimensionally Split Temperature Jumps

3 - Unsplit Velocity Jumps

4 - Dimensionally Split Velocity Jumps

When ShockMethod nonzero, this will create a “Mach” field in the output files.

Note: Method 1 has been used the most by the developer, and therefore is the primary method. Method 2 has been tested quite a bit, but the downsides of using a dimensionally split method are outlined in the above paper. Methods 3 and 4 are more experimental and will run, but results may vary.

Additional Shock Finding Parameters:

ShockTemperatureFloor - When calculating the mach number using temperature jumps, set the temperature floor in the calculation to this value.

StorePreShockFields - Optionally store the Pre-shock Density and Temperature during data output.