recoil_pressure_operation.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/exceptions.h>
#include <deal.II/base/vectorization.h>
 
#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/material_data.hpp>
#include <meltpooldg/core/scratch_data.hpp>
#include <meltpooldg/level_set/delta_approximation_phase_weighted.hpp>
#include <meltpooldg/phase_change/recoil_pressure_data.hpp>
 
#include <memory>
 
 
namespace MeltPoolDG::Evaporation
{
  template <typename number>
  class RecoilPressureModelBase
  {
  public:
    virtual ~RecoilPressureModelBase() = default;
 
    virtual number
    compute_recoil_pressure_coefficient(const number /*T*/) const
    {
      AssertThrow(false, dealii::ExcNotImplemented());
    };
 
    virtual dealii::VectorizedArray<number>
    compute_recoil_pressure_coefficient(const dealii::VectorizedArray<number> & /*T*/) const
    {
      AssertThrow(false, dealii::ExcNotImplemented());
    };
 
    virtual number
    compute_recoil_pressure_coefficient(const number /*T*/,
                                        const number /*m_dot*/,
                                        const number /*delta_coefficient*/) const
    {
      AssertThrow(false, dealii::ExcNotImplemented());
    };
 
    virtual dealii::VectorizedArray<number>
    compute_recoil_pressure_coefficient(
      const dealii::VectorizedArray<number> & /*T*/,
      const dealii::VectorizedArray<number> & /*m_dot*/,
      const dealii::VectorizedArray<number> & /*delta_coefficient*/) const
    {
      AssertThrow(false, dealii::ExcNotImplemented());
    };
  };
 
  template <typename number>
  class RecoilPressurePhenomenologicalModel : public RecoilPressureModelBase<number>
  {
    //@todo: to avoid compiler warnings regarding hidden overriden functions
    using RecoilPressureModelBase<number>::compute_recoil_pressure_coefficient;
 
  private:
    const RecoilPressureData<number> &recoil_data;
    const number                      boiling_temperature;
    const number                      molar_mass;
    const number                      latent_heat_evaporation;
    number                            activation_ramp_derivative;
 
  public:
    RecoilPressurePhenomenologicalModel(const RecoilPressureData<number> &recoil_data,
                                        const number                      boiling_temperature,
                                        const number                      molar_mass,
                                        const number                      latent_heat_evaporation);
 
    number
    compute_recoil_pressure_coefficient(const number T) const final;
 
    dealii::VectorizedArray<number>
    compute_recoil_pressure_coefficient(const dealii::VectorizedArray<number> &T) const final;
  };
 
  template <typename number>
  class RecoilPressureHybridModel : public RecoilPressureModelBase<number>
  {
    //@todo: to avoid compiler warnings regarding hidden overriden functions
    using RecoilPressureModelBase<number>::compute_recoil_pressure_coefficient;
 
  private:
    const RecoilPressureData<number>                 &recoil_data;
    const number                                      boiling_temperature;
    const number                                      density_coeff;
    const RecoilPressurePhenomenologicalModel<number> recoil_phenomenological;
 
  public:
    RecoilPressureHybridModel(const RecoilPressureData<number> &recoil_data,
                              const MaterialData<number>       &material_data);
 
    number
    compute_recoil_pressure_coefficient(const number T,
                                        const number m_dot,
                                        const number delta_coefficient) const final;
 
    dealii::VectorizedArray<number>
    compute_recoil_pressure_coefficient(
      const dealii::VectorizedArray<number> &T,
      const dealii::VectorizedArray<number> &m_dot,
      const dealii::VectorizedArray<number> &delta_coefficient) const final;
  };
 
 
  template <typename number>
  class RecoilPressureModelPressureAware : public RecoilPressureModelBase<number>
  {
    //@todo: to avoid compiler warnings regarding hidden overriden functions
    using RecoilPressureModelBase<number>::compute_recoil_pressure_coefficient;
 
  public:
    RecoilPressureModelPressureAware(
      const typename RecoilPressureData<number>::PressureAwareData &pressure_aware_data,
      const number                                                  boiling_temperature,
      const number                                                  molar_mass,
      const number                                                  latent_heat_evaporation);
 
    number
    compute_recoil_pressure_coefficient(const number T) const final;
 
    dealii::VectorizedArray<number>
    compute_recoil_pressure_coefficient(const dealii::VectorizedArray<number> &T) const final;
 
  private:
    const typename RecoilPressureData<number>::PressureAwareData pressure_aware_data;
 
    const number boiling_temperature;
 
    const number molar_mass;
 
    const number latent_heat_evaporation;
 
    std::vector<number> Kp;
 
    const number ambient_gas_pressure;
  };
 
 
  template <int dim, typename number>
  class RecoilPressureOperation
  {
  private:
    using VectorType      = dealii::LinearAlgebra::distributed::Vector<number>;
    using BlockVectorType = dealii::LinearAlgebra::distributed::BlockVector<number>;
 
    const ScratchData<dim, dim, number> &scratch_data;
    /*
     * indices for getting DoFHandler<dim> and Quadrature<dim> of relevant subproblems
     */
    const unsigned int flow_vel_dof_idx;
    const unsigned int flow_vel_quad_idx;
    const unsigned int flow_pressure_dof_idx;
    const unsigned int ls_dof_idx;
    const unsigned int heat_dof_idx;
 
    const bool                 do_level_set_pressure_gradient_interpolation;
    dealii::FullMatrix<number> ls_to_pressure_grad_interpolation_matrix;
 
    const RecoilPressureModelType model_type;
 
    std::unique_ptr<const RecoilPressureModelBase<number>> recoil_pressure_model;
 
    std::unique_ptr<const LevelSet::DeltaApproximationBase<number>> delta_phase_weighted;
 
    // For the one phase cut temperature, we need a dummy temperature value for the gas domain -
    // which doesn't have any temperature information. We use the value that is just below the
    // recoil pressure activation temperature, so the recoil pressure will reliably be zero in that
    // domain.
    const number dummy_temperature;
 
  public:
    RecoilPressureOperation(const ScratchData<dim, dim, number> &scratch_data_in,
                            const RecoilPressureData<number>    &recoil,
                            const MaterialData<number>          &material,
                            const unsigned int                   flow_vel_dof_idx_in,
                            const unsigned int                   flow_vel_quad_idx_in,
                            const unsigned int                   flow_pressure_dof_idx_in,
                            const unsigned int                   ls_dof_idx_in,
                            const unsigned int                   heat_dof_idx_in);
 
    void
    compute_recoil_pressure_force(VectorType        &force_rhs,
                                  const VectorType  &level_set_as_heaviside,
                                  const VectorType  &temperature,
                                  const VectorType  &evaporative_mass_flux,
                                  const unsigned int evapor_dof_idx,
                                  bool               zero_out = true) const;
  };
 
} // namespace MeltPoolDG::Evaporation