SpECTRE  v2026.04.01
Loading...
Searching...
No Matches
grmhd::GhValenciaDivClean::fd Namespace Reference

Finite difference functionality for the coupled Generalized Harmonic and ValenciaDivClean equations. More...

Namespaces

namespace  OptionTags
 Option tags for finite difference solver.
namespace  Tags
 Tags for finite difference solver

Classes

struct  BoundaryConditionGhostData
 Computes finite difference ghost data for external boundary conditions. More...
struct  FilterOptions
 Filtering/dissipation options. More...
class  MonotonisedCentralPrim
 Monotonised central reconstruction on the GRMHD primitive variables (see fd::reconstruction::monotonised_central() for details) and unlimited 3rd order (degree 2 polynomial) reconstruction on the metric variables. More...
class  PositivityPreservingAdaptiveOrderPrim
 Positivity-preserving adaptive order reconstruction. See fd::reconstruction::positivity_preserving_adaptive_order() for details. The rest mass density, electron fraction, and the pressure are kept positive. Use unlimited 5th order (degree 4 polynomial) reconstruction on the metric variables. More...
class  Reconstructor
 The base class from which all reconstruction schemes must inherit. More...
class  Wcns5zPrim
 Fifth order weighted nonlinear compact scheme reconstruction using the Z oscillation indicator. See fd::reconstruction::wcns5z() for details. More...

Typedefs

using tags_list_for_reconstruct
using tags_list_for_reconstruct_fd_neighbor

Functions

template<typename System>
void spacetime_derivatives (gsl::not_null< Variables< db::wrap_tags_in<::Tags::deriv, typename System::gradients_tags, tmpl::size_t< 3 >, Frame::Inertial > > * > result, const Variables< typename System::variables_tag::tags_list > &volume_evolved_variables, const DirectionalIdMap< 3, evolution::dg::subcell::GhostData > &all_ghost_data, const size_t &deriv_order, const Mesh< 3 > &volume_mesh, const InverseJacobian< DataVector, 3, Frame::ElementLogical, Frame::Inertial > &cell_centered_logical_to_inertial_inv_jacobian)
 Compute partial derivatives of the spacetime variables \(g_{ab}\), \(\Phi_{iab}\), and \(\Pi_{ab}\).
template<typename NeighborVariables, typename FirstGhTag>
void fill_neighbor_spacetime_variables (const gsl::not_null< DirectionMap< 3, gsl::span< const double > > * > ghost_cell_spacetime_vars, const DirectionalIdMap< 3, evolution::dg::subcell::GhostData > &all_ghost_data, const size_t number_of_gh_components)
 Helper function that takes spacetime variable data from a DirectionalIdMap containing all neighbor data and copies them to ghost_cell_spacetime_vars, a separate DirectionMap.
template<typename VariableTags>
void spacetime_kreiss_oliger_filter (gsl::not_null< Variables< VariableTags > * > result, const Variables< VariableTags > &volume_evolved_variables, const DirectionalIdMap< 3, evolution::dg::subcell::GhostData > &all_ghost_data, const Mesh< 3 > &volume_mesh, size_t order, double epsilon)
 Apply a Kreiss-Oliger filter to \(g_{ab}\), \(\Phi_{iab}\), and \(\Pi_{ab}\).
