Line data    Source code

       1           1 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : /// \file
       5             : /// Defines class ScalarWaveSystem.
       6             : 
       7             : #pragma once
       8             : 
       9             : #include <cstddef>
      10             : 
      11             : #include "DataStructures/VariablesTag.hpp"
      12             : #include "Evolution/Systems/ScalarWave/BoundaryConditions/BoundaryCondition.hpp"
      13             : #include "Evolution/Systems/ScalarWave/Characteristics.hpp"
      14             : #include "Evolution/Systems/ScalarWave/Equations.hpp"
      15             : #include "Evolution/Systems/ScalarWave/Tags.hpp"
      16             : #include "Evolution/Systems/ScalarWave/TimeDerivative.hpp"
      17             : #include "Utilities/TMPL.hpp"
      18             : 
      19             : /*!
      20             :  * \ingroup EvolutionSystemsGroup
      21             :  * \brief Items related to evolving the scalar wave equation.
      22             :  *
      23             :  * The equations of motion for the system augmented with constraint damping
      24             :  * terms are given by Eq. (15), (23) and (24) of \cite Holst2004wt (setting
      25             :  * background spacetime to Minkowskian):
      26             :  *
      27             :  * \f{align*}
      28             :  * \partial_t \psi =& -\Pi \\
      29             :  * \partial_t \Pi  =& -\partial^i \Phi_i \\
      30             :  * \partial_t \Phi_i =& -\partial_i \Pi + \gamma_2 (\partial_i \psi - \Phi_i)
      31             :  * \f}
      32             :  *
      33             :  * In our implementation here, to disable the constraint damping terms,
      34             :  * set \f$\gamma_2 = 0\f$.
      35             :  */
      36             : namespace ScalarWave {
      37             : 
      38             : template <size_t Dim>
      39           0 : struct System {
      40           0 :   using boundary_conditions_base = BoundaryConditions::BoundaryCondition<Dim>;
      41             : 
      42           0 :   static constexpr bool is_in_flux_conservative_form = false;
      43           0 :   static constexpr bool has_primitive_and_conservative_vars = false;
      44           0 :   static constexpr size_t volume_dim = Dim;
      45             : 
      46           0 :   using variables_tag =
      47             :       ::Tags::Variables<tmpl::list<Tags::Psi, Tags::Pi, Tags::Phi<Dim>>>;
      48           0 :   using flux_variables = tmpl::list<>;
      49           0 :   using gradient_variables = tmpl::list<Tags::Psi, Tags::Pi, Tags::Phi<Dim>>;
      50             : 
      51           0 :   using compute_volume_time_derivative_terms = TimeDerivative<Dim>;
      52             : 
      53           0 :   using compute_largest_characteristic_speed =
      54             :       Tags::ComputeLargestCharacteristicSpeed;
      55             : };
      56             : }  // namespace ScalarWave