www.mooseframework.org
AnisoHeatConduction.C
Go to the documentation of this file.
1 //* This file is part of the MOOSE framework
2 //* https://www.mooseframework.org
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 "AnisoHeatConduction.h"
11 
12 registerMooseObject("HeatTransferApp", AnisoHeatConduction);
13 
16 {
18  params.addClassDescription(
19  "Anisotropic HeatConduction kernel $\\nabla \\cdot -\\widetilde{k} \\nabla u$ "
20  "with weak form given by $(\\nabla \\psi_i, \\widetilde{k} \\nabla u)$.");
21  params.addParam<MaterialPropertyName>(
22  "thermal_conductivity",
23  "thermal_conductivity",
24  "Material property providing thermal conductivity of the material.");
25 
26  return params;
27 }
28 
30  : DerivativeMaterialInterface<Kernel>(parameters),
31  _k(getMaterialProperty<RankTwoTensor>("thermal_conductivity")),
32  _dk_dT(getMaterialPropertyDerivative<RankTwoTensor>("thermal_conductivity", _var.name()))
33 {
34 }
35 
36 Real
38 {
39  return _k[_qp] * _grad_u[_qp] * _grad_test[_i][_qp];
40 }
41 
42 Real
44 {
45  return _k[_qp] * _grad_phi[_j][_qp] * _grad_test[_i][_qp] +
46  _dk_dT[_qp] * _phi[_j][_qp] * _grad_u[_qp] * _grad_test[_i][_qp];
47 }
virtual Real computeQpJacobian() override
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
const MaterialProperty< RankTwoTensor > & _k
AnisoHeatConduction(const InputParameters &parameters)
registerMooseObject("HeatTransferApp", AnisoHeatConduction)
InputParameters validParams()
const std::string name
Definition: Setup.h:20
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
void addClassDescription(const std::string &doc_string)
virtual Real computeQpResidual() override
This kernel implements the Laplacian operator multiplied by a 2nd order tensor giving anisotropic (di...
const MaterialProperty< RankTwoTensor > & _dk_dT
static InputParameters validParams()