surface_tension_operation.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/lac/full_matrix.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/flow/surface_tension_data.hpp>
#include <meltpooldg/level_set/delta_approximation_phase_weighted.hpp>
 
#include <memory>
 
namespace MeltPoolDG::Flow
{
  template <int dim, typename number>
  class SurfaceTensionOperation
  {
  public:
    using VectorType      = dealii::LinearAlgebra::distributed::Vector<number>;
    using BlockVectorType = dealii::LinearAlgebra::distributed::BlockVector<number>;
 
    explicit SurfaceTensionOperation(const SurfaceTensionData<number>    &data_in,
                                     const ScratchData<dim, dim, number> &scratch_data,
                                     const VectorType                    &level_set_as_heaviside,
                                     const VectorType                    &solution_curvature,
                                     const unsigned int                   ls_dof_idx,
                                     const unsigned int                   curv_dof_idx,
                                     const unsigned int flow_vel_hanging_nodes_dof_idx,
                                     const unsigned int flow_pressure_hanging_nodes_dof_idx,
                                     const unsigned int flow_quad_idx);
 
    void
    register_temperature_and_normal_vector(const unsigned int     temp_dof_idx,
                                           const unsigned int     normal_dof_idx,
                                           const VectorType      *temperature,
                                           const BlockVectorType *solution_normal_vector);
    void
    register_solid_fraction(const unsigned int solid_dof_idx, const VectorType *solid);
 
    void
    compute_surface_tension(VectorType &force_rhs, const bool zero_out = true);
 
    number
    compute_time_step_limit(const number density_1, const number density_2);
 
  private:
    const SurfaceTensionData<number> &data;
 
    const ScratchData<dim, dim, number> &scratch_data;
 
    const VectorType &level_set_as_heaviside;
 
    const VectorType &solution_curvature;
 
    const unsigned int ls_dof_idx;
 
    const unsigned int curv_dof_idx;
 
    const unsigned int flow_vel_dof_idx;
 
    const unsigned int flow_pressure_hanging_nodes_dof_idx;
 
    const unsigned int flow_vel_quad_idx;
 
    const VectorType *temperature = nullptr;
 
    const BlockVectorType *solution_normal_vector = nullptr;
 
    const VectorType *solid = nullptr;
 
    unsigned int temp_dof_idx;
 
    unsigned int normal_dof_idx;
 
    unsigned int solid_dof_idx;
 
    const bool do_level_set_pressure_gradient_interpolation;
 
    dealii::FullMatrix<number> ls_to_pressure_grad_interpolation_matrix;
 
    std::unique_ptr<const LevelSet::DeltaApproximationBase<number>> delta_phase_weighted;
 
    const number alpha_residual;
 
    template <typename number_surface_tension_coeff>
    number_surface_tension_coeff
    local_compute_temperature_dependent_surface_tension_coefficient(
      const number_surface_tension_coeff &T);
  };
} // namespace MeltPoolDG::Flow