Kernel class for peridynamic heat conduction models. More...

#include <HeatConductionBPD.h>

Inheritance diagram for HeatConductionBPD:

Public Member Functions
	HeatConductionBPD (const InputParameters &parameters)

virtual void	computeResidual () override

virtual void	computeJacobian () override

virtual Real	computeQpResidual () override

Protected Member Functions
virtual void	computeLocalResidual () override
	Function to compute local contribution to the residual at the current nodes. More...

virtual void	computeLocalJacobian () override
	Function to compute local contribution to the diagonal Jacobian at the current nodes. More...

virtual void	computeNonlocalResidual ()
	Function to compute nonlocal contribution to the residual at the current nodes. More...

virtual void	prepare ()
	Function to precalculate data which will be used in the derived classes. More...

Protected Attributes
MooseVariable *	_bond_status_var
	Bond_status variable. More...

const bool	_use_full_jacobian
	Option to use full jacobian including nonlocal constribution or not. More...

RealGradient	_origin_vec_ij
	Vector for current bond under undefored configuration. More...

Real	_bond_status_ij
	Bond status of current bond/edge2. More...


const MaterialProperty< Real > &	_bond_heat_flow
	Bond based material properties. More...

const MaterialProperty< Real > &	_bond_dQdT


PeridynamicsMesh &	_pdmesh
	Parameters for peridynamic mesh information. More...

const unsigned int	_dim

const unsigned int	_nnodes

std::vector< Real >	_vols_ij

std::vector< Real >	_dg_vol_frac_ij

std::vector< Real >	_horiz_rad

Detailed Description

Kernel class for peridynamic heat conduction models.

Definition at line 22 of file HeatConductionBPD.h.

Constructor & Destructor Documentation

◆ HeatConductionBPD()

HeatConductionBPD::HeatConductionBPD ( const InputParameters & parameters )

Definition at line 25 of file HeatConductionBPD.C.

   : PeridynamicsKernelBase(parameters),
     _bond_heat_flow(getMaterialProperty<Real>("bond_heat_flow")),
     _bond_dQdT(getMaterialProperty<Real>("bond_dQdT"))
 {
 }

Member Function Documentation

◆ computeJacobian()

void PeridynamicsKernelBase::computeJacobian ( )

overridevirtualinherited

Definition at line 87 of file PeridynamicsKernelBase.C.

 {
   prepare();
  
   DenseMatrix<Number> & ke = _assembly.jacobianBlock(_var.number(), _var.number());
   _local_ke.resize(ke.m(), ke.n());
   _local_ke.zero();
  
   computeLocalJacobian();
  
   ke += _local_ke;
  
   if (_has_diag_save_in)
   {
     unsigned int rows = ke.m();
     DenseVector<Number> diag(rows);
     for (unsigned int i = 0; i < rows; ++i)
       diag(i) = _local_ke(i, i);
  
     Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
     for (unsigned int i = 0; i < _diag_save_in.size(); ++i)
       _diag_save_in[i]->sys().solution().add_vector(diag, _diag_save_in[i]->dofIndices());
   }
  
   _local_ke.zero();
  
   if (_use_full_jacobian)
     computeNonlocalJacobian();
 }

◆ computeLocalJacobian()

void HeatConductionBPD::computeLocalJacobian ( )

overrideprotectedvirtual

Function to compute local contribution to the diagonal Jacobian at the current nodes.

Reimplemented from PeridynamicsKernelBase.

Definition at line 40 of file HeatConductionBPD.C.

 {
   for (_i = 0; _i < _test.size(); ++_i)
     for (_j = 0; _j < _phi.size(); ++_j)
       _local_ke(_i, _j) += -1 * (_i == _j ? 1 : -1) * _bond_dQdT[0] * _bond_status_ij;
 }

◆ computeLocalResidual()

void HeatConductionBPD::computeLocalResidual ( )

overrideprotectedvirtual

Function to compute local contribution to the residual at the current nodes.

Implements PeridynamicsKernelBase.

Definition at line 33 of file HeatConductionBPD.C.

 {
   _local_re(0) = -_bond_heat_flow[0] * _bond_status_ij;
   _local_re(1) = -_local_re(0);
 }

◆ computeNonlocalResidual()

virtual void PeridynamicsKernelBase::computeNonlocalResidual ( )

inlineprotectedvirtualinherited

Function to compute nonlocal contribution to the residual at the current nodes.

Definition at line 43 of file PeridynamicsKernelBase.h.

43 {};

Referenced by PeridynamicsKernelBase::computeResidual().

◆ computeQpResidual()

virtual Real PeridynamicsKernelBase::computeQpResidual ( )

inlineoverridevirtualinherited

Definition at line 32 of file PeridynamicsKernelBase.h.

32 { return 0.0; }

◆ computeResidual()

void PeridynamicsKernelBase::computeResidual ( )

overridevirtualinherited

Definition at line 60 of file PeridynamicsKernelBase.C.

 {
   prepare();
  
   DenseVector<Number> & re = _assembly.residualBlock(_var.number());
   mooseAssert(re.size() == _nnodes,
               "PeridynamicsKernelBase is designed to only work with EDGE2 elements");
   _local_re.resize(re.size());
   _local_re.zero();
  
   computeLocalResidual();
  
   re += _local_re;
  
   if (_has_save_in)
   {
     Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
     for (unsigned int i = 0; i < _save_in.size(); ++i)
       _save_in[i]->sys().solution().add_vector(_local_re, _save_in[i]->dofIndices());
   }
  
   _local_re.zero();
  
   computeNonlocalResidual();
 }

