Calculates thermal conductivity and specific heat of the material. More...

#include <AnisoHeatConductionMaterial.h>

Inheritance diagram for AnisoHeatConductionMaterialTempl< is_ad >:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	AnisoHeatConductionMaterialTempl (const InputParameters &parameters)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

void	validateDerivativeMaterialPropertyBase (const std::string &base)

const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	initQpStatefulProperties () override

virtual void	computeQpProperties () override

Protected Attributes
const unsigned int	_dim

const Real	_ref_temp

const bool	_has_temp

const GenericVariableValue< is_ad > &	_T

const unsigned int	_T_var

const VariableName	_T_name

const std::string	_base_name

const RankTwoTensor	_user_provided_thermal_conductivity

GenericMaterialProperty< RankTwoTensor, is_ad > &	_thermal_conductivity

MaterialProperty< RankTwoTensor > &	_dthermal_conductivity_dT

const Function *const	_thermal_conductivity_temperature_coefficient_function

GenericMaterialProperty< Real, is_ad > &	_specific_heat

MaterialProperty< Real > &	_dspecific_heat_dT

const Function *const	_specific_heat_function

const MooseArray< ADPoint > *	_ad_q_point

Private Member Functions
auto	genericQPoints ()

auto	genericQPoints ()

auto	genericQPoints ()

bool	haveMaterialProperty (const std::string &prop_name)

std::vector< VariableName >	buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)

void	validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)

bool	isNotObjectVariable (const VariableName &name)

Private Attributes
FEProblemBase &	_dmi_fe_problem

Detailed Description

template<bool is_ad>
class AnisoHeatConductionMaterialTempl< is_ad >

Calculates thermal conductivity and specific heat of the material.

Definition at line 21 of file AnisoHeatConductionMaterial.h.

Constructor & Destructor Documentation

◆ AnisoHeatConductionMaterialTempl()

template<bool is_ad>

AnisoHeatConductionMaterialTempl< is_ad >::AnisoHeatConductionMaterialTempl ( const InputParameters & parameters )

Definition at line 44 of file AnisoHeatConductionMaterial.C.

  : DerivativeMaterialInterface<Material>(parameters),
    _dim(_subproblem.mesh().dimension()),
    _ref_temp(getParam<Real>("reference_temperature")),
 
    _has_temp(isCoupled("temperature")),
    _T(coupledGenericValue<is_ad>("temperature")),
    _T_var(coupled("temperature")),
    _T_name(coupledName("temperature", 0)),
    _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
 
    _user_provided_thermal_conductivity(getParam<std::vector<Real>>("thermal_conductivity")),
    _thermal_conductivity(
        declareGenericProperty<RankTwoTensor, is_ad>(_base_name + "thermal_conductivity")),
    _dthermal_conductivity_dT(
        declarePropertyDerivative<RankTwoTensor>(_base_name + "thermal_conductivity", _T_name)),
    _thermal_conductivity_temperature_coefficient_function(
        getParam<FunctionName>("thermal_conductivity_temperature_coefficient_function") != ""
            ? &getFunction("thermal_conductivity_temperature_coefficient_function")
            : nullptr),
 
    _specific_heat(declareGenericProperty<Real, is_ad>(_base_name + "specific_heat")),
    _dspecific_heat_dT(declarePropertyDerivative<Real>(_base_name + "specific_heat", _T_name)),
    _specific_heat_function(&getFunction("specific_heat")),
    _ad_q_point(is_ad ? &_assembly.adQPoints() : nullptr)
{
}

Member Function Documentation

◆ computeQpProperties()

template<bool is_ad>

void AnisoHeatConductionMaterialTempl< is_ad >::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 97 of file AnisoHeatConductionMaterial.C.

{
  const auto & temp_qp = _T[_qp];
  const auto & p = genericQPoints();
 
  if (_thermal_conductivity_temperature_coefficient_function)
  {
 
    _thermal_conductivity[_qp] =
        _user_provided_thermal_conductivity *
        (1.0 + _thermal_conductivity_temperature_coefficient_function->value(temp_qp, p) *
                   (temp_qp - _ref_temp));
    _dthermal_conductivity_dT[_qp] =
        _user_provided_thermal_conductivity *
        _thermal_conductivity_temperature_coefficient_function->timeDerivative(
            MetaPhysicL::raw_value(temp_qp), _q_point[_qp]) *
        MetaPhysicL::raw_value(temp_qp - _ref_temp);
  }
  else
    _thermal_conductivity[_qp] = _user_provided_thermal_conductivity;
 
  _specific_heat[_qp] = _specific_heat_function->value(temp_qp, p);
  _dspecific_heat_dT[_qp] =
      _specific_heat_function->timeDerivative(MetaPhysicL::raw_value(temp_qp), _q_point[_qp]);
}

◆ genericQPoints() [1/3]

template<bool is_ad>

auto AnisoHeatConductionMaterialTempl< is_ad >::genericQPoints ( )

private

◆ genericQPoints() [2/3]

auto AnisoHeatConductionMaterialTempl< true >::genericQPoints ( )

private

Definition at line 83 of file AnisoHeatConductionMaterial.C.

{
  return (*_ad_q_point)[_qp];
}

◆ genericQPoints() [3/3]

auto AnisoHeatConductionMaterialTempl< false >::genericQPoints ( )

private

Definition at line 90 of file AnisoHeatConductionMaterial.C.

{
  return _q_point[_qp];
}

◆ initQpStatefulProperties()

template<bool is_ad>

void AnisoHeatConductionMaterialTempl< is_ad >::initQpStatefulProperties ( )

overrideprotectedvirtual

Definition at line 75 of file AnisoHeatConductionMaterial.C.

{
  _thermal_conductivity[_qp] = _user_provided_thermal_conductivity;
  DerivativeMaterialInterface::initQpStatefulProperties();
}

◆ validParams()

template<bool is_ad>

InputParameters AnisoHeatConductionMaterialTempl< is_ad >::validParams ( )

static

Definition at line 23 of file AnisoHeatConductionMaterial.C.

{
  InputParameters params = Material::validParams();
 
  params.addCoupledVar("temperature", "Coupled variable for temperature.");
  params.addParam<Real>(
      "reference_temperature", 293.0, "Reference temperature for thermal conductivity in Kelvin.");
  params.addParam<std::string>("base_name", "Material property base name.");
  params.addRequiredParam<std::vector<Real>>("thermal_conductivity",
                                             "The thermal conductivity tensor values");
  params.addParam<FunctionName>(
      "thermal_conductivity_temperature_coefficient_function",
      "",
      "Temperature coefficient for thermal conductivity as a function of temperature.");
  params.addRequiredParam<FunctionName>("specific_heat",
                                        "Specific heat as a function of temperature.");
  params.addClassDescription("General-purpose material model for anisotropic heat conduction");
  return params;
}