Provides a heat source (J/m^3/s) from plastic deformation: (1 - porosity) * coeff * stress * plastic_strain_rate. More...

#include <PorousFlowPlasticHeatEnergy.h>

Inheritance diagram for PorousFlowPlasticHeatEnergy:

[legend]

Public Member Functions
	PorousFlowPlasticHeatEnergy (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 PorousFlowDictator &	_dictator
	PorousFlowDictator UserObject. More...

const bool	_strain_at_nearest_qp
	Whether the porosity uses the volumetric strain at the closest quadpoint. More...

const MaterialProperty< unsigned int > *const	_nearest_qp
	The nearest qp to the node. More...

const MaterialProperty< Real > &	_porosity
	Porosity at the nodes, but it can depend on grad(variables) which are actually evaluated at the qps. More...

const MaterialProperty< std::vector< Real > > &	_dporosity_dvar
	d(porosity)/d(PorousFlow variable) - these derivatives will be wrt variables at the nodes More...

const MaterialProperty< std::vector< RealGradient > > &	_dporosity_dgradvar
	d(porosity)/d(grad PorousFlow variable) - remember these derivatives will be wrt grad(vars) at qps More...

Real	_coeff
	coefficient of stress * plastic_strain_rate More...

const std::string	_base_name
	optional parameter that allows multiple mechanics models to be defined More...

const MaterialProperty< Real > &	_plastic_heat
	stress * plastic_strain_rate More...

const MaterialProperty< RankTwoTensor > &	_dplastic_heat_dstrain
	d(plastic_heat)/d(total_strain) More...

unsigned int	_ndisp
	umber of coupled displacement variables More...

std::vector< unsigned int >	_disp_var
	MOOSE variable number for the displacement variables. More...

Detailed Description

Provides a heat source (J/m^3/s) from plastic deformation: (1 - porosity) * coeff * stress * plastic_strain_rate.

Definition at line 24 of file PorousFlowPlasticHeatEnergy.h.

Constructor & Destructor Documentation

◆ PorousFlowPlasticHeatEnergy()

PorousFlowPlasticHeatEnergy::PorousFlowPlasticHeatEnergy ( const InputParameters & parameters )

Definition at line 36 of file PorousFlowPlasticHeatEnergy.C.

   : PlasticHeatEnergy(parameters),
     _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
     _strain_at_nearest_qp(getParam<bool>("strain_at_nearest_qp")),
     _nearest_qp(_strain_at_nearest_qp
                     ? &getMaterialProperty<unsigned int>("PorousFlow_nearestqp_nodal")
                     : nullptr),
     _porosity(getMaterialProperty<Real>("PorousFlow_porosity_nodal")),
     _dporosity_dvar(getMaterialProperty<std::vector<Real>>("dPorousFlow_porosity_nodal_dvar")),
     _dporosity_dgradvar(
         getMaterialProperty<std::vector<RealGradient>>("dPorousFlow_porosity_nodal_dgradvar"))
 {
 }

Member Function Documentation

◆ computeQpJacobian()

Real PorousFlowPlasticHeatEnergy::computeQpJacobian ( )

overrideprotectedvirtual

Reimplemented from PlasticHeatEnergy.

Definition at line 57 of file PorousFlowPlasticHeatEnergy.C.

 {
   return computeQpOffDiagJacobian(_var.number());
 }

◆ computeQpOffDiagJacobian()

Real PorousFlowPlasticHeatEnergy::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Reimplemented from PlasticHeatEnergy.

Definition at line 63 of file PorousFlowPlasticHeatEnergy.C.

 {
   // If the variable is not a PorousFlow variable, the Jacobian terms are 0
   if (_dictator.notPorousFlowVariable(jvar))
     return 0.0;
  
   const Real res_no_porosity = PlasticHeatEnergy::computeQpResidual();
   const Real jac_no_porosity = PlasticHeatEnergy::computeQpOffDiagJacobian(jvar);
  
   const unsigned pvar = _dictator.porousFlowVariableNum(jvar);
   const unsigned nearest_qp = (_strain_at_nearest_qp ? (*_nearest_qp)[_i] : _i);
  
   Real jac = (1.0 - _porosity[_i]) * jac_no_porosity -
              _dporosity_dgradvar[_i][pvar] * _grad_phi[_j][nearest_qp] * res_no_porosity;
   if (_i != _j)
     return jac;
  
   return jac - _dporosity_dvar[_i][pvar] * res_no_porosity;
 }

Referenced by computeQpJacobian().

◆ computeQpResidual()

Real PorousFlowPlasticHeatEnergy::computeQpResidual ( )

overrideprotectedvirtual

Reimplemented from PlasticHeatEnergy.

Definition at line 51 of file PorousFlowPlasticHeatEnergy.C.

 {
   return (1.0 - _porosity[_i]) * PlasticHeatEnergy::computeQpResidual();
 }

◆ validParams()

InputParameters PlasticHeatEnergy::validParams ( )

staticinherited

Definition at line 19 of file PlasticHeatEnergy.C.

 {
   InputParameters params = Kernel::validParams();
   params.addClassDescription("Plastic heat energy density = coeff * stress * plastic_strain_rate");
   params.addRequiredCoupledVar("displacements",
                                "The string of displacements suitable for the problem statement");
   params.addParam<std::string>("base_name", "Material property base name");
   params.addParam<Real>("coeff", 1.0, "Heat energy density = coeff * stress * plastic_strain_rate");
   return params;
 }