template<typename SpacetimeTagsToReconstruct, typename PrimTagsForReconstruction, typename PrimsTags, typename SpacetimeAndConsTags, size_t ThermodynamicDim, typename HydroReconstructor, typename SpacetimeReconstructor, typename ComputeGrmhdSpacetimeVarsFromReconstructedSpacetimeTags, typename PrimsTagsSentByNeighbor>
void reconstruct_prims_work (gsl::not_null< std::array< Variables< tags_list_for_reconstruct >, 3 > * > vars_on_lower_face, gsl::not_null< std::array< Variables< tags_list_for_reconstruct >, 3 > * > vars_on_upper_face, const HydroReconstructor &hydro_reconstructor, const SpacetimeReconstructor &spacetime_reconstructor, const ComputeGrmhdSpacetimeVarsFromReconstructedSpacetimeTags &spacetime_vars_for_grmhd, const Variables< PrimsTags > &volume_prims, const Variables< SpacetimeAndConsTags > &volume_spacetime_and_cons_vars, const EquationsOfState::EquationOfState< true, ThermodynamicDim > &eos, const Element< 3 > &element, const DirectionalIdMap< 3, Variables< PrimsTagsSentByNeighbor > > &neighbor_data, const Mesh< 3 > &subcell_mesh, size_t ghost_zone_size, bool compute_conservatives, bool reconstruct_density_times_temperature, const VariableFixing::FixToAtmosphere< 3 > *fix_to_atmosphere)
 Reconstructs \(\rho, p, Wv^i, B^i\), \(\Phi\), and the spacetime metric, then computes the Lorentz factor, upper spatial velocity, specific internal energy, and the conserved variables. All results are written into vars_on_lower_face and vars_on_upper_face.
template<typename SpacetimeTagsToReconstruct, typename PrimTagsForReconstruction, typename PrimsTagsSentByNeighbor, typename PrimsTags, size_t ThermodynamicDim, typename LowerHydroReconstructor, typename LowerSpacetimeReconstructor, typename UpperHydroReconstructor, typename UpperSpacetimeReconstructor, typename ComputeGrmhdSpacetimeVarsFromReconstructedSpacetimeTags>
void reconstruct_fd_neighbor_work (gsl::not_null< Variables< tags_list_for_reconstruct_fd_neighbor > * > vars_on_face, const LowerHydroReconstructor &reconstruct_lower_neighbor_hydro, const LowerSpacetimeReconstructor &reconstruct_lower_neighbor_spacetime, const UpperHydroReconstructor &reconstruct_upper_neighbor_hydro, const UpperSpacetimeReconstructor &reconstruct_upper_neighbor_spacetime, const ComputeGrmhdSpacetimeVarsFromReconstructedSpacetimeTags &spacetime_vars_for_grmhd, const Variables< PrimsTags > &subcell_volume_prims, const Variables< grmhd::GhValenciaDivClean::Tags::spacetime_reconstruction_tags > &subcell_volume_spacetime_vars, const EquationsOfState::EquationOfState< true, ThermodynamicDim > &eos, const Element< 3 > &element, const DirectionalIdMap< 3, evolution::dg::subcell::GhostData > &ghost_data, const Mesh< 3 > &subcell_mesh, const Direction< 3 > &direction_to_reconstruct, size_t ghost_zone_size, bool compute_conservatives, bool reconstruct_density_times_temperature, const VariableFixing::FixToAtmosphere< 3 > *fix_to_atmosphere)
 Reconstructs \(\rho, p, Wv^i, B^i\), \(\Phi\), the spacetime metric, \(\Phi_{iab}\), and \(\Pi_{ab}\), then computes the Lorentz factor, upper spatial velocity, specific internal energy, and the conserved variables. All results are written into vars_on_face.
template<typename System>
void register_derived_with_charm ()

Detailed Description

Finite difference functionality for the coupled Generalized Harmonic and ValenciaDivClean equations.

Typedef Documentation

◆ tags_list_for_reconstruct

