evaporation_model_knight.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/vectorization.h>
 
namespace MeltPoolDG::Evaporation
{
  template <typename number, typename number_2 = number>
  class EvaporationModelKnight
  {
  public:
    explicit EvaporationModelKnight(const number atmospheric_pressure,
                                    const number boiling_temperature_at_atmospheric_pressure,
                                    const number latent_heat_of_evaporation,
                                    const number specific_gas_constant,
                                    const number specific_heat_ratio_vapor);
 
    void
    reinit(const number_2    &T_liquid,
           const number_2    &Ma_gas,
           const unsigned int n_active_lanes = dealii::VectorizedArray<number>::size());
 
    [[nodiscard]] number_2
    get_evaporative_mass_flux() const
    {
      return evaporative_mass_flux;
    }
 
    [[nodiscard]] number_2
    get_temperature_jump() const
    {
      return temperature_jump;
    }
 
  private:
    const number atmospheric_pressure;
 
    const number boiling_temperature_at_atmospheric_pressure;
 
    const number latent_heat_of_evaporation;
 
    const number specific_gas_constant;
 
    const number specific_heat_ratio_vapor;
 
    const number temperature_constant;
 
    const number helper_mass_flux;
 
    const number helper_temperature_ratio;
 
    number_2 evaporative_mass_flux{};
 
    number_2 temperature_jump{};
 
    number_2
    compute_temperature_ratio(const number_2 &m_g) const;
 
    number_2
    compute_density_ratio(const number_2 &T_gas_over_T_sat,
                          const number_2 &m_g,
                          const number_2 &exp_m_g_2,
                          const number_2 &erfc_m_g) const;
 
    number_2
    compute_factor_beta(const number_2 &T_gas_over_T_sat,
                        const number_2 &rho_gas_over_rho_sat,
                        const number_2 &m_g,
                        const number_2 &exp_m_g_2) const;
 
    void
    apply_boiling_threshold(const number_2 &T_sat);
  };
} // namespace MeltPoolDG::Evaporation