Compute the curvature source terms for the flux-balanced grey M1 radiation transport. More...

#include <Sources.hpp>

Public Types
using	return_tags = tmpl::list< ::Tags::Source< Tags::TildeE< Frame::Inertial, NeutrinoSpecies > >..., ::Tags::Source< Tags::TildeS< Frame::Inertial, NeutrinoSpecies > >... >

using	argument_tags = tmpl::list< Tags::TildeE< Frame::Inertial, NeutrinoSpecies >..., Tags::TildeS< Frame::Inertial, NeutrinoSpecies >..., Tags::TildeP< Frame::Inertial, NeutrinoSpecies >..., gr::Tags::Lapse< DataVector >, ::Tags::deriv< gr::Tags::Lapse< DataVector >, tmpl::size_t< 3 >, Frame::Inertial >, ::Tags::deriv< gr::Tags::Shift< DataVector, 3 >, tmpl::size_t< 3 >, Frame::Inertial >, ::Tags::deriv< gr::Tags::SpatialMetric< DataVector, 3 >, tmpl::size_t< 3 >, Frame::Inertial >, gr::Tags::InverseSpatialMetric< DataVector, 3 >, gr::Tags::ExtrinsicCurvature< DataVector, 3 >, gr::Tags::SpatialMetric< DataVector, 3 >, Tags::GreyEmissivity< NeutrinoSpecies >..., Tags::GreyAbsorptionOpacity< NeutrinoSpecies >..., Tags::GreyScatteringOpacity< NeutrinoSpecies >..., Tags::TildeJ< NeutrinoSpecies >..., Tags::TildeHNormal< NeutrinoSpecies >..., Tags::TildeHSpatial< Frame::Inertial, NeutrinoSpecies >..., hydro::Tags::SpatialVelocity< DataVector, 3 >, hydro::Tags::LorentzFactor< DataVector >, gr::Tags::SqrtDetSpatialMetric< DataVector > >

Static Public Member Functions
static void	apply (const gsl::not_null< typename Tags::TildeE< Frame::Inertial, NeutrinoSpecies >::type * >... sources_tilde_e, const gsl::not_null< typename Tags::TildeS< Frame::Inertial, NeutrinoSpecies >::type * >... sources_tilde_s, 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 > &d_lapse, const tnsr::iJ< DataVector, 3 > &d_shift, const tnsr::ijj< DataVector, 3 > &d_spatial_metric, const tnsr::II< DataVector, 3 > &inv_spatial_metric, const tnsr::ii< DataVector, 3 > &extrinsic_curvature, const tnsr::ii< DataVector, 3 > &spatial_metric, const Scalar< DataVector > &emissivity, const Scalar< DataVector > &absorption_opacity, const Scalar< DataVector > &scattering_opacity, const typename Tags::TildeJ< NeutrinoSpecies >::type &... tilde_j, const typename Tags::TildeHNormal< NeutrinoSpecies >::type &... tilde_h_normal, const typename Tags::TildeHSpatial< Frame::Inertial, NeutrinoSpecies >::type &... tilde_h_spatial, const tnsr::I< DataVector, 3 > &spatial_velocity, const Scalar< DataVector > &lorentz, const Scalar< DataVector > &sqrt_det_spatial_metric)

Detailed Description

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

Compute the curvature source terms for the flux-balanced grey M1 radiation transport.

A flux-balanced system has the generic form:

\[ \partial_t U_i + \partial_m F^m(U_i) = S(U_i) \]

where \(F^a()\) denotes the flux of a conserved variable \(U_i\) and \(S()\) denotes the source term for the conserved variable.

For the grey M1 formalism (neglecting coupling to the fluid):

\begin{align*} S({\tilde E}) &= \alpha \tilde P^{ij} K_{ij} - \tilde S^i \partial_i \alpha,\\ S({\tilde S_i}) &= -\tilde E \partial_i \alpha + \tilde S_k \partial_i \beta^k + \frac{1}{2} \alpha \tilde P^{jk} \partial_i \gamma_{jk}, \end{align*}

where \({\tilde E}\), \({\tilde S_i}\), \({\tilde P}^{ij}\) are the densitized energy, momentum, and pressure tensor of the neutrinos/photons, \(K_{ij}\) is the extrinsic curvature, and \(\alpha\), \(\beta^i\), \(\gamma_{ij}\) are the lapse, shift and 3-metric.

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

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

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

Public Types

Static Public Member Functions

Detailed Description