using grmhd::GhValenciaDivClean::fd::tags_list_for_reconstruct
Initial value:
tmpl::push_front<grmhd::ValenciaDivClean::fd::tags_list_for_reconstruct,
gr::Tags::SpacetimeMetric<DataVector, 3>,
gh::Tags::Pi<DataVector, 3>, gh::Tags::Phi<DataVector, 3>>
grmhd::ValenciaDivClean::fd::tags_list_for_reconstruct
tmpl::list< grmhd::ValenciaDivClean::Tags::TildeD, grmhd::ValenciaDivClean::Tags::TildeYe, grmhd::ValenciaDivClean::Tags::TildeTau, grmhd::ValenciaDivClean::Tags::TildeS< Frame::Inertial >, grmhd::ValenciaDivClean::Tags::TildeB< Frame::Inertial >, grmhd::ValenciaDivClean::Tags::TildePhi, hydro::Tags::RestMassDensity< DataVector >, hydro::Tags::ElectronFraction< DataVector >, hydro::Tags::SpecificInternalEnergy< DataVector >, hydro::Tags::SpatialVelocity< DataVector, 3 >, hydro::Tags::MagneticField< DataVector, 3 >, hydro::Tags::DivergenceCleaningField< DataVector >, hydro::Tags::LorentzFactor< DataVector >, hydro::Tags::Pressure< DataVector >, hydro::Tags::Temperature< DataVector >, hydro::Tags::LorentzFactorTimesSpatialVelocity< DataVector, 3 >, ::Tags::Flux< grmhd::ValenciaDivClean::Tags::TildeD, tmpl::size_t< 3 >, Frame::Inertial >, ::Tags::Flux< grmhd::ValenciaDivClean::Tags::TildeYe, tmpl::size_t< 3 >, Frame::Inertial >, ::Tags::Flux< grmhd::ValenciaDivClean::Tags::TildeTau, tmpl::size_t< 3 >, Frame::Inertial >, ::Tags::Flux< grmhd::ValenciaDivClean::Tags::TildeS< Frame::Inertial >, tmpl::size_t< 3 >, Frame::Inertial >, ::Tags::Flux< grmhd::ValenciaDivClean::Tags::TildeB< Frame::Inertial >, tmpl::size_t< 3 >, Frame::Inertial >, ::Tags::Flux< grmhd::ValenciaDivClean::Tags::TildePhi, tmpl::size_t< 3 >, Frame::Inertial >, gr::Tags::Lapse< DataVector >, gr::Tags::Shift< DataVector, 3 >, hydro::Tags::SpatialVelocityOneForm< DataVector, 3, Frame::Inertial >, gr::Tags::SpatialMetric< DataVector, 3 >, gr::Tags::SqrtDetSpatialMetric< DataVector >, gr::Tags::InverseSpatialMetric< DataVector, 3 >, evolution::dg::Actions::detail::NormalVector< 3 > > tags_list_for_reconstruct
The list of tags used by reconstruct_prims_work and reconstruct_fd_neighbor_work
Definition ReconstructWork.hpp:50
gh::Tags::Phi
Auxiliary variable which is analytically the spatial derivative of the spacetime metric.
Definition Tags.hpp:38
gh::Tags::Pi
Conjugate momentum to the spacetime metric.
Definition Tags.hpp:27
gr::Tags::SpacetimeMetric
Definition Tags.hpp:16

◆ tags_list_for_reconstruct_fd_neighbor

using grmhd::GhValenciaDivClean::fd::tags_list_for_reconstruct_fd_neighbor
Initial value:
tmpl::append<
tmpl::front<detail::tags_list_for_reconstruct_split_lapse>,
tmpl::push_front<tmpl::back<detail::tags_list_for_reconstruct_split_lapse>,
gh::Tags::ConstraintGamma1, gh::Tags::ConstraintGamma2,
gr::Tags::Lapse<DataVector>>>
gh::Tags::ConstraintGamma1
Constraint dammping parameter for the generalized harmonic system (cf. lindblom2005qh).
Definition ConstraintDampingTags.hpp:62
gh::Tags::ConstraintGamma2
Constraint dammping parameter for the generalized harmonic system (cf. lindblom2005qh).
Definition ConstraintDampingTags.hpp:70
gr::Tags::Lapse
Definition Tags.hpp:65

Function Documentation

◆ fill_neighbor_spacetime_variables()

template<typename NeighborVariables, typename FirstGhTag>
void grmhd::GhValenciaDivClean::fd::fill_neighbor_spacetime_variables ( const gsl::not_null< DirectionMap< 3, gsl::span< const double > > * > ghost_cell_spacetime_vars,
const DirectionalIdMap< 3, evolution::dg::subcell::GhostData > & all_ghost_data,
const size_t number_of_gh_components )

