laser_heat_source_projection_based.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/function.h>
#include <deal.II/base/point.h>
 
#include <deal.II/lac/la_parallel_block_vector.h>
#include <deal.II/lac/la_parallel_vector.h>
 
#include <meltpooldg/core/scratch_data.hpp>
#include <meltpooldg/heat/laser_data.hpp>
#include <meltpooldg/level_set/delta_approximation_phase_weighted.hpp>
#include <meltpooldg/level_set/delta_approximation_phase_weighted_data.hpp>
 
#include <memory>
 
namespace MeltPoolDG::Heat
{
  template <int dim, typename number>
  class LaserHeatSourceProjectionBased
  {
    using VectorType      = dealii::LinearAlgebra::distributed::Vector<number>;
    using BlockVectorType = dealii::LinearAlgebra::distributed::BlockVector<number>;
 
  public:
    LaserHeatSourceProjectionBased(
      const LaserData<number>                                   &laser_data_in,
      const std::shared_ptr<const dealii::Function<dim, number>> intensity_profile_in,
      const bool variable_properties_over_interface_in,
      const LevelSet::DeltaApproximationPhaseWeightedData<number>
        &delta_approximation_phase_weighted_data);
 
    void
    compute_interfacial_heat_source(VectorType                          &heat_source_vector,
                                    const ScratchData<dim, dim, number> &scratch_data,
                                    const unsigned int                   heat_dof_idx,
                                    const VectorType                    &level_set_heaviside,
                                    const unsigned int                   ls_dof_idx,
                                    const bool                           zero_out       = true,
                                    const BlockVectorType               *normal_vector  = nullptr,
                                    const unsigned int                   normal_dof_idx = 0) const;
 
    void
    compute_interfacial_heat_source_sharp(VectorType                          &heat_rhs,
                                          const ScratchData<dim, dim, number> &scratch_data,
                                          const unsigned int                   heat_dof_idx,
                                          const VectorType                    &level_set_heaviside,
                                          const unsigned int                   ls_dof_idx,
                                          const bool                           zero_out = true,
                                          const BlockVectorType *normal_vector          = nullptr,
                                          const unsigned int     normal_dof_idx         = 0) const;
 
    void
    compute_interfacial_heat_source_sharp_conforming(
      VectorType                          &heat_rhs,
      const ScratchData<dim, dim, number> &scratch_data,
      const unsigned int                   heat_dof_idx,
      const unsigned int                   heat_quad_idx,
      const VectorType                    &level_set_heaviside,
      const unsigned int                   ls_dof_idx,
      const bool                           zero_out,
      const BlockVectorType               *normal_vector,
      const unsigned int                   normal_dof_idx) const;
 
  private:
    number
    local_compute_interfacial_heat_source(const dealii::Point<dim>             &p,
                                          const dealii::Tensor<1, dim, number> &normal_vector,
                                          const number                          delta_value,
                                          const number                          heaviside) const;
 
    const LaserData<number> &laser_data;
 
    const std::shared_ptr<const dealii::Function<dim, number>> intensity_profile;
 
    const dealii::Tensor<1, dim, number> laser_direction;
 
    const bool variable_properties_over_interface;
 
    std::unique_ptr<const LevelSet::DeltaApproximationBase<number>> delta_phase_weighted;
  };
} // namespace MeltPoolDG::Heat