single_track_melt.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/function.h>
#include <deal.II/base/parameter_handler.h>
#include <deal.II/base/point.h>
 
#include <deal.II/distributed/tria.h>
 
#include <deal.II/grid/grid_in.h>
 
#include <memory>
#include <string>
 
#include "../../../mp-heat-transfer/heat_transfer_case.hpp"
#include "../../melt_pool_case.hpp"
 
namespace MeltPoolDG::Simulation::SingleTrackMelt
{
  using namespace dealii;
 
  template <int dim>
  class GasInletVelocity : public Function<dim>
  {
  public:
    GasInletVelocity(const double inlet_velocity);
 
    double
    value(const Point<dim> &p, const unsigned int component = 0) const override;
 
    const double inlet_velocity;
  };
 
  template <int dim>
  class InitialSignedDistance : public Function<dim>
  {
  public:
    InitialSignedDistance(const double y_interface);
 
    double
    value(const Point<dim> &p, const unsigned int component = 0) const override;
 
  private:
    const double y_interface;
  };
 
  template <int dim, typename Number, typename CaseClass>
  class SimulationSingleTrackMelt : public CaseClass
  {
  public:
    SimulationSingleTrackMelt(std::string parameter_file, const MPI_Comm mpi_communicator);
 
    bool
    add_case_specific_parameters(ParameterHandler &prm) override;
 
    void
    create_spatial_discretization() override;
 
    void
    set_boundary_conditions() override;
 
    void
    set_field_conditions() override;
 
  private:
    std::string mesh_file_name = "single-track-melt.msh";
 
    double T_initial = 298.15;
 
    double inlet_temperature = 298.15;
 
    double inlet_velocity = 100.0;
 
    double outlet_pressure = 0.0;
 
    double y_interface = 0.0;
 
    mutable unsigned int n_time_step = 0;
 
    dealii::Point<dim> local_refinement_1_bottom_left;
 
    dealii::Point<dim> local_refinement_1_top_right;
 
    dealii::Point<dim> local_refinement_2_bottom_left;
 
    dealii::Point<dim> local_refinement_2_top_right;
 
    unsigned int n_additional_refinement_1 = 0;
 
    unsigned int n_additional_refinement_2 = 0;
  };
} // namespace MeltPoolDG::Simulation::SingleTrackMelt