boundary_condition_functions.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <deal.II/base/exception_macros.h>
#include <deal.II/base/function.h>
 
#include <meltpooldg/compressible_flow/boundary_conditions.hpp>
#include <meltpooldg/compressible_flow/data_types.hpp>
#include <meltpooldg/utilities/better_enum.hpp>
 
#include <memory>
 
namespace MeltPoolDG::CompressibleFlow
{
  BETTER_ENUM(RampUpType, char, none, linear, exponential, cosine);
 
  template <int dim, typename number>
  class ConservativeVariablesFunction : public dealii::Function<dim, number>
  {
  public:
    ConservativeVariablesFunction(
      const number                                         initial_time,
      const std::shared_ptr<dealii::Function<dim, number>> density,
      const std::shared_ptr<dealii::Function<dim, number>> velocity,
      const std::shared_ptr<dealii::Function<dim, number>> inner_energy,
      const std::shared_ptr<dealii::Function<dim, number>> mass_fractions);
 
    ConservativeVariablesFunction(
      const number                                         initial_time,
      const std::shared_ptr<dealii::Function<dim, number>> density,
      const std::shared_ptr<dealii::Function<dim, number>> velocity,
      const std::shared_ptr<dealii::Function<dim, number>> inner_energy);
 
    ConservativeVariablesFunction(const number       initial_time,
                                  const std::string &density_function_string,
                                  const std::string &velocity_function_string,
                                  const std::string &inner_energy_function_string,
                                  const std::string &mass_fraction_functions_string,
                                  const unsigned     n_mass_fractions);
 
 
    ConservativeVariablesFunction(const number       initial_time,
                                  const std::string &density_function_string,
                                  const std::string &velocity_function_string,
                                  const std::string &inner_energy_function_string);
 
    void
    set_velocity_ramp_up(const number duration, const RampUpType type)
    {
      velocity_ramp_up.duration = duration;
      velocity_ramp_up.type     = type;
    }
 
    void
    set_time(const number new_time) override;
 
    number
    value(const dealii::Point<dim, number> &loc, const unsigned int component) const override;
 
  private:
    const std::shared_ptr<dealii::Function<dim, number>> density;
 
    const std::shared_ptr<dealii::Function<dim, number>> velocity;
 
    const std::shared_ptr<dealii::Function<dim, number>> inner_energy;
 
    const std::shared_ptr<dealii::Function<dim, number>> mass_fractions;
 
    struct RampUpParameters
    {
      number duration{0.};
 
      RampUpType type{RampUpType::none};
    };
 
    RampUpParameters velocity_ramp_up;
  };
 
 
  BETTER_ENUM(FreeJetProfile, char, cosine, constant)
 
  
  template <int dim, typename number>
  class FreeJetVelocityFunction : public dealii::Function<dim, number>
  {
  public:
    FreeJetVelocityFunction(const number                     peak_inflow_velocity,
                            const FreeJetProfile             free_jet_profile,
                            const dealii::Point<dim, number> free_jet_center,
                            const number                     free_jet_diameter);
 
    number
    value(const dealii::Point<dim, number> &location, const unsigned int component) const final;
 
  private:
    const FreeJetProfile free_jet_profile;
 
    const dealii::Point<dim, number> free_jet_center;
 
    const number free_jet_radius;
 
    const number peak_inflow_velocity;
  };
 
  template <int dim, typename number>
  class FreeJetInflow : public dealii::Function<dim, number>
  {
  public:
    FreeJetInflow(const number                     initial_time,
                  const number                     density,
                  const number                     inner_energy,
                  const number                     peak_inflow_velocity,
                  const FreeJetProfile             free_jet_profile,
                  const dealii::Point<dim, number> free_jet_center,
                  const number                     free_jet_diameter);
 
    FreeJetInflow(const number                     initial_time,
                  const number                     density,
                  const number                     inner_energy,
                  const std::vector<number>       &species_mass_fractions,
                  const number                     peak_inflow_velocity,
                  const FreeJetProfile             free_jet_profile,
                  const dealii::Point<dim, number> free_jet_center,
                  const number                     free_jet_diameter);
 
    void
    set_velocity_ramp_up(const number duration, const RampUpType type);
 
    void
    set_time(const number new_time) override;
 
    number
    value(const dealii::Point<dim, number> &location, const unsigned int component) const final;
 
  private:
    ConservativeVariablesFunction<dim, number> conservative_variables;
 
    unsigned n_mass_fractions = 1;
 
    struct
    {
      number jet_hole_radius;
 
      dealii::Point<dim, number> jet_hole_center;
    } free_jet_parameters;
 
    bool
    point_inside_jet_hole(const dealii::Point<dim, number> &location) const;
  };
 
} // namespace MeltPoolDG::CompressibleFlow