KKSACBulkBase child class for the phase concentration term \( - \sum_j \frac{dF_1}{dc_1} \frac{dh_j}{d\eta_i} (c_j) \) in the the Allen-Cahn bulk residual. More...

#include <KKSMultiACBulkC.h>

Inheritance diagram for KKSMultiACBulkC:

[legend]

Public Member Functions
	KKSMultiACBulkC (const InputParameters &parameters)

virtual void	initialSetup ()

Static Public Member Functions
static InputParameters	validParams ()

Protected Types
enum	PFFunctionType

Protected Member Functions
virtual Real	computeDFDOP (PFFunctionType type)

virtual Real	computeQpOffDiagJacobian (unsigned int jvar)

virtual Real	precomputeQpResidual ()

virtual Real	precomputeQpJacobian ()

virtual Real	computeDFDOP (PFFunctionType type)=0

Protected Attributes
MaterialPropertyName	_c1_name
	Names of phase concentration variables. More...

std::vector< const VariableValue * >	_cjs

std::vector< unsigned int >	_cjs_var

const MaterialProperty< Real > &	_prop_dF1dc1
	Derivative of the free energy function \( \frac d{dc_1} F_1 \). More...

const MaterialProperty< Real > &	_prop_d2F1dc12
	Second derivative of the free energy function \( \frac {d^2}{dc_1^2} F_1 \). More...

std::vector< const MaterialProperty< Real > * >	_prop_d2F1dc1darg
	Mixed partial derivatives of the free energy function wrt c1 and any other coupled variables \( \frac {d^2}{dc_1 dv} F_1 \). More...

unsigned int	_nvar
	Number of coupled variables. More...

VariableName	_etai_name
	name of order parameter that derivatives are taken wrt (needed to retrieve the derivative material properties) More...

unsigned int	_etai_var
	index of order parameter that derivatives are taken wrt More...

std::vector< MaterialPropertyName >	_Fj_names
	Names of free energy functions for each phase \( F_j \). More...

unsigned int	_num_j

std::vector< const MaterialProperty< Real > * >	_prop_Fj
	Values of the free energy functions for each phase \( F_j \). More...

std::vector< std::vector< const MaterialProperty< Real > * > >	_prop_dFjdarg
	Derivatives of the free energy functions (needed for off-diagonal Jacobians) More...

std::vector< MaterialPropertyName >	_hj_names
	switching function names More...

std::vector< const MaterialProperty< Real > * >	_prop_hj
	Values of the switching functions for each phase \( h_j \). More...

std::vector< const MaterialProperty< Real > * >	_prop_dhjdetai
	Derivatives of the switching functions wrt the order parameter for this kernel. More...

std::vector< const MaterialProperty< Real > * >	_prop_d2hjdetai2
	Second derivatives of the switching functions wrt the order parameter for this kernel. More...

std::vector< std::vector< const MaterialProperty< Real > * > >	_prop_d2hjdetaidarg
	Second derivatives of the switching functions (needed for off-diagonal Jacobians) More...

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

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

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

Detailed Description

KKSACBulkBase child class for the phase concentration term \( - \sum_j \frac{dF_1}{dc_1} \frac{dh_j}{d\eta_i} (c_j) \) in the the Allen-Cahn bulk residual.

The non-linear variable for this Kernel is the order parameter 'eta_i'.

Definition at line 27 of file KKSMultiACBulkC.h.

Member Enumeration Documentation

◆ PFFunctionType

enum ACBulk::PFFunctionType

protectedinherited

Definition at line 37 of file ACBulk.h.

   {
     Jacobian,
     Residual
   };

Constructor & Destructor Documentation

◆ KKSMultiACBulkC()

KKSMultiACBulkC::KKSMultiACBulkC ( const InputParameters & parameters )

Definition at line 26 of file KKSMultiACBulkC.C.

   : KKSMultiACBulkBase(parameters),
     _c1_name(getVar("cj_names", 0)
                  ->name()), // Can use any dFj/dcj since they are equal so pick first cj in the list
     _cjs(_num_j),
     _cjs_var(_num_j),
     _prop_dF1dc1(getMaterialPropertyDerivative<Real>(_Fj_names[0],
                                                      _c1_name)), // Use first Fj in list for dFj/dcj
     _prop_d2F1dc12(getMaterialPropertyDerivative<Real>(_Fj_names[0], _c1_name, _c1_name))
 {
   if (_num_j != coupledComponents("cj_names"))
     paramError("cj_names", "Need to pass in as many cj_names as Fj_names");
  
   // Load concentration variables into the arrays
   for (unsigned int i = 0; i < _num_j; ++i)
   {
     _cjs[i] = &coupledValue("cj_names", i);
     _cjs_var[i] = coupled("cj_names", i);
   }
  
   // Resize to number of coupled variables (_nvar from KKSMultiACBulkBase constructor)
   _prop_d2F1dc1darg.resize(_nvar);
  
   // Iterate over all coupled variables
   for (unsigned int i = 0; i < _nvar; ++i)
   {
     MooseVariableFEBase * cvar = _coupled_moose_vars[i];
  
     // get second partial derivatives wrt c1 and other coupled variable
     _prop_d2F1dc1darg[i] =
         &getMaterialPropertyDerivative<Real>(_Fj_names[0], _c1_name, cvar->name());
   }
 }

