Compute the Allen-Cahn interface term with the weak form residual \( \left(\nabla (L \psi), 1/2 {\partial \kappa} \over {\partial \nabla \eta_{\alpha i}} \sum \{(\nabla \eta_{\beta j})^2 \} \right) \). More...

#include <ACInterface2DMultiPhase1.h>

Inheritance diagram for ACInterface2DMultiPhase1:

[legend]

Public Member Functions
	ACInterface2DMultiPhase1 (const InputParameters &parameters)

virtual void	initialSetup ()

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

virtual Real	computeQpOffDiagJacobian (unsigned int jvar)

Real	sumSquareGradEta ()

RealGradient	gradL ()

RealGradient	gradKappa ()

RealGradient	nablaLPsi ()
	the \( \nabla(L\psi) \) term More...

RealGradient	kappaNablaLPsi ()
	the \( \kappa\nabla(L\psi) \) term More...

Protected Attributes
const MaterialProperty< RealGradient > &	_dkappadgrad_etaa
	Interfacial parameter. More...

const MaterialProperty< RealTensorValue > &	_d2kappadgrad_etaa

unsigned int	_num_etas
	Order parameters. More...

std::vector< const VariableValue * >	_eta

std::vector< const VariableGradient * >	_grad_eta

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

const MaterialProperty< Real > &	_kappa
	Interfacial parameter. More...

const bool	_variable_L
	flag set if L is a function of non-linear variables in args More...

const MaterialProperty< Real > &	_dkappadop
	kappa derivative w.r.t. order parameter More...

const unsigned int	_nvar
	number of coupled variables More...

std::vector< const MaterialProperty< Real > * >	_dkappadarg
	kappa derivative w.r.t. other coupled variables More...

std::vector< const VariableGradient * >	_gradarg
	Gradients for all coupled variables. More...


const MaterialProperty< Real > &	_dLdop
	Mobility derivatives w.r.t. order parameter. More...

const MaterialProperty< Real > &	_d2Ldop2


std::vector< const MaterialProperty< Real > * >	_dLdarg
	Mobility derivative w.r.t. other coupled variables. More...

std::vector< const MaterialProperty< Real > * >	_d2Ldargdop

std::vector< std::vector< const MaterialProperty< Real > * > >	_d2Ldarg2

Detailed Description

Compute the Allen-Cahn interface term with the weak form residual \( \left(\nabla (L \psi), 1/2 {\partial \kappa} \over {\partial \nabla \eta_{\alpha i}} \sum \{(\nabla \eta_{\beta j})^2 \} \right) \).

Definition at line 27 of file ACInterface2DMultiPhase1.h.

Constructor & Destructor Documentation

◆ ACInterface2DMultiPhase1()

ACInterface2DMultiPhase1::ACInterface2DMultiPhase1 ( const InputParameters & parameters )

Definition at line 34 of file ACInterface2DMultiPhase1.C.

   : ACInterface(parameters),
     _dkappadgrad_etaa(getMaterialProperty<RealGradient>("dkappadgrad_etaa_name")),
     _d2kappadgrad_etaa(getMaterialProperty<RealTensorValue>("d2kappadgrad_etaa_name")),
     _num_etas(coupledComponents("etas")),
     _eta(_num_etas),
     _grad_eta(_num_etas)
 {
 }

Member Function Documentation

◆ computeQpJacobian()

Real ACInterface2DMultiPhase1::computeQpJacobian ( )

protectedvirtual

Reimplemented from ACInterface.

