fixed_melt_pool.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 <meltpooldg/level_set/level_set_type.hpp>
#include <meltpooldg/utilities/numbers.hpp>
 
#include <string>
 
#include "../../heat_transfer_case.hpp"
 
namespace MeltPoolDG::Simulation::FixedMeltPool
{
  template <int dim, typename number>
  class FixedMeltPoolGeometry : public dealii::Function<dim, number>
  {
  public:
    FixedMeltPoolGeometry(const LevelSet::LevelSetType level_set_type, const number eps = 0.0);
 
    number
    value(const dealii::Point<dim> &p, const unsigned int /*component*/) const override;
 
  private:
    number
    signed_distance_level(const dealii::Point<dim> &p) const;
 
    number
    signed_distance_pool(const dealii::Point<dim> &p) const;
 
    number
    signed_distance_beads(const dealii::Point<dim> &p) const;
 
    bool
    is_inside_melt_pool(const dealii::Point<dim> &p) const;
 
    const number y_level = 10e-6;
 
    const dealii::Point<dim> centre{0.0, 0.0};
    const number             centre_radius = 50e-6;
    const number             bead_radius   = 10e-6;
 
    const number melt_pool_radius = centre_radius + bead_radius;
    const number bead_center_x    = centre[0] + melt_pool_radius;
 
    const LevelSet::LevelSetType level_set_type;
    const number                 eps;
  };
 
  template <int dim, typename number>
  class SimulationFixedMeltPool : public Heat::HeatTransferCase<dim, number>
  {
  private:
    const number a = 100e-6;
 
    std::string grid_file;
 
    const number T_hat = 500.0;
 
    number interface_thickness = numbers::invalid_double;
 
  public:
    SimulationFixedMeltPool(std::string parameter_file, const MPI_Comm mpi_communicator);
 
    bool
    add_case_specific_parameters(dealii::ParameterHandler &prm) override;
 
    void
    create_spatial_discretization() override;
 
    void
    set_boundary_conditions() override;
 
    void
    set_field_conditions() override;
  };
 
} // namespace MeltPoolDG::Simulation::FixedMeltPool