Isotropic diffusion kernel that takes a diffusion coefficient of type Real. More...

#include <PolycrystalMatDiffusion.h>

Inheritance diagram for PolycrystalMatDiffusion:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

typedef std::vector< int >	JvarMap

Public Member Functions
	PolycrystalMatDiffusion (const InputParameters &parameters)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

void	validateDerivativeMaterialPropertyBase (const std::string &base)

virtual void	computeOffDiagJacobian (unsigned int jvar) override

unsigned int	mapJvarToCvar (unsigned int jvar)

int	mapJvarToCvar (unsigned int jvar, const JvarMap &jvar_map)

bool	mapJvarToCvar (unsigned int jvar, unsigned int &cvar)

const JvarMap &	getJvarMap ()

const JvarMap &	getParameterJvarMap (std::string parameter_name)

const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

virtual void	initialSetup () override

virtual RealGradient	precomputeQpJacobian () override

virtual RealGradient	precomputeQpCJacobian ()

virtual GenericRealVectorValue< is_ad >	precomputeQpResidual () override

Protected Attributes
const unsigned int	_op_num
	Number of grain order parameters for cases where there is dependence on grain OP gradients. More...

const MaterialProperty< typename GradientType< Real >::type > &	_dDdgradc
	diffusion coefficient derivatives w.r.t. variables that have explicit dependence on gradients More...

std::vector< const MaterialProperty< typename GradientType< Real >::type > *>	_dDdgradeta

unsigned int	_surface_op_var
	For solid-pore systems, mame of the order parameter identifies the solid-pore surface. More...

const bool	_add_grain_op_gradients
	Variable to allow user to control whether grain OP gradient contributions are added to Jacobian. More...

const MaterialProperty< T > &	_ddiffusivity_dc

std::vector< const MaterialProperty< T > * >	_ddiffusivity_darg

const bool	_is_coupled

unsigned int	_v_var

const GenericMaterialProperty< T, is_ad > &	_diffusivity

const GenericVariableGradient< is_ad > &	_grad_v

	usingGenericKernelGradMembers

const unsigned int	_n_args

Detailed Description

Isotropic diffusion kernel that takes a diffusion coefficient of type Real.

Extends MatDiffusion to allow improved numerical performance in polycrystalline microstructures by adding Jacobian contributions for cases where diffusivity depends on gradients of order parameters.

Definition at line 20 of file PolycrystalMatDiffusion.h.

Constructor & Destructor Documentation

◆ PolycrystalMatDiffusion()

PolycrystalMatDiffusion::PolycrystalMatDiffusion ( const InputParameters & parameters )

Definition at line 23 of file PolycrystalMatDiffusion.C.

   : PolycrystalMatDiffusionBase<Real>(parameters)
 {
 }

Member Function Documentation

◆ computeQpOffDiagJacobian()

Real PolycrystalMatDiffusionBase< Real >::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtualinherited

Reimplemented from MatDiffusionBase< Real >.

Definition at line 98 of file PolycrystalMatDiffusionBase.h.

 {
   Real sum = MatDiffusionBase<T>::computeQpOffDiagJacobian(jvar);
   if (jvar == _surface_op_var)
     sum += _dDdgradc[this->_qp] * this->_grad_phi[this->_j][this->_qp] * this->_grad_v[this->_qp] *
            this->_grad_test[this->_i][this->_qp];
 
   if (_add_grain_op_gradients)
   {
     for (unsigned int k = 0; k < _op_num; ++k)
       if (jvar == this->coupled("grain_op_vars", k))
         sum += (*_dDdgradeta[k])[this->_qp] * this->_grad_phi[this->_j][this->_qp] *
                this->_grad_v[this->_qp] * this->_grad_test[this->_i][this->_qp];
   }
 
   return sum;
 }

◆ validParams()

InputParameters PolycrystalMatDiffusion::validParams ( )

static

Definition at line 15 of file PolycrystalMatDiffusion.C.

 {
   InputParameters params = PolycrystalMatDiffusionBase<Real>::validParams();
   params.addClassDescription(
       "Diffusion equation Kernel that takes an isotropic Diffusivity from a material property");
   return params;
 }