evaporation_tools.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/vectorization.h>
 
#include <cmath>
 
namespace MeltPoolDG::Evaporation
{
  template <typename number, typename number_2>
  inline number_2
  compute_saturated_gas_pressure(const number_2 &T_sat,
                                 const number    boiling_temperature,
                                 const number    ambient_gas_pressure,
                                 const number    temperature_constant)
  {
    const number T_v = boiling_temperature;
    return ambient_gas_pressure * std::exp(-temperature_constant * (1. / T_sat - 1. / T_v));
  }
 
  template <typename number>
  inline dealii::VectorizedArray<number>
  compute_linear_scaling_coeff(const dealii::VectorizedArray<number> &T,
                               const number                           T_ac,
                               const number                           T_v)
  {
    return compare_and_apply_mask<dealii::SIMDComparison::less_than>(
      T,
      T_ac,
      0.0,
      compare_and_apply_mask<dealii::SIMDComparison::greater_than_or_equal>(
        T, T_v, 1., (T - T_ac) / (T_v - T_ac)));
  }
 
  template <typename number>
  inline number
  compute_linear_scaling_coeff(const number T, const number T_ac, const number T_v)
  {
    return (T < T_ac) ? 0.0 : (T >= T_v) ? 1. : (T - T_ac) / (T_v - T_ac);
  }
} // namespace MeltPoolDG::Evaporation