Definition at line 64 of file ACInterface2DMultiPhase1.C.

 {
   // dsum is the derivative \f$ \frac\partial{\partial \eta} \left( \nabla (L\psi) \right) \f$
   RealGradient dsum = _dLdop[_qp] * _phi[_j][_qp] * _grad_test[_i][_qp];
  
   // compute the derivative of the gradient of the mobility
   if (_variable_L)
   {
     RealGradient dgradL =
         _grad_phi[_j][_qp] * _dLdop[_qp] + _grad_u[_qp] * _phi[_j][_qp] * _d2Ldop2[_qp];
  
     for (unsigned int i = 0; i < _nvar; ++i)
       dgradL += (*_gradarg[i])[_qp] * _phi[_j][_qp] * (*_d2Ldargdop[i])[_qp];
  
     dsum += dgradL * _test[_i][_qp];
   }
   Real jac1 = 0.5 * dsum * _dkappadgrad_etaa[_qp] * sumSquareGradEta();
   Real jac2 =
       0.5 * nablaLPsi() * (_d2kappadgrad_etaa[_qp] * _grad_phi[_j][_qp]) * sumSquareGradEta();
   Real jac3 = nablaLPsi() * _dkappadgrad_etaa[_qp] * _grad_u[_qp] * _grad_phi[_j][_qp];
   return jac1 + jac2 + jac3;
 }

◆ computeQpOffDiagJacobian()

Real ACInterface2DMultiPhase1::computeQpOffDiagJacobian ( unsigned int jvar )

protectedvirtual

Reimplemented from ACInterface.

Definition at line 88 of file ACInterface2DMultiPhase1.C.

 {
   // get the coupled variable jvar is referring to
   const unsigned int cvar = mapJvarToCvar(jvar);
  
   // dsum is the derivative \f$ \frac\partial{\partial \eta} \left( \nabla (L\psi) \right) \f$
   RealGradient dsum = (*_dLdarg[cvar])[_qp] * _phi[_j][_qp] * _grad_test[_i][_qp];
  
   // compute the derivative of the gradient of the mobility
   if (_variable_L)
   {
     RealGradient dgradL = _grad_phi[_j][_qp] * (*_dLdarg[cvar])[_qp] +
                           _grad_u[_qp] * _phi[_j][_qp] * (*_d2Ldargdop[cvar])[_qp];
  
     for (unsigned int i = 0; i < _nvar; ++i)
       dgradL += (*_gradarg[i])[_qp] * _phi[_j][_qp] * (*_d2Ldarg2[cvar][i])[_qp];
  
     dsum += dgradL * _test[_i][_qp];
   }
  
   // the gradient of coupled variable etab
   _grad_eta[0] = &coupledGradient("etas", 0);
  
   Real jac1 = 0.5 * dsum * _dkappadgrad_etaa[_qp] * sumSquareGradEta();
   Real jac2 =
       -0.5 * nablaLPsi() * (_d2kappadgrad_etaa[_qp] * _grad_phi[_j][_qp]) * sumSquareGradEta();
   Real jac3 = nablaLPsi() * _dkappadgrad_etaa[_qp] * (*_grad_eta[0])[_qp] * _grad_phi[_j][_qp];
   return jac1 + jac2 + jac3;
 }

◆ computeQpResidual()

Real ACInterface2DMultiPhase1::computeQpResidual ( )

protectedvirtual

Reimplemented from ACInterface.

Definition at line 58 of file ACInterface2DMultiPhase1.C.

 {
   return 0.5 * nablaLPsi() * _dkappadgrad_etaa[_qp] * sumSquareGradEta();
 }

◆ gradKappa()

RealGradient ACInterface::gradKappa ( )

protectedinherited

◆ gradL()

RealGradient ACInterface::gradL ( )

protectedinherited

Definition at line 77 of file ACInterface.C.

 {
   RealGradient g = _grad_u[_qp] * _dLdop[_qp];
   for (unsigned int i = 0; i < _nvar; ++i)
     g += (*_gradarg[i])[_qp] * (*_dLdarg[i])[_qp];
   return g;
 }

Referenced by ACInterface::computeQpJacobian(), ACInterface::computeQpOffDiagJacobian(), and ACInterface::nablaLPsi().

◆ initialSetup()

void ACInterface::initialSetup ( )

virtualinherited

Definition at line 70 of file ACInterface.C.

 {
   validateCoupling<Real>("mob_name");
   validateCoupling<Real>("kappa_name");
 }

◆ kappaNablaLPsi()

RealGradient ACInterface::kappaNablaLPsi ( )

protectedinherited

the \( \kappa\nabla(L\psi) \) term

