Initial conditions for relativistic MHD Riemann problems. More...

#include <RiemannProblem.hpp>

Classes
struct	AdiabaticIndex

struct	Lapse

struct	LeftMagneticField

struct	LeftPressure

struct	LeftRestMassDensity

struct	LeftSpatialVelocity

struct	RightMagneticField

struct	RightPressure

struct	RightRestMassDensity

struct	RightSpatialVelocity

struct	ShiftX

Public Types
using	equation_of_state_type = EquationsOfState::IdealFluid< true >

using	options = tmpl::list< AdiabaticIndex, LeftRestMassDensity, RightRestMassDensity, LeftPressure, RightPressure, LeftSpatialVelocity, RightSpatialVelocity, LeftMagneticField, RightMagneticField, Lapse, ShiftX >

Public Types inherited from grmhd::AnalyticDataBase
template<typename DataType >
using	tags = tmpl::push_back< typename gr::AnalyticSolution< 3 >::template tags< DataType >, hydro::Tags::RestMassDensity< DataType >, hydro::Tags::ElectronFraction< DataType >, hydro::Tags::SpecificInternalEnergy< DataType >, hydro::Tags::Temperature< DataType >, hydro::Tags::Pressure< DataType >, hydro::Tags::SpatialVelocity< DataType, 3 >, hydro::Tags::MagneticField< DataType, 3 >, hydro::Tags::DivergenceCleaningField< DataType >, hydro::Tags::LorentzFactor< DataType >, hydro::Tags::SpecificEnthalpy< DataType > >

Public Member Functions
	RiemannProblem (const RiemannProblem &)=default

RiemannProblem &	operator= (const RiemannProblem &)=default

	RiemannProblem (RiemannProblem &&)=default

RiemannProblem &	operator= (RiemannProblem &&)=default

	RiemannProblem (double adiabatic_index, double left_rest_mass_density, double right_rest_mass_density, double left_pressure, double right_pressure, const std::array< double, 3 > &left_spatial_velocity, const std::array< double, 3 > &right_spatial_velocity, const std::array< double, 3 > &left_magnetic_field, const std::array< double, 3 > &right_magnetic_field, double lapse, double shift)

auto	get_clone () const -> std::unique_ptr< evolution::initial_data::InitialData > override

template<typename DataType , typename Tag1 , typename Tag2 , typename... Tags>
tuples::TaggedTuple< Tag1, Tag2, Tags... >	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< Tag1, Tag2, Tags... >) const
	Retrieve a collection of hydrodynamic variables at position x.

template<typename DataType , typename Tag >
tuples::TaggedTuple< Tag >	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< Tag >) const
	Retrieve the metric variables.

const EquationsOfState::IdealFluid< true > &	equation_of_state () const

void	pup (PUP::er &) override


template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::RestMassDensity< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::RestMassDensity< DataType > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::ElectronFraction< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::ElectronFraction< DataType > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::SpecificInternalEnergy< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::SpecificInternalEnergy< DataType > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::Pressure< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Pressure< DataType > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::SpatialVelocity< DataType, 3 > >) const -> tuples::TaggedTuple< hydro::Tags::SpatialVelocity< DataType, 3 > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::MagneticField< DataType, 3 > >) const -> tuples::TaggedTuple< hydro::Tags::MagneticField< DataType, 3 > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::DivergenceCleaningField< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::DivergenceCleaningField< DataType > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::LorentzFactor< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::LorentzFactor< DataType > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::SpecificEnthalpy< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::SpecificEnthalpy< DataType > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< gr::Tags::Lapse< DataType > >) const -> tuples::TaggedTuple< gr::Tags::Lapse< DataType > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< gr::Tags::Shift< DataType, 3 > >) const -> tuples::TaggedTuple< gr::Tags::Shift< DataType, 3 > >
	Retrieve the GRMHD variables at a given position.