Member Function Documentation

◆ computeDFDOP() [1/2]

Real KKSMultiACBulkC::computeDFDOP ( PFFunctionType type )

protectedvirtual

Definition at line 61 of file KKSMultiACBulkC.C.

 {
   Real sum = 0.0;
  
   switch (type)
   {
     case Residual:
       for (unsigned int n = 0; n < _num_j; ++n)
         sum += (*_prop_dhjdetai[n])[_qp] * (*_cjs[n])[_qp];
  
       return -_prop_dF1dc1[_qp] * sum;
  
     case Jacobian:
       // For when this kernel is used in the Lagrange multiplier equation
       // In that case the Lagrange multiplier is the nonlinear variable
       if (_etai_var != _var.number())
         return 0.0;
  
       // For when eta_i is the nonlinear variable
       for (unsigned int n = 0; n < _num_j; ++n)
         sum += (*_prop_d2hjdetai2[n])[_qp] * (*_cjs[n])[_qp];
  
       return -_phi[_j][_qp] * _prop_dF1dc1[_qp] * sum;
   }
  
   mooseError("Invalid type passed in");
 }

◆ computeDFDOP() [2/2]

virtual Real ACBulk< Real >::computeDFDOP ( PFFunctionType type )

protectedpure virtualinherited

◆ computeQpOffDiagJacobian()

Real KKSMultiACBulkC::computeQpOffDiagJacobian ( unsigned int jvar )

protectedvirtual

Reimplemented from ACBulk< Real >.

Definition at line 90 of file KKSMultiACBulkC.C.

 {
   // first get dependence of mobility _L on other variables using parent class
   // member function
   Real res = ACBulk<Real>::computeQpOffDiagJacobian(jvar);
  
   Real sum = 0.0;
   // Then add dependence of KKSACBulkC on other variables
   // Treat cj variables specially, as they appear in the residual
   if (jvar == _cjs_var[0])
   {
     for (unsigned int n = 0; n < _num_j; ++n)
       sum += (*_prop_dhjdetai[n])[_qp] * (*_cjs[n])[_qp];
  
     res -= _L[_qp] * (sum * _prop_d2F1dc12[_qp] + _prop_dF1dc1[_qp] * (*_prop_dhjdetai[0])[_qp]) *
            _phi[_j][_qp] * _test[_i][_qp];
     return res;
   }
  
   for (unsigned int i = 1; i < _num_j; ++i)
   {
     if (jvar == _cjs_var[i])
     {
       res -=
           _L[_qp] * _prop_dF1dc1[_qp] * (*_prop_dhjdetai[i])[_qp] * _phi[_j][_qp] * _test[_i][_qp];
       return res;
     }
   }
  
   //  for all other vars get the coupled variable jvar is referring to
   const unsigned int cvar = mapJvarToCvar(jvar);
  
   for (unsigned int n = 0; n < _num_j; ++n)
     sum += _prop_dF1dc1[_qp] * (*_prop_d2hjdetaidarg[n][cvar])[_qp] * (*_cjs[n])[_qp] +
            (*_prop_d2F1dc1darg[cvar])[_qp] * (*_prop_dhjdetai[n])[_qp] * (*_cjs[n])[_qp];
  
   res -= _L[_qp] * sum * _phi[_j][_qp] * _test[_i][_qp];
  
   return res;
 }

◆ initialSetup()

void KKSMultiACBulkBase::initialSetup ( )

virtualinherited

Reimplemented from ACBulk< Real >.

Definition at line 76 of file KKSMultiACBulkBase.C.

 {
   ACBulk<Real>::initialSetup();
  
   for (unsigned int n = 0; n < _num_j; ++n)
   {
     validateNonlinearCoupling<Real>(_Fj_names[n]);
     validateNonlinearCoupling<Real>(_hj_names[n]);
   }
 }

◆ precomputeQpJacobian()

Real ACBulk< Real >::precomputeQpJacobian ( )

protectedvirtualinherited

