evaporation_model_saturated_vapor_pressure.hpp

Go to the documentation of this file.

#pragma once
#include <meltpooldg/phase_change/evaporation_model_base.hpp>
#include <meltpooldg/phase_change/recoil_pressure_data.hpp>
 
namespace MeltPoolDG::Evaporation
{
  template <typename number>
  class EvaporationModelSaturatedVaporPressure : public EvaporationModelBase<number>
  {
  private:
    const RecoilPressureData<number> recoil_data;
    const number                     boiling_temperature;
    const number                     sticking_constant;
    const number                     molar_mass;
    const number                     latent_heat_evaporation;
 
    // Model constant computed as molar_mass/(2*pi*molar_gas_constant)
    const number Cm;
 
  public:
    EvaporationModelSaturatedVaporPressure(const RecoilPressureData<number> &recoil_data,
                                           const number                      boiling_temperature,
                                           const number                      molar_mass,
                                           const number latent_heat_evaporation);
 
    /*
     * The evaporative mass flux is computed as
     *
     *                                    ______
     *                                   /  M
     * .                                /-------
     * m(T) = 0.82 · c_s  · p_sat(T) · √  2πR T
     *
     * where c_s is the sticking coefficient, p_sat(T) is the saturated vapor
     * pressure, M the molar mass, R the molar gas constant and T the temperature.
     */
    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