RadiationTransport::M1Grey::BoundaryCorrections::Rusanov< tmpl::list< NeutrinoSpecies... > > Class Template Reference

A Rusanov/local Lax-Friedrichs Riemann solver. More...

#include <Rusanov.hpp>

Public Types
using	options = tmpl::list<>
using	dg_package_field_tags
using	dg_package_data_temporary_tags = tmpl::list<>
using	dg_package_data_primitive_tags = tmpl::list<>
using	dg_package_data_volume_tags = tmpl::list<>
using	dg_boundary_terms_volume_tags = tmpl::list<>

Public Member Functions
	Rusanov (const Rusanov &)=default
Rusanov &	operator= (const Rusanov &)=default
	Rusanov (Rusanov &&)=default
Rusanov &	operator= (Rusanov &&)=default
void	pup (PUP::er &p) override
std::unique_ptr< BoundaryCorrection >	get_clone () const override
double	dg_package_data (const gsl::not_null< typename Tags::TildeE< Frame::Inertial, NeutrinoSpecies >::type * >... packaged_tilde_e, const gsl::not_null< typename Tags::TildeS< Frame::Inertial, NeutrinoSpecies >::type * >... packaged_tilde_s, const gsl::not_null< typename Tags::TildeE< Frame::Inertial, NeutrinoSpecies >::type * >... packaged_normal_dot_flux_tilde_e, const gsl::not_null< typename Tags::TildeS< Frame::Inertial, NeutrinoSpecies >::type * >... packaged_normal_dot_flux_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::Flux< Tags::TildeE< Frame::Inertial, NeutrinoSpecies >, tmpl::size_t< 3 >, Frame::Inertial >::type &... flux_tilde_e, const typename ::Tags::Flux< Tags::TildeS< Frame::Inertial, NeutrinoSpecies >, tmpl::size_t< 3 >, Frame::Inertial >::type &... flux_tilde_s, const tnsr::i< DataVector, 3, Frame::Inertial > &normal_covector, const tnsr::I< DataVector, 3, Frame::Inertial > &normal_vector, const std::optional< tnsr::I< DataVector, 3, Frame::Inertial > > &mesh_velocity, const std::optional< Scalar< DataVector > > &normal_dot_mesh_velocity) const
void	dg_boundary_terms (const gsl::not_null< typename Tags::TildeE< Frame::Inertial, NeutrinoSpecies >::type * >... boundary_correction_tilde_e, const gsl::not_null< typename Tags::TildeS< Frame::Inertial, NeutrinoSpecies >::type * >... boundary_correction_tilde_s, const typename Tags::TildeE< Frame::Inertial, NeutrinoSpecies >::type &... tilde_e_int, const typename Tags::TildeS< Frame::Inertial, NeutrinoSpecies >::type &... tilde_s_int, const typename Tags::TildeE< Frame::Inertial, NeutrinoSpecies >::type &... normal_dot_flux_tilde_e_int, const typename Tags::TildeS< Frame::Inertial, NeutrinoSpecies >::type &... normal_dot_flux_tilde_s_int, const typename Tags::TildeE< Frame::Inertial, NeutrinoSpecies >::type &... tilde_e_ext, const typename Tags::TildeS< Frame::Inertial, NeutrinoSpecies >::type &... tilde_s_ext, const typename Tags::TildeE< Frame::Inertial, NeutrinoSpecies >::type &... normal_dot_flux_tilde_e_ext, const typename Tags::TildeS< Frame::Inertial, NeutrinoSpecies >::type &... normal_dot_flux_tilde_s_ext, dg::Formulation dg_formulation) const
Public Member Functions inherited from evolution::BoundaryCorrection
	BoundaryCorrection (const BoundaryCorrection &)=default
BoundaryCorrection &	operator= (const BoundaryCorrection &)=default
	BoundaryCorrection (BoundaryCorrection &&)=default
BoundaryCorrection &	operator= (BoundaryCorrection &&)=default
	BoundaryCorrection (CkMigrateMessage *msg)

Static Public Attributes
static constexpr Options::String	help

Detailed Description

template<typename... NeutrinoSpecies>
class RadiationTransport::M1Grey::BoundaryCorrections::Rusanov< tmpl::list< NeutrinoSpecies... > >

A Rusanov/local Lax-Friedrichs Riemann solver.

Let \(U\) be the state vector of evolved variables, \(F^i\) the corresponding fluxes, and \(n_i\) be the outward directed unit normal to the interface. Denoting \(F := n_i F^i\), the Rusanov boundary correction is

\begin{align*}G_\text{Rusanov} = \frac{F_\text{int} - F_\text{ext}}{2} - \frac{\text{max}\left(\{|\lambda_\text{int}|\}, \{|\lambda_\text{ext}|\}\right)}{2} \left(U_\text{ext} - U_\text{int}\right), \end{align*}

where "int" and "ext" stand for interior and exterior, and \(\{|\lambda|\}\) is the set of characteristic/signal speeds. The minus sign in front of the \(F_{\text{ext}}\) is necessary because the outward directed normal of the neighboring element has the opposite sign, i.e. \(n_i^{\text{ext}}=-n_i^{\text{int}}\).

For radiation, the max characteristic/signal speed is 1 (the speed of light). Note that the characteristic speeds of this system are not yet fully implemented, see the warning in the documentation for the characteristics.

Note

In the strong form the dg_boundary_terms function returns \(G - F_\text{int}\)

Member Typedef Documentation

dg_package_field_tags

template<typename... NeutrinoSpecies>

using RadiationTransport::M1Grey::BoundaryCorrections::Rusanov< tmpl::list< NeutrinoSpecies... > >::dg_package_field_tags

Initial value:

 tmpl::list<
      Tags::TildeE<Frame::Inertial, NeutrinoSpecies>...,
      Tags::TildeS<Frame::Inertial, NeutrinoSpecies>...,
      ::Tags::NormalDotFlux<Tags::TildeE<Frame::Inertial, NeutrinoSpecies>>...,
      ::Tags::NormalDotFlux<Tags::TildeS<Frame::Inertial, NeutrinoSpecies>>...>

Member Function Documentation

get_clone()

template<typename... NeutrinoSpecies>

std::unique_ptr< BoundaryCorrection > RadiationTransport::M1Grey::BoundaryCorrections::Rusanov< tmpl::list< NeutrinoSpecies... > >::get_clone ( ) const

inlineoverridevirtual

Implements evolution::BoundaryCorrection.

Member Data Documentation

help

template<typename... NeutrinoSpecies>

Options::String RadiationTransport::M1Grey::BoundaryCorrections::Rusanov< tmpl::list< NeutrinoSpecies... > >::help

staticconstexpr

Initial value:

= {
      "Computes the Rusanov or local Lax-Friedrichs boundary correction term "
      "for the M1Grey radiation transport system."}

---

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

• src/Evolution/Systems/RadiationTransport/M1Grey/BoundaryCorrections/Rusanov.hpp