A class for defining the time derivative of the heat equation. More...

#include <HeatCapacityConductionTimeDerivative.h>

Inheritance diagram for HeatCapacityConductionTimeDerivative:

Public Member Functions
	HeatCapacityConductionTimeDerivative (const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Protected Attributes

const MaterialProperty< Real > &	_heat_capacity
	Heat capacity and its derivatives with respect to temperature and other coupled variables. More...

const MaterialProperty< Real > &	_d_heat_capacity_dT

std::vector< const MaterialProperty< Real > * >	_d_heat_capacity_dargs

Detailed Description

A class for defining the time derivative of the heat equation.

By default this Kernel computes: \( C_p * \frac{\partial T}{\partial t}, \) where \( C_p \) is material property for the "heat_capacity".

Definition at line 30 of file HeatCapacityConductionTimeDerivative.h.

Constructor & Destructor Documentation

◆ HeatCapacityConductionTimeDerivative()

HeatCapacityConductionTimeDerivative::HeatCapacityConductionTimeDerivative ( const InputParameters & parameters )

Definition at line 33 of file HeatCapacityConductionTimeDerivative.C.

   : DerivativeMaterialInterface<JvarMapKernelInterface<TimeDerivative>>(parameters),
     _heat_capacity(getMaterialProperty<Real>("heat_capacity")),
     _d_heat_capacity_dT(getMaterialPropertyDerivative<Real>("heat_capacity", _var.name()))
 {
   // get number of coupled variables
   unsigned int nvar = _coupled_moose_vars.size();
  
   // reserve space for derivatives
   _d_heat_capacity_dargs.resize(nvar);
  
   // iterate over all coupled variables
   for (unsigned int i = 0; i < nvar; ++i)
     _d_heat_capacity_dargs[i] =
         &getMaterialPropertyDerivative<Real>("heat_capacity", _coupled_moose_vars[i]->name());
 }

Member Function Documentation

◆ computeQpJacobian()

Real HeatCapacityConductionTimeDerivative::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 58 of file HeatCapacityConductionTimeDerivative.C.

 {
   // on-diagonal Jacobian with all terms that may depend on the kernel variable
   return _heat_capacity[_qp] * TimeDerivative::computeQpJacobian() +
          _d_heat_capacity_dT[_qp] * _phi[_j][_qp] * TimeDerivative::computeQpResidual();
 }

◆ computeQpOffDiagJacobian()

Real HeatCapacityConductionTimeDerivative::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 66 of file HeatCapacityConductionTimeDerivative.C.

 {
   // get the coupled variable jvar is referring to
   const unsigned int cvar = mapJvarToCvar(jvar);
  
   // off-diagonal contribution with terms that depend on coupled variables
   return (*_d_heat_capacity_dargs[cvar])[_qp] * _phi[_j][_qp] * TimeDerivative::computeQpResidual();
 }

◆ computeQpResidual()

Real HeatCapacityConductionTimeDerivative::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 52 of file HeatCapacityConductionTimeDerivative.C.

 {
   return _heat_capacity[_qp] * TimeDerivative::computeQpResidual();
 }

◆ validParams()

InputParameters HeatCapacityConductionTimeDerivative::validParams ( )

static

Definition at line 17 of file HeatCapacityConductionTimeDerivative.C.

 {
   InputParameters params = TimeDerivative::validParams();
   params.addClassDescription("Time derivative term $C_p \\frac{\\partial T}{\\partial t}$ of "
                              "the heat equation with the heat capacity $C_p$ as an argument.");
  
   // Density may be changing with deformation, so we must integrate
   // over current volume by setting the use_displaced_mesh flag.
   params.set<bool>("use_displaced_mesh") = true;
  
   params.addParam<MaterialPropertyName>(
       "heat_capacity", "heat_capacity", "Property name of the heat capacity material property");
   params.addCoupledVar("args", "Vector of additional arguments of the heat capacity");
   return params;
 }