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/DataBox/Prefixes.hpp"
       9              : #include "DataStructures/DataBox/Tag.hpp"
      10              : #include "DataStructures/Tensor/TypeAliases.hpp"
      11              : #include "Domain/Tags.hpp"
      12              : #include "Evolution/DiscontinuousGalerkin/TimeDerivativeDecisions.hpp"
      13              : #include "Evolution/Systems/NewtonianEuler/Sources/Source.hpp"
      14              : #include "Evolution/Systems/NewtonianEuler/Tags.hpp"
      15              : #include "PointwiseFunctions/Hydro/Tags.hpp"
      16              : #include "Time/Tags/Time.hpp"
      17              : #include "Utilities/TMPL.hpp"
      18              : 
      19              : /// \cond
      20              : namespace gsl {
      21              : template <typename T>
      22              : class not_null;
      23              : }  // namespace gsl
      24              : class DataVector;
      25              : /// \endcond
      26              : 
      27              : namespace NewtonianEuler {
      28              : namespace detail {
      29              : template <size_t Dim>
      30              : void fluxes_impl(
      31              :     gsl::not_null<tnsr::I<DataVector, Dim>*> mass_density_cons_flux,
      32              :     gsl::not_null<tnsr::IJ<DataVector, Dim>*> momentum_density_flux,
      33              :     gsl::not_null<tnsr::I<DataVector, Dim>*> energy_density_flux,
      34              :     gsl::not_null<Scalar<DataVector>*> enthalpy_density,
      35              :     const tnsr::I<DataVector, Dim>& momentum_density,
      36              :     const Scalar<DataVector>& energy_density,
      37              :     const tnsr::I<DataVector, Dim>& velocity,
      38              :     const Scalar<DataVector>& pressure);
      39              : }  // namespace detail
      40              : 
      41              : /*!
      42              :  * \brief Compute the time derivative of the conserved variables for the
      43              :  * Newtonian Euler system
      44              :  */
      45              : template <size_t Dim>
      46            1 : struct TimeDerivativeTerms {
      47              :  private:
      48            0 :   struct EnthalpyDensity : db::SimpleTag {
      49            0 :     using type = Scalar<DataVector>;
      50              :   };
      51              : 
      52              :  public:
      53            0 :   using temporary_tags = tmpl::list<EnthalpyDensity>;
      54            0 :   using argument_tags =
      55              :       tmpl::list<Tags::MassDensityCons, Tags::MomentumDensity<Dim>,
      56              :                  Tags::EnergyDensity,
      57              :                  hydro::Tags::SpatialVelocity<DataVector, Dim>,
      58              :                  hydro::Tags::Pressure<DataVector>,
      59              :                  hydro::Tags::SpecificInternalEnergy<DataVector>,
      60              :                  hydro::Tags::EquationOfState<false, 2>,
      61              :                  domain::Tags::Coordinates<Dim, Frame::Inertial>, ::Tags::Time,
      62              :                  NewtonianEuler::Tags::SourceTerm<Dim>>;
      63              : 
      64            0 :   static evolution::dg::TimeDerivativeDecisions<Dim> apply(
      65              :       // Time derivatives returned by reference. All the tags in the
      66              :       // variables_tag in the system struct.
      67              :       const gsl::not_null<Scalar<DataVector>*> non_flux_terms_dt_mass_density,
      68              :       const gsl::not_null<tnsr::I<DataVector, Dim>*>
      69              :           non_flux_terms_dt_momentum_density,
      70              :       const gsl::not_null<Scalar<DataVector>*> non_flux_terms_dt_energy_density,
      71              : 
      72              :       // Fluxes returned by reference. Listed in the system struct as
      73              :       // flux_variables.
      74              :       const gsl::not_null<tnsr::I<DataVector, Dim>*> mass_density_cons_flux,
      75              :       const gsl::not_null<tnsr::IJ<DataVector, Dim>*> momentum_density_flux,
      76              :       const gsl::not_null<tnsr::I<DataVector, Dim>*> energy_density_flux,
      77              : 
      78              :       // Temporaries returned by reference. Listed in temporary_tags above.
      79              :       const gsl::not_null<Scalar<DataVector>*> enthalpy_density,
      80              : 
      81              :       // Arguments listed in argument_tags above
      82              :       const Scalar<DataVector>& mass_density_cons,
      83              :       const tnsr::I<DataVector, Dim>& momentum_density,
      84              :       const Scalar<DataVector>& energy_density,
      85              :       const tnsr::I<DataVector, Dim>& velocity,
      86              :       const Scalar<DataVector>& pressure,
      87              :       const Scalar<DataVector>& specific_internal_energy,
      88              :       const EquationsOfState::EquationOfState<false, 2>& eos,
      89              :       const tnsr::I<DataVector, Dim>& coords, const double time,
      90              :       const Sources::Source<Dim>& source) {
      91              :     detail::fluxes_impl(mass_density_cons_flux, momentum_density_flux,
      92              :                         energy_density_flux, enthalpy_density, momentum_density,
      93              :                         energy_density, velocity, pressure);
      94              : 
      95              :     get(*non_flux_terms_dt_mass_density) = 0.0;
      96              :     for (size_t i = 0; i < Dim; ++i) {
      97              :       non_flux_terms_dt_momentum_density->get(i) = 0.0;
      98              :     }
      99              :     get(*non_flux_terms_dt_energy_density) = 0.0;
     100              :     const auto eos_2d = eos.promote_to_2d_eos();
     101              :     source(non_flux_terms_dt_mass_density, non_flux_terms_dt_momentum_density,
     102              :            non_flux_terms_dt_energy_density, mass_density_cons,
     103              :            momentum_density, energy_density, velocity, pressure,
     104              :            specific_internal_energy, *eos_2d, coords, time);
     105              :     return {true};
     106              :   }
     107              : };
     108              : 
     109              : }  // namespace NewtonianEuler