Definition at line 98 of file ACInterface.C.

 {
   return _kappa[_qp] * nablaLPsi();
 }

Referenced by ACInterface::computeQpJacobian(), and ACInterface::computeQpResidual().

◆ nablaLPsi()

RealGradient ACInterface::nablaLPsi ( )

protectedinherited

the \( \nabla(L\psi) \) term

Definition at line 86 of file ACInterface.C.

 {
   // sum is the product rule gradient \f$ \nabla (L\psi) \f$
   RealGradient sum = _L[_qp] * _grad_test[_i][_qp];
  
   if (_variable_L)
     sum += gradL() * _test[_i][_qp];
  
   return sum;
 }

Referenced by computeQpJacobian(), ACInterface2DMultiPhase2::computeQpJacobian(), ACInterface2DMultiPhase2::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), computeQpResidual(), and ACInterface::kappaNablaLPsi().

◆ sumSquareGradEta()

Real ACInterface2DMultiPhase1::sumSquareGradEta ( )

protected

Definition at line 45 of file ACInterface2DMultiPhase1.C.

 {
   // get the sum of square of gradients of all order parameters
   Real SumSquareGradOp = _grad_u[_qp] * _grad_u[_qp];
   for (unsigned int i = 0; i < _num_etas; ++i)
   {
     _grad_eta[i] = &coupledGradient("etas", i);
     SumSquareGradOp += (*_grad_eta[i])[_qp] * (*_grad_eta[i])[_qp];
   }
   return SumSquareGradOp;
 }

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

Member Data Documentation

◆ _d2kappadgrad_etaa

const MaterialProperty<RealTensorValue>& ACInterface2DMultiPhase1::_d2kappadgrad_etaa

protected

Definition at line 41 of file ACInterface2DMultiPhase1.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ _d2Ldarg2

std::vector<std::vector<const MaterialProperty<Real> *> > ACInterface::_d2Ldarg2

protectedinherited

Definition at line 67 of file ACInterface.h.

Referenced by ACInterface::ACInterface(), ACInterface::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), and ACInterface2DMultiPhase2::computeQpOffDiagJacobian().

◆ _d2Ldargdop

std::vector<const MaterialProperty<Real> *> ACInterface::_d2Ldargdop

protectedinherited

Definition at line 66 of file ACInterface.h.

Referenced by ACInterface::ACInterface(), ACInterface::computeQpJacobian(), ACInterface2DMultiPhase2::computeQpJacobian(), computeQpJacobian(), ACInterface::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), and ACInterface2DMultiPhase2::computeQpOffDiagJacobian().

◆ _d2Ldop2

const MaterialProperty<Real>& ACInterface::_d2Ldop2

protectedinherited

Definition at line 55 of file ACInterface.h.

Referenced by ACInterface::computeQpJacobian(), computeQpJacobian(), and ACInterface2DMultiPhase2::computeQpJacobian().

◆ _dkappadarg

std::vector<const MaterialProperty<Real> *> ACInterface::_dkappadarg

protectedinherited

kappa derivative w.r.t. other coupled variables

Definition at line 71 of file ACInterface.h.

Referenced by ACInterface::ACInterface(), and ACInterface::computeQpOffDiagJacobian().

◆ _dkappadgrad_etaa

const MaterialProperty<RealGradient>& ACInterface2DMultiPhase1::_dkappadgrad_etaa

protected

Interfacial parameter.

Definition at line 40 of file ACInterface2DMultiPhase1.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

◆ _dkappadop

const MaterialProperty<Real>& ACInterface::_dkappadop

protectedinherited

kappa derivative w.r.t. order parameter

Definition at line 59 of file ACInterface.h.

Referenced by ACInterface::computeQpJacobian(), and ACInterface::computeQpOffDiagJacobian().

◆ _dLdarg

std::vector<const MaterialProperty<Real> *> ACInterface::_dLdarg

protectedinherited

Mobility derivative w.r.t. other coupled variables.

Definition at line 65 of file ACInterface.h.

Referenced by ACInterface::ACInterface(), ACInterface::computeQpOffDiagJacobian(), ACInterface2DMultiPhase2::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), and ACInterface::gradL().

