parameterized_thermal_material.hpp

Go to the documentation of this file.

// Copyright (c) Lawrence Livermore National Security, LLC and
// other Smith Project Developers. See the top-level LICENSE file for
// details.
//
// SPDX-License-Identifier: (BSD-3-Clause)
 
#pragma once
 
#include "smith/physics/materials/thermal_material.hpp"
 
namespace smith::heat_transfer {
 
class ParameterizedLinearIsotropicConductor {
 public:
  ParameterizedLinearIsotropicConductor(double density = 1.0, double specific_heat_capacity = 1.0,
                                        double conductivity_offset = 0.0)
      : density_(density), specific_heat_capacity_(specific_heat_capacity), conductivity_offset_(conductivity_offset)
  {
    SLIC_ERROR_ROOT_IF(conductivity_offset_ < 0.0,
                       "Conductivity must be positive in the linear isotropic conductor material model.");
 
    SLIC_ERROR_ROOT_IF(density_ < 0.0, "Density must be positive in the linear isotropic conductor material model.");
 
    SLIC_ERROR_ROOT_IF(specific_heat_capacity_ < 0.0,
                       "Specific heat capacity must be positive in the linear isotropic conductor material model.");
  }
 
  template <typename T1, typename T2, typename T3, typename T4>
  SMITH_HOST_DEVICE auto operator()(const T1& /* x */, const T2& /* temperature */, const T3& temperature_gradient,
                                    const T4& parameter) const
  {
    return smith::tuple{density_ * specific_heat_capacity_,
                        -1.0 * (conductivity_offset_ + get<0>(parameter)) * temperature_gradient};
  }
 
  static constexpr int numParameters() { return 1; }
 
 private:
  double density_;
 
  double specific_heat_capacity_;
 
  double conductivity_offset_;
};
 
struct ParameterizedIsotropicConductorWithLinearConductivityVsTemperature {
  ParameterizedIsotropicConductorWithLinearConductivityVsTemperature(double density = 1.0,
                                                                     double specific_heat_capacity = 1.0,
                                                                     double conductivity_offset = 1.0,
                                                                     double d_conductivity_d_temperature = 0.0)
      : density_(density),
        specific_heat_capacity_(specific_heat_capacity),
        conductivity_offset_(conductivity_offset),
        d_conductivity_d_temperature_(d_conductivity_d_temperature)
  {
    SLIC_ERROR_ROOT_IF(density_ < 0.0, "Density must be positive in the linear isotropic conductor material model.");
 
    SLIC_ERROR_ROOT_IF(specific_heat_capacity_ < 0.0,
                       "Specific heat capacity must be positive in the linear isotropic conductor material model.");
  }
 
  template <typename T1, typename T2, typename T3, typename T4>
  SMITH_HOST_DEVICE auto operator()(const T1& /* x */, const T2& temperature, const T3& temperature_gradient,
                                    const T4& parameter) const
  {
    const auto currentConductivity =
        conductivity_offset_ + get<0>(parameter) + d_conductivity_d_temperature_ * temperature;
#if defined(SMITH_USE_CUDA) || defined(SMITH_USE_HIP)
    assert(smith::get_value(currentConductivity) >= 0.0);
#else
    SLIC_ERROR_ROOT_IF(
        smith::get_value(currentConductivity) < 0.0,
        "Conductivity in the IsotropicConductorWithLinearConductivityVsTemperature model has gone negative.");
#endif
    return smith::tuple{density_ * specific_heat_capacity_, -1.0 * currentConductivity * temperature_gradient};
  }
 
  static constexpr int numParameters() { return 1; }
 
 private:
  double density_;
 
  double specific_heat_capacity_;
 
  double conductivity_offset_;
 
  double d_conductivity_d_temperature_;
};
 
struct ParameterizedSource {
  double source_offset_ = 0.0;
 
  template <typename T1, typename T2, typename T3, typename T4>
  SMITH_HOST_DEVICE auto operator()(const T1& /* x */, const double /* time */, const T2& /* temperature */,
                                    const T3& /* temperature_gradient */, const T4& parameter) const
  {
    return source_offset_ + parameter;
  }
 
  static constexpr int numParameters() { return 1; }
};
 
struct ParameterizedFlux {
  double flux_offset_ = 0.0;
 
  template <typename T1, typename T2, typename T3, typename T4>
  SMITH_HOST_DEVICE auto operator()(const T1& /* x */, const T2& /* normal */, const double /* time */,
                                    const T3& /* temperature */, const T4& parameter) const
  {
    return flux_offset_ + parameter;
  }
 
  static constexpr int numParameters() { return 1; }
};
 
}  // namespace smith::heat_transfer