The fluxes of the conservative variables in the M1 scheme. More...

#include <Fluxes.hpp>

Public Types
using	return_tags = implementation defined

using	argument_tags = implementation defined

Static Public Member Functions
static void	apply (const gsl::not_null< typename ::Tags::Flux< Tags::TildeE< Frame::Inertial, NeutrinoSpecies >, tmpl::size_t< 3 >, Frame::Inertial >::type * >... tilde_e_flux, const gsl::not_null< typename ::Tags::Flux< Tags::TildeS< Frame::Inertial, NeutrinoSpecies >, tmpl::size_t< 3 >, Frame::Inertial >::type * >... tilde_s_flux, const typename Tags::TildeE< Frame::Inertial, NeutrinoSpecies >::type &... tilde_e, const typename Tags::TildeS< Frame::Inertial, NeutrinoSpecies >::type &... tilde_s, const typename Tags::TildeP< Frame::Inertial, NeutrinoSpecies >::type &... tilde_p, const Scalar< DataVector > &lapse, const tnsr::I< DataVector, 3, Frame::Inertial > &shift, const tnsr::ii< DataVector, 3, Frame::Inertial > &spatial_metric, const tnsr::II< DataVector, 3, Frame::Inertial > &inv_spatial_metric)

Detailed Description

template<typename... NeutrinoSpecies>
struct RadiationTransport::M1Grey::ComputeFluxes< NeutrinoSpecies >

The fluxes of the conservative variables in the M1 scheme.

$\begin{aligned} F^{i} (\tilde{E}) = & α γ^{i j} {\tilde{S}}_{j} - β^{j} \tilde{E} \\ F^{i} ({\tilde{S}}_{j}) = & α {\tilde{P}}^{i k} γ_{k j} - β^{i} {\tilde{S}}_{j} \end{aligned}$

where the conserved variables $\tilde{E}$ , ${\tilde{S}}_{i}$ , are a generalized mass-energy density and momentum density. Furthermore, ${\tilde{P}}^{i j}$ is the pressure tensor density of the radiation field, $α$ is the lapse, $β^{i}$ is the shift, $γ_{i j}$ the 3-metric, and $γ^{i j}$ its inverse.

In the main function, we loop over all neutrino species, and then call the actual implementation of the fluxes.

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

src/Evolution/Systems/RadiationTransport/M1Grey/Fluxes.hpp

Public Types

Static Public Member Functions

Detailed Description