Definition at line 103 of file ACBulk.h.

 {
   // Get free energy derivative and Jacobian
   Real dFdop = computeDFDOP(Residual);
  
   Real JdFdop = computeDFDOP(Jacobian);
  
   // Set Jacobian value using product rule
   return _L[_qp] * JdFdop + _dLdop[_qp] * _phi[_j][_qp] * dFdop;
 }

◆ precomputeQpResidual()

Real ACBulk< Real >::precomputeQpResidual ( )

protectedvirtualinherited

Definition at line 92 of file ACBulk.h.

 {
   // Get free energy derivative from function
   Real dFdop = computeDFDOP(Residual);
  
   // Set residual
   return _L[_qp] * dFdop;
 }

◆ validParams()

InputParameters ACBulk< Real >::validParams ( )

staticinherited

Definition at line 74 of file ACBulk.h.

 {
   InputParameters params = ::validParams<KernelValue>();
   params.addClassDescription("Allen-Cahn base Kernel");
   params.addParam<MaterialPropertyName>("mob_name", "L", "The mobility used with the kernel");
   params.addCoupledVar("args", "Vector of arguments of the mobility");
   return params;
 }

Member Data Documentation

◆ _c1_name

MaterialPropertyName KKSMultiACBulkC::_c1_name

protected

Names of phase concentration variables.

Definition at line 37 of file KKSMultiACBulkC.h.

Referenced by KKSMultiACBulkC().

◆ _cjs

std::vector<const VariableValue *> KKSMultiACBulkC::_cjs

protected

Definition at line 38 of file KKSMultiACBulkC.h.

Referenced by computeDFDOP(), computeQpOffDiagJacobian(), and KKSMultiACBulkC().

◆ _cjs_var

std::vector<unsigned int> KKSMultiACBulkC::_cjs_var

protected

Definition at line 39 of file KKSMultiACBulkC.h.

Referenced by computeQpOffDiagJacobian(), and KKSMultiACBulkC().

◆ _dLdarg

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

protectedinherited

Mobility derivative w.r.t coupled variables.

Definition at line 52 of file ACBulk.h.

◆ _dLdop

const MaterialProperty<Real >& ACBulk< Real >::_dLdop

protectedinherited

Mobility derivative w.r.t. order parameter.

Definition at line 49 of file ACBulk.h.

◆ _etai_name

VariableName KKSMultiACBulkBase::_etai_name

protectedinherited

name of order parameter that derivatives are taken wrt (needed to retrieve the derivative material properties)

Definition at line 39 of file KKSMultiACBulkBase.h.

Referenced by KKSMultiACBulkBase::KKSMultiACBulkBase().

◆ _etai_var

unsigned int KKSMultiACBulkBase::_etai_var

protectedinherited

index of order parameter that derivatives are taken wrt

Definition at line 42 of file KKSMultiACBulkBase.h.

Referenced by computeDFDOP(), KKSMultiACBulkF::computeDFDOP(), and KKSMultiACBulkF::computeQpOffDiagJacobian().

◆ _Fj_names

std::vector<MaterialPropertyName> KKSMultiACBulkBase::_Fj_names

protectedinherited

Names of free energy functions for each phase \( F_j \).

Definition at line 45 of file KKSMultiACBulkBase.h.

Referenced by KKSMultiACBulkBase::initialSetup(), KKSMultiACBulkBase::KKSMultiACBulkBase(), and KKSMultiACBulkC().

◆ _hj_names

std::vector<MaterialPropertyName> KKSMultiACBulkBase::_hj_names

protectedinherited

switching function names

Definition at line 55 of file KKSMultiACBulkBase.h.

Referenced by KKSMultiACBulkBase::initialSetup(), and KKSMultiACBulkBase::KKSMultiACBulkBase().

◆ _L

const MaterialProperty<Real >& ACBulk< Real >::_L

protectedinherited

Mobility.

Definition at line 46 of file ACBulk.h.

◆ _num_j

unsigned int KKSMultiACBulkBase::_num_j

protectedinherited

Definition at line 46 of file KKSMultiACBulkBase.h.

Referenced by computeDFDOP(), KKSMultiACBulkF::computeDFDOP(), KKSMultiACBulkF::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), KKSMultiACBulkBase::initialSetup(), KKSMultiACBulkBase::KKSMultiACBulkBase(), and KKSMultiACBulkC().

◆ _nvar

unsigned int KKSMultiACBulkBase::_nvar

protectedinherited

Number of coupled variables.

Definition at line 36 of file KKSMultiACBulkBase.h.

Referenced by KKSMultiACBulkBase::KKSMultiACBulkBase(), and KKSMultiACBulkC().

◆ _prop_d2F1dc12

const MaterialProperty<Real>& KKSMultiACBulkC::_prop_d2F1dc12

protected

Second derivative of the free energy function \( \frac {d^2}{dc_1^2} F_1 \).