template<typename DataType >
auto	variables (const tnsr::I< DataType, 3 > &x, tmpl::list< hydro::Tags::Temperature< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Temperature< DataType > >
	Retrieve the GRMHD variables at a given position.

virtual auto	get_clone () const -> std::unique_ptr< InitialData >=0

Public Member Functions inherited from hydro::TemperatureInitialization< RiemannProblem >
auto	variables (const tnsr::I< DataType, Dim > &x, tmpl::list< hydro::Tags::Temperature< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Temperature< DataType > >

auto	variables (const tnsr::I< DataType, Dim > &x, const double t, tmpl::list< hydro::Tags::Temperature< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Temperature< DataType > >

auto	variables (ExtraVars &extra_variables, const tnsr::I< DataType, Dim > &x, Args &... extra_args, tmpl::list< hydro::Tags::Temperature< DataType > >) const -> tuples::TaggedTuple< hydro::Tags::Temperature< DataType > >

Static Public Attributes
static constexpr Options::String	help

Static Public Attributes inherited from grmhd::AnalyticDataBase
static constexpr size_t	volume_dim = 3_st

Friends
bool	operator== (const RiemannProblem &lhs, const RiemannProblem &rhs)

bool	operator!= (const RiemannProblem &lhs, const RiemannProblem &rhs)

Detailed Description

Initial conditions for relativistic MHD Riemann problems.

The standard problems were first collected in [10]. A complete Riemann solver for the RMHD equations is presented in [76] and can be downloaded from https://www.brunogiacomazzo.org/?page_id=395

The domain is from \([-0.5,0.5]^3\) with periodic boundary conditions in \(y\) and \(z\). The problems are usually run at a resolution of 1600 finite difference/finite volume grid points and the initial discontinuity is located at \(x=0\).

While the problems were originally run in Minkowski space, changing the lapse \(\alpha\) to a non-unity constant value and/or the shift \(\beta^x\) to a non-zero constant value allows testing some of the metric terms in the evolution equations in a fairly simple setup.

Below are the initial conditions for the 5 different Balsara Riemann problems. Please note that RP5 has a different final time than the rest, and that RP1 has a different adiabatic index than the rest.

RP1:

AdiabaticIndex: 2.0
LeftDensity: 1.0
LeftPressure: 1.0
LeftVelocity: [0.0, 0.0, 0.0]
LeftMagneticField: [0.5, 1.0, 0.0]
RightDensity: 0.125
RightPressure: 0.1
RightVelocity: [0.0, 0.0, 0.0]
RightMagneticField: [0.5, -1.0, 0.0]
Lapse: 1.0
ShiftX: 0.0
Final time: 0.4

RP2:

AdiabaticIndex: 1.66666666666666666
LeftDensity: 1.0
LeftPressure: 30.0
LeftVelocity: [0.0, 0.0, 0.0]
LeftMagneticField: [5.0, 6.0, 6.0]
RightDensity: 1.0
RightPressure: 1.0
RightVelocity: [0.0, 0.0, 0.0]
RightMagneticField: [5.0, 0.7, 0.7]
Lapse: 1.0
ShiftX: 0.0
Final time: 0.4

RP3:

AdiabaticIndex: 1.66666666666666666
LeftDensity: 1.0
LeftPressure: 1000.0
LeftVelocity: [0.0, 0.0, 0.0]
LeftMagneticField: [10.0, 7.0, 7.0]
RightDensity: 1.0
RightPressure: 0.1
RightVelocity: [0.0, 0.0, 0.0]
RightMagneticField: [10.0, 0.7, 0.7]
Lapse: 1.0
ShiftX: 0.0
Final time: 0.4

RP4:

AdiabaticIndex: 1.66666666666666666
LeftDensity: 1.0
LeftPressure: 0.1
LeftVelocity: [0.999, 0.0, 0.0]
LeftMagneticField: [10.0, 7.0, 7.0]
RightDensity: 1.0
RightPressure: 0.1
RightVelocity: [-0.999, 0.0, 0.0]
RightMagneticField: [10.0, -7.0, -7.0]
Lapse: 1.0
ShiftX: 0.0
Final time: 0.4

RP5:

AdiabaticIndex: 1.66666666666666666
LeftDensity: 1.08
LeftPressure: 0.95
LeftVelocity: [0.4, 0.3, 0.2]
LeftMagneticField: [2.0, 0.3, 0.3]
RightDensity: 1.0
RightPressure: 1.0
RightVelocity: [-0.45, -0.2, 0.2]
RightMagneticField: [2.0, -0.7, 0.5]
Lapse: 1.0
ShiftX: 0.0
Final time: 0.55

Member Function Documentation

◆ get_clone()

auto grmhd::AnalyticData::RiemannProblem::get_clone ( ) const -> std::unique_ptr< evolution::initial_data::InitialData >

overridevirtual

Implements evolution::initial_data::InitialData.

Member Data Documentation

◆ help

constexpr Options::String grmhd::AnalyticData::RiemannProblem::help

staticconstexpr

Initial value:

= {
      "Analytic initial data for a GRMHD Riemann problem. The fluid variables "
      "are set homogeneously on either half of the domain left and right of "
      "x=0."}

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

src/PointwiseFunctions/AnalyticData/GrMhd/RiemannProblem.hpp

Classes

Public Types

Public Member Functions

Static Public Attributes

Friends

Detailed Description

Member Function Documentation

◆ get_clone()

Member Data Documentation

◆ help