Material class for bond based peridynamic heat conduction model based on regular spatial discretization. More...

#include <ThermalConstantHorizonMaterialBPD.h>

Inheritance diagram for ThermalConstantHorizonMaterialBPD:

Public Member Functions
	ThermalConstantHorizonMaterialBPD (const InputParameters &parameters)

Protected Member Functions
virtual void	computePeridynamicsParams (const Real ave_thermal_conductivity) override
	Function to compute micro-conductivity. More...

virtual void	computeProperties () override

void	setupMeshRelatedData ()
	Function to setup mesh related data to be used in this class. More...

Protected Attributes
const MaterialProperty< Real > &	_thermal_conductivity
	Thermal conductivity. More...

Real	_Kij
	Micro-conductivity. More...


MooseVariable *	_temp_var
	Temperature variable and values. More...

std::vector< Real >	_temp


MaterialProperty< Real > &	_bond_heat_flow
	Material properties to be stored. More...

MaterialProperty< Real > &	_bond_dQdT


PeridynamicsMesh &	_pdmesh
	Mesh related information for material points of current bond/element. More...

const unsigned int	_dim

const unsigned int	_nnodes

std::vector< Real >	_horiz_rad

std::vector< Real >	_node_vol

std::vector< Real >	_horiz_vol

RealGradient	_origin_vec

Real	_origin_length

Detailed Description

Material class for bond based peridynamic heat conduction model based on regular spatial discretization.

Definition at line 23 of file ThermalConstantHorizonMaterialBPD.h.

Constructor & Destructor Documentation

◆ ThermalConstantHorizonMaterialBPD()

ThermalConstantHorizonMaterialBPD::ThermalConstantHorizonMaterialBPD ( const InputParameters & parameters )

Definition at line 27 of file ThermalConstantHorizonMaterialBPD.C.

   : ThermalMaterialBaseBPD(parameters)
 {
 }

Member Function Documentation

◆ computePeridynamicsParams()

void ThermalConstantHorizonMaterialBPD::computePeridynamicsParams ( const Real ave_thermal_conductivity )

overrideprotectedvirtual

Function to compute micro-conductivity.

Parameters

ave_thermal_conductivity average thermal conductivity for the current element

Implements ThermalMaterialBaseBPD.

Definition at line 34 of file ThermalConstantHorizonMaterialBPD.C.

 {
   _Kij = 3.0 * ave_thermal_conductivity / M_PI *
          (1.0 / (MathUtils::pow(_horiz_rad[0], _dim + 1)) +
           1.0 / (MathUtils::pow(_horiz_rad[1], _dim + 1)));
 }

◆ computeProperties()

void ThermalMaterialBaseBPD::computeProperties ( )

overrideprotectedvirtualinherited

Definition at line 42 of file ThermalMaterialBaseBPD.C.

 {
   setupMeshRelatedData(); // function from base class
  
   Real ave_thermal_conductivity = 0.0;
   for (unsigned int qp = 0; qp < _qrule->n_points(); ++qp)
     ave_thermal_conductivity += _thermal_conductivity[qp] * _JxW[qp] * _coord[qp];
   ave_thermal_conductivity /= _assembly.elemVolume();
  
   // nodal temperature
   _temp[0] = _temp_var->getNodalValue(*_current_elem->node_ptr(0));
   _temp[1] = _temp_var->getNodalValue(*_current_elem->node_ptr(1));
  
   // compute peridynamic micro-conductivity: _Kij
   computePeridynamicsParams(ave_thermal_conductivity);
  
   for (_qp = 0; _qp < _nnodes; ++_qp)
   {
     // residual term
     _bond_heat_flow[_qp] =
         _Kij * (_temp[1] - _temp[0]) / _origin_length * _node_vol[0] * _node_vol[1];
  
     // derivative of the residual term
     _bond_dQdT[_qp] = -_Kij / _origin_length * _node_vol[0] * _node_vol[1];
   }
 }

◆ setupMeshRelatedData()

void PeridynamicsMaterialBase::setupMeshRelatedData ( )

protectedinherited

Function to setup mesh related data to be used in this class.

Definition at line 34 of file PeridynamicsMaterialBase.C.

 {
   for (unsigned int i = 0; i < _nnodes; ++i)
   {
     _horiz_rad[i] = _pdmesh.getHorizon(_current_elem->node_id(i));
     _node_vol[i] = _pdmesh.getPDNodeVolume(_current_elem->node_id(i));
     _horiz_vol[i] = _pdmesh.getHorizVolume(_current_elem->node_id(i));
   }
  
   _origin_vec = _pdmesh.getPDNodeCoord(_current_elem->node_id(1)) -
                 _pdmesh.getPDNodeCoord(_current_elem->node_id(0));
   _origin_length = _origin_vec.norm();
 }