◆ prepare()

void PeridynamicsKernelBase::prepare ( )

protectedvirtualinherited

Function to precalculate data which will be used in the derived classes.

Definition at line 43 of file PeridynamicsKernelBase.C.

 {
   for (unsigned int nd = 0; nd < _nnodes; ++nd)
   {
     _vols_ij[nd] = _pdmesh.getPDNodeVolume(_current_elem->node_id(nd));
     dof_id_type id_ji_in_ij =
         _pdmesh.getNeighborIndex(_current_elem->node_id(nd), _current_elem->node_id(1 - nd));
     _dg_vol_frac_ij[nd] = _pdmesh.getDefGradVolFraction(_current_elem->node_id(nd), id_ji_in_ij);
     _horiz_rad[nd] = _pdmesh.getHorizon(_current_elem->node_id(nd));
   }
  
   _origin_vec_ij = _pdmesh.getPDNodeCoord(_current_elem->node_id(1)) -
                    _pdmesh.getPDNodeCoord(_current_elem->node_id(0));
   _bond_status_ij = _bond_status_var->getElementalValue(_current_elem);
 }

Referenced by PeridynamicsKernelBase::computeJacobian(), PeridynamicsKernelBase::computeResidual(), and MechanicsBasePD::prepare().

Member Data Documentation

◆ _bond_dQdT

const MaterialProperty<Real>& HeatConductionBPD::_bond_dQdT

protected

Definition at line 33 of file HeatConductionBPD.h.

Referenced by computeLocalJacobian().

◆ _bond_heat_flow

const MaterialProperty<Real>& HeatConductionBPD::_bond_heat_flow

protected

Bond based material properties.

Definition at line 32 of file HeatConductionBPD.h.

Referenced by computeLocalResidual().

◆ _bond_status_ij

Real PeridynamicsKernelBase::_bond_status_ij

protectedinherited

Bond status of current bond/edge2.

Definition at line 74 of file PeridynamicsKernelBase.h.

Referenced by computeLocalJacobian(), computeLocalResidual(), and PeridynamicsKernelBase::prepare().

◆ _bond_status_var

MooseVariable* PeridynamicsKernelBase::_bond_status_var

protectedinherited

Bond_status variable.

Definition at line 56 of file PeridynamicsKernelBase.h.

Referenced by PeridynamicsKernelBase::prepare().

◆ _dg_vol_frac_ij

std::vector<Real> PeridynamicsKernelBase::_dg_vol_frac_ij

protectedinherited

Definition at line 66 of file PeridynamicsKernelBase.h.

Referenced by PeridynamicsKernelBase::prepare().

◆ _dim

const unsigned int PeridynamicsKernelBase::_dim

protectedinherited

Definition at line 63 of file PeridynamicsKernelBase.h.

◆ _horiz_rad

std::vector<Real> PeridynamicsKernelBase::_horiz_rad

protectedinherited

Definition at line 67 of file PeridynamicsKernelBase.h.

Referenced by PeridynamicsKernelBase::prepare().

◆ _nnodes

const unsigned int PeridynamicsKernelBase::_nnodes

protectedinherited

Definition at line 64 of file PeridynamicsKernelBase.h.

Referenced by HeatSourceBPD::computeLocalResidual(), PeridynamicsKernelBase::computeResidual(), and PeridynamicsKernelBase::prepare().

◆ _origin_vec_ij

RealGradient PeridynamicsKernelBase::_origin_vec_ij

protectedinherited

Vector for current bond under undefored configuration.

Definition at line 71 of file PeridynamicsKernelBase.h.

Referenced by PeridynamicsKernelBase::prepare().

◆ _pdmesh

PeridynamicsMesh& PeridynamicsKernelBase::_pdmesh

protectedinherited

Parameters for peridynamic mesh information.

Definition at line 62 of file PeridynamicsKernelBase.h.

Referenced by HeatSourceBPD::computeLocalResidual(), and PeridynamicsKernelBase::prepare().

◆ _use_full_jacobian

const bool PeridynamicsKernelBase::_use_full_jacobian

protectedinherited

Option to use full jacobian including nonlocal constribution or not.

Definition at line 59 of file PeridynamicsKernelBase.h.

Referenced by PeridynamicsKernelBase::computeJacobian().

◆ _vols_ij

std::vector<Real> PeridynamicsKernelBase::_vols_ij

protectedinherited

Definition at line 65 of file PeridynamicsKernelBase.h.

Referenced by HeatSourceBPD::computeLocalResidual(), and PeridynamicsKernelBase::prepare().

The documentation for this class was generated from the following files:

peridynamics/include/kernels/HeatConductionBPD.h
peridynamics/src/kernels/HeatConductionBPD.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HeatConductionBPD()

Member Function Documentation

◆ computeJacobian()

◆ computeLocalJacobian()

◆ computeLocalResidual()

◆ computeNonlocalResidual()

◆ computeQpResidual()

◆ computeResidual()

◆ prepare()

Member Data Documentation

◆ _bond_dQdT

◆ _bond_heat_flow

◆ _bond_status_ij

◆ _bond_status_var

◆ _dg_vol_frac_ij

◆ _dim

◆ _horiz_rad

◆ _nnodes

◆ _origin_vec_ij

◆ _pdmesh

◆ _use_full_jacobian

◆ _vols_ij