ThermalConductivity.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "ThermalConductivity.h"

registerMooseObject("HeatTransferApp", ThermalConductivity);

InputParameters
ThermalConductivity::validParams()
{
  InputParameters params = SideAverageValue::validParams();
  params.addRequiredParam<Real>("dx", "Length between sides of sample in length_scale");
  params.addRequiredParam<PostprocessorName>(
      "flux", "Heat flux out of 'cold' boundary in solution units, should always be positive");
  params.addRequiredParam<PostprocessorName>("T_hot", "Temperature on 'hot' boundary in K");
  params.addParam<Real>("length_scale", 1e-8, "Length scale of the solution, default is 1e-8");
  params.addParam<Real>("k0", 0.0, "Initial value of the thermal conductivity");
  params.addClassDescription("Computes the effective thermal conductivity averaged on a boundary.");
  return params;
}

ThermalConductivity::ThermalConductivity(const InputParameters & parameters)
  : SideAverageValue(parameters),
    _dx(getParam<Real>("dx")),
    _flux(getPostprocessorValue("flux")),
    _T_hot(getPostprocessorValue("T_hot")),
    _length_scale(getParam<Real>("length_scale")),
    _k0(getParam<Real>("k0")),
    _step_zero(declareRestartableData<bool>("step_zero", true)),
    _value(0.0)
{
}

void
ThermalConductivity::finalize()
{
  SideAverageValue::finalize();

  const Real T_cold = _integral_value / _volume;
  Real Th_cond = 0.0;
  if (_t_step >= 1)
    _step_zero = false;

  // Calculate effective thermal conductivity in W/(length_scale-K)
  if (std::abs(_T_hot - T_cold) > 1.0e-20)
    Th_cond = std::abs(_flux) * _dx / std::abs(_T_hot - T_cold);

  if (_step_zero)
    _value = _k0;
  else
    _value = Th_cond / _length_scale; // In W/(m-K)
}

Real
ThermalConductivity::getValue() const
{
  return _value;
}