types.hh

Go to the documentation of this file.

/* SPDX-FileCopyrightText: Copyright (c) 2016 - 2025, the adamantine authors.
 * SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
 */
 
#ifndef TYPES_HH
#define TYPES_HH
 
#include <Kokkos_NumericTraits.hpp>
 
#include <array>
#include <string>
#include <unordered_map>
 
namespace dealii
{
namespace LinearAlgebra
{
}
 
namespace LA = LinearAlgebra;
} // namespace dealii
 
namespace adamantine
{
 
static std::unordered_map<std::string, double> g_unit_scaling_factor{
    {"millimeter", 1e-3},        {"centimeter", 1e-2},
    {"inch", 2.54e-2},           {"meter", 1.},
    {"milliwatt", 1e-3},         {"watt", 1.},
    {"millimeter/second", 1e-3}, {"centimer/second", 1e-2},
    {"meter/second", 1.}};
 
enum class StateProperty
{
  density,
  specific_heat,
  thermal_conductivity_x,
  thermal_conductivity_y,
  thermal_conductivity_z,
  emissivity,
  radiation_heat_transfer_coef,
  convection_heat_transfer_coef,
  // Mechanical properties only make sense for the solid state. They do not make
  // sense for the powder and the liquid state
  lame_first_parameter,
  lame_second_parameter,
  thermal_expansion_coef,
  // Density is used both by the thermal and the mechanical simulation. We need
  // to duplicate the property because the mechanical code cannot run on the
  // GPU.
  density_s,
  plastic_modulus,
  isotropic_hardening,
  elastic_limit,
  SIZE,
  SIZE_MECHANICAL = 7
};
 
static unsigned int constexpr g_n_state_properties =
    static_cast<unsigned int>(StateProperty::SIZE);
 
static unsigned int constexpr g_n_mechanical_state_properties =
    static_cast<unsigned int>(StateProperty::SIZE_MECHANICAL);
 
static unsigned int constexpr g_n_thermal_state_properties =
    g_n_state_properties - g_n_mechanical_state_properties;
 
enum class Property
{
  liquidus,
  solidus,
  latent_heat,
  radiation_temperature_infty,
  convection_temperature_infty,
  SIZE
};
 
static unsigned int constexpr g_n_properties =
    static_cast<unsigned int>(Property::SIZE);
 
static std::array<std::string, 3> const material_state_names = {
    {"solid", "liquid", "powder"}};
 
static std::array<std::string, 5> const property_names = {
    {"liquidus", "solidus", "latent_heat", "radiation_temperature_infty",
     "convection_temperature_infty"}};
 
static std::array<std::string, 15> const state_property_names = {
    {"density", "specific_heat", "thermal_conductivity_x",
     "thermal_conductivity_y", "thermal_conductivity_z", "emissivity",
     "radiation_heat_transfer_coef", "convection_heat_transfer_coef",
     "lame_first_parameter", "lame_second_parameter", "thermal_expansion_coef",
     "density", "plastic_modulus", "isotropic_hardening", "elastic_limit"}};
 
enum Timing
{
  main,
  refine,
  add_material_search,
  add_material_activate,
  da_experimental_data,
  da_dof_mapping,
  da_covariance_sparsity,
  da_obs_covariance,
  da_update_ensemble,
  evol_time,
  evol_time_eval_th_ph,
  evol_time_update_bound_mat_prop,
  output,
  n_timers
};
 
struct Constant
{
  static double constexpr stefan_boltzmann = 5.670374419e-8;
};
 
template <int dim>
struct axis;
 
// dim == 2 specialization
template <>
struct axis<2>
{
  static unsigned int constexpr x = 0;
  static unsigned int constexpr y =
      Kokkos::Experimental::finite_max_v<unsigned int>;
  static unsigned int constexpr z = 1;
};
 
// dim == 3 specialization
template <>
struct axis<3>
{
  static unsigned int constexpr x = 0;
  static unsigned int constexpr y = 1;
  static unsigned int constexpr z = 2;
};
 
} // namespace adamantine
 
#endif