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/Tensor/TypeAliases.hpp"
      10             : #include "Evolution/Systems/NewtonianEuler/TagsDeclarations.hpp"
      11             : #include "PointwiseFunctions/Hydro/Tags.hpp"
      12             : #include "Utilities/TMPL.hpp"
      13             : 
      14             : /// \cond
      15             : namespace gsl {
      16             : template <typename T>
      17             : class not_null;
      18             : }  // namespace gsl
      19             : 
      20             : class DataVector;
      21             : /// \endcond
      22             : 
      23             : namespace NewtonianEuler {
      24             : 
      25             : /*!
      26             :  * \brief Compute the fluxes of the conservative variables of the
      27             :  * Newtonian Euler system
      28             :  *
      29             :  * The fluxes are \f$(\text{Dim} + 2)\f$ vectors of
      30             :  * dimension \f$\text{Dim}\f$. Denoting the flux of the conservative
      31             :  * variable \f$u\f$ as \f$F(u)\f$, one has
      32             :  *
      33             :  * \f{align*}
      34             :  * F^i(\rho) &= S^i\\
      35             :  * F^i(S^j) &= S^i v^j + \delta^{ij}p\\
      36             :  * F^i(e) &= (e + p)v^i
      37             :  * \f}
      38             :  *
      39             :  * where \f$S^i\f$ is the momentum density, \f$e\f$ is the energy density,
      40             :  * \f$v^i\f$ is the velocity, \f$p\f$ is the pressure, and \f$\delta^{ij}\f$
      41             :  * is the Kronecker delta. This form of the fluxes combines conservative and
      42             :  * primitive variables (while the velocity appears explicitly, the pressure
      43             :  * implicitly depends, for instance, on the mass density and the specific
      44             :  * internal energy), so the conversion from one variable set to the other
      45             :  * must be known prior to calling this function.
      46             :  */
      47             : template <size_t Dim>
      48           1 : struct ComputeFluxes {
      49           0 :   using return_tags = tmpl::list<
      50             :       ::Tags::Flux<Tags::MassDensityCons, tmpl::size_t<Dim>, Frame::Inertial>,
      51             :       ::Tags::Flux<Tags::MomentumDensity<Dim>, tmpl::size_t<Dim>,
      52             :                    Frame::Inertial>,
      53             :       ::Tags::Flux<Tags::EnergyDensity, tmpl::size_t<Dim>, Frame::Inertial>>;
      54             : 
      55           0 :   using argument_tags =
      56             :       tmpl::list<Tags::MomentumDensity<Dim>, Tags::EnergyDensity,
      57             :                  hydro::Tags::SpatialVelocity<DataVector, Dim>,
      58             :                  hydro::Tags::Pressure<DataVector>>;
      59             : 
      60           0 :   static void apply(
      61             :       gsl::not_null<tnsr::I<DataVector, Dim>*> mass_density_cons_flux,
      62             :       gsl::not_null<tnsr::IJ<DataVector, Dim>*> momentum_density_flux,
      63             :       gsl::not_null<tnsr::I<DataVector, Dim>*> energy_density_flux,
      64             :       const tnsr::I<DataVector, Dim>& momentum_density,
      65             :       const Scalar<DataVector>& energy_density,
      66             :       const tnsr::I<DataVector, Dim>& velocity,
      67             :       const Scalar<DataVector>& pressure);
      68             : };
      69             : 
      70             : }  // namespace NewtonianEuler