Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <complex>
       7             : #include <cstddef>
       8             : #include <limits>
       9             : #include <memory>
      10             : #include <vector>
      11             : 
      12             : #include "DataStructures/SpinWeighted.hpp"
      13             : #include "Evolution/Systems/Cce/AnalyticSolutions/SphericalMetricData.hpp"
      14             : #include "Evolution/Systems/Cce/AnalyticSolutions/WorldtubeData.hpp"
      15             : #include "Evolution/Systems/Cce/Tags.hpp"
      16             : #include "Options/String.hpp"
      17             : #include "Utilities/Gsl.hpp"
      18             : #include "Utilities/Literals.hpp"
      19             : #include "Utilities/Serialization/CharmPupable.hpp"
      20             : #include "Utilities/TMPL.hpp"
      21             : 
      22             : /// \cond
      23             : class DataVector;
      24             : class ComplexDataVector;
      25             : /// \endcond
      26             : 
      27             : namespace Cce::Solutions {
      28             : 
      29             : /*!
      30             :  * \brief Computes the analytic data for the rotating Schwarzschild solution
      31             :  * described in \cite Barkett2019uae, section VI.C.
      32             :  *
      33             :  * \details This is a comparatively simple test which simply determines the
      34             :  * Schwarzschild metric in transformed coordinates given by \f$\phi\rightarrow
      35             :  * \phi + \omega u\f$, where \f$u\f$ is the retarded time.
      36             :  */
      37           1 : struct RotatingSchwarzschild : public SphericalMetricData {
      38           0 :   struct ExtractionRadius {
      39           0 :     using type = double;
      40           0 :     static constexpr Options::String help{
      41             :         "The extraction radius of the spherical solution"};
      42           0 :     static type lower_bound() { return 0.0; }
      43             :   };
      44           0 :   struct Mass {
      45           0 :     using type = double;
      46           0 :     static constexpr Options::String help{
      47             :         "The mass of the Schwarzschild black hole"};
      48           0 :     static type lower_bound() { return 0.0; }
      49             :   };
      50           0 :   struct Frequency {
      51           0 :     using type = double;
      52           0 :     static constexpr Options::String help{
      53             :         "The frequency of the coordinate rotation."};
      54           0 :     static type lower_bound() { return 0.0; }
      55             :   };
      56             : 
      57           0 :   using options = tmpl::list<ExtractionRadius, Mass, Frequency>;
      58             : 
      59           0 :   static constexpr Options::String help = {
      60             :       "Analytic solution representing a Schwarzschild black hole in a rotating "
      61             :       "frame"};
      62             : 
      63           0 :   WRAPPED_PUPable_decl_template(RotatingSchwarzschild);  // NOLINT
      64             : 
      65           0 :   explicit RotatingSchwarzschild(CkMigrateMessage* msg)
      66             :       : SphericalMetricData(msg) {}
      67             : 
      68             :   // clang doesn't manage to use = default correctly in this case
      69             :   // NOLINTNEXTLINE(modernize-use-equals-default)
      70           0 :   RotatingSchwarzschild() {}
      71             : 
      72           0 :   RotatingSchwarzschild(double extraction_radius, double mass,
      73             :                         double frequency);
      74             : 
      75           0 :   ~RotatingSchwarzschild() override = default;
      76             : 
      77           0 :   std::unique_ptr<WorldtubeData> get_clone() const override;
      78             : 
      79           0 :   void pup(PUP::er& p) override;
      80             : 
      81             :  protected:
      82             :   /// A no-op as the rotating Schwarzschild solution does not have substantial
      83             :   /// shared computation to prepare before the separate component calculations.
      84           1 :   void prepare_solution(const size_t /*output_l_max*/,
      85             :                         const double /*time*/) const override {}
      86             : 
      87             :   /*!
      88             :    * \brief Compute the spherical coordinate metric from the closed-form
      89             :    * rotating Schwarzschild metric.
      90             :    *
      91             :    * \details The rotating Schwarzschild takes the coordinate form
      92             :    * \cite Barkett2019uae,
      93             :    *
      94             :    * \f{align}{
      95             :    * ds^2 = -\left(1 - \frac{2 M}{r} - \omega^2 r^2 \sin^2 \theta\right) dt^2
      96             :    * + \frac{1}{1 - \frac{2 M}{r}} dr^2
      97             :    * + 2 \omega r^2 \sin^2 \theta dt d\phi + r^2 d\Omega^2
      98             :    * \f}
      99             :    */
     100           1 :   void spherical_metric(
     101             :       gsl::not_null<
     102             :           tnsr::aa<DataVector, 3, ::Frame::Spherical<::Frame::Inertial>>*>
     103             :           spherical_metric,
     104             :       size_t l_max, double time) const override;
     105             : 
     106             :   /*!
     107             :    * \brief Compute the radial derivative of the spherical coordinate metric
     108             :    * from the closed-form rotating Schwarzschild metric.
     109             :    *
     110             :    * \details The rotating Schwarzschild takes the coordinate form
     111             :    * \cite Barkett2019uae,
     112             :    *
     113             :    * \f{align}{
     114             :    * \partial_r g_{a b} dx^a dx^b =& -\left(\frac{2 M}{r^2} - 2 \omega^2 r
     115             :    * \sin^2 \theta\right) dt^2 - \frac{2 M}{(r - 2 M)^2} dr^2
     116             :    * + 4 \omega r \sin^2 \theta dt d\phi + 2 r d\Omega^2
     117             :    * \f}
     118             :    */
     119           1 :   void dr_spherical_metric(
     120             :       gsl::not_null<
     121             :           tnsr::aa<DataVector, 3, ::Frame::Spherical<::Frame::Inertial>>*>
     122             :           dr_spherical_metric,
     123             :       size_t l_max, double time) const override;
     124             : 
     125             :   /*!
     126             :    * \brief The time derivative of the spherical coordinate metric in the
     127             :    * rotating Schwarzschild metric vanishes.
     128             :    */
     129           1 :   void dt_spherical_metric(
     130             :       gsl::not_null<
     131             :           tnsr::aa<DataVector, 3, ::Frame::Spherical<::Frame::Inertial>>*>
     132             :           dt_spherical_metric,
     133             :       size_t l_max, double time) const override;
     134             : 
     135           0 :   using WorldtubeData::variables_impl;
     136             : 
     137           1 :   using SphericalMetricData::variables_impl;
     138             : 
     139             :   /// The News vanishes in the rotating Schwarzschild metric
     140           1 :   void variables_impl(
     141             :       gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, -2>>*> news,
     142             :       size_t l_max, double time,
     143             :       tmpl::type_<Tags::News> /*meta*/) const override;
     144             : 
     145           0 :   double frequency_ = std::numeric_limits<double>::signaling_NaN();
     146           0 :   double mass_ = std::numeric_limits<double>::signaling_NaN();
     147             : };
     148             : }  // namespace Cce::Solutions