evaporation_model_recoil_pressure.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/vectorization.h>
 
#include <meltpooldg/phase_change/evaporation_model_base.hpp>
#include <meltpooldg/phase_change/recoil_pressure_data.hpp>
#include <meltpooldg/phase_change/recoil_pressure_operation.hpp>
 
namespace MeltPoolDG::Evaporation
{
  template <typename number>
  class EvaporationModelRecoilPressure : public EvaporationModelBase<number>
  {
  private:
    const RecoilPressurePhenomenologicalModel<number> recoil_model;
 
    // according to Meier 2020
    const number sticking_constant;
    const number Cm;                   // molar_mass/(2*pi*molar_gas_constant)
    const number temperature_constant; // molar_latent_heat_of_evaporation / molar_gas_constant
 
  public:
    EvaporationModelRecoilPressure(const RecoilPressureData<number> &recoil_data,
                                   const number                      boiling_temperature,
                                   const number                      molar_mass,
                                   const number                      latent_heat_evaporation);
 
    number
    local_compute_evaporative_mass_flux(const number T) const final;
 
    dealii::VectorizedArray<number>
    local_compute_evaporative_mass_flux_vec(const dealii::VectorizedArray<number> &T) const final;
 
    number
    local_compute_evaporative_mass_flux_derivative(const number T) const final;
 
    dealii::VectorizedArray<number>
    local_compute_evaporative_mass_flux_vec_derivative(
      const dealii::VectorizedArray<number> &T) const final;
  };
} // namespace MeltPoolDG::Evaporation