case_utils.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/function.h>
#include <deal.II/base/mpi.h>
#include <deal.II/base/parameter_handler.h>
#include <deal.II/base/point.h>
 
#include <deal.II/grid/tria.h>
 
#include <meltpooldg/compressible_flow/boundary_condition_functions.hpp>
#include <meltpooldg/compressible_flow/boundary_conditions.hpp>
#include <meltpooldg/compressible_flow/data_types.hpp>
 
#include <array>
#include <memory>
#include <string>
#include <vector>
 
namespace MeltPoolDG::CompressibleFlow
{
  template <int dim, typename number>
  struct InputDefinedFreeJetInflow
  {
    std::shared_ptr<dealii::Function<dim, number>>
    create_free_jet_inflow_boundary(const number     start_time,
                                    const number     velocity_ramp_up_duration,
                                    const RampUpType velocity_ramp_up_type,
                                    const unsigned   n_species = 1) const;
    void
    add_parameters(dealii::ParameterHandler &prm);
 
    number jet_hole_diameter{0.};
 
    number jet_hole_density{0.};
 
    number jet_hole_inner_energy{0.};
 
    std::vector<number> jet_hole_mass_fractions;
 
    dealii::Point<dim, number> jet_hole_center;
 
    number jet_peak_velocity{0.};
 
    FreeJetProfile free_jet_profile{FreeJetProfile::constant};
  };
 
 
  template <int dim, typename number>
  struct InputDefinedBoundaryCondition
  {
    BoundaryConditionType type = BoundaryConditionType::slip_wall;
 
    std::string density;
 
    std::string velocity;
 
    std::string pressure;
 
    std::string energy;
 
    std::string species_mass_fractions;
 
    struct RampUpParameters
    {
      number duration{0.};
      RampUpType type{RampUpType::none};
    };
 
    RampUpParameters inflow_velocity_ramp_up;
 
    InputDefinedFreeJetInflow<dim, number> free_jet_inflow;
 
    std::shared_ptr<dealii::Function<dim>>
    create_boundary_function(const number start_time, const int n_species = 1) const;
 
    void
    add_user_boundary_condition_parameters(dealii::ParameterHandler &prm,
                                           const std::string        &boundary_location);
  };
 
  template <int dim, typename number>
  void
  add_hyper_rectangle_custom_boundary_condition_parameters(
    dealii::ParameterHandler                                        &prm,
    std::array<InputDefinedBoundaryCondition<dim, number>, 2 * dim> &boundary_conditions);
 
  template <int dim, typename number>
  struct InputDefinedInitialCondition
  {
    std::string density;
 
    std::string velocity;
 
    std::string inner_energy;
 
    std::string species_mass_fractions;
 
    std::shared_ptr<dealii::Function<dim, number>>
    create_initial_condition_function(const number initial_time, const int n_species = 1) const;
 
    void
    add_parameters(dealii::ParameterHandler &prm);
  };
 
  template <int dim, typename number>
  struct InputDefinedSubdividedHyperRectangleDomain
  {
    void
    create_triangulation(std::shared_ptr<dealii::Triangulation<dim>> &triangulation,
                         const MPI_Comm                              &mpi_communicator,
                         const unsigned                               n_global_refinements);
 
    void
    add_parameters(dealii::ParameterHandler &prm);
 
    std::vector<number> domain_dimensions;
 
    std::vector<unsigned> domain_base_discretization;
 
    std::array<std::vector<number>, 2> region_of_interest_corner;
 
    unsigned region_of_interest_refinement_times = 0;
 
    std::vector<bool> periodic_boundaries = std::vector<bool>(dim, false);
  };
 
} // namespace MeltPoolDG::CompressibleFlow