initial_state_four_state_contrast.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// Copyright (C) 2023 - 2026 by the ryujin authors
//
 
#pragma once
 
#include <ryujin/base/compile_time_options.h>
#include <ryujin/typed/initial_state.h>
 
namespace ryujin
{
  namespace EulerInitialStates
  {
    template <typename Description, int dim, typename Number>
    class FourStateContrast : public InitialState<Description, dim, Number>
    {
    public:
      using HyperbolicSystem = typename Description::HyperbolicSystem;
      using View = typename HyperbolicSystem::template View<dim, Number>;
      using state_type = typename View::state_type;
      using state_type_2d =
          typename HyperbolicSystem::template View<2, Number>::state_type;
 
      FourStateContrast(const HyperbolicSystem &hyperbolic_system,
                        const std::string &subsection)
          : InitialState<Description, dim, Number>("four state contrast",
                                                   subsection)
          , hyperbolic_system_(hyperbolic_system)
      {
        primitive_bottom_left_[0] = 1.4;
        primitive_bottom_left_[1] = 0.;
        primitive_bottom_left_[2] = 0.;
        if constexpr (View::have_energy_equation)
          primitive_bottom_left_[3] = 1.;
        this->add_parameter("primitive state bottom left",
                            primitive_bottom_left_,
                            "Primitive state [rho, u, v, p] on bottom left (or "
                            "[rho, u] for the barotropic Euler module)");
 
        primitive_bottom_right_[0] = 1.4;
        primitive_bottom_right_[1] = 0.;
        primitive_bottom_right_[2] = 0.;
        if constexpr (View::have_energy_equation)
          primitive_bottom_right_[3] = 1.;
        this->add_parameter("primitive state bottom right",
                            primitive_bottom_right_,
                            "Primitive state [rho, u, v, p] on bottom right "
                            "(or [rho, u] for the barotropic Euler module)");
 
        primitive_top_left_[0] = 1.4;
        primitive_top_left_[1] = 0.;
        primitive_top_left_[2] = 0.;
        if constexpr (View::have_energy_equation)
          primitive_top_left_[3] = 1.;
        this->add_parameter("primitive state top left",
                            primitive_top_left_,
                            "Primitive state [rho, u, v, p] on top left (or "
                            "[rho, u] for the barotropic Euler module)");
 
        primitive_top_right_[0] = 1.4;
        primitive_top_right_[1] = 0.;
        primitive_top_right_[2] = 0.;
        if constexpr (View::have_energy_equation)
          primitive_top_right_[3] = 1.;
        this->add_parameter("primitive state top right",
                            primitive_top_right_,
                            "Primitive state [rho, u, v, p] on top right (or "
                            "[rho, u] for the barotropic Euler module)");
 
        const auto convert_states = [&]() {
          const auto view = hyperbolic_system_.template view<dim, Number>();
          if constexpr (dim != 1) {
            state_bottom_left_ =
                view.from_initial_state(primitive_bottom_left_);
            state_bottom_right_ =
                view.from_initial_state(primitive_bottom_right_);
            state_top_left_ = view.from_initial_state(primitive_top_left_);
            state_top_right_ = view.from_initial_state(primitive_top_right_);
          }
        };
        this->parse_parameters_call_back.connect(convert_states);
        convert_states();
      }
 
      state_type compute(const dealii::Point<dim> &point, Number /*t*/) final
      {
        if constexpr (dim == 1) {
          AssertThrow(false, dealii::ExcNotImplemented());
          __builtin_trap();
 
        } else {
 
          const auto top = point[0] >= 0. ? state_top_right_ : state_top_left_;
          const auto bottom =
              point[0] >= 0. ? state_bottom_right_ : state_bottom_left_;
          return (point[1] >= 0. ? top : bottom);
        }
      }
 
    private:
      const HyperbolicSystem &hyperbolic_system_;
 
      state_type_2d primitive_bottom_left_;
      state_type_2d primitive_bottom_right_;
      state_type_2d primitive_top_left_;
      state_type_2d primitive_top_right_;
 
      state_type state_bottom_left_;
      state_type state_bottom_right_;
      state_type state_top_left_;
      state_type state_top_right_;
    };
  } // namespace EulerInitialStates
} // namespace ryujin