Line data Source code
1 0 : // Distributed under the MIT License.
2 : // See LICENSE.txt for details.
3 :
4 : #pragma once
5 :
6 : #include <memory>
7 : #include <optional>
8 : #include <pup.h>
9 : #include <string>
10 : #include <type_traits>
11 : #include <unordered_map>
12 :
13 : #include "DataStructures/DataBox/Prefixes.hpp"
14 : #include "DataStructures/DataVector.hpp"
15 : #include "DataStructures/Tensor/Tensor.hpp"
16 : #include "DataStructures/Variables.hpp"
17 : #include "Domain/CoordinateMaps/CoordinateMap.hpp"
18 : #include "Domain/CoordinateMaps/Tags.hpp"
19 : #include "Domain/ElementMap.hpp"
20 : #include "Domain/FunctionsOfTime/FunctionOfTime.hpp"
21 : #include "Domain/FunctionsOfTime/Tags.hpp"
22 : #include "Domain/Structure/Direction.hpp"
23 : #include "Domain/Structure/DirectionMap.hpp"
24 : #include "Domain/Tags.hpp"
25 : #include "Evolution/BoundaryConditions/Type.hpp"
26 : #include "Evolution/DgSubcell/GhostZoneLogicalCoordinates.hpp"
27 : #include "Evolution/DgSubcell/Tags/Coordinates.hpp"
28 : #include "Evolution/DgSubcell/Tags/Mesh.hpp"
29 : #include "Evolution/Systems/Burgers/BoundaryConditions/BoundaryCondition.hpp"
30 : #include "Evolution/Systems/Burgers/FiniteDifference/Factory.hpp"
31 : #include "Evolution/Systems/Burgers/FiniteDifference/Reconstructor.hpp"
32 : #include "Evolution/Systems/Burgers/FiniteDifference/Tags.hpp"
33 : #include "Evolution/Systems/Burgers/Tags.hpp"
34 : #include "Evolution/TypeTraits.hpp"
35 : #include "NumericalAlgorithms/Spectral/Mesh.hpp"
36 : #include "Options/String.hpp"
37 : #include "PointwiseFunctions/AnalyticData/Tags.hpp"
38 : #include "PointwiseFunctions/AnalyticSolutions/AnalyticSolution.hpp"
39 : #include "PointwiseFunctions/InitialDataUtilities/InitialData.hpp"
40 : #include "Utilities/Gsl.hpp"
41 : #include "Utilities/Serialization/CharmPupable.hpp"
42 : #include "Utilities/TMPL.hpp"
43 :
44 : /// \cond
45 : namespace Tags {
46 : struct Time;
47 : } // namespace Tags
48 : /// \endcond
49 :
50 : namespace Burgers::BoundaryConditions {
51 : /*!
52 : * \brief Sets Dirichlet boundary conditions using the analytic solution or
53 : * analytic data.
54 : */
55 1 : class DirichletAnalytic final : public BoundaryCondition {
56 : private:
57 0 : using flux_tag =
58 : ::Tags::Flux<Burgers::Tags::U, tmpl::size_t<1>, Frame::Inertial>;
59 :
60 : public:
61 : /// \brief What analytic solution/data to prescribe.
62 1 : struct AnalyticPrescription {
63 0 : static constexpr Options::String help =
64 : "What analytic solution/data to prescribe.";
65 0 : using type = std::unique_ptr<evolution::initial_data::InitialData>;
66 : };
67 :
68 0 : using options = tmpl::list<AnalyticPrescription>;
69 0 : static constexpr Options::String help{
70 : "DirichletAnalytic boundary conditions setting the value of U to "
71 : "the analytic solution or analytic data."};
72 :
73 0 : DirichletAnalytic() = default;
74 0 : DirichletAnalytic(DirichletAnalytic&&) = default;
75 0 : DirichletAnalytic& operator=(DirichletAnalytic&&) = default;
76 0 : DirichletAnalytic(const DirichletAnalytic&);
77 0 : DirichletAnalytic& operator=(const DirichletAnalytic&);
78 0 : ~DirichletAnalytic() override = default;
79 :
80 0 : explicit DirichletAnalytic(
81 : std::unique_ptr<evolution::initial_data::InitialData>
82 : analytic_prescription);
83 :
84 0 : explicit DirichletAnalytic(CkMigrateMessage* msg);
85 :
86 0 : WRAPPED_PUPable_decl_base_template(
87 : domain::BoundaryConditions::BoundaryCondition, DirichletAnalytic);
88 :
89 0 : auto get_clone() const -> std::unique_ptr<
90 : domain::BoundaryConditions::BoundaryCondition> override;
91 :
92 0 : static constexpr evolution::BoundaryConditions::Type bc_type =
93 : evolution::BoundaryConditions::Type::Ghost;
94 :
95 0 : void pup(PUP::er& p) override;
96 :
97 0 : using dg_interior_evolved_variables_tags = tmpl::list<>;
98 0 : using dg_interior_temporary_tags =
99 : tmpl::list<domain::Tags::Coordinates<1, Frame::Inertial>>;
100 0 : using dg_gridless_tags = tmpl::list<::Tags::Time>;
101 :
102 0 : std::optional<std::string> dg_ghost(
103 : const gsl::not_null<Scalar<DataVector>*> u,
104 : const gsl::not_null<tnsr::I<DataVector, 1, Frame::Inertial>*> flux_u,
105 : const std::optional<
106 : tnsr::I<DataVector, 1, Frame::Inertial>>& /*face_mesh_velocity*/,
107 : const tnsr::i<DataVector, 1, Frame::Inertial>& /*normal_covector*/,
108 : const tnsr::I<DataVector, 1, Frame::Inertial>& coords,
109 : [[maybe_unused]] const double time) const;
110 :
111 0 : using fd_interior_evolved_variables_tags = tmpl::list<>;
112 0 : using fd_interior_temporary_tags =
113 : tmpl::list<evolution::dg::subcell::Tags::Mesh<1>>;
114 0 : using fd_gridless_tags =
115 : tmpl::list<::Tags::Time, domain::Tags::FunctionsOfTime,
116 : domain::Tags::ElementMap<1, Frame::Grid>,
117 : domain::CoordinateMaps::Tags::CoordinateMap<1, Frame::Grid,
118 : Frame::Inertial>,
119 : fd::Tags::Reconstructor>;
120 :
121 0 : void fd_ghost(const gsl::not_null<Scalar<DataVector>*> u,
122 : const Direction<1>& direction, const Mesh<1> subcell_mesh,
123 : const double time,
124 : const std::unordered_map<
125 : std::string,
126 : std::unique_ptr<domain::FunctionsOfTime::FunctionOfTime>>&
127 : functions_of_time,
128 : const ElementMap<1, Frame::Grid>& logical_to_grid_map,
129 : const domain::CoordinateMapBase<Frame::Grid, Frame::Inertial,
130 : 1>& grid_to_inertial_map,
131 : const fd::Reconstructor& reconstructor) const;
132 :
133 : private:
134 0 : static void flux_impl(
135 : gsl::not_null<tnsr::I<DataVector, 1, Frame::Inertial>*> flux,
136 : const Scalar<DataVector>& u_analytic);
137 :
138 0 : std::unique_ptr<evolution::initial_data::InitialData> analytic_prescription_;
139 : };
140 : } // namespace Burgers::BoundaryConditions
|
