This kernel computes the off-diagonal jacobian of elastic energy in AllenCahn respect to displacements. More...

#include <AllenCahnElasticEnergyOffDiag.h>

Inheritance diagram for AllenCahnElasticEnergyOffDiag:

Public Member Functions
	AllenCahnElasticEnergyOffDiag (const InputParameters &parameters)

Protected Member Functions
Real	computeQpResidual () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Protected Attributes
const MaterialProperty< Real > &	_L
	Mobility. More...

const MaterialProperty< RankTwoTensor > &	_d2Fdcdstrain
	Free energy material properties and derivatives. More...


const unsigned int	_ndisp
	Displacement variables used for off-diagonal Jacobian. More...

std::vector< unsigned int >	_disp_var

Detailed Description

This kernel computes the off-diagonal jacobian of elastic energy in AllenCahn respect to displacements.

Definition at line 29 of file AllenCahnElasticEnergyOffDiag.h.

Constructor & Destructor Documentation

◆ AllenCahnElasticEnergyOffDiag()

AllenCahnElasticEnergyOffDiag::AllenCahnElasticEnergyOffDiag ( const InputParameters & parameters )

Definition at line 31 of file AllenCahnElasticEnergyOffDiag.C.

   : DerivativeMaterialInterface<JvarMapKernelInterface<Kernel>>(parameters),
     _L(getMaterialProperty<Real>("mob_name")),
     _ndisp(coupledComponents("displacements")),
     _disp_var(_ndisp),
     _d2Fdcdstrain(getMaterialProperty<RankTwoTensor>("d2Fdcdstrain"))
 {
   for (unsigned int i = 0; i < _ndisp; ++i)
     _disp_var[i] = coupled("displacements", i);
 }

Member Function Documentation

◆ computeQpOffDiagJacobian()

Real AllenCahnElasticEnergyOffDiag::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 43 of file AllenCahnElasticEnergyOffDiag.C.

 {
   for (unsigned int c_comp = 0; c_comp < _ndisp; ++c_comp)
     if (jvar == _disp_var[c_comp])
     {
       const Real dxddFdc = _L[_qp] * _test[_i][_qp];
       const Real d2Fdcdstrain_comp =
           (_d2Fdcdstrain[_qp].column(c_comp) + _d2Fdcdstrain[_qp].row(c_comp)) / 2.0 *
           _grad_phi[_j][_qp];
       return dxddFdc * d2Fdcdstrain_comp;
     }
  
   return 0.0;
 }