Kernel class to implement hear source term for peridynamic heat conduction models. More...

#include <HeatSourceBPD.h>

Inheritance diagram for HeatSourceBPD:

Public Member Functions
	HeatSourceBPD (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	computeNonlocalResidual ()
	Function to compute nonlocal contribution to the residual at the current nodes. More...

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

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

Protected Attributes
const Function &	_power_density
	Power density function. More...

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...


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 to implement hear source term for peridynamic heat conduction models.

Definition at line 22 of file HeatSourceBPD.h.

Constructor & Destructor Documentation

◆ HeatSourceBPD()

HeatSourceBPD::HeatSourceBPD ( const InputParameters & parameters )

Definition at line 29 of file HeatSourceBPD.C.

   : PeridynamicsKernelBase(parameters), _power_density(getFunction("power_density"))
 {
 }

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()

virtual void PeridynamicsKernelBase::computeLocalJacobian ( )

inlineprotectedvirtualinherited

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

Reimplemented in HeatConductionBPD.

Definition at line 48 of file PeridynamicsKernelBase.h.

48 {};

Referenced by PeridynamicsKernelBase::computeJacobian().

◆ computeLocalResidual()

void HeatSourceBPD::computeLocalResidual ( )

overrideprotectedvirtual

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

Implements PeridynamicsKernelBase.

Definition at line 35 of file HeatSourceBPD.C.

 {
   Point p;
   for (unsigned int i = 0; i < _nnodes; ++i)
   {
     std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(i));
     _local_re(i) = -_power_density.value(_t, p) * _vols_ij[i] / neighbors.size();
   }
 }

◆ 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_status_ij

Real PeridynamicsKernelBase::_bond_status_ij

protectedinherited

Bond status of current bond/edge2.

Definition at line 74 of file PeridynamicsKernelBase.h.

Referenced by HeatConductionBPD::computeLocalJacobian(), HeatConductionBPD::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 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 computeLocalResidual(), and PeridynamicsKernelBase::prepare().

◆ _power_density

const Function& HeatSourceBPD::_power_density

protected

Power density function.

Definition at line 31 of file HeatSourceBPD.h.

Referenced by computeLocalResidual().

◆ _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 computeLocalResidual(), and PeridynamicsKernelBase::prepare().

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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HeatSourceBPD()

Member Function Documentation

◆ computeJacobian()

◆ computeLocalJacobian()

◆ computeLocalResidual()

◆ computeNonlocalResidual()

◆ computeQpResidual()

◆ computeResidual()

◆ prepare()

Member Data Documentation

◆ _bond_status_ij

◆ _bond_status_var

◆ _dg_vol_frac_ij

◆ _dim

◆ _horiz_rad

◆ _nnodes

◆ _origin_vec_ij

◆ _pdmesh

◆ _power_density

◆ _use_full_jacobian

◆ _vols_ij