obstacle_forces.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <deal.II/base/tensor.h>
 
#include <deal.II/particles/particle_accessor.h>
 
#include "meltpooldg/particles/particle_accessor.hpp"
 
#include <memory>
#include <utility>
 
namespace MeltPoolDG
{
  template <int dim, typename number, typename ObstacleType>
  class ObstacleField;
 
  template <int dim, typename number, typename ObstacleType>
  struct ObstacleLoad
  {
  private:
    struct ObstacleLoadConcept
    {
      virtual ~ObstacleLoadConcept() = default;
 
      virtual void
      add_load_to_obstacles(ObstacleField<dim, number, ObstacleType> &obstacle_field) const = 0;
    };
 
    template <typename ObstacleLoadType>
    struct ObstacleLoadModel final : public ObstacleLoadConcept
    {
      explicit ObstacleLoadModel(ObstacleLoadType &&obstacle_load)
        : obstacle_load(std::move(obstacle_load))
      {}
 
      void
      add_load_to_obstacles(ObstacleField<dim, number, ObstacleType> &obstacle_field) const override
      {
        obstacle_load.add_load_to_obstacles(obstacle_field);
      }
 
    private:
      const ObstacleLoadType obstacle_load;
    };
 
    std::unique_ptr<ObstacleLoadConcept> obstacle_load_pimpl;
 
  public:
    template <typename ObstacleLoadType>
    explicit ObstacleLoad(ObstacleLoadType &&obstacle_load)
      : obstacle_load_pimpl(
          std::make_unique<ObstacleLoadModel<ObstacleLoadType>>(std::move(obstacle_load)))
    {}
 
    void
    add_load_to_obstacles(ObstacleField<dim, number, ObstacleType> &obstacle_field) const
    {
      obstacle_load_pimpl->add_load_to_obstacles(obstacle_field);
    }
  };
 
  template <int dim, typename number, typename ObstacleType>
  struct ObstacleGravitationalForce
  {
  public:
    explicit ObstacleGravitationalForce(const number gravitational_constant)
      : gravitational_constant(gravitational_constant)
    {}
 
    void
    add_load_to_obstacles(ObstacleField<dim, number, ObstacleType> &obstacle_field) const
    {
      for (DEMParticleAccessor<dim, number> obstacle :
           obstacle_field.locally_owned_particle_range())
        {
          dealii::Tensor<1, dim, number> force;
          force[dim - 1] = -gravitational_constant * obstacle.mass();
          obstacle.add_force(force);
        }
    }
 
  private:
    const number gravitational_constant;
  };
} // namespace MeltPoolDG