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 = implementation defined

using	argument_tags = implementation defined

Static Public Member Functions
static std::tuple< int, evolution::dg::subcell::RdmpTciData >	apply (const Scalar< DataVector > &subcell_tilde_d, const Scalar< DataVector > &subcell_tilde_ye, const Scalar< DataVector > &subcell_tilde_tau, const tnsr::I< DataVector, 3, Frame::Inertial > &subcell_tilde_b, const Scalar< DataVector > &subcell_rest_mass_density, const Scalar< DataVector > &subcell_pressure, bool vars_needed_fixing, const Mesh< 3 > &dg_mesh, const Mesh< 3 > &subcell_mesh, const evolution::dg::subcell::RdmpTciData &past_rdmp_tci_data, const TciOptions &tci_options, const evolution::dg::subcell::SubcellOptions &subcell_options, double persson_exponent, bool need_rdmp_data_only)

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:

List of checks
Description	TCI status
if `min(tilde_d)` is less than `tci_options.minimum_rest_mass_density_times_lorentz_factor`, then we remain on FD.	`+1`
if `min(tilde_ye)` is less than `tci_options.minimum_rest_mass_density_times_lorentz_factor` times `tci_options.minimum_ye`, then we remain on FD.	`+2`
if `min(tilde_tau)` is less than `tci_options.minimum_tilde_tau`, then we remain on FD.	`+3`
if `grmhd::ValenciaDivClean::Tags::VariablesNeededFixing` is `true` and the maximum of rest mass density on FD grid is greater than `tci_options.atmosphere_density`, then we remain on FD.	`+4`
apply the Persson TCI to $\tilde{D}$ , ${\tilde{Y}}_{e}$ , and pressure if the maximum of rest mass density on FD grid is greater than `tci_options.atmosphere_density`.	`+5`
apply the Persson TCI to ${\tilde{Y}}_{e}$ if the maximum of rest mass density on FD grid is greater than `tci_options.atmosphere_density`.	`+6`
apply the Persson TCI to pressure if the maximum of rest mass density on FD grid is greater than `tci_options.atmosphere_density`.	`+7`
apply the RDMP TCI to `TildeD`	`+8`
apply the RDMP TCI to `TildeYe`	`+9`
apply the RDMP TCI to `TildeTau`	`+10`
apply the RDMP TCI to `TildeB`	`+11`
apply the Persson TCI to the magnitude of ${\tilde{B}}^{n + 1}$ if its magnitude is greater than `tci_options.magnetic_field_cutoff`.	`+12`

The second column of the table above denotes the value of an integer stored as the first element of the returned std::tuple, which indicates the particular kind of check that failed. For example, if the tenth check (RDMP TCI to TildeTau) fails and cell is marked as troubled, an integer with value +10 is stored in the first slot of the returned tuple. Note that this integer is marking only the first check to fail, since checks are done in a particular sequence as listed above. If all checks are passed and cell is not troubled, it is returned with the value 0.

Note: We adopt positive integers to mark TCI status from FD grid returned by TciOnFdGrid class. Negative integers are reserved for TCIs on DG grid; see TciOnDgGrid and its documentation.

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

src/Evolution/Systems/GrMhd/ValenciaDivClean/Subcell/TciOnFdGrid.hpp

Public Types

Static Public Member Functions

Detailed Description