material_data.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/parameter_handler.h>
 
#include <meltpooldg/utilities/enum.hpp>
#include <meltpooldg/utilities/numbers.hpp>
 
#include <string>
 
namespace MeltPoolDG
{
  BETTER_ENUM(
    MaterialTemplate,
    char,
    none,            // all material parameters must be specified
    stainless_steel, // melt and solid material parameters will be set to stainless steel values
    Ti64,            // melt and solid material parameters will be set to Ti-6Al-4V values
    Ti64Benchmark    // melt and solid material parameters will be set to Ti-6Al-4V values for the
                     // benchmark paper together with PolyMTL.
  )
  BETTER_ENUM(SolidLiquidPropertiesTransitionType,
              char,
              mushy_zone, // solid and liquid properties are interpolated between the solidus and
                          // liquidus temperature
              sharp // solid and liquid properties jump at the melting point, which is set via the
                    // solidus temperature
  )
  BETTER_ENUM(TwoPhaseFluidPropertiesTransitionType,
              char,
              sharp,  // properties jump at heaviside = 0.5
              smooth, // properties are smeared between the phases proportional to the heaviside
              consistent_with_evaporation //  same as "smooth", but the density is interpolated
                                          //  proportional by the harmonic mean.
                                          // consistent with the evaporation formulation
  )
 
  
  BETTER_ENUM(ApparentCapacityType,
              char,
 
              poly4_bell,
 
              constant,
 
              qlq)
 
  template <typename number>
  struct MaterialPhaseData
  {
    number thermal_conductivity   = 0.0;
    number specific_heat_capacity = 0.0;
    number density                = 0.0;
    number dynamic_viscosity      = 0.0;
 
    void
    add_parameters(dealii::ParameterHandler &prm, const std::string &phase_name);
  };
 
  template <typename number>
  struct MaterialData
  {
  private:
    MaterialTemplate material_template = MaterialTemplate::none;
 
  public:
    MaterialPhaseData<number> gas;
 
    MaterialPhaseData<number> liquid;
 
    MaterialPhaseData<number> solid;
 
    number solidus_temperature        = 0.0;
    number liquidus_temperature       = 0.0;
    number latent_heat_of_fusion      = 0.0;
    number boiling_temperature        = 0.0;
    number latent_heat_of_evaporation = 0.0;
    number molar_mass                 = 0.0;
 
    number specific_enthalpy_reference_temperature = numbers::invalid_double;
 
    ApparentCapacityType apparent_capacity_type = ApparentCapacityType::qlq;
 
    SolidLiquidPropertiesTransitionType solid_liquid_properties_transition_type =
      SolidLiquidPropertiesTransitionType::mushy_zone;
    TwoPhaseFluidPropertiesTransitionType two_phase_fluid_properties_transition_type =
      TwoPhaseFluidPropertiesTransitionType::smooth;
 
    void
    add_parameters(dealii::ParameterHandler &prm);
 
    void
    check_parameters_heat_transfer(const bool do_two_phase, const bool do_solidification) const;
 
    static MaterialData<number>
    create_stainless_steel_material_data();
 
    static MaterialData<number>
    create_Ti64_material_data();
 
    static MaterialData<number>
    create_Ti64_benchmark_material_data();
  };
} // namespace MeltPoolDG