Line data Source code
1 0 : // Distributed under the MIT License.
2 : // See LICENSE.txt for details.
3 :
4 : #pragma once
5 :
6 : #include <cstddef>
7 :
8 : #include "DataStructures/VariablesTag.hpp"
9 : #include "Evolution/Imex/GuessResult.hpp"
10 : #include "Evolution/Imex/Protocols/ImexSystem.hpp"
11 : #include "Evolution/Systems/RadiationTransport/M1Grey/BoundaryConditions/BoundaryCondition.hpp"
12 : #include "Evolution/Systems/RadiationTransport/M1Grey/Characteristics.hpp"
13 : #include "Evolution/Systems/RadiationTransport/M1Grey/M1Closure.hpp"
14 : #include "Evolution/Systems/RadiationTransport/M1Grey/M1HydroCoupling.hpp"
15 : #include "Evolution/Systems/RadiationTransport/M1Grey/Tags.hpp"
16 : #include "Evolution/Systems/RadiationTransport/M1Grey/TimeDerivativeTerms.hpp"
17 : #include "PointwiseFunctions/GeneralRelativity/Tags.hpp"
18 : #include "PointwiseFunctions/Hydro/Tags.hpp"
19 : #include "Utilities/ProtocolHelpers.hpp"
20 : #include "Utilities/TMPL.hpp"
21 :
22 : /// \ingroup EvolutionSystemsGroup
23 : /// \brief Items related to general relativistic radiation transport
24 : namespace RadiationTransport {
25 : /// The M1 scheme for radiation transport
26 : ///
27 : /// References:
28 : /// - Post-merger evolution of a neutron star-black hole binary with
29 : /// neutrino transport, \cite Foucart2015vpa
30 : namespace M1Grey {
31 :
32 : template <typename NeutrinoSpeciesList>
33 0 : struct System;
34 :
35 : template <typename... NeutrinoSpecies>
36 : struct System<tmpl::list<NeutrinoSpecies...>>
37 : : tt::ConformsTo<imex::protocols::ImexSystem> {
38 : static std::string name() { return "M1Grey"; }
39 :
40 : static constexpr bool is_in_flux_conservative_form = true;
41 : static constexpr bool has_primitive_and_conservative_vars = false;
42 : static constexpr size_t volume_dim = 3;
43 : // If coupling to hydro, we'll want 3D equations of state
44 : // i.e. P(rho,T,Ye)... but this is not implemented yet.
45 : // For early tests of M1, we'll ignore coupling to the fluid
46 : // and provide analytical expressions for its 4-velocity / LorentzFactor
47 :
48 : using boundary_conditions_base =
49 : BoundaryConditions::BoundaryCondition<tmpl::list<NeutrinoSpecies...>>;
50 :
51 : using variables_tag = ::Tags::Variables<
52 : tmpl::list<Tags::TildeE<Frame::Inertial, NeutrinoSpecies>...,
53 : Tags::TildeS<Frame::Inertial, NeutrinoSpecies>...>>;
54 : using flux_variables =
55 : tmpl::list<Tags::TildeE<Frame::Inertial, NeutrinoSpecies>...,
56 : Tags::TildeS<Frame::Inertial, NeutrinoSpecies>...>;
57 : using gradient_variables = tmpl::list<>;
58 : using primitive_variables_tag = ::Tags::Variables<tmpl::list<
59 : Tags::ClosureFactor<NeutrinoSpecies>...,
60 : Tags::TildeP<Frame::Inertial, NeutrinoSpecies>...,
61 : Tags::TildeJ<NeutrinoSpecies>..., Tags::TildeHNormal<NeutrinoSpecies>...,
62 : Tags::TildeHSpatial<Frame::Inertial, NeutrinoSpecies>...,
63 : Tags::M1HydroCouplingNormal<NeutrinoSpecies>...,
64 : Tags::M1HydroCouplingSpatial<Frame::Inertial, NeutrinoSpecies>...>>;
65 : // gr::tags_for_hydro contains all these tags plus SqrtDetSpatialMetric,
66 : // so it can be used when adding M1 coupling to hydro
67 : using spacetime_variables_tag = ::Tags::Variables<
68 : tmpl::list<gr::Tags::Lapse<DataVector>, gr::Tags::Shift<DataVector, 3>,
69 : gr::Tags::SpatialMetric<DataVector, 3>,
70 : gr::Tags::InverseSpatialMetric<DataVector, 3>,
71 : gr::Tags::SqrtDetSpatialMetric<DataVector>,
72 : ::Tags::deriv<gr::Tags::Lapse<DataVector>, tmpl::size_t<3>,
73 : Frame::Inertial>,
74 : ::Tags::deriv<gr::Tags::Shift<DataVector, 3>, tmpl::size_t<3>,
75 : Frame::Inertial>,
76 : ::Tags::deriv<gr::Tags::SpatialMetric<DataVector, 3>,
77 : tmpl::size_t<3>, Frame::Inertial>,
78 : gr::Tags::ExtrinsicCurvature<DataVector, 3>>>;
79 : using hydro_variables_tag = ::Tags::Variables<
80 : tmpl::list<hydro::Tags::LorentzFactor<DataVector>,
81 : hydro::Tags::SpatialVelocity<DataVector, 3>,
82 : Tags::GreyEmissivity<NeutrinoSpecies>...,
83 : Tags::GreyAbsorptionOpacity<NeutrinoSpecies>...,
84 : Tags::GreyScatteringOpacity<NeutrinoSpecies>...>>;
85 :
86 : using compute_volume_time_derivative_terms =
87 : TimeDerivativeTerms<NeutrinoSpecies...>;
88 :
89 : using inverse_spatial_metric_tag =
90 : gr::Tags::InverseSpatialMetric<DataVector, 3>;
91 :
92 : template <typename Species>
93 0 : struct ImplicitSector : tt::ConformsTo<imex::protocols::ImplicitSector> {
94 0 : using tensors = tmpl::list<Tags::TildeE<Frame::Inertial, Species>,
95 : Tags::TildeS<Frame::Inertial, Species>>;
96 : // Need initial_guess
97 0 : using initial_guess = imex::GuessExplicitResult;
98 :
99 0 : struct M1Solve {
100 : // tags in addition to sector tensors to be made available in databox
101 0 : using tags_from_evolution = tmpl::list<
102 : // spacetime
103 : gr::Tags::Lapse<DataVector>, gr::Tags::SpatialMetric<DataVector, 3>,
104 : gr::Tags::SqrtDetSpatialMetric<DataVector>,
105 : gr::Tags::InverseSpatialMetric<DataVector, 3>,
106 : // neutrino
107 : Tags::GreyEmissivity<NeutrinoSpecies>...,
108 : Tags::GreyAbsorptionOpacity<NeutrinoSpecies>...,
109 : Tags::GreyScatteringOpacity<NeutrinoSpecies>...,
110 : // hydro
111 : hydro::Tags::LorentzFactor<DataVector>,
112 : hydro::Tags::SpatialVelocity<DataVector, 3>>;
113 :
114 0 : using simple_tags = tmpl::list<::Tags::Variables<tmpl::list<
115 : Tags::ClosureFactor<NeutrinoSpecies>...,
116 : Tags::TildeP<Frame::Inertial, NeutrinoSpecies>...,
117 : Tags::TildeJ<NeutrinoSpecies>...,
118 : Tags::TildeHNormal<NeutrinoSpecies>...,
119 : Tags::TildeHSpatial<Frame::Inertial, NeutrinoSpecies>...>>>;
120 0 : using compute_tags = tmpl::list<>;
121 :
122 0 : using source_prep = tmpl::list<ComputeM1Closure<tmpl::list<Species>>>;
123 0 : using jacobian_prep = tmpl::list<ComputeM1Closure<tmpl::list<Species>>>;
124 :
125 0 : using source = ComputeM1HydroCoupling<tmpl::list<Species>>;
126 0 : using jacobian = ComputeM1HydroCouplingJacobian<Species>;
127 : };
128 :
129 0 : using solve_attempts = tmpl::list<M1Solve>;
130 : };
131 :
132 : using implicit_sectors = tmpl::list<ImplicitSector<NeutrinoSpecies>...>;
133 : };
134 : } // namespace M1Grey
135 : } // namespace RadiationTransport
|