Helper function that takes spacetime variable data from a DirectionalIdMap containing all neighbor data and copies them to ghost_cell_spacetime_vars, a separate DirectionMap.

Template Parameters
NeighborVariablesa Variables type containing all tags stored in the neighbor data
FirstGhTagthe first spacetime variables tag in NeighborVariables. Note: it is assumed that the spacetime variables are consecutively stored in the Variables
Parameters
ghost_cell_spacetime_varsa DirectionMap to be filled with spans storing the spacetime data
all_ghost_dataa DirectionalIdMap containing all neighbor data
number_of_gh_componentsthe number of independent components for all spacetime variables in NeighborVariables

◆ reconstruct_fd_neighbor_work()

template<typename SpacetimeTagsToReconstruct, typename PrimTagsForReconstruction, typename PrimsTagsSentByNeighbor, typename PrimsTags, size_t ThermodynamicDim, typename LowerHydroReconstructor, typename LowerSpacetimeReconstructor, typename UpperHydroReconstructor, typename UpperSpacetimeReconstructor, typename ComputeGrmhdSpacetimeVarsFromReconstructedSpacetimeTags>
void grmhd::GhValenciaDivClean::fd::reconstruct_fd_neighbor_work ( gsl::not_null< Variables< tags_list_for_reconstruct_fd_neighbor > * > vars_on_face,
const LowerHydroReconstructor & reconstruct_lower_neighbor_hydro,
const LowerSpacetimeReconstructor & reconstruct_lower_neighbor_spacetime,
const UpperHydroReconstructor & reconstruct_upper_neighbor_hydro,
const UpperSpacetimeReconstructor & reconstruct_upper_neighbor_spacetime,
const ComputeGrmhdSpacetimeVarsFromReconstructedSpacetimeTags & spacetime_vars_for_grmhd,
const Variables< PrimsTags > & subcell_volume_prims,
const Variables< grmhd::GhValenciaDivClean::Tags::spacetime_reconstruction_tags > & subcell_volume_spacetime_vars,
const EquationsOfState::EquationOfState< true, ThermodynamicDim > & eos,
const Element< 3 > & element,
const DirectionalIdMap< 3, evolution::dg::subcell::GhostData > & ghost_data,
const Mesh< 3 > & subcell_mesh,
const Direction< 3 > & direction_to_reconstruct,
size_t ghost_zone_size,
bool compute_conservatives,
bool reconstruct_density_times_temperature,
const VariableFixing::FixToAtmosphere< 3 > * fix_to_atmosphere )

Reconstructs \(\rho, p, Wv^i, B^i\), \(\Phi\), the spacetime metric, \(\Phi_{iab}\), and \(\Pi_{ab}\), then computes the Lorentz factor, upper spatial velocity, specific internal energy, and the conserved variables. All results are written into vars_on_face.

This is used on DG elements to reconstruct their subcell neighbors' solution on the shared faces.

◆ spacetime_derivatives()

template<typename System>
void grmhd::GhValenciaDivClean::fd::spacetime_derivatives ( gsl::not_null< Variables< db::wrap_tags_in<::Tags::deriv, typename System::gradients_tags, tmpl::size_t< 3 >, Frame::Inertial > > * > result,
const Variables< typename System::variables_tag::tags_list > & volume_evolved_variables,
const DirectionalIdMap< 3, evolution::dg::subcell::GhostData > & all_ghost_data,
const size_t & deriv_order,
const Mesh< 3 > & volume_mesh,
const InverseJacobian< DataVector, 3, Frame::ElementLogical, Frame::Inertial > & cell_centered_logical_to_inertial_inv_jacobian )

Compute partial derivatives of the spacetime variables \(g_{ab}\), \(\Phi_{iab}\), and \(\Pi_{ab}\).

The derivatives are computed using FD of order deriv_order