obstacle_data.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/parameter_handler.h>
 
#include <meltpooldg/particles/cohesive_forces_data.hpp>
#include <meltpooldg/particles/contact_forces_data.hpp>
 
 
namespace MeltPoolDG
{
  template <typename number>
  struct ObstacleData
  {
    bool stationary_obstacles = false;
 
    std::string obstacle_state_input_file;
 
    struct
    {
      bool do_refine_obstacles = false;
 
      number inner_fractional_distance_to_surface = 1.;
 
      number outer_fractional_distance_to_surface = 1.;
    } amr;
 
    SphericalParticleCohesiveForceData<number> cohesive_forces;
 
    SphericalParticleContactData<number> contact_forces;
 
    number gravity_constant = 0.0;
 
    bool do_rayleigh_time_stepping = false;
 
    number rayleigh_time_stepping_safety_factor = 0.2;
 
    void
    add_parameters(dealii::ParameterHandler &prm)
    {
      prm.enter_subsection("obstacles");
      {
        prm.add_parameter("stationary",
                          stationary_obstacles,
                          "Set to true if obstacles shall be stationary.");
        prm.add_parameter("obstacle state input file",
                          obstacle_state_input_file,
                          "File in which the obstacle initial state data is stored.");
        prm.add_parameter(
          "gravity constant",
          gravity_constant,
          "Gravity constant data defining the gravitational acceleration acting on all particles.");
        prm.add_parameter(
          "rayleigh time stepping",
          do_rayleigh_time_stepping,
          "If true, the time step size is determined from the critical Rayleigh time step.");
        prm.add_parameter("rayleigh time stepping safety factor",
                          rayleigh_time_stepping_safety_factor,
                          "Fraction of the critical Rayleigh time step used as time step size.");
        prm.enter_subsection("amr");
        {
          prm.add_parameter(
            "do refine",
            amr.do_refine_obstacles,
            "If true and AMR is enabled, regions around obstacles are adaptively refined "
            "according to the specified refinement parameters.");
          prm.add_parameter(
            "inner surface distance",
            amr.inner_fractional_distance_to_surface,
            "Fraction of a reference length (e.g., particle radius) used to determine the "
            "distance from the obstacle surface within which cells inside the obstacles "
            "are refined.");
          prm.add_parameter(
            "outer surface distance",
            amr.outer_fractional_distance_to_surface,
            "Fraction of a reference length (e.g., particle radius) used to determine the "
            "distance from the obstacle surface within which cells outside the obstacles "
            "are refined.");
        }
        prm.leave_subsection();
 
        cohesive_forces.add_parameters(prm);
        contact_forces.add_parameters(prm);
      }
      prm.leave_subsection();
    }
  };
} // namespace MeltPoolDG