Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "ThermalSensitivity.h"
      11             : 
      12             : registerMooseObject("HeatTransferApp", ThermalSensitivity);
      13             : 
      14             : InputParameters
      15          85 : ThermalSensitivity::validParams()
      16             : {
      17          85 :   InputParameters params = Material::validParams();
      18          85 :   params.addClassDescription("Computes cost sensitivity needed for multimaterial SIMP method.");
      19         170 :   params.addRequiredCoupledVar("design_density", "Design density variable name.");
      20         170 :   params.addRequiredCoupledVar("temperature", "temperature");
      21         170 :   params.addRequiredParam<MaterialPropertyName>("thermal_conductivity",
      22             :                                                 "DerivativeParsedMaterial for cost of materials.");
      23          85 :   return params;
      24           0 : }
      25             : 
      26          66 : ThermalSensitivity::ThermalSensitivity(const InputParameters & parameters)
      27             :   : DerivativeMaterialInterface<Material>(parameters),
      28          66 :     _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
      29          66 :     _sensitivity(declareProperty<Real>(_base_name + "thermal_sensitivity")),
      30          66 :     _design_density(coupledValue("design_density")),
      31          66 :     _design_density_name(coupledName("design_density", 0)),
      32          66 :     _grad_temperature(coupledGradient("temperature")),
      33          66 :     _thermal_conductivity(
      34         132 :         getMaterialPropertyByName<Real>(getParam<MaterialPropertyName>("thermal_conductivity"))),
      35         132 :     _dTdp(getMaterialPropertyDerivativeByName<Real>(
      36          66 :         getParam<MaterialPropertyName>("thermal_conductivity"), _design_density_name))
      37             : {
      38          66 : }
      39             : 
      40             : void
      41      524272 : ThermalSensitivity::computeQpProperties()
      42             : {
      43             :   const Real thermal_compliance =
      44      524272 :       0.5 * _thermal_conductivity[_qp] * _grad_temperature[_qp] * _grad_temperature[_qp];
      45      524272 :   _sensitivity[_qp] = -_dTdp[_qp] * thermal_compliance / _thermal_conductivity[_qp];
      46      524272 : }