Definition at line 44 of file KKSMultiACBulkC.h.

Referenced by computeQpOffDiagJacobian().

◆ _prop_d2F1dc1darg

std::vector<const MaterialProperty<Real> *> KKSMultiACBulkC::_prop_d2F1dc1darg

protected

Mixed partial derivatives of the free energy function wrt c1 and any other coupled variables \( \frac {d^2}{dc_1 dv} F_1 \).

Definition at line 47 of file KKSMultiACBulkC.h.

Referenced by computeQpOffDiagJacobian(), and KKSMultiACBulkC().

◆ _prop_d2hjdetai2

std::vector<const MaterialProperty<Real> *> KKSMultiACBulkBase::_prop_d2hjdetai2

protectedinherited

Second derivatives of the switching functions wrt the order parameter for this kernel.

Definition at line 64 of file KKSMultiACBulkBase.h.

Referenced by computeDFDOP(), KKSMultiACBulkF::computeDFDOP(), and KKSMultiACBulkBase::KKSMultiACBulkBase().

◆ _prop_d2hjdetaidarg

std::vector<std::vector<const MaterialProperty<Real> *> > KKSMultiACBulkBase::_prop_d2hjdetaidarg

protectedinherited

Second derivatives of the switching functions (needed for off-diagonal Jacobians)

Definition at line 67 of file KKSMultiACBulkBase.h.

Referenced by computeQpOffDiagJacobian(), KKSMultiACBulkF::computeQpOffDiagJacobian(), and KKSMultiACBulkBase::KKSMultiACBulkBase().

◆ _prop_dF1dc1

const MaterialProperty<Real>& KKSMultiACBulkC::_prop_dF1dc1

protected

Derivative of the free energy function \( \frac d{dc_1} F_1 \).

Definition at line 42 of file KKSMultiACBulkC.h.

Referenced by computeDFDOP(), and computeQpOffDiagJacobian().

◆ _prop_dFjdarg

std::vector<std::vector<const MaterialProperty<Real> *> > KKSMultiACBulkBase::_prop_dFjdarg

protectedinherited

Derivatives of the free energy functions (needed for off-diagonal Jacobians)

Definition at line 52 of file KKSMultiACBulkBase.h.

Referenced by KKSMultiACBulkF::computeQpOffDiagJacobian(), and KKSMultiACBulkBase::KKSMultiACBulkBase().

◆ _prop_dhjdetai

std::vector<const MaterialProperty<Real> *> KKSMultiACBulkBase::_prop_dhjdetai

protectedinherited

Derivatives of the switching functions wrt the order parameter for this kernel.

Definition at line 61 of file KKSMultiACBulkBase.h.

Referenced by computeDFDOP(), KKSMultiACBulkF::computeDFDOP(), computeQpOffDiagJacobian(), KKSMultiACBulkF::computeQpOffDiagJacobian(), and KKSMultiACBulkBase::KKSMultiACBulkBase().

◆ _prop_Fj

std::vector<const MaterialProperty<Real> *> KKSMultiACBulkBase::_prop_Fj

protectedinherited

Values of the free energy functions for each phase \( F_j \).

Definition at line 49 of file KKSMultiACBulkBase.h.

Referenced by KKSMultiACBulkF::computeDFDOP(), KKSMultiACBulkF::computeQpOffDiagJacobian(), and KKSMultiACBulkBase::KKSMultiACBulkBase().

◆ _prop_hj

std::vector<const MaterialProperty<Real> *> KKSMultiACBulkBase::_prop_hj

protectedinherited

Values of the switching functions for each phase \( h_j \).

Definition at line 58 of file KKSMultiACBulkBase.h.

Referenced by KKSMultiACBulkBase::KKSMultiACBulkBase().

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

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

Public Member Functions

Static Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ PFFunctionType

Constructor & Destructor Documentation

◆ KKSMultiACBulkC()

Member Function Documentation

◆ computeDFDOP() [1/2]

◆ computeDFDOP() [2/2]

◆ computeQpOffDiagJacobian()

◆ initialSetup()

◆ precomputeQpJacobian()

◆ precomputeQpResidual()

◆ validParams()

Member Data Documentation

◆ _c1_name

◆ _cjs

◆ _cjs_var

◆ _dLdarg

◆ _dLdop

◆ _etai_name

◆ _etai_var

◆ _Fj_names

◆ _hj_names

◆ _L

◆ _num_j

◆ _nvar

◆ _prop_d2F1dc12

◆ _prop_d2F1dc1darg

◆ _prop_d2hjdetai2

◆ _prop_d2hjdetaidarg

◆ _prop_dF1dc1

◆ _prop_dFjdarg

◆ _prop_dhjdetai

◆ _prop_Fj

◆ _prop_hj