evaporative_cooling.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/vectorization.h>
 
#include <meltpooldg/core/material_data.hpp>
#include <meltpooldg/phase_change/evaporation_data.hpp>
#include <meltpooldg/phase_change/evaporation_model_base.hpp>
 
#include <memory>
 
 
namespace MeltPoolDG::Evaporation
{
  template <typename number>
  class EvaporativeCooling
  {
  public:
    EvaporativeCooling(const EvaporationData<number> &evapor_data,
                       const MaterialData<number>    &material_data,
                       const bool                     setup_internal_mass_flux_operator = false);
 
 
    template <typename ValueType>
    inline ValueType
    compute_evaporative_cooling(const ValueType                  &mass_flux,
                                [[maybe_unused]] const ValueType &temperature) const;
 
 
    inline number
    compute_evaporative_cooling(const number temperature) const;
 
    inline dealii::VectorizedArray<number>
    compute_evaporative_cooling(const dealii::VectorizedArray<number> &temperature) const;
 
 
    template <typename ValueType>
    inline ValueType
    compute_evaporative_cooling_derivative_constant_mass_flux(
      [[maybe_unused]] const ValueType &mass_flux) const;
 
 
    inline number
    compute_evaporative_cooling_derivative_with_temperature_dependent_mass_flux(
      const number temperature) const;
 
    inline dealii::VectorizedArray<number>
    compute_evaporative_cooling_derivative_with_temperature_dependent_mass_flux(
      const dealii::VectorizedArray<number> &temperature) const;
 
  private:
    template <typename ValueType>
    inline ValueType
    compute_phenomenological_specific_enthalpy(const ValueType &temperature) const;
 
    const bool                                    do_phenomenological_recoil_pressure;
    const number                                  latent_heat_of_evaporation;
    const number                                  specific_heat_capacity;
    const number                                  specific_enthalpy_reference_temperature;
    const number                                  boiling_temperature;
    number                                        activation_temperature;
    bool                                          ramp_enabled;
    number                                        activation_ramp_derivative = 0;
    std::unique_ptr<EvaporationModelBase<number>> mass_flux_operator;
  };
} // namespace MeltPoolDG::Evaporation