ThermalCompositeSiCProperties.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 "ThermalCompositeSiCProperties.h"
#include "libmesh/utility.h"

registerMooseObject("SolidPropertiesApp", ThermalCompositeSiCProperties);

InputParameters
ThermalCompositeSiCProperties::validParams()
{
  InputParameters params = ThermalSolidProperties::validParams();

  params.addRangeCheckedParam<Real>("density", 3216.0, "density > 0.0", "(Constant) density");
  params.addClassDescription("Composite silicon carbide thermal properties.");
  return params;
}

ThermalCompositeSiCProperties::ThermalCompositeSiCProperties(const InputParameters & parameters)
  : ThermalSolidProperties(parameters),
    _rho_const(getParam<Real>("density")),
    _c1(925.65),
    _c2(0.3772),
    _c3(7.9259e-5),
    _c4(3.1946e7)
{
}

Real
ThermalCompositeSiCProperties::cp_from_T(const Real & T) const
{
  return _c1 + _c2 * T - _c3 * Utility::pow<2>(T) - _c4 / Utility::pow<2>(T);
}

void
ThermalCompositeSiCProperties::cp_from_T(const Real & T, Real & cp, Real & dcp_dT) const
{
  cp = cp_from_T(T);
  dcp_dT = _c2 - 2.0 * _c3 * T + 2.0 * _c4 / Utility::pow<3>(T);
}

Real
ThermalCompositeSiCProperties::cp_integral(const Real & T) const
{
  return _c1 * T + 0.5 * _c2 * Utility::pow<2>(T) - _c3 / 3.0 * Utility::pow<3>(T) + _c4 / T;
}

Real
ThermalCompositeSiCProperties::k_from_T(const Real & T) const
{
  return -1.71e-11 * Utility::pow<4>(T) + 7.35e-8 * Utility::pow<3>(T) -
         1.10e-4 * Utility::pow<2>(T) + 0.061 * T + 7.97;
}

void
ThermalCompositeSiCProperties::k_from_T(const Real & T, Real & k, Real & dk_dT) const
{
  k = k_from_T(T);
  dk_dT = -6.84e-11 * Utility::pow<3>(T) + 2.205e-7 * Utility::pow<2>(T) - 2.2e-4 * T + 0.061;
}

Real
ThermalCompositeSiCProperties::rho_from_T(const Real & /* T */) const
{
  return _rho_const;
}

void
ThermalCompositeSiCProperties::rho_from_T(const Real & T, Real & rho, Real & drho_dT) const
{
  rho = rho_from_T(T);
  drho_dT = 0.0;
}