◆ _dLdop

const MaterialProperty<Real>& ACInterface::_dLdop

protectedinherited

Mobility derivatives w.r.t. order parameter.

Definition at line 54 of file ACInterface.h.

Referenced by ACInterface::computeQpJacobian(), ACInterface2DMultiPhase2::computeQpJacobian(), computeQpJacobian(), and ACInterface::gradL().

◆ _eta

std::vector<const VariableValue *> ACInterface2DMultiPhase1::_eta

protected

Definition at line 45 of file ACInterface2DMultiPhase1.h.

◆ _grad_eta

std::vector<const VariableGradient *> ACInterface2DMultiPhase1::_grad_eta

protected

Definition at line 46 of file ACInterface2DMultiPhase1.h.

Referenced by computeQpOffDiagJacobian(), and sumSquareGradEta().

◆ _gradarg

std::vector<const VariableGradient *> ACInterface::_gradarg

protectedinherited

Gradients for all coupled variables.

Definition at line 74 of file ACInterface.h.

Referenced by ACInterface::ACInterface(), ACInterface::computeQpJacobian(), ACInterface2DMultiPhase2::computeQpJacobian(), computeQpJacobian(), ACInterface::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), ACInterface2DMultiPhase2::computeQpOffDiagJacobian(), and ACInterface::gradL().

◆ _kappa

const MaterialProperty<Real>& ACInterface::_kappa

protectedinherited

Interfacial parameter.

Definition at line 48 of file ACInterface.h.

Referenced by ACInterface::computeQpJacobian(), ACInterface2DMultiPhase2::computeQpJacobian(), ACInterface::computeQpOffDiagJacobian(), ACInterface2DMultiPhase2::computeQpOffDiagJacobian(), and ACInterface::kappaNablaLPsi().

◆ _L

const MaterialProperty<Real>& ACInterface::_L

protectedinherited

Mobility.

Definition at line 46 of file ACInterface.h.

Referenced by ACInterface::computeQpJacobian(), ACInterface::computeQpOffDiagJacobian(), and ACInterface::nablaLPsi().

◆ _num_etas

unsigned int ACInterface2DMultiPhase1::_num_etas

protected

Order parameters.

Definition at line 44 of file ACInterface2DMultiPhase1.h.

Referenced by sumSquareGradEta().

◆ _nvar

const unsigned int ACInterface::_nvar

protectedinherited

number of coupled variables

Definition at line 62 of file ACInterface.h.

Referenced by ACInterface::ACInterface(), ACInterface::computeQpJacobian(), ACInterface2DMultiPhase2::computeQpJacobian(), computeQpJacobian(), ACInterface::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), ACInterface2DMultiPhase2::computeQpOffDiagJacobian(), and ACInterface::gradL().

◆ _variable_L

const bool ACInterface::_variable_L

protectedinherited

flag set if L is a function of non-linear variables in args

Definition at line 51 of file ACInterface.h.

Referenced by ACInterface::computeQpJacobian(), computeQpJacobian(), ACInterface2DMultiPhase2::computeQpJacobian(), ACInterface::computeQpOffDiagJacobian(), ACInterface2DMultiPhase2::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), and ACInterface::nablaLPsi().

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

phase_field/include/kernels/ACInterface2DMultiPhase1.h
phase_field/src/kernels/ACInterface2DMultiPhase1.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ACInterface2DMultiPhase1()

Member Function Documentation

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ gradKappa()

◆ gradL()

◆ initialSetup()

◆ kappaNablaLPsi()

◆ nablaLPsi()

◆ sumSquareGradEta()

Member Data Documentation

◆ _d2kappadgrad_etaa

◆ _d2Ldarg2

◆ _d2Ldargdop

◆ _d2Ldop2

◆ _dkappadarg

◆ _dkappadgrad_etaa

◆ _dkappadop

◆ _dLdarg

◆ _dLdop

◆ _eta

◆ _grad_eta

◆ _gradarg

◆ _kappa

◆ _L

◆ _num_etas

◆ _nvar

◆ _variable_L