SpECTRE  v2021.10.04
grmhd::ValenciaDivClean::subcell::TciOnFdGrid Struct Reference

The troubled-cell indicator run on the FD grid to check if the corresponding DG solution is admissible. More...

#include <TciOnFdGrid.hpp>

Public Types

using return_tags = tmpl::list<>
 
using argument_tags = tmpl::list< evolution::dg::subcell::Tags::Inactive< grmhd::ValenciaDivClean::Tags::TildeD >, evolution::dg::subcell::Tags::Inactive< grmhd::ValenciaDivClean::Tags::TildeTau >, evolution::dg::subcell::Tags::Inactive< grmhd::ValenciaDivClean::Tags::TildeB<> >, grmhd::ValenciaDivClean::Tags::VariablesNeededFixing, domain::Tags::Mesh< 3 >, Tags::TciOptions >
 

Static Public Member Functions

static bool apply (const Scalar< DataVector > &tilde_d, const Scalar< DataVector > &tilde_tau, const tnsr::I< DataVector, 3, Frame::Inertial > &tilde_b, bool vars_needed_fixing, const Mesh< 3 > &dg_mesh, const TciOptions &tci_options, double persson_exponent)
 

Detailed Description

The troubled-cell indicator run on the FD grid to check if the corresponding DG solution is admissible.

The following checks are done in the order they are listed:

  • if grmhd::ValenciaDivClean::Tags::VariablesNeededFixing is true then we remain on FD. (Note: this could be relaxed in the future if we need to allow switching from FD to DG in the atmosphere and the current approach isn't working.)
  • if min(tilde_d) is less than tci_options.minimum_rest_mass_density_times_lorentz_factor or if min(tilde_tau) is less than tci_options.minimum_tilde_tau then the we remain on FD.
  • apply the Persson TCI to \(\tilde{D}\) and \(\tilde{\tau}\)

The documentation for this struct was generated from the following file: