Line data Source code
1 0 : // Distributed under the MIT License.
2 : // See LICENSE.txt for details.
3 :
4 : #pragma once
5 :
6 : #include <array>
7 : #include <cstddef>
8 : #include <utility>
9 :
10 : #include "DataStructures/DataBox/Prefixes.hpp"
11 : #include "DataStructures/Tensor/TypeAliases.hpp"
12 : #include "Domain/Structure/DirectionalIdMap.hpp"
13 : #include "Evolution/Systems/GrMhd/ValenciaDivClean/TagsDeclarations.hpp"
14 : #include "Evolution/VariableFixing/FixToAtmosphere.hpp"
15 : #include "PointwiseFunctions/GeneralRelativity/TagsDeclarations.hpp"
16 : #include "PointwiseFunctions/Hydro/TagsDeclarations.hpp"
17 : #include "Utilities/TMPL.hpp"
18 :
19 : /// \cond
20 : class DataVector;
21 : template <typename TagsList>
22 : class Variables;
23 : namespace gsl {
24 : template <typename>
25 : class not_null;
26 : } // namespace gsl
27 : template <size_t Dim>
28 : class Direction;
29 : template <size_t Dim>
30 : class ElementId;
31 : template <size_t Dim>
32 : class Element;
33 : template <size_t Dim>
34 : class Mesh;
35 : namespace EquationsOfState {
36 : template <bool IsRelativistic, size_t ThermodynamicDim>
37 : class EquationOfState;
38 : } // namespace EquationsOfState
39 : namespace evolution::dg::subcell {
40 : class GhostData;
41 : } // namespace evolution::dg::subcell
42 : namespace evolution::dg::Actions::detail {
43 : template <size_t Dim>
44 : struct NormalVector;
45 : } // namespace evolution::dg::Actions::detail
46 : /// \endcond
47 :
48 : namespace grmhd::ValenciaDivClean::fd {
49 : /// \brief The list of tags used by `reconstruct_prims_work` and
50 : /// `reconstruct_fd_neighbor_work`
51 1 : using tags_list_for_reconstruct = tmpl::list<
52 : grmhd::ValenciaDivClean::Tags::TildeD,
53 : grmhd::ValenciaDivClean::Tags::TildeYe,
54 : grmhd::ValenciaDivClean::Tags::TildeTau,
55 : grmhd::ValenciaDivClean::Tags::TildeS<Frame::Inertial>,
56 : grmhd::ValenciaDivClean::Tags::TildeB<Frame::Inertial>,
57 : grmhd::ValenciaDivClean::Tags::TildePhi,
58 : hydro::Tags::RestMassDensity<DataVector>,
59 : hydro::Tags::ElectronFraction<DataVector>,
60 : hydro::Tags::SpecificInternalEnergy<DataVector>,
61 : hydro::Tags::SpatialVelocity<DataVector, 3>,
62 : hydro::Tags::MagneticField<DataVector, 3>,
63 : hydro::Tags::DivergenceCleaningField<DataVector>,
64 : hydro::Tags::LorentzFactor<DataVector>, hydro::Tags::Pressure<DataVector>,
65 : hydro::Tags::Temperature<DataVector>,
66 : hydro::Tags::LorentzFactorTimesSpatialVelocity<DataVector, 3>,
67 : ::Tags::Flux<grmhd::ValenciaDivClean::Tags::TildeD, tmpl::size_t<3>,
68 : Frame::Inertial>,
69 : ::Tags::Flux<grmhd::ValenciaDivClean::Tags::TildeYe, tmpl::size_t<3>,
70 : Frame::Inertial>,
71 : ::Tags::Flux<grmhd::ValenciaDivClean::Tags::TildeTau, tmpl::size_t<3>,
72 : Frame::Inertial>,
73 : ::Tags::Flux<grmhd::ValenciaDivClean::Tags::TildeS<Frame::Inertial>,
74 : tmpl::size_t<3>, Frame::Inertial>,
75 : ::Tags::Flux<grmhd::ValenciaDivClean::Tags::TildeB<Frame::Inertial>,
76 : tmpl::size_t<3>, Frame::Inertial>,
77 : ::Tags::Flux<grmhd::ValenciaDivClean::Tags::TildePhi, tmpl::size_t<3>,
78 : Frame::Inertial>,
79 : gr::Tags::Lapse<DataVector>, gr::Tags::Shift<DataVector, 3>,
80 : hydro::Tags::SpatialVelocityOneForm<DataVector, 3, Frame::Inertial>,
81 : gr::Tags::SpatialMetric<DataVector, 3>,
82 : gr::Tags::SqrtDetSpatialMetric<DataVector>,
83 : gr::Tags::InverseSpatialMetric<DataVector, 3>,
84 : evolution::dg::Actions::detail::NormalVector<3>>;
85 :
86 : /*!
87 : * \brief Reconstructs the `PrimTagsForReconstruction` (usually
88 : * \f$\rho, p, Wv^i, B^i\f$, and \f$\Phi\f$), and if `compute_conservatives` is
89 : * true computes the Lorentz factor, upper spatial velocity, specific internal
90 : * energy, and the conserved variables.
91 : *
92 : * All results are written into `vars_on_lower_face` and `vars_on_upper_face`.
93 : *
94 : * The reason the `PrimTagsForReconstruction` can be specified separately is
95 : * because some variables might need separate reconstruction methods from
96 : * others, e.g. to guarantee the reconstructed solution is positive.
97 : */
98 : template <typename PrimTagsForReconstruction, typename PrimsTagsVolume,
99 : size_t ThermodynamicDim, typename F, typename PrimsTagsSentByNeighbor>
100 1 : void reconstruct_prims_work(
101 : gsl::not_null<std::array<Variables<tags_list_for_reconstruct>, 3>*>
102 : vars_on_lower_face,
103 : gsl::not_null<std::array<Variables<tags_list_for_reconstruct>, 3>*>
104 : vars_on_upper_face,
105 : const F& reconstruct, const Variables<PrimsTagsVolume>& volume_prims,
106 : const EquationsOfState::EquationOfState<true, ThermodynamicDim>& eos,
107 : const Element<3>& element,
108 : const DirectionalIdMap<3, Variables<PrimsTagsSentByNeighbor>>&
109 : neighbor_data,
110 : const Mesh<3>& subcell_mesh, size_t ghost_zone_size,
111 : bool compute_conservatives, bool reconstruct_density_times_temperature,
112 : const VariableFixing::FixToAtmosphere<3>* fix_to_atmosphere);
113 :
114 : /*!
115 : * \brief Reconstructs the `PrimTagsForReconstruction` and if
116 : * `compute_conservatives` is `true` computes the Lorentz factor, upper spatial
117 : * velocity, specific internal energy, and the conserved variables.
118 : *
119 : * All results are written into `vars_on_face`.
120 : *
121 : * This is used on DG elements to reconstruct their subcell neighbors' solution
122 : * on the shared faces.
123 : *
124 : * The reason the `PrimTagsForReconstruction` can be specified separately is
125 : * because some variables might need separate reconstruction methods from
126 : * others, e.g. to guarantee the reconstructed solution is positiv
127 : */
128 : template <typename PrimTagsForReconstruction, typename PrimsTagsSentByNeighbor,
129 : typename TagsList, typename PrimsTags, size_t ThermodynamicDim,
130 : typename F0, typename F1>
131 1 : void reconstruct_fd_neighbor_work(
132 : gsl::not_null<Variables<tags_list_for_reconstruct>*> vars_on_face,
133 : const F0& reconstruct_lower_neighbor, const F1& reconstruct_upper_neighbor,
134 : const Variables<PrimsTags>& subcell_volume_prims,
135 : const EquationsOfState::EquationOfState<true, ThermodynamicDim>& eos,
136 : const Element<3>& element,
137 : const DirectionalIdMap<3, evolution::dg::subcell::GhostData>& ghost_data,
138 : const Mesh<3>& subcell_mesh, const Direction<3>& direction_to_reconstruct,
139 : size_t ghost_zone_size, bool compute_conservatives,
140 : bool reconstruct_density_times_temperature,
141 : const VariableFixing::FixToAtmosphere<3>* fix_to_atmosphere);
142 : } // namespace grmhd::ValenciaDivClean::fd
|
